A trivariate chi-squared distribution derived from the complex Wishart distribution

The joint density for a particular trivariate chi-squared distribution given by the diagonal elements of a complex Wishart matrix is derived. This distribution has applications in the processing of multilook synthetic aperture radar data. The expression for the density is in the form of an infinite series that converges rapidly and is simple and fast to compute. The expression is shown to reduce to known forms for a number of special cases and is validated by simulation. The characteristic function is also derived and used to relate joint moments of the trivariate distribution to the parameters of the density function.     

On Wishart distribution: Some extensions

This paper proposes a unified approach that enables the Wishart distribution to be studied simultaneously in the real, complex, quaternion and octonion cases under elliptical models. In particular, the matrix multivariate elliptical distribution, the noncentral generalised Wishart distribution, the joint density of the eigenvalues and the distribution of the maximum eigenvalue are obtained for real normed division algebras.     

Approximating by the Wishart distribution

Approximations of density functions are considered in the multivariate case. The results are presented with the help of matrix derivatives, powers of Kronecker products and Taylor expansions of functions with matrix argument. In particular, an approximation by the Wishart distribution is discussed. It is shown that in many situations the distributions should be centred. The results are applied to the approximation of the distribution of the sample covariance matrix and to the distribution of the non-central Wishart distribution.     

Use of the hyperbolic differential operator to find a scalar statistic which has constant regression on the mean of a sample of Wishart matrices

Scalar polynomial statistics are found which have constant regression on the mean of a sample of Wishart matrices. The method used is to differentiate the characteristic function associated with the Wishart distribution, thus expressing the constant regression condition as a differential equation which is satisfied by the Wishart characteristic function. In this respect, use is made of the hyperbolic differential operator.     

Noncentral matrix quadratic forms of the skew elliptical variables

In this paper, the noncentral matrix quadratic forms of the skew elliptical variables are studied. A family of the matrix variate noncentral generalized Dirichlet distributions is introduced as the extension of the noncentral Wishart distributions, the Dirichlet distributions and the noncentral generalized Dirichlet distributions. Main distributional properties are investigated. These include probability density and closure property under linear transformation and marginalization, the joint distribution of the sub-matrices of the matrix quadratic forms in the skew elliptical variables and the moment generating functions and Bartlett's decomposition of the matrix quadratic forms in the skew normal variables. Two versions of the noncentral Cochran's Theorem for the matrix variate skew normal distributions are obtained, providing sufficient and necessary conditions for the quadratic forms in the skew normal variables to have the matrix variate noncentral generalized Dirichlet distributions. Applications include the properties of the least squares estimation in multivariate linear model and the robustness property of the Wilk's likelihood ratio statistic in the family of the matrix variate skew elliptical distributions.     

Wishart ratios with dependent structure: New members of the bimatrix beta type IV

In multivariate statistics under normality, the problems of interest are random covariance matrices (known as Wishart matrices) and “ratios” of Wishart matrices that arise in multivariate analysis of variance (MANOVA) (see 24). The bimatrix variate beta type IV distribution (also known in the literature as bimatrix variate generalised beta; matrix variate generalization of a bivariate beta type I) arises from “ratios” of Wishart matrices. In this paper, we add a further independent Wishart random variate to the “denominator” of one of the ratios; this results in deriving the exact expression for the density function of the bimatrix variate extended beta type IV distribution. The latter leads to the proposal of the bimatrix variate extended F distribution. Some interesting characteristics of these newly introduced bimatrix distributions are explored. Lastly, we focus on the bivariate extended beta type IV distribution (that is an extension of bivariate Jones’ beta) with emphasis on P(X₁<X₂) where X₁ is the random stress variate and X₂ is the random strength variate.     

On the structure of the Wishart distribution

In this paper it is shown that every nonnegative definite symmetric random matrix with independent diagonal elements and at least one nondegenerate nondiagonal element has a noninfinitely divisible distribution. Using this result it is established that every Wishart distribution W_p(k, Σ, M) with both p and rank (Σ) ≥ 2 is noninfinitely divisible. The paper also establishes that any Wishart matrix having distribution W_p(k, Σ, 0) has the joint distribution of its elements in the rth row and rth column to be infinitely divisible for every r = 1,2,…,p.     

On the evaluation of some distributions that arise in simultaneous tests for the equality of the latent roots of the covariance matrix

In this paper, the authors consider the evaluation of the distribution functions of the ratios of the intermediate roots to the trace of the real Wishart matrix as well as the ratios of the individual roots to the trace of the complex Wishart matrix. In addition, the authors consider the evaluation of the distribution functions of the ratios of the extreme roots of the Wishart matrix in the real and complex cases. Some applications and tables of the above distributions are also given.     

The distribution of the characteristic roots ofS 1 S 2 −1 under violations in the complex case and power comparisons of four tests

