Sampling nested Archimedean copulas

We give algorithms for sampling from non-exchangeable Archimedean copulas created by the nesting of Archimedean copula generators, where in the most general algorithm the generators may be nested to an arbitrary depth. These algorithms are based on mixture representations of these copulas using Laplace transforms. While in principle the approach applies to all nested Archimedean copulas, in practice the approach is restricted to certain cases where we are able to sample distributions with given Laplace transforms. Precise instructions are given for the case when all generators are taken from the Gumbel parametric family or the Clayton family; the Gumbel case in particular proves very easy to simulate.     

Estimation and inference of the joint conditional distribution for multivariate longitudinal data using nonparametric copulas

In this paper we study estimating the joint conditional distributions of multivariate longitudinal outcomes using regression models and copulas. For the estimation of marginal models, we consider a class of time-varying transformation models and combine the two marginal models using nonparametric empirical copulas. Our models and estimation method can be applied in many situations where the conditional mean-based models are not good enough. Empirical copulas combined with time-varying transformation models may allow quite flexible modelling for the joint conditional distributions for multivariate longitudinal data. We derive the asymptotic properties for the copula-based estimators of the joint conditional distribution functions. For illustration we apply our estimation method to an epidemiological study of childhood growth and blood pressure.     

Constructing copula functions with weighted geometric means

The weighted arithmetic mean of two copulas is a copula. In some cases, geometric and harmonic means also provide copulas. There are copulas specially appropriate to be combined by using weighted geometric means. With this method of construction we combine Farlie–Gumbel–Morgentern and Ali–Mikhail–Haq copulas to obtain families of copulas which can be expressed in terms of double power series. The Gumbel–Barnett copula is also considered and a new copula is proposed, which arises as the first order approximation of the weighted geometric mean of two copulas. Invariance of two multivariate distributions (Cuadras–Augé and Johnson–Kotz) by weighted geometric and arithmetic means is also studied.     

Copules archimédiennes et families de lois bidimensionnelles dont les marges sont données

Every bivariate distribution function with continuous marginals can be represented in terms of a unique copula, that is, in terms of a distribution function on the unit square with uniform marginals. This paper is concerned with a special class of copulas called Archimedean, which includes the uniform representation of many standard bivariate distributions. Conditions are given under which these copulas are stochastically ordered and pointwise limits of sequences of Archimedean copulas are examined. We also provide two new one-parameter families of bivariate distributions which include as limiting cases the Frechet bounds and the independence distribution.     

Examination and visualisation of the simplifying assumption for vine copulas in three dimensions

Vine copulas are a highly flexible class of dependence models, which are based on the decomposition of the density into bivariate building blocks. For applications one usually makes the simplifying assumption that copulas of conditional distributions are independent of the variables on which they are conditioned. However this assumption has been criticised for being too restrictive. We examine both simplified and non‐simplified vine copulas in three dimensions and investigate conceptual differences. We show and compare contour surfaces of three‐dimensional vine copula models, which prove to be much more informative than the contour lines of the bivariate marginals. Our investigation shows that non‐simplified vine copulas can exhibit arbitrarily irregular shapes, whereas simplified vine copulas appear to be smooth extrapolations of their bivariate margins to three dimensions. In addition to a variety of constructed examples, we also investigate a three‐dimensional subset of the well‐known uranium data set and visually detect the fact that a non‐simplified vine copula is necessary to capture its complex dependence structure.     

A bivariate extension of the beta generated distribution derived from copulas

In this paper, we introduce a new class of bivariate distributions whose marginals are beta-generated distributions. Copulas are employed to construct this bivariate extension of the beta-generated distributions. It is shown that when Archimedean copulas and convex beta generators are used in generating bivariate distributions, the copulas of the resulting distributions also belong to the Archimedean family. The dependence of the proposed bivariate distributions is examined. Simulation results for beta generators and an application to financial risk management are presented.     

CONDITIONAL EXPECTATION FORMULAE FOR COPULAS

