Representing Sparse Gaussian DAGs as Sparse R-Vines Allowing for Non-Gaussian Dependence

Modeling dependence in high-dimensional systems has become an increasingly important topic. Most approaches rely on the assumption of a multivariate Gaussian distribution such as statistical models on directed acyclic graphs (DAGs). They are based on modeling conditional independencies and are scalable to high dimensions. In contrast, vine copula models accommodate more elaborate features like tail dependence and asymmetry, as well as independent modeling of the marginals. This flexibility comes however at the cost of exponentially increasing complexity for model selection and estimation. We show a novel connection between DAGs with limited number of parents and truncated vine copulas under sufficient conditions. This motivates a more general procedure exploiting the fast model selection and estimation of sparse DAGs while allowing for non-Gaussian dependence using vine copulas. By numerical examples in hundreds of dimensions, we demonstrate that our approach outperforms the standard method for vine structure selection. Supplementary material for this article is available online.     

Dependence Properties of Conditional Distributions of some Copula Models

For multivariate data from an observational study, inferences of interest can include conditional probabilities or quantiles for one variable given other variables. For statistical modeling, one could fit a parametric multivariate model, such as a vine copula, to the data and then use the model-based conditional distributions for further inference. Some results are derived for properties of conditional distributions under different positive dependence assumptions for some copula-based models. The multivariate version of the stochastically increasing ordering of conditional distributions is introduced for this purpose. Results are explained in the context of multivariate Gaussian distributions, as properties for Gaussian distributions can help to understand the properties of copula extensions based on vines.     

A practical non-parametric copula algorithm for system reliability with correlations

System reliability analysis involving correlated random variables is challenging because the failure probability cannot be uniquely determined under the given probability information. This paper proposes a system reliability evaluation method based on non-parametric copulas. The approximated joint probability distribution satisfying the constraints specified by correlations has the maximal relative entropy with respect to the joint probability distribution of independent random variables. Thus the reliability evaluation is unbiased from the perspective of information theory. The estimation of the non-parametric copula parameters from Pearson linear correlation, Spearman rank correlation, and Kendall rank correlation are provided, respectively. The approximated maximum entropy distribution is then integrated with the first and second order system reliability method. Four examples are adopted to illustrate the accuracy and efficiency of the proposed method. It is found that traditional system reliability method encodes excessive dependence information for correlated random variables and the estimated failure probability can be significantly biased.     

Modeling multi-country mortality dependence and its application in pricing survivor index swaps—A dynamic copula approach

This paper introduces mortality dependence in multi-country mortality modeling using a dynamic copula approach. Specifically, we use time-varying copula models to capture the mortality dependence structure across countries, examining both symmetric and asymmetric dependence structures. In addition, to capture the phenomenon of a heavy tail for the multi-country mortality index, we consider not only the setting of Gaussian innovations but also non-Gaussian innovations under the Lee–Carter framework model. As tests of the goodness of fit of different dynamic copula models, the pattern of mortality dependence, and the distribution of the innovations, we used empirical mortality data from Finland, France, the Netherlands, and Sweden. To understand the effect of mortality dependence on longevity derivatives, we also built a valuation framework for pricing a survivor index swap, then investigated the fair swap rates of a survivor swap numerically. We demonstrate that failing to consider the dynamic copula mortality model and non-Gaussian innovations would lead to serious underestimations of the swap rates and loss reserves.     

Modeling dependence structure among European markets and among Asian-Pacific markets: a regime switching regular vine copula approach

This paper investigates the structure of dependence among twelve European markets and among twelve Asian-Pacific markets. The dynamic of the dependence structure is described by a two-state regime switching model. The dependence structure during a bull phase is modelled by the Gaussian copula, while dependence during a bear phase is modelled by the regular vine copula. We analyze the regular vine structure in the second regime precisely. We perform a simplification procedure using a likelihood-ratio test and discuss the substitution of general regular vines by canonical vines or drawable vines. The analysis confirms the two-state nature of financial markets in addition to asymmetric and heavy-tailed dependences. Additionally, the European market has proven to be more strongly connected than the Asian-Pacific market, and European dependences are deeper in terms of conditional dependences. The results can be used by international investors by taking into account differences of both analyzed regions. Additionally, the analysis may help with the crisis prediction. The shift time to the market phase describing crisis times occurs significantly before the crisis itself.     

