On Pearson-Kotz Dirichlet distributions

In this paper, we discuss some basic distributional and asymptotic properties of the Pearson-Kotz Dirichlet multivariate distributions. These distributions, which appear as the limit of conditional Dirichlet random vectors, possess many appealing properties and are interesting from theoretical as well as applied points of view. We illustrate an application concerning the approximation of the joint conditional excess distribution of elliptically symmetric random vectors.     

High-dimensional generation of Bernoulli random vectors

The objective of this paper is to explore different modeling strategies to generate high-dimensional Bernoulli vectors. We discuss the multivariate Bernoulli (MB) distribution, probe its properties and examine three models for generating random vectors. A latent multivariate normal model whose bivariate distributions are approximated with Plackett distributions with univariate normal distributions is presented. A conditional mean model is examined where the conditional probability of success depends on previous history of successes. A mixture of beta distributions is also presented that expresses the probability of the MB vector as a product of correlated binary random variables. Each method has a domain of effectiveness. The latent model offers unpatterned correlation structures while the conditional mean and the mixture model provide computational feasibility for high-dimensional generation of MB vectors.     

A multivariate nonparametric test of independence

A new nonparametric approach to the problem of testing the joint independence of two or more random vectors in arbitrary dimension is developed based on a measure of association determined by interpoint distances. The population independence coefficient takes values between 0 and 1, and equals zero if and only if the vectors are independent. We show that the corresponding statistic has a finite limit distribution if and only if the two random vectors are independent; thus we have a consistent test for independence. The coefficient is an increasing function of the absolute value of product moment correlation in the bivariate normal case, and coincides with the absolute value of correlation in the Bernoulli case. A simple modification of the statistic is affine invariant. The independence coefficient and the proposed statistic both have a natural extension to testing the independence of several random vectors. Empirical performance of the test is illustrated via a comparative Monte Carlo study.     

A family of kurtosis orderings for multivariate distributions

In this paper, a family of kurtosis orderings for multivariate distributions is proposed and studied. Each ordering characterizes in an affine invariant sense the movement of probability mass from the “shoulders” of a distribution to either the center or the tails or both. All even moments of the Mahalanobis distance of a random vector from its mean (if exists) preserve a subfamily of the orderings. For elliptically symmetric distributions, each ordering determines the distributions up to affine equivalence. As applications, the orderings are used to study elliptically symmetric distributions. Ordering results are established for three important families of elliptically symmetric distributions: Kotz type distributions, Pearson Type VII distributions, and Pearson Type II distributions.     

Generalized skew-elliptical distributions and their quadratic forms

This paper introduces generalized skew-elliptical distributions (GSE), which include the multivariate skew-normal, skew-t, skew-Cauchy, and skew-elliptical distributions as special cases. GSE are weighted elliptical distributions but the distribution of any even function in GSE random vectors does not depend on the weight function. In particular, this holds for quadratic forms in GSE random vectors. This property is beneficial for inference from non-random samples. We illustrate the latter point on a data set of Australian athletes.     

Multivariate analysis of variance with fewer observations than the dimension

In this article, we consider the problem of testing a linear hypothesis in a multivariate linear regression model which includes the case of testing the equality of mean vectors of several multivariate normal populations with common covariance matrix Σ, the so-called multivariate analysis of variance or MANOVA problem. However, we have fewer observations than the dimension of the random vectors. Two tests are proposed and their asymptotic distributions under the hypothesis as well as under the alternatives are given under some mild conditions. A theoretical comparison of these powers is made.     

Extreme behavior of multivariate phase-type distributions

This paper investigates the limiting distributions of the componentwise maxima and minima of suitably normalized iid multivariate phase-type random vectors. In the case of maxima, a large parametric class of multivariate extreme value (MEV) distributions is obtained. The flexibility of this new class is exemplified in the bivariate setup. For minima, it is shown that the dependence structure of the Marshall-Olkin class arises in the limit.     

渐近正态随机变量函数的极限分布