Not only are copula functions joint distribution functions in their own right, they also provide a link between multivariate distributions and their lower‐dimensional marginal distributions. Copulas have a structure that allows us to characterize all possible multivariate distributions, and therefore they have the potential to be a very useful statistical tool. Although copulas can be traced back to 1959, there is still much scope for new results, as most of the early work was theoretical rather than practical. We focus on simple practical tools based on conditional expectation, because such tools are not widely available. When dealing with data sets in which the dependence throughout the sample is variable, we suggest that copula‐based regression curves may be more accurate predictors of specific outcomes than linear models. We derive simple conditional expectation formulae in terms of copulas and apply them to a combination of simulated and real data.     

Model-based clustering of Gaussian copulas for mixed data

Clustering of mixed data is important yet challenging due to a shortage of conventional distributions for such data. In this article, we propose a mixture model of Gaussian copulas for clustering mixed data. Indeed copulas, and Gaussian copulas in particular, are powerful tools for easily modeling the distribution of multivariate variables. This model clusters data sets with continuous, integer, and ordinal variables (all having a cumulative distribution function) by considering the intra-component dependencies in a similar way to the Gaussian mixture. Indeed, each component of the Gaussian copula mixture produces a correlation coefficient for each pair of variables and its univariate margins follow standard distributions (Gaussian, Poisson, and ordered multinomial) depending on the nature of the variable (continuous, integer, or ordinal). As an interesting by-product, this model generalizes many well-known approaches and provides tools for visualization based on its parameters. The Bayesian inference is achieved with a Metropolis-within-Gibbs sampler. The numerical experiments, on simulated and real data, illustrate the benefits of the proposed model: flexible and meaningful parameterization combined with visualization features.     

A flexible and tractable class of one-factor copulas

Copulas are a useful tool to model multivariate distributions. While there exist various families of bivariate copulas, the construction of flexible and yet tractable copulas suitable for high-dimensional applications is much more challenging. This is even more true if one is concerned with the analysis of extreme values. In this paper, we construct a class of one-factor copulas and a family of extreme-value copulas well suited for high-dimensional applications and exhibiting a good balance between tractability and flexibility. The inference for these copulas is performed by using a least-squares estimator based on dependence coefficients. The modeling capabilities of the copulas are illustrated on simulated and real datasets.     

Gluing Copulas

We present a new way of constructing n-copulas, by scaling and gluing finitely many n-copulas. Gluing for bivariate copulas produces a copula that coincides with the independence copula on some grid of horizontal and vertical sections. Examples illustrate how gluing can be applied to build complicated copulas from simple ones. Finally, we investigate the analytical as well as statistical properties of the copulas obtained by gluing, in particular, the behavior of Spearman's ρ and Kendall's τ.     

Bivariate Tail Dependence and the Generation of Multivariate Extreme Value Distributions

We define, in a probabilistic way, a parametric family of multivariate extreme value distributions. We derive its copula, which is a mixture of several complete dependent copulas and total independent copulas, and the bivariate tail dependence and extremal coefficients. Based on the obtained results for these coefficients, we propose a method to build multivariate extreme value distributions with prescribed tail/extremal coefficients. We illustrate the results with examples.     

Copulas with Truncation-Invariance Property

For a truncation-invariant copula, truncation does not change the dependence structure as well as all nonparametric measures of association such as Kendall's tau and Spearman's rho. In this article, we show that the products of algebraically independent Archimedean multivariate Clayton copulas and standard uniform distributions are the only truncation-invariant copulas.     

Bayesian bivariate survival analysis using the power variance function copula

Copula models have become increasingly popular for modelling the dependence structure in multivariate survival data. The two-parameter Archimedean family of Power Variance Function (PVF) copulas includes the Clayton, Positive Stable (Gumbel) and Inverse Gaussian copulas as special or limiting cases, thus offers a unified approach to fitting these important copulas. Two-stage frequentist procedures for estimating the marginal distributions and the PVF copula have been suggested by Andersen (Lifetime Data Anal 11:333–350, 2005), Massonnet et al. (J Stat Plann Inference 139(11):3865–3877, 2009) and Prenen et al. (J R Stat Soc Ser B 79(2):483–505, 2017) which first estimate the marginal distributions and conditional on these in a second step to estimate the PVF copula parameters. Here we explore an one-stage Bayesian approach that simultaneously estimates the marginal and the PVF copula parameters. For the marginal distributions, we consider both parametric as well as semiparametric models. We propose a new method to simulate uniform pairs with PVF dependence structure based on conditional sampling for copulas and on numerical approximation to solve a target equation. In a simulation study, small sample properties of the Bayesian estimators are explored. We illustrate the usefulness of the methodology using data on times to appendectomy for adult twins in the Australian NH&MRC Twin registry. Parameters of the marginal distributions and the PVF copula are simultaneously estimated in a parametric as well as a semiparametric approach where the marginal distributions are modelled using Weibull and piecewise exponential distributions, respectively.     

