S S symmetry
Article Multivariate Skew t-Distribution: Asymptotics for Parameter Estimators and Extension to Skew t-Copula
Tõnu Kollo , Meelis Käärik and Anne Selart *
Institute of Mathematics and Statistics, University of Tartu, Narva mnt 18, 51009 Tartu, Estonia; [email protected] (T.K.); [email protected] (M.K.) * Correspondence: [email protected]
Abstract: Symmetric elliptical distributions have been intensively used in data modeling and ro- bustness studies. The area of applications was considerably widened after transforming elliptical distributions into the skew elliptical ones that preserve several good properties of the corresponding symmetric distributions and increase possibilities of data modeling. We consider three-parameter p-variate skew t-distribution where p-vector µ is the location parameter, Σ : p × p is the positive definite scale parameter, p-vector α is the skewness or shape parameter, and the number of degrees of freedom ν is fixed. Special attention is paid to the two-parameter distribution when µ = 0 that is useful for construction of the skew t-copula. Expressions of the parameters are presented through the moments and parameter estimates are found by the method of moments. Asymptotic normality is established for the estimators of Σ and α. Convergence to the asymptotic distributions is examined in simulation experiments. Keywords: asymptotic normality; inverse chi-distribution; multivariate cumulants; multivariate Citation: Kollo, T.; Käärik, M.; moments; skew normal distribution; skew t-copula; skew t-distribution Selart, A. Multivariate Skew t-Distribution: Asymptotics for MSC: 62E20; 62H10; 62H12 Parameter Estimators and Extension to Skew t-Copula. Symmetry 2021, 13, 1059. https://doi.org/10.3390/ sym13061059 1. Introduction The idea of constructing multivariate skewed distributions from symmetric elliptical Academic Editor: Nicola Maria distributions has attracted many statisticians over last 25 years. This generalization pre- Rinaldo Loperfido serves several good properties of corresponding symmetric distributions and in the same time increases possibilities of data modeling. In 1996, Azzalini and Dalla Valle introduced Received: 20 May 2021 multivariate skew normal distribution [1]. Their idea of transforming multivariate normal Accepted: 8 June 2021 Published: 11 June 2021 distribution into a skewed distribution has been carried over to several other distributions and many generalizations have been made since 1996 (see for instance, Theodossiou [2], t Publisher’s Note: MDPI stays neutral Arellano-Valle and Genton [3]). From these developments multivariate skew -distribution with regard to jurisdictional claims in has become most frequently used [4]. This distribution became an important tool in many published maps and institutional affil- data modeling applications because of heavier tail area than the normal and skew normal iations. distributions (see, for instance, Lee and McLachan [5], Marchenko [6], Parisi and Liseo [7], Fernandez and Steel [8], Jones [9]). A review of applications in finance and insurance mathematics is given in Adcock et al. [10]. Multivariate skew t-distribution has also been successfully used as a data model in environmental studies and modelling coastal flood- ings (Ghizzolini et al. [11], Thompson and Shen [12]). Several other applications of skew Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. elliptical distributions can be found in the collective monograph edited by Genton [13] This article is an open access article and in the book by Azzalini and Capitanio [14] which includes a comprehensive list of distributed under the terms and references on the topic. Another useful property that makes multivariate t-distribution conditions of the Creative Commons and skew t-distribution attractive in financial and actuarial applications is the possibility Attribution (CC BY) license (https:// to take tail dependence between marginals into account. Tail dependence of multivariate creativecommons.org/licenses/by/ t- and skew t-distributions has been intensively studied in recent years, see Fung and 4.0/). Seneta [15] or Joe and Li [16], for example. In Kollo et al. [17] it is shown that the upper
Symmetry 2021, 13, 1059. https://doi.org/10.3390/sym13061059 https://www.mdpi.com/journal/symmetry Symmetry 2021, 13, 1059 2 of 22
tail dependence coefficient of the skew t-distribution can be several times bigger than the corresponding value of multivariate t-distribution. Moments and inferential aspects re- lated to the multivariate skew t-distribution have been also actively studied in Padoan [18], Galarza et al. [19], Zhou and He [20], Lee and McLachlan [5], DiCiccio and Monti [21], for instance. A construction based on the skew t-distribution—skew t-copula—has made this distribution even more attractive in applications where correlated differently distributed marginals are joined into a multivariate distribution. Demarta and McNeil introduced skew t-copula based on multivariate generalized hyperbolic distribution [22], Kollo and Pettere [23] defined skew t-copula on the basis of the multivariate skew t-distribution following Azzalini and Capitanio [4]. The paper is organized as follows. In Section2, we introduce necessary notation and notions from matrix algebra. In Section3, skew t-distribution is defined through skew normal and inverse χ-distribution and the first four central moments are found via the moments of skew normal and inverse χ-distribution. The asymptotic normality for parameter estimators in the two-parameter case is also established in this section. In Section4, obtained analytical results are illustrated by simulations for different scenarios and the speed of convergence to asymptotic distribution is also addressed. In Section5, multivariate skew t-copula is defined and some properties are discussed. Technical proofs are presented in AppendixA. A comparison of density contours of Gaussian, skew normal and skew t-copula is presented in Supplementary Materials.
2. Notation and Preliminaries To derive asymptotic normal distribution for the parameter estimators of skew t- distribution, we use matrix technique. Necessary notions and notation are presented in what follows. Let us first denote the moments of a random vector x in the following matrix representation:
0 0 0 0 mk(x) = E x ⊗ x ⊗ x ⊗ x ⊗ ... = E xx ⊗ xx ⊗ ... . | {z } k times
The corresponding central moment is
mk(x) = mk(x − Ex).
We use notions vec -operator, commutation matrix, Kronecker product, matrix derivative and their properties. For deeper insight into this technique we refer to Magnus and Neudecker [24], or Kollo and von Rosen [25]. If A is a p × q matrix, A = (a1, ... , aq), the pq column vector vec A is formed by stacking the columns ai under each other; that is a1 vec A = .... aq
We denote a p × p identity matrix by Ip and a pq × pq commutation matrix by Kp,q. Commutation matrix Kp,q is an orthogonal pq × pq partitioned matrix consisting of q × p blocks with the property 0 Kp,qvec A = vec A where A is a p × q-matrix. The Kronecker product A ⊗ B of a p × q matrix A and an r × s matrix B is a pr × qs partitioned matrix consisting of the r × s blocks
A ⊗ B = [aijB], i = 1, . . . , p; j = 1, . . . , q.
The following properties of the notions above are repeatedly used. It is assumed that the dimensions of the matrices match with the operations. 1. (A ⊗ B)(C ⊗ D) = (AC) ⊗ (BD); Symmetry 2021, 13, 1059 3 of 22
2. ab0 = a ⊗ b0 = b0 ⊗ a; 3. vec (ABC) = (C0 ⊗ A)vec B; 4. Kp,1 = Ip, K1,q = Iq; 5. for A : p × q and B : r × s, A ⊗ B = Kp,r(B ⊗ A)Ks,q. Further, when deriving asymptotic distributions of the estimators of Σ and α in Section3 , we use the matrix derivative where partial derivatives ∂ykl are ordered by the ∂xij Kronecker product (see Magnus and Neudecker, [24])
dY d = ⊗ vec Y, dX dvec 0X with X : p × q and Y = Y(X) : r × s. It follows from the definition that
dY dvec Y dvec Y dY = = = . dX dvec 0X dX dvec 0X The main properties of the matrix derivative are listed below. dX 1. dX = Ipq; d(Y+Z) dY dZ 2. dX = dX + dX , Z = Z(X) : r × s; dAXB 0 3. dX = B ⊗ A; dZ dZ dY 4. If Z = Z(Y), Y = Y(X) then dX = dY dX ; 5. If W = YZ, Z = Z(X) : s × t then dW 0 dY dZ dX = (Z ⊗ Ir) dX + (It ⊗ Y) dX ; dX−1 0 −1 −1 6. If X : p × p then dX = −(X ) ⊗ X ; p×p d|X| −1 0 7. If X ∈ R then dX = |X|vec (X ) , where | · | denotes the determinant; 8. If W = Y ⊗ Z, Z = Z(X) : m × n, then dW h dY dZ i dX = (Is ⊗ Kr,n ⊗ Im) (Irs ⊗ vec Z) dX + (vec Y ⊗ Imn) dX . 9. If y = y(X), Z = Z(X) : m × n, then d(yZ) dy dZ dX = vec Z dX + y dX .
3. Multivariate Skew tν-Distribution 3.1. Definition Multivariate skew t-distribution has been defined in different ways, see Kotz and Nadarajah [26] (Chapter 5). Kshirsagar [27] was the first to introduce a noncentral multi- variate t-distribution, but the density expression was complicated and the distribution was not used in applications. From more recent definitions we refer to Fang et al. [28], but as the marginals in their definition have different degrees of freedom, the distribution had problems with practical use. Gupta [29] defined the multivariate skew t-distribution in a similar way like Azzalini and Capitanio [4], but their definition had several advantages and has become widely used in applications (for example, it is possible to use explicit expression of the moment generating function of skew normal distribution for calculating moments of skew t-distribution). We follow the approach by Azzalini and Capitanio [4]. A p-vector x has multivariate skew t-distribution with ν degrees of freedom if it is the ratio of a p-variate skew normal vector and a chi-distributed random variable with ν degrees of freedom. We have slightly modified their definition of the skew normal vector and use the covariance matrix as a parameter instead of the correlation matrix. Let y be a p-variate skew normal vector with density
( ) = ( ) ( 0 ) f y 2 fNp(0,Σ) y Φ α y , (1)
( ) where fNp(0,Σ) y is the density function (density) of the p-variate normal distribution Np(0, Σ) and Φ(·) is the cumulative distribution function (distribution function) of N(0, 1). Symmetry 2021, 13, 1059 4 of 22
Denote the distribution of vector y with density (1) by y ∼ SNp(Σ, α), where Σ : p × p is the scale parameter and α ∈ Rp is the shape or skewing parameter.
Definition 1. Let y be skew normally distributed with density (1) and let Z be χ-distributed with ν degrees of freedom, independent from y. Then √ y x = µ + ν (2) Z
is skew tν-distributed, x ∼ Stν(µ, Σ, α).
The density function of x ∼ Stν(µ, Σ, α) is of the form (a modification of Azzalini and Capitanio [4]) r 0 p + ν f (x) = 2 f (x)F + α (x − µ) , x ν,p p ν Q + ν 0 −1 where Q = (x − µ) Σ (x − µ), Fk is the distribution function of the univariate t-distribution with k degrees of freedom and fν,p is the density function of the p-variate t-distribution with ν degrees of freedom and scale matrix Σ
− ν+p ν+p " # 2 Γ (x − µ)0Σ−1(x − µ) f (x) = 2 1 + . ν,p p/2 ν 1/2 (νπ) Γ 2 |Σ| ν
3.2. Moments By Formula (2), a skew t-distributed x is constructed using independent skew normal random vector y and χ-distributed Z. Because of the independence of y and Z, the moments of x can be found as products of the moments of the skew normal and inverse χ-distribution. The k-th moment, central moment and cumulant of a skew-normal random vector y are denoted by mk(y), mk(y) and ck(y), respectively. Explicit expressions of the raw moments, central moments and cumulants of y can be obtained by differentiating the moment generating function
0 1 t0Σt α Σt M(t) = 2e 2 Φ √ 1 + α0Σα
and the cumulant generating function
C(t) = ln M(t).
The first four cumulants for a skew-normal random vector y are (Kollo et al. [30]) r 2 2 c (y) = Ey = δ, c (y) = Dy = Σ − δδ0, (3) 1 π 2 π r 2 4 c (y) = m (y) = − 1 δ ⊗ δ0 ⊗ δ, (4) 3 3 π π 8 3 c (y) = 1 − δ ⊗ δ0 ⊗ δ ⊗ δ0, 4 π π where δ equals Σα δ = √ . (5) 1 + α0Σα From general relation between the fourth order cumulants and central moments we have (see, e.g., Kollo and von Rosen [25] (p. 187))
0 m4(y) = c4(y) + (Ip2 + Kp,p)(Dy ⊗ Dy) + vec (Dy)vec (Dy). (6) Symmetry 2021, 13, 1059 5 of 22
After some algebra, we get the expression for α from (5):
Σ−1δ α = p . (7) 1 − δ0Σ−1δ The quantity under square root in (7) influences the stability of parameter estimates and will be used later in a related discussion. Thus, we add a notation for it:
η = 1 − δ0Σ−1δ. (8) 0 Note also that each component δi in vector δ = (δ1, ... , δp) is directly connected to the corresponding marginal distribution, whereas this property does not hold for vector 0 α = (α1, ... , αp) . For a more comprehensive discussion on parametrizations of skew- normal vector see Käärik et al. [31]. To derive the central moments of the multivariate skew tν-distribution, we also need the moments of the skew normally distributed vector y
m2(y) = Σ; (9) r 2 m (y) = (vec Σδ0 + δ ⊗ Σ + Σ ⊗ δ − δδ0 ⊗ δ). (10) 3 π 1 Recall that a random variable Vν = Z is inverse χν-distributed if its probability density function is (Lee [32] (p. 378))
1 1 f (x) = x−ν−1 exp − . ν/2−1 ν 2 2 Γ 2 2x
The moments of the inverse χν-distribution can be obtained by integration.
Lemma 1. The moments of the inverse χν-distribution are:
ν−1 Γ 2 EVν = √ , ν > 1; ν 2Γ 2 Γ ν−k k 2 EVν = , ν > k. k/2 ν 2 Γ 2
Proof. By standard integration technique, we can derive the expectation:
Z ∞ 1 1 EV = x x−ν−1 exp − dx ν ν/2−1 ν 2 0 2 Γ 2 2x − − 1 Z ∞ Γ ν 1 1 Γ ν 1 = √2 x−ν exp − dx = 2√ . ν ν−1 − − 2 ν Γ 0 2 1 ν 1 2x Γ 2 2 2 2Γ 2 2 Following the same approach, the k-th order moment can be derived as follows Symmetry 2021, 13, 1059 6 of 22
Z ∞ 1 1 EVk = xk x−ν−1 exp − dx ν ν/2−1 ν 2 0 2 Γ 2 2x ν−1 −1 ν−1 2 2 Γ Z ∞ 1 1 = 2 x−(ν−1)−1xk−1 exp − dx ν −1 ν ν−1 − − 2 2 2 Γ 0 2 1 ν 1 2x 2 2 Γ 2 k−1 Γ ν−k = = 2 > ∏ EVν−i k , ν k. 2 ν i=0 2 Γ 2
The expressions of the second, third and fourth moment are special cases of this result:
1 EV2 = , ν > 2; ν ν − 2 EV EV3 = ν , ν > 3; ν ν − 3 1 EV4 = , ν > 4. ν (ν − 2)(ν − 4) After some elementary algebra we get expressions of the corresponding central moments
7 − 2ν m (V ) = EV + 2(EV )2 ; 3 ν ν (ν − 3)(ν − 2) ν
1 ν − 5 m (V ) = + 2 (EV )2 − 3(EV )4. 4 ν (ν − 2)(ν − 4) (ν − 2)(ν − 3) ν ν
To keep the notation simpler, we denote V = Vν. In order to derive the asymptotic normality for parameter estimators of Stν, we need the first 4 central moments of the distribution. Let us start with general expressions for central moments m3(x) and m4(x).
Lemma 2. Let x be a p-dimensional random vector with finite fourth order moments. Then the third and fourth central moment of x can be expressed as follows:
0 0 m3(x) = m3(x) − (Ip2 + Kp,p)(m2(x) ⊗ Ex) − vec m2(x)Ex + 2ExEx ⊗ Ex, 0 0 0 m4(x) = m4(x) − C(x)(Ip2 + Kp,p) − (Ip2 + Kp,p)(C(x)) − 3ExEx ⊗ ExEx 0 + (Ip2 + Kp,p)A(x) + A(x)Kp,p + B(x) + (A(x) + B(x)) ,
where
0 0 0 0 A(x) = Ex ⊗ m2(x) ⊗ Ex; B(x) = Ex ⊗ vec m2(x) ⊗ Ex ; C(x) = m3(x) ⊗ Ex .
Proof of the lemma is presented in Appendix A.1. Applying this lemma to skew t-distributed x, the next result follows.
Lemma 3. Let x ∼ Stν(µ, Σ, α). Then its mean Ex, covariance matrix Dx, and central moments m3(x) and m4(x) can be expressed through the moments of inverse χν-distributed V and y ∼ SNp(Σ, α) in the following form: √ Ex = µ + νEVEy; (11)
0 0 Dx = νm2(V)m2(y) − ExEx = ν(m2(V)Dy + DVEyEy ); (12) Symmetry 2021, 13, 1059 7 of 22
3/2 3 3 3 0 m3(x) = ν {EV m3(y) − 2(EV − (EV) )(EyEy ⊗ Ey) 3 2 0 + (EV − EVEV )[(Ip2 + Kp,p)(m2(y) ⊗ Ey) + vec m2(y)Ey ]}; (13)
2 4 4 3 0 m4(x) = ν {EV m4(y) + (EV − EVEV )(C(y)(Ip2 + Kp,p) + (Ip2 + Kp,p)C(y) ) + 3(EV4 − (EV)4)EyEy0 ⊗ EyEy0 4 2 2 0 − (EV − EV (EV) )[(Ip2 + Kp,p)A(y) + A(y)Kp,p + B(y) + (A(y) + B(y)) ]}. (14)
Proof is straightforward and follows from Lemma2 after substituting in the required moments of skew normal distribution from Formulas (9) and (10) and the moments of inverse χ-distribution given in Lemma1.
3.3. Point Estimates Let us now focus on the case when µ = 0 and estimate the parameters of the two- parameter distribution Stν(Σ, α) using the method of moments. Let us have a sample of size n, n > p, from Stν(Σ, α) with sample mean x and unbiased sample covariance matrix S. Now, the estimates for parameters Σ and α through the moments of V and y can be obtained from Expressions (11) and (12):