基于渐近正态随机变量,导出随机变量函数极限分布的两个一般性理论结果.作为应用,证明了渐近正态随机变量一系列具体函数的极限分布,其中包括泊松随机变量平方根的渐近正态性,以及随机变量部分和在正则化常数是随机变量情况下的渐近正态性.     

一般形式下斜正态随机向量的矩

本文给出一般形式下斜正态随机向量及其平方型的矩公式. 作为应用, 计算出了斜正态随机向量的多元偏度和峰度.     

On Affine Equivariant Multivariate Quantiles

An extension of univariate quantiles in the multivariate set-up has been proposed and studied. The proposed approach is affine equivariant, and it is based on an adaptive transformation retransformation procedure. Behadur type linear representations of the proposed quantiles are established and consequently asymptotic distributions are also derived. As applications of these multivariate quantiles, we develop some affine equivariant quantile contour plots which can be used to study the geometry of the data cloud as well as the underlying probability distribution and to detect outliers. These quantiles can also be used to construct affine invariant versions of multivariate Q-Q plots which are useful in checking how well a given multivariate probability distribution fits the data and for comparing the distributions of two data sets. We illustrate these applications with some simulated and real data sets. We also indicate a way of extending the notion of univariate L-estimates and trimmed means in the multivariate set-up using these affine equivariant quantiles.     

Pseudo-inverse multivariate/matrix-variate distributions

The Moore-Penrose inverse of a singular or nonsquare matrix is not only existent but also unique. In this paper, we derive the Jacobian of the transformation from such a matrix to the transpose of its Moore-Penrose inverse. Using this Jacobian, we investigate the distribution of the Moore-Penrose inverse of a random matrix and propose the notion of pseudo-inverse multivariate/matrix-variate distributions. For arbitrary multivariate or matrix-variate distributions, we can develop the corresponding pseudo-inverse distributions. In particular, we present pseudo-inverse multivariate normal distributions, pseudo-inverse Dirichlet distributions, pseudo-inverse matrix-variate normal distributions and pseudo-inverse Wishart distributions.     

An alternative method for computing mean and covariance matrix of some multivariate distributions

Computing the mean and covariance matrix of some multivariate distributions, in particular, multivariate normal distribution and Wishart distribution are considered in this article. It involves a matrix transformation of the normal random vector into a random vector whose components are independent normal random variables, and then integrating univariate integrals for computing the mean and covariance matrix of a multivariate normal distribution. Moment generating function technique is used for computing the mean and covariances between the elements of a Wishart matrix. In this article, an alternative method that uses matrix differentiation and differentiation of the determinant of a matrix is presented. This method does not involve any integration.     

Nonparametric tests of independence between random vectors

A nonparametric test of the mutual independence between many numerical random vectors is proposed. This test is based on a characterization of mutual independence defined from probabilities of half-spaces in a combinatorial formula of Möbius. As such, it is a natural generalization of tests of independence between univariate random variables using the empirical distribution function. If the number of vectors is p and there are n observations, the test is defined from a collection of processes R_n,A, where A is a subset of {1,…,p} of cardinality |A|>1, which are asymptotically independent and Gaussian. Without the assumption that each vector is one-dimensional with a continuous cumulative distribution function, any test of independence cannot be distribution free. The critical values of the proposed test are thus computed with the bootstrap which is shown to be consistent. Another similar test, with the same asymptotic properties, for the serial independence of a multivariate stationary sequence is also proposed. The proposed test works when some or all of the marginal distributions are singular with respect to Lebesgue measure. Moreover, in singular cases described in Section 4, the test inherits useful invariance properties from the general affine invariance property.     

Testing for affine equivalence of elliptically symmetric distributions

Let X and Y be d-dimensional random vectors having elliptically symmetric distributions. Call X and Y affinely equivalent if Y has the same distribution as AX+b for some nonsingular d×d-matrix A and some . This paper studies a class of affine invariant tests for affine equivalence under certain moment restrictions. The test statistics are measures of discrepancy between the empirical distributions of the norm of suitably standardized data.     

Joint limit distributions of exceedances point processes and partial sums of gaussian vector sequence