基于“藤”结构的高维动态Copula的构建

高维化和动态化是当前Copula理论研究和应用的两个重要方向.采用图形建模工具中"藤"的层叠结构,以二元动态Copula取代原有二元静态Copula作为"藤"的节点,将高维Copula建模中"藤"的方法与动态Copula相结合,构造了"动态藤Copula".实证表明,高维动态藤Copula较相应的高维静态藤Copula对数据的概率模型的似然率更高.     

Modeling dependence based on mixture copulas and its application in risk management

This paper is concerned with the statistical modeling of the dependence structure of multivariate financial data using the copula, and the application of copula functions in VaR valuation. After the introduction of the pure copula method and the maximum and minimum mixture copula method, authors present a new algorithm based on the more generalized mixture copula functions and the dependence measure, and apply the method to the portfolio of Shanghai stock composite index and Shenzhen stock component index. Comparing with the results from various methods, one can find that the mixture copula method is better than the pure Gaussian copula method and the maximum and minimum mixture copula method on different VaR level.     

Robust pair‐copula based forecasts of realized volatility

A useful application for copula functions is modeling the dynamics in the conditional moments of a time series. Using copulas, one can go beyond the traditional linear ARMA (p,q) modeling, which is solely based on the behavior of the autocorrelation function, and capture the entire dependence structure linking consecutive observations. This type of serial dependence is best represented by a canonical vine decomposition, and we illustrate this idea in the context of emerging stock markets, modeling linear and nonlinear temporal dependences of Brazilian series of realized volatilities. However, the analysis of intraday data collected from e‐markets poses some specific challenges. The large amount of real‐time information calls for heavy data manipulation, which may result in gross errors. Atypical points in high‐frequency intraday transaction prices may contaminate the series of daily realized volatilities, thus affecting classical statistical inference and leading to poor predictions. Therefore, in this paper, we propose to robustly estimate pair‐copula models using the weighted minimum distance and the weighted maximum likelihood estimates (WMLE). The excellent performance of these robust estimates for pair‐copula models are assessed through a comprehensive set of simulations, from which the WMLE emerged as the best option for members of the elliptical copula family. We evaluate and compare alternative volatility forecasts and show that the robustly estimated canonical vine‐based forecasts outperform the competitors. Copyright © 2013 John Wiley & Sons, Ltd.     

A copula entropy approach to correlation measurement at the country level

The entropy optimization approach has widely been applied in finance for a long time, notably in the areas of market simulation, risk measurement, and financial asset pricing. In this paper, we propose copula entropy models with two and three variables to measure dependence in stock markets, which extend the copula theory and are based on Jaynes’s information criterion. Both of them are usually applied under the non-Gaussian distribution assumption. Comparing with the linear correlation coefficient and the mutual information, the strengths and advantages of the copula entropy approach are revealed and confirmed. We also propose an algorithm for the copula entropy approach to obtain the numerical results. With the experimental data analysis at the country level and the economic circle theory in international economy, the validity of the proposed approach is approved; evidently, it captures the non-linear correlation, multi-dimensional correlation, and correlation comparisons without common variables. We would like to make it clear that correlation illustrates dependence, but dependence is not synonymous with correlation. Copulas can capture some special types of dependence, such as tail dependence and asymmetric dependence, which other conventional probability distributions, such as the normal p.d.f. and the Student’s t p.d.f., cannot.     

Copula selection for graphical models in continuous Estimation of Distribution Algorithms