Summary The joint density function of the latent roots ofS ₁ S ₂ ⁻¹ under violations is obtained whereS ₁ has a complex non-central Wishart distributionW _c (p,n ₁,Σ ₁,Ω) andS ₂, an independent complex central Wishart,W _c (p,n ₂,Σ ₂, 0). The density and moments of Hotelling's trace are also derived under violations. Further, the non-null distributions of the following four criteria in the two-roots case are studied for tests of three hypotheses: Hotelling's trace, Pillai's trace, Wilks' criterion and Roy's largest root. In addition, tabulations of powers are carried out and power comparisons for tests of each of three hypotheses based on the four criteria are made in the complex case extending such work of Pillai and Jayachandran in the classical Gaussian case. The findings in the complex Gaussian are generally similar to those in the classical.     

Estimation of Wishart mean matrices under simple tree ordering

For Wishart density functions, we study the risk dominance problems of the restricted maximum likelihood estimators of mean matrices with respect to the Kullback-Leibler loss function over restricted parameter space under the simple tree ordering set. The results are directly applied to the estimation of covariance matrices for the completely balanced multivariate multi-way random effects models without interactions.     

Importance Sampling for p-value Computations in Multivariate Tests

Abstract

An importance sampling procedure is developed to approximate the distribution of an arbitrary function of the eigenvalues for a matrix beta random matrix or a Wishart random matrix. The procedure is easily implemented and provides confidence intervals for the p-values of many of the commonly used test statistics in multivariate analysis. An adaptive procedure allows for the control of either absolute error or relative error in this p-value estimation through the choice of importance sample size.     

Minimax and admissible minimax estimators of the mean of a multivariate normal distribution for unknown covariance matrix

Let X be a p-variate (p ≥ 3) vector normally distributed with mean μ and covariance Σ, and let A be a p × p random matrix distributed independent of X, according to the Wishart distribution W(n, Σ). For estimating μ, we consider estimators of the form δ = δ(X, A). We obtain families of Bayes, minimax and admissible minimax estimators with respect to the quadratic loss function (δ ? μ)′ Σ^?1(δ ? μ) where Σ is unknown. This paper extends previous results of the author [1], given for the case in which the covariance matrix of the distribution is of the form σ²I, where σ is known.     

Properties of the singular, inverse and generalized inverse partitioned Wishart distributions

In this paper we discuss the distributions and independency properties of several generalizations of the Wishart distribution. First, an analog to Muirhead [R.J. Muirhead, Aspects of Multivariate Statistical Theory, Wiley, New York, 1982] Theorem 3.2.10 for the partitioned matrix is established in the case of arbitrary partitioning for singular and inverse Wishart distributions. Second, the density of is derived in the case of singular, non-central singular, inverse and generalized inverse Wishart distributions. The importance of the derived results is illustrated with an example from portfolio theory.     

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.     

Expressions for some hypergeometric functions of matrix argument with applications

Reasonably simple expressions are given for some hypergeometric functions when the size of the argument matrix or matrices is two. Applications of these expressions in connection with the distributions of the latent roots of a 2 × 2 Wishart matrix are also given.     

Nonnull distribution of likelihood ratio criterion for reality of covariance matrix

In this paper the distribution of the likelihood ratio test for testing the reality of the covariance matrix of a complex multivariate normal distribution is investigated. Some simplifications in the noncentral distribution are made and the noncentral distribution is derived for the special case where the rank of the noncentrality matrix is two. In the null case exact expressions for the distribution are given up to p = 6, and percentage points are tabulated. These percentage points were compared with percentage points derived from an asymptotic expansion of the distribution, and the accuracy of the approximation was found to be sufficient for several practical situations.     

Asymptotic nonnull distributions for tests for reality of a covariance matrix

In this paper asymptotic nonnull distributions are derived for two statistics used in testing for the reality of the covariance matrix in a complex Gaussian distribution.     

On the projections of the Dirichlet probability distribution on symmetric cones

In this paper, we introduce the Riesz-Dirichlet distribution on a symmetric cone as an extension of the Dirichlet distribution defined by the Wishart distribution. We also show that some projections of these distributions related to the Pierce decomposition are also Dirichlet.     

Asymptotic distribution of Wishart matrix for block-wise dispersion of population eigenvalues

This paper deals with the asymptotic distribution of Wishart matrix and its application to the estimation of the population matrix parameter when the population eigenvalues are block-wise infinitely dispersed. We show that the appropriately normalized eigenvectors and eigenvalues asymptotically generate two Wishart matrices and one normally distributed random matrix, which are mutually independent. For a family of orthogonally equivariant estimators, we calculate the asymptotic risks with respect to the entropy or the quadratic loss function and derive the asymptotically best estimator among the family. We numerically show (1) the convergence in both the distributions and the risks are quick enough for a practical use, (2) the asymptotically best estimator is robust against the deviation of the population eigenvalues from the block-wise infinite dispersion.     

A note on a Mar?enko-Pastur type theorem for time series

In this note we develop an extension of the Mar?enko-Pastur theorem to time series model with temporal correlations. The limiting spectral distribution (LSD) of the sample covariance matrix is characterised by an explicit equation for its Stieltjes transform depending on the spectral density of the time series. A numerical algorithm is then given to compute the density functions of these LSD’s.