Construction of asymmetric multivariate copulas

In this paper we introduce two methods for the construction of asymmetric multivariate copulas. The first is connected with products of copulas. The second approach generalises the Archimedean copulas. The resulting copulas are asymmetric and may have more than two parameters in contrast to most of the parametric families of copulas described in the literature. We study the properties of the proposed families of copulas such as the dependence of two components (Kendall’s tau, tail dependence), marginal distributions and the generation of random variates.     

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.     

Constructing hierarchical Archimedean copulas with Lévy subordinators

A probabilistic interpretation for hierarchical Archimedean copulas based on Lévy subordinators is given. Independent exponential random variables are divided by group-specific Lévy subordinators which are evaluated at a common random time. The resulting random vector has a hierarchical Archimedean survival copula. This approach suggests an efficient sampling algorithm and allows one to easily construct several new parametric families of hierarchical Archimedean copulas.     

From Archimedean to Liouville copulas

We use a recent characterization of the d-dimensional Archimedean copulas as the survival copulas of d-dimensional simplex distributions (McNeil and Nešlehová (2009) [1]) to construct new Archimedean copula families, and to examine the relationship between their dependence properties and the radial parts of the corresponding simplex distributions. In particular, a new formula for Kendall’s tau is derived and a new dependence ordering for non-negative random variables is introduced which generalises the Laplace transform order. We then generalise the Archimedean copulas to obtain Liouville copulas, which are the survival copulas of Liouville distributions and which are non-exchangeable in general. We derive a formula for Kendall’s tau of Liouville copulas in terms of the radial parts of the corresponding Liouville distributions.     

Nonparametric rank-based tests of bivariate extreme-value dependence

A new class of tests of extreme-value dependence for bivariate copulas is proposed. It is based on the process comparing the empirical copula with a natural nonparametric rank-based estimator of the unknown copula under extreme-value dependence. A multiplier technique is used to compute approximate p-values for several candidate test statistics. Extensive Monte Carlo experiments were carried out to compare the resulting procedures with the tests of extreme-value dependence recently studied in Ben Ghorbal et al. (2009) [1] and Kojadinovic and Yan (2010) [19]. The finite-sample performance study of the tests is complemented by local power calculations.     

Characterization of dependence of multidimensional Lévy processes using Lévy copulas

This paper suggests Lévy copulas in order to characterize the dependence among components of multidimensional Lévy processes. This concept parallels the notion of a copula on the level of Lévy measures. As for random vectors, a version of Sklar's theorem states that the law of a general multivariate Lévy process is obtained by combining arbitrary univariate Lévy processes with an arbitrary Lévy copula. We construct parametric families of Lévy copulas and prove a limit theorem, which indicates how to obtain the Lévy copula of a multivariate Lévy process X from the ordinary copula of the random vector X_t for small t.     

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.     

On rank correlation measures for non-continuous random variables

For continuous random variables, many dependence concepts and measures of association can be expressed in terms of the corresponding copula only and are thus independent of the marginal distributions. This interrelationship generally fails as soon as there are discontinuities in the marginal distribution functions. In this paper, we consider an alternative transformation of an arbitrary random variable to a uniformly distributed one. Using this technique, the class of all possible copulas in the general case is investigated. In particular, we show that one of its members—the standard extension copula introduced by Schweizer and Sklar—captures the dependence structures in an analogous way the unique copula does in the continuous case. Furthermore, we consider measures of concordance between arbitrary random variables and obtain generalizations of Kendall's tau and Spearman's rho that correspond to the sample version of these quantities for empirical distributions.     

Empirical likelihood based confidence intervals for copulas

Copula as an effective way of modeling dependence has become more or less a standard tool in risk management, and a wide range of applications of copula models appear in the literature of economics, econometrics, insurance, finance, etc. How to estimate and test a copula plays an important role in practice, and both parametric and nonparametric methods have been studied in the literature. In this paper, we focus on interval estimation and propose an empirical likelihood based confidence interval for a copula. A simulation study and a real data analysis are conducted to compare the finite sample behavior of the proposed empirical likelihood method with the bootstrap method based on either the empirical copula estimator or the kernel smoothing copula estimator.     

On the distribution of the (un)bounded sum of random variables