In this paper, we study the joint limit distributions of point processes of exceedances and partial sums of multivariate Gaussian sequences and show that the point processes and partial sums are asymptotically independent under some mild conditions. As a result, for a sequence of standardized stationary Gaussian vectors, we obtain that the point process of exceedances formed by the sequence (centered at the sample mean) converges in distribution to a Poisson process and it is asymptotically independent of the partial sums. The asymptotic joint limit distributions of order statistics and partial sums are also investigated under different conditions.     

On a robust and efficient maximum depth estimator

The best breakdown point robustness is one of the most outstanding features of the univariate median. For this robustness property, the median, however, has to pay the price of a low efficiency at normal and other light-tailed models. Affine equivariant multivariate analogues of the univariate median with high breakdown points were constructed in the past two decades. For the high breakdown robustness, most of them also have to sacrifice their efficiency at normal and other models, nevertheless. The affine equivariant maximum depth estimator proposed and studied in this paper turns out to be an exception. Like the univariate median, it also possesses a highest breakdown point among all its multivariate competitors. Unlike the univariate median, it is also highly efficient relative to the sample mean at normal and various other distributions, overcoming the vital low-efficiency shortcoming of the univariate and other multivariate generalized medians. The paper also studies the asymptotics of the estimator and establishes its limit distribution without symmetry and other strong assumptions that are typically imposed on the underlying distribution. This work was supported by Natural Science Foundation of USA (Grant Nos. DMS-0071976, DMS-0234078) and by the Southwestern University of Finance and Economics Third Period Construction Item Funds of the 211 Project (Grant No. 211D3T06)     

A robust and efficient adaptive reweighted estimator of multivariate location and scatter

This article proposes a reweighted estimator of multivariate location and scatter, with weights adaptively computed from the data. Its breakdown point and asymptotic behavior under elliptical distributions are established. This adaptive estimator is able to attain simultaneously the maximum possible breakdown point for affine equivariant estimators and full asymptotic efficiency at the multivariate normal distribution. For the special case of hard-rejection weights and the MCD as initial estimator, it is shown to be more efficient than its non-adaptive counterpart for a broad range of heavy-tailed elliptical distributions. A Monte Carlo study shows that the adaptive estimator is as robust as its non-adaptive relative for several types of bias-inducing contaminations, while it is remarkably more efficient under normality for sample sizes as small as 200.     

Limit theorems for conditional distributions with regard for large deviations

Possible limit laws are studied for the multivariate conditional distribution of a subset of components of the sum of independent identically distributed random vectors under the condition that other components belong to the domain of large deviations. It is assumed that the considered distribution is absolutely continuous and belongs to the domain of attraction of the normal law but possesses “heavy tails.” The approach suggested is based on the local theorem for large deviations. Institute of Mathematics, Ukrainian Academy of Sciences, Kiev; Kopernik Institute, Torun, Poland. Translated from Ukrainskii Matematicheskii Zhurnal, Vol. 51, No. 8, pp. 1054–1064, August, 1999.     

Multivariate quadratic forms of random vectors

We obtain the distribution of the sum of n random vectors and the distribution of their quadratic forms: their densities are expanded in series of Hermite and Laguerre polynomials. We do not suppose that these vectors are independent. In particular, we apply these results to multivariate quadratic forms of Gaussian vectors. We obtain also their densities expanded in Mac Laurin series or in the form of an integral. By this last result, we introduce a new method of computation which can be much simpler than the previously known techniques. In particular, we introduce a new method in the very classical univariate case. We remark that we do not assume the independence of normal variables.     

On a new multivariate IFR ageing notion based on the standard construction

Many criteria of ageing for random variables or vectors have been proposed in the literature over many years. For instance, a random variable is increasing in failure rate (IFR) if, and only if, it can be ordered with an exponentially distributed random variable in the classical univariate convex transform order. A new multivariate generalization of the convex transform order has recently been proposed in the literature. In this work, we propose a new multivariate IFR notion for multivariate distributions based on comparisons in this new order with a properly defined exponentially distributed random vector. Properties, applications, and illustrations of this new notion are given as well. Copyright © 2016 John Wiley & Sons, Ltd.