Dynamic stochastic copula models: estimation,inference and applications

We propose a new dynamic copula model in which the parameter characterizing dependence follows an autoregressive process. As this model class includes the Gaussian copula with stochastic correlation process, it can be viewed as a generalization of multivariate stochastic volatility models. Despite the complexity of the model, the decoupling of marginals and dependence parameters facilitates estimation. We propose estimation in two steps, where first the parameters of the marginal distributions are estimated, and then those of the copula. Parameters of the latent processes (volatilities and dependence) are estimated using efficient importance sampling. We discuss goodness‐of‐fit tests and ways to forecast the dependence parameter. For two bivariate stock index series, we show that the proposed model outperforms standard competing models. Copyright © 2010 John Wiley & Sons, Ltd.     

Modelling dependence using skew t copulas: Bayesian inference and applications

We construct a copula from the skew t distribution of Sahu et al. ( 2003 ). This copula can capture asymmetric and extreme dependence between variables, and is one of the few copulas that can do so and still be used in high dimensions effectively. However, it is difficult to estimate the copula model by maximum likelihood when the multivariate dimension is high, or when some or all of the marginal distributions are discrete‐valued, or when the parameters in the marginal distributions and copula are estimated jointly. We therefore propose a Bayesian approach that overcomes all these problems. The computations are undertaken using a Markov chain Monte Carlo simulation method which exploits the conditionally Gaussian representation of the skew t distribution. We employ the approach in two contemporary econometric studies. The first is the modelling of regional spot prices in the Australian electricity market. Here, we observe complex non‐Gaussian margins and nonlinear inter‐regional dependence. Accurate characterization of this dependence is important for the study of market integration and risk management purposes. The second is the modelling of ordinal exposure measures for 15 major websites. Dependence between websites is important when measuring the impact of multi‐site advertising campaigns. In both cases the skew t copula substantially outperforms symmetric elliptical copula alternatives, demonstrating that the skew t copula is a powerful modelling tool when coupled with Bayesian inference. Copyright © 2010 John Wiley & Sons, Ltd.     

Inversion copulas from nonlinear state space models with an application to inflation forecasting

We propose the construction of copulas through the inversion of nonlinear state space models. These copulas allow for new time series models that have the same serial dependence structure as a state space model, but with an arbitrary marginal distribution, and flexible density forecasts. We examine the time series properties of the copulas, outline serial dependence measures, and estimate the models using likelihood-based methods. Copulas constructed from three example state space models are considered: a stochastic volatility model with an unobserved component, a Markov switching autoregression, and a Gaussian linear unobserved component model. We show that all three inversion copulas with flexible margins improve the fit and density forecasts of quarterly U.S. broad inflation and electricity inflation.     

Copula-based multivariate GARCH model with uncorrelated dependent errors

Multivariate GARCH (MGARCH) models are usually estimated under multivariate normality. In this paper, for non-elliptically distributed financial returns, we propose copula-based multivariate GARCH (C-MGARCH) model with uncorrelated dependent errors, which are generated through a linear combination of dependent random variables. The dependence structure is controlled by a copula function. Our new C-MGARCH model nests a conventional MGARCH model as a special case. The aim of this paper is to model MGARCH for non-normal multivariate distributions using copulas. We model the conditional correlation (by MGARCH) and the remaining dependence (by a copula) separately and simultaneously. We apply this idea to three MGARCH models, namely, the dynamic conditional correlation (DCC) model of Engle [Engle, R.F., 2002. Dynamic conditional correlation: A simple class of multivariate generalized autoregressive conditional heteroskedasticity models. Journal of Business and Economic Statistics 20, 339–350], the varying correlation (VC) model of Tse and Tsui [Tse, Y.K., Tsui, A.K., 2002. A multivariate generalized autoregressive conditional heteroscedasticity model with time-varying correlations. Journal of Business and Economic Statistics 20, 351–362], and the BEKK model of Engle and Kroner [Engle, R.F., Kroner, K.F., 1995. Multivariate simultaneous generalized ARCH. Econometric Theory 11, 122–150]. Empirical analysis with three foreign exchange rates indicates that the C-MGARCH models outperform DCC, VC, and BEKK in terms of in-sample model selection and out-of-sample multivariate density forecast, and in terms of these criteria the choice of copula functions is more important than the choice of the volatility models.     

Estimation of copula-based semiparametric time series models