This paper presents the use of graphical models and copula functions in Estimation of Distribution Algorithms (EDAs) for solving multivariate optimization problems. It is shown in this work how the incorporation of copula functions and graphical models for modeling the dependencies among variables provides some theoretical advantages over traditional EDAs. By means of copula functions and two well known graphical models, this paper presents a novel approach for defining new EDAs. Either dependence is modeled by a copula function chosen from a predefined set of six functions that aim to cover a wide range of inter-relations. It is also shown how the use of mutual information in the learning of graphical models implies a natural way of employing copula entropies. The experimental results on separable and non-separable functions show that the two new EDAs, which adopt copula functions to model dependencies, perform better than their original version with Gaussian variables.     

Tail order and intermediate tail dependence of multivariate copulas

In order to study copula families that have tail patterns and tail asymmetry different from multivariate Gaussian and t copulas, we introduce the concepts of tail order and tail order functions. These provide an integrated way to study both tail dependence and intermediate tail dependence. Some fundamental properties of tail order and tail order functions are obtained. For the multivariate Archimedean copula, we relate the tail heaviness of a positive random variable to the tail behavior of the Archimedean copula constructed from the Laplace transform of the random variable, and extend the results of Charpentier and Segers [7] [A. Charpentier, J. Segers, Tails of multivariate Archimedean copulas, Journal of Multivariate Analysis 100 (7) (2009) 1521–1537] for upper tails of Archimedean copulas. In addition, a new one-parameter Archimedean copula family based on the Laplace transform of the inverse Gamma distribution is proposed; it possesses patterns of upper and lower tails not seen in commonly used copula families. Finally, tail orders are studied for copulas constructed from mixtures of max-infinitely divisible copulas.     

Tail negative dependence and its applications for aggregate loss modeling

Tail order of copulas can be used to describe the strength of dependence in the tails of a joint distribution. When the value of tail order is larger than the dimension, it may lead to tail negative dependence. First, we prove results on conditions that lead to tail negative dependence for Archimedean copulas. Using the conditions, we construct new parametric copula families that possess upper tail negative dependence. Among them, a copula based on a scale mixture with a generalized gamma random variable (GGS copula) is useful for modeling asymmetric tail negative dependence. We propose mixed copula regression based on the GGS copula for aggregate loss modeling of a medical expenditure panel survey dataset. For this dataset, we find that there exists upper tail negative dependence between loss frequency and loss severity, and the introduction of tail negative dependence structures significantly improves the aggregate loss modeling.     

Second order tail asymptotics for the sum of dependent, tail-independent regularly varying risks

In this paper we consider dependent random variables with common regularly varying marginal distribution. Under the assumption that these random variables are tail-independent, it is well known that the tail of the sum behaves like in the independence case. Under some conditions on the marginal distributions and the dependence structure (including Gaussian copula’s and certain Archimedean copulas) we provide the second-order asymptotic behavior of the tail of the sum.     

Tail Dependence from a Distributional Point of View

The dependence structure in the tails of bivariate random variables is studied by means of appropriate copulae. Weak convergence results show that these copulae are natural dependence structures for joint tail events. The results obtained apply to particular types of copulae such as archimedean copulae and the Gaussian copula. Further, connections to multivariate extreme value theory are investigated and a two-dimensional Pickands–Balkema–de Haan Theorem type is derived. Finally, a counterexample showing that the tail dependence coefficients do not completely determine the dependence structure of bivariate rare events is provided.     

Monte Carlo evaluation of biological variation: Random generation of correlated non-Gaussian model parameters

When modelling the behaviour of horticultural products, demonstrating large sources of biological variation, we often run into the issue of non-Gaussian distributed model parameters. This work presents an algorithm to reproduce such correlated non-Gaussian model parameters for use with Monte Carlo simulations. The algorithm works around the problem of non-Gaussian distributions by transforming the observed non-Gaussian probability distributions using a proposed SKN-distribution function before applying the covariance decomposition algorithm to generate Gaussian random co-varying parameter sets. The proposed SKN-distribution function is based on the standard Gaussian distribution function and can exhibit different degrees of both skewness and kurtosis. This technique is demonstrated using a case study on modelling the ripening of tomato fruit evaluating the propagation of biological variation with time.     

On approximating max-stable processes and constructing extremal copula functions

