Perlman and Wellner’s Circular and Transformed Circular Copulas are Particular Beta and t Copulas

: All but one of the copulas in a recent paper in Symmetry by Perlman and Wellner can be identiﬁed as particular members of either the beta or t families of elliptical copulas


Introduction
In an interesting paper in this journal, Perlman and Wellner ( [1]; henceforth PW) explored the copulas arising from uniform distributions on the unit ball in R d and, when d = 2, the copulas arising from certain transformations thereof. Copulas are probability distributions whose marginal distributions are all uniformly distributed; they have a major role to play in multivariate statistical analysis. Amongst the more prominent examples of copulas in the statistical literature are "elliptical copulas" which are those based on marginal transformations to uniformity of distributions whose densities have elliptically symmetric contours; these, of course, include distributions whose densities have spherically symmetric contours as special cases. Prominent amongst multivariate elliptical/spherical distributions/copulas are multivariate t distributions/copulas and multivariate symmetric beta distributions/copulas; in elliptical distribution form, the latter were introduced as multivariate Pearson Type II distributions ( [2,3]), and a scaled version of the former are called multivariate Pearson Type VII distributions. See Fang et al. [4].
It turns out that all the copulas in PW can be identified as members of one or other of these families except for their "spherical copula" (PW, Section 4). In two cases (Section 2 below), this is just a not-especially-helpful renaming of the copulas. In the third (Section 3 below), the link is perhaps surprising and worthy of some explication.

Immediate Links with Beta Copulas
For clarity, suppose that d = 2. The bivariate spherically symmetric beta distribution has the simple density function for parameter β > 0. Its marginal densities are both symmetric beta distributions on [−1, 1] with parameter b = β + 1/2 where, for example, Clearly, these marginals are uniform if β = 1/2 in which case the bivariate density is which is density (6) of PW. That is, PW's "circular copula" is the bivariate spherically symmetric beta (or Pearson Type II) copula with parameter β = 1/2 in the above parametrisation. (This identification says nothing about its uniqueness, as proved by PW, nor, like the other identifications I make in this note, does it help with obtaining its distribution function-often just called its "copula"-which is a major focus of PW's work.) Now introduce the elliptically symmetric beta distribution associated with the above spherically symmetric beta distribution. In the usual way (e.g., [4]), it arises by linearly transforming the random variables that follow the spherically symmetric beta distribution by premultiplying their vector (X, Y ) T by a square root of the matrix where −1 < ρ < 1 turns out to be the correlation coefficient. The resulting density, in a standard parametrisation, is When β = 1/2 and ρ = sin γ, this can be seen to equate to (32) of PW. That is, PW's "one-parameter family of elliptical copulas" is the bivariate elliptically symmetric beta (or Pearson Type II) copula with parameter β = 1/2 in the above parametrisation. Algebraically, the non-standard derivation of this distribution by PW (Section 5) arises from explicitly working directly with the Cholesky square root of Σ, essentially their (29).

Less Immediate Links with t Copulas
In their Section 6, PW observe that if {X, Y } are the random variables distributed according to the "circular copula", then the non-linearly transformed variables {U, V } where are marginally uniformly distributed and hence jointly distributed according to a copula. In their Proposition 6.1, PW show this copula to have density I now observe that this is the bivariate spherically symmetric t copula with degrees of freedom 2. (I have not seen this simple formula for, in short, the t 2 copula density published anywhere before.) To verify this claim, start from the density of the bivariate spherically symmetric t distribution, again simple and given by f t,s (x, y; ν) = 1 2π Here, the parameter ν > 0, for reasons concerned with sampling derivations of such distributions, is known as the degrees of freedom (henceforth d.f.). As is well known, its marginal densities are both ("Student") t distributions on ν d.f. The t ν copula density on [0, 1] × [0, 1] can then be derived from f t,s via the formula (e.g., Nelsen [5]) , 0 ≤ w, z ≤ 1.
Here, f ν and F ν are the density and distribution functions of the univariate t ν distribution. A further linear transformation from [0, 1] × [0, 1] to [−1, 1] × [−1, 1], as in PW, then sees the t ν copula density become In general, F ν involves the incomplete beta function but, as is well known in the case of ν = 1 (the Cauchy case) and under-appreciated in the case of ν = 2 (Jones [6]), in each of those two cases it reduces to a simple invertible formula. In particular, .
Using these formulae and that of f t,s (x, y; 2) in c † t (w, z; 2) results in c t (u, v; 2) = c P W (u, v) by straightforward algebraic manipulations.
A more interesting route to the same result takes advantage of general relationships between univariate and bivariate spherically symmetric beta and t distributions. First, let B be a random variable following the symmetric beta distribution on [−1, 1] with parameter b = ν/2 and let T be a random variable following the t distribution with parameter ν. The following relationship is a version of one that is long established (Cacoullos [7,8]) but still relatively little known among statisticians: Second, let {B 1 , B 2 } be random variables following the bivariate spherically symmetric beta distribution with density f b,s (x, y; ν/2), i.e., with β = ν/2, and let {T 1 , T 2 } be random variables following the bivariate spherically symmetric t distribution with density f t,s (x, y; ν). Then, the relationship between {B 1 , B 2 } and {T 1 , T 2 }, which is also known and is a rather easy consequence of the polar representations of these bivariate spherically symmetric distributions, is that Recalling PW's notation that {B 1 , B 2 } = {X, Y }, it follows that the random variables follow the bivariate spherically symmetric t 2 distribution. But marginally transforming {T 1 , T 2 } back to beta marginals (in this case, uniform marginals since b = 1) via B i = T i / 2 + T 2 i , i = 1, 2, yields a distribution with the same copula as {T 1 , T 2 }-the t 2 copula. And As PW note, their construction of {U, V } "extends readily to generate a [d-dimensional] copula". I will now outline briefly how the above argument generalises to show that this copula is the d-dimensional spherically symmetric t copula with degrees of freedom 2. Start from {X 1 , ..., X d } being uniformly distributed on the unit ball, which still corresponds to a spherically symmetric beta distribution with β = 1. PW's elegant conditional argument (p. 594) carries through to the random variables following the uniform distribution marginally and hence a copula jointly. The d-dimensional relationship between spherically symmetric beta and t distributions can be readily shown to be Then, the argument based on the same marginal transformations back to the beta (uniform) hold as for d = 2. In d dimensions, the t 2 copula density on [−1, 1] d actually turns out to be c t (u 1 , ..., u d ) = Γ((d/2) + 1) (It is not clear how much, if anything, this identification would have helped PW with their task of obtaining its distribution function.)

Conclusions
Most of the copulas in PW ( [1]) have been identified as either beta or t copulas. The simple explicit form of what has been shown in Section 3 to be the spherically symmetric t 2 copula density is an interesting byproduct of PW's work.