This paper studies the estimation of a class of copula-based semiparametric stationary Markov models. These models are characterized by nonparametric marginal distributions and parametric copula functions, while the copulas capture all the scale-free temporal dependence of the processes. Simple estimators of the marginal distribution and the copula parameter are provided, and their asymptotic properties are established under easily verifiable conditions. These results are used to obtain root-n consistent and asymptotically normal estimators of important features of the transition distribution such as the (nonlinear) conditional moments and conditional quantiles. The semiparametric conditional quantile estimators are automatically monotonic across quantiles, which is attractive for portfolio conditional value-at-risk calculations.     

Non-Gaussian models for CoVaR estimation

In this paper we show how to obtain estimates of CoVaR based on models that take into consideration some stylized facts about multivariate financial time series of equity log returns: heavy tails, negative skew, asymmetric dependence, and volatility clustering. While the volatility clustering effect is captured by AR-GARCH dynamics of the Glosten-Jagannathan-Runkle (GJR) type, the other stylized facts are explained by non-Gaussian multivariate models and copula functions. We compare the different models in the period from January 2007 to March 2020. Our empirical study conducted on a sample of listed banks in the euro area confirms that, in measuring CoVaR, it is important to capture the time-varying dynamics of the volatility. Additionally, a correct assessment of the heaviness of the tails and of the dependence structure is needed in the evaluation of this systemic risk measure.     

Improving forecasting performance using covariate-dependent copula models

Copulas provide an attractive approach to the construction of multivariate distributions with flexible marginal distributions and different forms of dependences. Of particular importance in many areas is the possibility of forecasting the tail-dependences explicitly. Most of the available approaches are only able to estimate tail-dependences and correlations via nuisance parameters, and cannot be used for either interpretation or forecasting. We propose a general Bayesian approach for modeling and forecasting tail-dependences and correlations as explicit functions of covariates, with the aim of improving the copula forecasting performance. The proposed covariate-dependent copula model also allows for Bayesian variable selection from among the covariates of the marginal models, as well as the copula density. The copulas that we study include the Joe-Clayton copula, the Clayton copula, the Gumbel copula and the Student’s $t$-copula. Posterior inference is carried out using an efficient MCMC simulation method. Our approach is applied to both simulated data and the S&P 100 and S&P 600 stock indices. The forecasting performance of the proposed approach is compared with those of other modeling strategies based on log predictive scores. A value-at-risk evaluation is also performed for the model comparisons.     

Modelling soccer matches using bivariate discrete distributions with general dependence structure

In this paper copulas are used to generate bivariate discrete distributions. These distributions are fitted to soccer data from the English Premier League. An interesting aspect of these data is that the primary variable of interest, the discrete pair shots-for and shots-against, exhibits negative dependence; thus, in particular, we apply bivariate Poisson-related distributions that allow such dependence. The paper focuses on Archimedian copulas, for which the dependence structure is fully determined by a one-dimensional projection that is invariant under marginal transformations. Diagnostic plots for copula fit based on this projection are adapted to deal with discrete variables. Covariates relating to within-match contributions such as numbers of passes and tackles are introduced to explain variability in shot outcomes. The results of this analysis would appear to support the notion that playing the 'beautiful game' is an effective strategy—more passes and crosses contribute to more effective play and more shots on the goal.     

Modeling non-normal corporate bond yield spreads by copula

This research focuses on modeling for how corporate bond yield spreads are affected by explanatory variables such as equity volatility, interest rate volatility, r, slope, rating, liquidity, coupon rate, and maturity. The existing literature assumes normality and linearity in the analysis, which is not the case in our sample. Thus, through a powerful and flexible copula approach, we study the dependence at the mean of the joint distribution by using the Gaussian copula marginal regression method and the dependence structure at the tails by using various copula functions. To our knowledge, this is the first application of the copula marginal regression model to bond market data. In addition, we employ several copula functions to test for the tail dependence between yield spreads and other explanatory variables. We find stronger tail dependence in the joint upper tail for the relation between equity volatility and yield spreads, among others. This result indicates the positive effect of equity volatility on yield spreads in the upper tail is greater than that in the low tail. This finding should be useful to practitioners, such as investors. By relying on better-fitting, more meaningful statistical models, this paper contributes to the extant literature on how corporate bond yield spreads are determined.     

The t Copula and Related Copulas