We propose a general treatment of random variables aggregation accounting for the dependence among variables and bounded or unbounded support of their sum. The approach is based on the extension to the concept of convolution to dependent variables, involving copula functions. We show that some classes of copula functions (such as Marshall-Olkin and elliptical) cannot be used to represent the dependence structure of two variables whose sum is bounded, while Archimedean copulas can be applied only if the generator becomes linear beyond some point. As for the application, we study the problem of capital allocation between risks when the sum of losses is bounded.     

On the construction of copulas and quasi-copulas with given diagonal sections

We study a method, which we call a copula (or quasi-copula) diagonal splice, for creating new functions by joining portions of two copulas (or quasi-copulas) with a common diagonal section. The diagonal splice of two quasi-copulas is always a quasi-copula, and we find a necessary and sufficient condition for the diagonal splice of two copulas to be a copula. Applications of this method include the construction of absolutely continuous asymmetric copulas with a prescribed diagonal section, and determining the best-possible upper bound on the set of copulas with a particular type of diagonal section. Several examples illustrate our results.     

Estimating copula densities, using model selection techniques

Recently a new way of modeling dependence has been introduced considering a sequence of parametric copula models, covering more and more dependency aspects and thus giving a closer approximation to the true copula density. The method uses contamination families based on Legendre polynomials. It has been shown that in general after a few steps accurate approximations are obtained. In this paper selection of the adequate number of steps is considered, and estimation of the unknown parameters within the chosen contamination family is established, thus obtaining an estimator of the unknown copula density. There should be a balance between the complexity of the model and the number of parameters to be estimated. High complexity gives a low model error, but a large stochastic or estimation error, while a very simple model gives a small stochastic error, but a large model error. Techniques from model selection are applied, thus letting the data tell us which aspects are important enough to capture into the model. Natural and simple estimators of the involved Fourier coefficients complete the procedure. Theoretical results show that the expected quadratic error is reduced by the selection rule to the same order of magnitude as in a classical parametric problem. The method is applied on a real data set, illustrating that the new method describes the data set very well: the error involved in the classical Gaussian copula density is reduced with no fewer than 50%.     

Dependence structure of conditional Archimedean copulas

In this article, copulas associated to multivariate conditional distributions in an Archimedean model are characterized. It is shown that this popular class of dependence structures is closed under the operation of conditioning, but that the associated conditional copula has a different analytical form in general. It is also demonstrated that the extremal copula for conditional Archimedean distributions is no longer the Fréchet upper bound, but rather a member of the Clayton family. Properties of these conditional distributions as well as conditional versions of tail dependence indices are also considered.     

Parameter estimation of a bivariate compound Poisson process

In this article, we review the concept of a Lévy copula to describe the dependence structure of a bivariate compound Poisson process. In this first statistical approach we consider a parametric model for the Lévy copula and estimate the parameters of the full dependent model based on a maximum likelihood approach. This approach ensures that the estimated model remains in the class of multivariate compound Poisson processes. A simulation study investigates the small sample behaviour of the MLEs, where we also suggest a new simulation algorithm. Finally, we apply our method to Danish fire insurance data.     

A New Graphical Tool for Copula Selection

The selection of copulas is an important aspect of dependence modeling. In many practical applications, only a limited number of copulas is tested, and the modeling applications usually are restricted to the bivariate case. One explanation is the fact that no graphical copula tool exists that allows us to assess the goodness-of-fit of a large set of (possible higher-dimensional) copula functions at once. This article seeks to overcome this problem by developing a new graphical tool for the copula selection, based on a statistical analysis technique called “principal coordinate analysis.” The advantage is three-fold. First, when projecting the empirical copula of a modeling application on a two-dimensional (2D) copula space, it allows us to visualize the fit of a whole collection of multivariate copulas at once. Second, the visual tool allows us to identify “search” directions for potential fit improvements (e.g., through the use of copula transforms). Finally, the tool makes it also possible to give a 2D visual overview of a large number of known copula families, leading to a better understanding and a more efficient use of the different copula families. The robustness of the new graphical tool is investigated by means of a small simulation study, and the practical use of the tool is demonstrated for two 2D and two 3D (three-dimensional) fitting examples. MATLAB code through the examples is available online in the supplementary materials.     

On testing equality of pairwise rank correlations in a multivariate random vector