There is an infinite number of parameters in the definition of multivariate maxima of moving maxima (M4) processes, which poses challenges in statistical applications where workable models are preferred. This paper establishes sufficient conditions under which an M4 process with infinite number of parameters may be approximated by an M4 process with finite number of parameters. In statistical inferences, the paper focuses on a family of sectional multivariate extreme value copula (SMEVC) functions which is derived from the joint distribution functions of M4 processes. A new non-standard parameter estimation procedure is introduced, which is based on order statistics of ratios of (transformed) marginal unit Fréchet random variables, and is shown via simulation to be more efficient than a semi-parametric estimation procedure. In real data analysis, empirical results show that SMEVCs are more flexible for modeling various dependence structures, and perform better than the widely used Gumbel-Hougaard copulas.     

A generalized beta copula with applications in modeling multivariate long-tailed data

This work proposes a new copula class that we call the MGB2 copula. The new copula originates from extracting the dependence function of the multivariate GB2 distribution (MGB2) whose marginals follow the univariate generalized beta distribution of the second kind (GB2). The MGB2 copula can capture non-elliptical and asymmetric dependencies among marginal coordinates and provides a simple formulation for multi-dimensional applications. This new class features positive tail dependence in the upper tail and tail independence in the lower tail. Furthermore, it includes some well-known copula classes, such as the Gaussian copula, as special or limiting cases.To illustrate the usefulness of the MGB2 copula, we build a trivariate MGB2 copula model of bodily injury liability closed claims. Extended GB2 distributions are chosen to accommodate the right-skewness and the long-tailedness of the outcome variables. For the regression component, location parameters with continuous predictors are introduced using a nonlinear additive function. For comparison purposes, we also consider the Gumbel and t copulas, alternatives that capture the upper tail dependence. The paper introduces a conditional plot graphical tool for assessing the validation of the MGB2 copula. Quantitative and graphical assessment of the goodness of fit demonstrate the advantages of the MGB2 copula over the other copulas.     

Robust portfolio optimization with copulas

Conditional Value at Risk (CVaR) is widely used in portfolio optimization as a measure of risk. CVaR is clearly dependent on the underlying probability distribution of the portfolio. We show how copulas can be introduced to any problem that involves distributions and how they can provide solutions for the modeling of the portfolio. We use this to provide the copula formulation of the CVaR of a portfolio. Given the critical dependence of CVaR on the underlying distribution, we use a robust framework to extend our approach to Worst Case CVaR (WCVaR). WCVaR is achieved through the use of rival copulas. These rival copulas have the advantage of exploiting a variety of dependence structures, symmetric and not. We compare our model against two other models, Gaussian CVaR and Worst Case Markowitz. Our empirical analysis shows that WCVaR can asses the risk more adequately than the two competitive models during periods of crisis.     

Dependence modeling in non-life insurance using the Bernstein copula

This paper illustrates the modeling of dependence structures of non-life insurance risks using the Bernstein copula. We conduct a goodness-of-fit analysis and compare the Bernstein copula with other widely used copulas. Then, we illustrate the use of the Bernstein copula in a value-at-risk and tail-value-at-risk simulation study. For both analyses we utilize German claims data on storm, flood, and water damage insurance for calibration. Our results highlight the advantages of the Bernstein copula, including its flexibility in mapping inhomogeneous dependence structures and its easy use in a simulation context due to its representation as mixture of independent Beta densities. Practitioners and regulators working toward appropriate modeling of dependences in a risk management and solvency context can benefit from our results.     

Tail dependence functions and vine copulas

Tail dependence and conditional tail dependence functions describe, respectively, the tail probabilities and conditional tail probabilities of a copula at various relative scales. The properties as well as the interplay of these two functions are established based upon their homogeneous structures. The extremal dependence of a copula, as described by its extreme value copulas, is shown to be completely determined by its tail dependence functions. For a vine copula built from a set of bivariate copulas, its tail dependence function can be expressed recursively by the tail dependence and conditional tail dependence functions of lower-dimensional margins. The effect of tail dependence of bivariate linking copulas on that of a vine copula is also investigated.