The t copula and its properties are described with a focus on issues related to the dependence of extreme values. The Gaussian mixture representation of a multivariate t distribution is used as a starting point to construct two new copulas, the skewed t copula and the grouped t copula, which allow more heterogeneity in the modelling of dependent observations. Extreme value considerations are used to derive two further new copulas: the t extreme value copula is the limiting copula of componentwise maxima of t distributed random vectors; the t lower tail copula is the limiting copula of bivariate observations from a t distribution that are conditioned to lie below some joint threshold that is progressively lowered. Both these copulas may be approximated for practical purposes by simpler, better-known copulas, these being the Gumbel and Clayton copulas respectively.     

Estimation and model selection of semiparametric multivariate survival functions under general censorship

We study estimation and model selection of semiparametric models of multivariate survival functions for censored data, which are characterized by possibly misspecified parametric copulas and nonparametric marginal survivals. We obtain the consistency and root-n$n$ asymptotic normality of a two-step copula estimator to the pseudo-true copula parameter value according to KLIC, and provide a simple consistent estimator of its asymptotic variance, allowing for a first-step nonparametric estimation of the marginal survivals. We establish the asymptotic distribution of the penalized pseudo-likelihood ratio statistic for comparing multiple semiparametric multivariate survival functions subject to copula misspecification and general censorship. An empirical application is provided.     

Financial dependence analysis: applications of vine copulas

This paper features the application of a novel and recently developed method of statistical and mathematical analysis to the assessment of financial risk, namely regular vine copulas. Dependence modelling using copulas is a popular tool in financial applications but is usually applied to pairs of securities. Vine copulas offer greater flexibility and permit the modelling of complex dependence patterns using the rich variety of bivariate copulas that can be arranged and analysed in a tree structure to facilitate the analysis of multiple dependencies. We apply regular vine copula analysis to a sample of stocks comprising the Dow Jones index to assess their interdependencies and to assess how their correlations change in different economic circumstances using three different sample periods around Global Financial Crisis (GFC).: pre‐GFC (January 2005 to July 2007), GFC (July 2007 to September 2009) and post‐GFC periods (September 2009 to December 2011). The empirical results suggest that the dependencies change in a complex manner, and there is evidence of greater reliance on the Student‐t copula in the copula choice within the tree structures for the GFC period, which is consistent with the existence of larger tails in the distributions of returns for this period. One of the attractions of this approach to risk modelling is the flexibility in the choice of distributions used to model co‐dependencies. The practical application of regular vine metrics is demonstrated via an example of the calculation of the Value at Risk of a portfolio of stocks.     

Comparing probabilistic forecasts of the daily minimum and maximum temperature

Understanding changes in the frequency, severity, and seasonality of daily temperature extremes is important for public policy decisions regarding heat waves and cold snaps. A heat wave is sometimes defined in terms of both the daily minimum and maximum temperature, which necessitates the generation of forecasts of their joint distribution. In this paper, we develop time series models with the aim of providing insight and producing forecasts of the joint distribution that can challenge the accuracy of forecasts based on ensemble predictions from a numerical weather prediction model. We use ensemble model output statistics to recalibrate the raw ensemble predictions for the marginal distributions, with ensemble copula coupling used to capture the dependency between the marginal distributions. In terms of time series modelling, we consider a bivariate VARMA-MGARCH model. We use daily Spanish data recorded over a 65-year period, and find that, for the 5-year out-of-sample period, the recalibrated ensemble predictions outperform the time series models in terms of forecast accuracy.     

Estimation of multivariate models for time series of possibly different lengths

We consider the problem of estimating parametric multivariate density models when unequal amounts of data are available on each variable. We focus in particular on the case that the unknown parameter vector may be partitioned into elements relating only to a marginal distribution and elements relating to the copula. In such a case we propose using a multi‐stage maximum likelihood estimator (MSMLE) based on all available data rather than the usual one‐stage maximum likelihood estimator (1SMLE) based only on the overlapping data. We provide conditions under which the MSMLE is not less asymptotically efficient than the 1SMLE, and we examine the small sample efficiency of the estimators via simulations. The analysis in this paper is motivated by a model of the joint distribution of daily Japanese yen–US dollar and euro–US dollar exchange rates. We find significant evidence of time variation in the conditional copula of these exchange rates, and evidence of greater dependence during extreme events than under the normal distribution. Copyright © 2006 John Wiley & Sons, Ltd.     

A Unifying Framework for Marginalised Random‐Intercept Models of Correlated Binary Outcomes