Spearman’s rank-correlation coefficient (also called Spearman’s rho) represents one of the best-known measures to quantify the degree of dependence between two random variables. As a copula-based dependence measure, it is invariant with respect to the distribution’s univariate marginal distribution functions. In this paper, we consider statistical tests for the hypothesis that all pairwise Spearman’s rank correlation coefficients in a multivariate random vector are equal. The tests are nonparametric and their asymptotic distributions are derived based on the asymptotic behavior of the empirical copula process. Only weak assumptions on the distribution function, such as continuity of the marginal distributions and continuous partial differentiability of the copula, are required for obtaining the results. A nonparametric bootstrap method is suggested for either estimating unknown parameters of the test statistics or for determining the associated critical values. We present a simulation study in order to investigate the power of the proposed tests. The results are compared to a classical parametric test for equal pairwise Pearson’s correlation coefficients in a multivariate random vector. The general setting also allows the derivation of a test for stochastic independence based on Spearman’s rho.     

Crisis and risk dependencies

The knowledge of the multivariate stochastic dependence between the returns of asset classes is of importance for many finance applications, such as asset allocation or risk management. By means of goodness-of-fit tests, we analyze for a multitude of portfolios consisting of different asset classes whether the stochastic dependence between the portfolios’ constituents can be adequately described by multivariate versions of some standard parametric copula functions. Furthermore, we test whether the stochastic dependence between the returns of different asset classes has changed during the recent financial crisis. The main findings are: First, whether a specific copula assumption can be rejected or not, crucially depends on the asset class and the time period considered. Second, different goodness-of-fit tests for copulas can yield very different results and these differences can vary for different asset classes and for different tested copulas. Third, even when using various goodness-of-fit tests for copulas, it is not always possible to differentiate between various copula assumptions. Fourth, during the financial crisis, copula assumptions are more frequently rejected. However, the results also raise some concerns over the suitability of goodness-of-fit tests for copulas as a diagnostic tool for identifying stressed risk dependencies.     

Estimation of nonstrict Archimedean copulas and its application to quantum networks

Bivariate nonstrict Archimedean copulas form a subclass of Archimedean copulas and are able to model the dependence structure of random variables that do not take on low quantiles simultaneously; i.e. their domain includes a set, the so‐called zero set, with positive Lebesgue measure but zero probability mass. Standard methods to fit a parametric Archimedean copula, e.g. classical maximum likelihood estimation, are either getting computationally more involved or even fail when dealing with this subclass. We propose an alternative method for estimating the parameter of a nonstrict Archimedean copula that is based on the zero set and the functional form of its boundary curve. This estimator is fast to compute and can be applied to absolutely continuous copulas but also allows singular components. In a simulation study, we compare its performance to that of the standard estimators. Finally, the estimator is applied when modeling the dependence structure of quantities describing the quality of transmission in a quantum network, and it is shown how this model can be used effectively to detect potential intruders in this network. Copyright © 2014 John Wiley & Sons, Ltd.     

Lévy-frailty copulas

A parametric family of n-dimensional extreme-value copulas of Marshall–Olkin type is introduced. Members of this class arise as survival copulas in Lévy-frailty models. The underlying probabilistic construction introduces dependence to initially independent exponential random variables by means of first-passage times of a Lévy subordinator. Jumps of the subordinator correspond to a singular component of the copula. Additionally, a characterization of completely monotone sequences via the introduced family of copulas is derived. An alternative characterization is given by Hausdorff’s moment problem in terms of random variables with compact support. The resulting correspondence between random variables, Lévy subordinators, and copulas is studied and illustrated with several examples. Finally, it is used to provide a general methodology for sampling the copula in many cases. The new class is shown to share some properties with Archimedean copulas regarding construction and analytical form. Finally, the parametric form allows us to compute different measures of dependence and the Pickands representation.     

Tails of multivariate Archimedean copulas

A complete and user-friendly directory of tails of Archimedean copulas is presented which can be used in the selection and construction of appropriate models with desired properties. The results are synthesized in the form of a decision tree: Given the values of some readily computable characteristics of the Archimedean generator, the upper and lower tails of the copula are classified into one of three classes each, one corresponding to asymptotic dependence and the other two to asymptotic independence. For a long list of single-parameter families, the relevant tail quantities are computed so that the corresponding classes in the decision tree can easily be determined. In addition, new models with tailor-made upper and lower tails can be constructed via a number of transformation methods. The frequently occurring category of asymptotic independence turns out to conceal a surprisingly rich variety of tail dependence structures.