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.     

Local asymptotic normality in a stationary model for spatial extremes

De Haan and Pereira (2006) [6] provided models for spatial extremes in the case of stationarity, which depend on just one parameter β>0 measuring tail dependence, and they proposed different estimators for this parameter. We supplement this framework by establishing local asymptotic normality (LAN) of a corresponding point process of exceedances above a high multivariate threshold. Standard arguments from LAN theory then provide the asymptotic minimum variance within the class of regular estimators of β. It turns out that the relative frequency of exceedances is a regular estimator sequence with asymptotic minimum variance, if the underlying observations follow a multivariate extreme value distribution or a multivariate generalized Pareto distribution.     

On Mardia's and Song's measures of kurtosis in elliptical distributions

The main objective of this paper is the calculation and the comparative study of two general measures of multivariate kurtosis, namely Mardia's measure β_2,p and Song's measure S(f). In this context, general formulas for the said measures are derived for the broad family of the elliptically contoured symmetric distributions and also for specific members of this family, like the multivariate t-distribution, the multivariate Pearson type II, the multivariate Pearson type VII, the multivariate symmetric Kotz type distribution and the uniform distribution in the unit sphere. Analytic expressions for computing Shannon and Rényi entropies are obtained under the elliptic family. The behaviour of Mardia's and Song's measures, their similarities and differences, possible interpretations and uses in practice are investigated by comparing them in specific members of the elliptic family of multivariate distributions. An empirical estimator of Song's measure is moreover proposed and its asymptotic distribution is investigated under the elliptic family of multivariate distributions.     

Approximation and estimation of very small probabilities of multivariate extreme events

This article discusses modelling of the tail of a multivariate distribution function by means of a large deviation principle (LDP), and its application to the estimation of the probability p _n of a multivariate extreme event from a sample of n iid random vectors, with \(p_{n}\in [n^{-\tau _{2}},n^{-\tau _{1}}]\) for some t ₁>1 and t ₂>t ₁. One way to view the classical tail limits is as limits of probability ratios. In contrast, the tail LDP provides asymptotic bounds or limits for log-probability ratios. After standardising the marginals to standard exponential, tail dependence is represented by a homogeneous rate function I. Furthermore, the tail LDP can be extended to represent both dependence and marginals, the latter implying marginal log-Generalised Weibull tail limits. A connection is established between the tail LDP and residual tail dependence (or hidden regular variation) and a recent extension of it. Under a smoothness assumption, they are implied by the tail LDP. Based on the tail LDP, a simple estimator for very small probabilities of extreme events is formulated. It avoids estimation of I by making use of its homogeneity. Strong consistency in the sense of convergence of log-probability ratios is proven. Simulations and an application illustrate the difference between the classical approach and the LDP-based approach.     

The matrix-t distribution and its applications in predictive inference

The predictive distributions of the future responses and regression matrix under the multivariate elliptically contoured distributions are derived using structural approach. The predictive distributions are obtained as matrix-t which are identical to those obtained under matrix normal and matrix-t distributions. This gives inference robustness with respect to departures from the reference case of independent sampling from the matrix normal or dependent but uncorrelated sampling from matrix-t distributions. Some successful applications of matrix-t distribution in the field of spatial prediction have been addressed.     

The multivariate Behrens-Fisher distribution

The main purpose of this paper is the study of the multivariate Behrens-Fisher distribution. It is defined as the convolution of two independent multivariate Student t distributions. Some representations of this distribution as the mixture of known distributions are shown. An important result presented in the paper is the elliptical condition of this distribution in the special case of proportional scale matrices of the Student t distributions in the defining convolution. For the bivariate Behrens-Fisher problem, the authors propose a non-informative prior distribution leading to highest posterior density (H.P.D.) regions for the difference of the mean vectors whose coverage probability matches the frequentist coverage probability more accurately than that obtained using the independence-Jeffreys prior distribution, even with small samples.     

Orthant tail dependence of multivariate extreme value distributions