A stochastic representation and sampling algorithm for nested Archimedean copulas

A general sampling algorithm for nested Archimedean copulas was recently suggested. It is given in two different forms, a recursive or an explicit one. The explicit form allows for a simpler version of the algorithm which is numerically more stable and faster since less function evaluations are required. The algorithm can also be given in general form, not being restricted to a particular nesting such as fully nested Archimedean copulas. Further, several examples are given.     

Gaussian copula distributions for mixed data,with application in discrimination

The construction of a joint model for mixed discrete and continuous random variables that accounts for their associations is an important statistical problem in many practical applications. In this paper, we use copulas to construct a class of joint distributions of mixed discrete and continuous random variables. In particular, we employ the Gaussian copula to generate joint distributions for mixed variables. Examples include the robit-normal and probit-normal-exponential distributions, the first for modelling the distribution of mixed binary-continuous data and the second for a mixture of continuous, binary and trichotomous variables. The new class of joint distributions is general enough to include many mixed-data models currently available. We study properties of the distributions and outline likelihood estimation; a small simulation study is used to investigate the finite-sample properties of estimates obtained by full and pairwise likelihood methods. Finally, we present an application to discriminant analysis of multiple correlated binary and continuous data from a study involving advanced breast cancer patients.     

Duality of the Multivariate Distributions of Marshall–Olkin Type and Tail Dependence

A dual class of the multivariate distributions of Marshall–Olkin type is introduced, and their copulas are presented and utilized to derive explicit expressions of the distributional tail dependencies, which describe the amount of dependence in the upper-orthant tail or lower-orthant tail of a multivariate distribution and can be used in the study of dependence among extreme values. A sufficient condition under which tail dependencies of two such distributions can be compared are obtained. Some examples are also presented to illustrate our results.     

Mixing and moments properties of a non-stationary copula-based Markov process

Abstract

We provide conditions under which a non-stationary copula-based Markov process is geometric β-mixing and geometric ρ-mixing. Our results generalize some results of Beare who considers the stationary case. As a particular case we introduce a stochastic process, that we call convolution-based Markov process, whose construction is obtained by using the C-convolution operator which allows the increments to be dependent. Within this subclass of processes we characterize a modified version of the standard random walk where copulas and marginal distributions involved are in the same elliptical family. We study mixing and moments properties to identify the differences compared to the standard case.     

Estimating non-simplified vine copulas using penalized splines

Vine copulas (or pair-copula constructions) have become an important tool for high-dimensional dependence modeling. Typically, so-called simplified vine copula models are estimated where bivariate conditional copulas are approximated by bivariate unconditional copulas. We present the first nonparametric estimator of a non-simplified vine copula that allows for varying conditional copulas using penalized hierarchical B-splines. Throughout the vine copula, we test for the simplifying assumption in each edge, establishing a data-driven non-simplified vine copula estimator. To overcome the curse of dimensionality, we approximate conditional copulas with more than one conditioning argument by a conditional copula with the first principal component as conditioning argument. An extensive simulation study is conducted, showing a substantial improvement in the out-of-sample Kullback–Leibler divergence if the null hypothesis of a simplified vine copula can be rejected. We apply our method to the famous uranium data and present a classification of an eye state data set, demonstrating the potential benefit that can be achieved when conditional copulas are modeled.     

Various measures of dependence of a new asymmetric generalized Farlie–Gumbel–Morgenstern copulas

ABSTRACT

In this paper, we discuss an asymmetric extension of Farlie–Gumbel–Morgenstern copulas studied by several authors and obtain the range of the parameter. We derive an expression for regression function and the properties of these copulas are studied in detail. Also, explicit expressions for various measures of association are obtained. These measures are numerically compared for some particular parametric values of the copulas.