We demonstrate that many current approaches for marginal modelling of correlated binary outcomes produce likelihoods that are equivalent to the copula‐based models herein. These general copula models of underlying latent threshold random variables yield likelihood‐based models for marginal fixed effects estimation and interpretation in the analysis of correlated binary data with exchangeable correlation structures. Moreover, we propose a nomenclature and set of model relationships that substantially elucidates the complex area of marginalised random‐intercept models for binary data. A diverse collection of didactic mathematical and numerical examples are given to illustrate concepts. Copyright © 2013 John Wiley & Sons, Ltd.     

Default models based on scale mixtures of Marshall-Olkin copulas: properties and applications

We present a unification of the Archimedean and the Lévy-frailty copula model for portfolio default models. The new default model exhibits a copula known as scale mixture of Marshall-Olkin copulas and an investigation of the dependence structure reveals that desirable properties of both original models are combined. This allows for a wider range of dependence patterns, while the analytical tractability is retained. Furthermore, simultaneous defaults and default clustering are incorporated. In addition, a hierarchical extension is presented which allows for a heterogeneous dependence structure. Finally, the model is applied to the pricing of CDO contracts. For this purpose, an efficient Laplace transform inversion approach is developed. Supporting a separation of marginal default probabilities and dependence structure, the model can be calibrated to CDS contracts in a first step. In a second step, the calibration of several parametric families to CDO contracts demonstrates a good fitting quality, which further emphasizes the suitability of the approach.     

Nonparametric inference on multivariate versions of Blomqvist’s beta and related measures of tail dependence

We consider nonparametric estimation of multivariate versions of Blomqvist’s beta, also known as the medial correlation coefficient. For a two-dimensional population, the sample version of Blomqvist’s beta describes the proportion of data which fall into the first or third quadrant of a two-way contingency table with cutting points being the sample medians. Asymptotic normality and strong consistency of the estimators are established by means of the empirical copula process, imposing weak conditions on the copula. Though the asymptotic variance takes a complicated form, we are able to derive explicit formulas for large families of copulas. For the copulas of elliptically contoured distributions we obtain a variance stabilizing transformation which is similar to Fisher’s z-transformation. This allows for an explicit construction of asymptotic confidence bands used for hypothesis testing and eases the analysis of asymptotic efficiency. The computational complexity of estimating Blomqvist’s beta corresponds to the sample size n, which is lower than the complexity of most competing dependence measures.      

Interpoint Distance Test of Homogeneity for Multivariate Mixture Models

Finite mixtures offer a rich class of distributions for modelling of a variety of random phenomena in numerous fields of study. Using the sample interpoint distances (IPDs), we propose the IPD‐test statistic for testing the hypothesis of homogeneity of mixture of multivariate power series distribution or multivariate normal distribution. We derive the distribution of the IPDs that are drawn from a finite mixture of the multivariate power series distribution and multivariate normal distribution. Based on the empirical distribution of the IPDs, we construct a bootstrap test of homogeneity for other multivariate finite mixture models. The IPD test is applied to mixture models for matrix‐valued distributions and a test of homogeneity for Wishart mixture is presented. Numerical comparisons show that IPD test has accurate type I errors and is more powerful in most multivariate cases than the expectation–maximization (EM) test and modified likelihood ratio test.     

Volatility linkages in the spot and futures market in Australia: a copula approach

To better characterize the dependence structure of the joint returns distribution, we propose to blend copula functions with Asymmetric GARCH (AGARCH) models, which are allowed for generalized error distribution. We model the copula’s marginals by the AGARCH processes that can differentiate between the impacts of positive and negative shocks on the returns volatility while taking the large kurtosis of the returns into account. An application of the procedure is elaborated on the All Ordinaries Index and its corresponding Share Price Index on future contracts in Australia. The findings reveal that the two spot and future markets show a strong right tail dependence but no left tail dependence. This provides a very useful knowledge for the risk management and hedging in futures markets.     

Nonparametric estimation of structural change points in volatility models for time series

We propose a hybrid estimation procedure that combines the least squares and nonparametric methods to estimate change points of volatility in time series models. Its main advantage is that it does not require any specific form of marginal or transitional densities of the process. We also establish the asymptotic properties of the estimators when the regression and conditional volatility functions are not known. The proposed tests for change points of volatility are shown to be consistent and more powerful than the nonparametric ones in the literature. Finally, we provide simulations and empirical results using the Hong Kong stock market index (HSI) series.