The orthant tail dependence describes the relative deviation of upper- (or lower-) orthant tail probabilities of a random vector from similar orthant tail probabilities of a subset of its components, and can be used in the study of dependence among extreme values. Using the conditional approach, this paper examines the extremal dependence properties of multivariate extreme value distributions and their scale mixtures, and derives the explicit expressions of orthant tail dependence parameters for these distributions. Properties of the tail dependence parameters, including their relations with other extremal dependence measures used in the literature, are discussed. Various examples involving multivariate exponential, multivariate logistic distributions and copulas of Archimedean type are presented to illustrate the results.     

Tail dependence between order statistics

In this work, we introduce the s,k-extremal coefficients for studying the tail dependence between the s-th lower and k-th upper order statistics of a normalized random vector. If its margins have tail dependence then so do their order statistics, with the strength of bivariate tail dependence decreasing as two order statistics become farther apart. Some general properties are derived for these dependence measures which can be expressed via copulas of random vectors. Its relations with other extremal dependence measures used in the literature are discussed, such as multivariate tail dependence coefficients, the coefficient η of tail dependence, coefficients based on tail dependence functions, the extremal coefficient ?, the multivariate extremal index and an extremal coefficient for min-stable distributions. Several examples are presented to illustrate the results, including multivariate exponential and multivariate Gumbel distributions widely used in applications.     

On a new multivariate two-sample test

In this paper we propose a new test for the multivariate two-sample problem. The test statistic is the difference of the sum of all the Euclidean interpoint distances between the random variables from the two different samples and one-half of the two corresponding sums of distances of the variables within the same sample. The asymptotic null distribution of the test statistic is derived using the projection method and shown to be the limit of the bootstrap distribution. A simulation study includes the comparison of univariate and multivariate normal distributions for location and dispersion alternatives. For normal location alternatives the new test is shown to have power similar to that of the t- and T²-Test.     

Peaks-over-threshold stability of multivariate generalized Pareto distributions

It is well-known that the univariate generalized Pareto distributions (GPD) are characterized by their peaks-over-threshold (POT) stability. We extend this result to multivariate GPDs.It is also shown that this POT stability is asymptotically shared by distributions which are in a certain neighborhood of a multivariate GPD. A multivariate extreme value distribution is a typical example.The usefulness of the results is demonstrated by various applications. We immediately obtain, for example, that the excess distribution of a linear portfolio with positive weights a_i, i≤d, is independent of the weights, if (U₁,…,U_d) follows a multivariate GPD with identical univariate polynomial or Pareto margins, which was established by Macke [On the distribution of linear combinations of multivariate EVD and GPD distributed random vectors with an application to the expected shortfall of portfolios, Diploma Thesis, University of Würzburg, 2004, (in German)] and Falk and Michel [Testing for tail independence in extreme value models. Ann. Inst. Statist. Math. 58 (2006) 261-290]. This implies, for instance, that the expected shortfall as a measure of risk fails in this case.     

Multivariate nonparametric tests in a randomized complete block design

In this paper multivariate extensions of the Friedman and Page tests for the comparison of several treatments are introduced. Related unadjusted and adjusted treatment effect estimates for the multivariate response variable are also found and their properties discussed. The test statistics and estimates are analogous to the traditional univariate methods. In test constructions, the univariate ranks are replaced by multivariate spatial ranks (J. Nonparam. Statist. 5 (1995) 201). Asymptotic theory is developed to provide approximations for the limiting distributions of the test statistics and estimates. Limiting efficiencies of the tests and treatment effect estimates are found in the multivariate normal and t distribution cases. The tests are rotation invariant only, but affine invariant versions can be easily constructed. The theory is illustrated by an example.     

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.     

Predictivistic characterizations of multivariate student-t models

De Finetti style theorems characterize models (predictive distributions) as mixtures of the likelihood function and the prior distribution, beginning from some judgment of invariance about observable quantities. The likelihood function generally has its functional form identified from invariance assumptions only. However, we need additional conditions on observable quantities (typically, assumptions on conditional expectations) to identify the prior distribution. In this paper, we consider some well-known invariance assumptions and establish additional conditions on observable quantities in order to obtain a predictivistic characterization of the multivariate and matrix-variate Student-t distributions as well as for the Student-t linear model. As a byproduct, a characterization for the Pearson type II distribution is provided.     

Finite sample tail behavior of multivariate location estimators

A finite sample performance measure of multivariate location estimators is introduced based on “tail behavior”. The tail performance of multivariate “monotone” location estimators and the halfspace depth based “non-monotone” location estimators including the Tukey halfspace median and multivariate L-estimators is investigated. The connections among the finite sample performance measure, the finite sample breakdown point, and the halfspace depth are revealed. It turns out that estimators with high breakdown point or halfspace depth have “appealing” tail performance. The tail performance of the halfspace median is very appealing and also robust against underlying population distributions, while the tail performance of the sample mean is very sensitive to underlying population distributions. These findings provide new insights into the notions of the halfspace depth and breakdown point and identify the important role of tail behavior as a quantitative measure of robustness in the multivariate location setting.     

Tail Dependence Comparison of Survival Marshall–Olkin Copulas

The multivariate tail dependence describes 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. We derive an explicit expression of tail dependence of multivariate survival Marshall–Olkin copulas, and obtain a sufficient condition under which tail dependencies of two survival Marshall–Olkin copulas can be compared. Some examples are also presented to illustrate our results.      

More on the inadmissibility of step-up

Cohen and Sackrowitz [Characterization of Bayes procedures for multiple endpoint problems and inadmissibility of the step-up procedure, Ann. Statist. 33 (2005) 145-158] proved that the step-up multiple testing procedure is inadmissible for a multivariate normal model with unknown mean vector and known intraclass covariance matrix. The hypotheses tested are each mean is zero vs. each mean is positive. The risk function is a 2×1 vector where one component is average size and the other component is one minus average power. In this paper, we extend the inadmissibility result to several different models, to two-sided alternatives, and to other risk functions. The models include one-parameter exponential families, independent t-variables, independent χ²-variables, t-tests arising from the analysis of variance, and t-tests arising from testing treatments against a control. The additional risk functions are linear combinations where one component is the false discovery rate (FDR).     

A new representation of Student's t as a function of independent t's,with a generalization to the matrix t

A new synthetic representation of Student's t statistic is given which exhibits the random variable as a nonlinear function of two independent t variables. A generalization is given to the matricvariate t with an interpretation in terms of “natural conjugate” inference for multivariate normal sampling.     

A multivariate extension of Sarhan and Balakrishnan’s bivariate distribution and its ageing and dependence properties

A new class of bivariate distributions (NBD) was recently introduced by Sarhan and Balakrishnan [A.M. Sarhan, N. Balakrishnan, A new class of bivariate distributions and its mixture, J. Multivariate Anal. 98 (2007) 1508-1527]. In this note, we give the joint survival function of a multivariate extension of the NBD, which is not an absolutely continuous multivariate distribution, and its marginal and extreme order statistics distributions are also derived. The multivariate ageing and dependence properties of the proposed n-dimensional distribution are also discussed, and then we analyze the stochastic ageing of its marginals and its minimum and maximum order statistics.     

On the Student's t-distribution and the t-statistic

In this paper we provide rather weak conditions on a distribution which would guarantee that the t-statistic of a random vector of order n follows the t-distribution with n-1 degrees of freedom. The results sharpen the earlier conclusions of Mauldon [Characterizing properties of statistical distributions, Quart. J. Math. 2(7) (1956) 155-160] and the more recent advances due to Bondesson [When is the t-statistic t-distributed, Sankhyā, Ser. A 45 (1983) 338-345]. The basic tool involved in the derivations is the vertical density representation originally suggested by Troutt [A theorem on the density of the density ordinate and an alternative interpretation of the Box-Muller method, Statistics 22(3) (1991) 463-466; Vertical density representation and a further remark on the Box-Muller method, Statistics 24 (1993) 81-83]. Several illustrative examples are presented.     

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.