On some generalized wishart expectations

In a recent paper Muirhead (1986) derived certain useful identities involving expectations taken with respect to the Wishart distribution. This note generalizes the above results by taking expectations with respect to a generalized version of the Wishart distribution, considered by Sutradhar and Ali (1989), based on a multivariate tdistribution.     

The Density Functions of the Singular Quaternion Normal Matrix and the Singular Quaternion Wishart Matrix

In this article, we give the density functions of the singular quaternion normal matrix and the singular quaternion Wishart matrix. Furthermore, we also give the density functions of certain singular quaternion β-matrix and the singular quaternion F-matrix in terms of the density function of the singular quaternion Wishart matrix and hypergeometric functions of quaternion matrix argument.     

Zonal Polynomials and Hypergeometric Functions of Quaternion Matrix Argument

We define zonal polynomials of quaternion matrix argument and deduce some impor-tant formulae of zonal polynomials and hypergeometric functions of quaternion matrix argument. As an application, we give the distributions of the largest and smallest eigenvalues of a quaternion central Wishart matrix W ～ ?W(n, Σ), respectively.     

SIMULATIONS AND DERIVED APPROXIMATIONS FOR THE MEANS AND STANDARD DEVIATIONS OF THE CHARACTERISTIC ROOTS OF A WISHART MATRIX

Monte Carlo simulations were done to estimate the means and standard deviations of the characteristic roots of a Wishart matrix which can be used in computing tests of hypotheses concerning multiplicative terms in balanced linear-bilinear (multiplicative) models for an m × n table of data. In this report we extend the previous results (Mandel, 1971; Cornelius, 1980) to r ≤ 199, c ≤ 149 or r ≤ 149, c ≤ 199, where r and c are row and column degrees of freedom, respectively, of the two-way array of residuals (with total degrees of freedom rc) after fitting the linear effects. For 187 combinations of r and c at intervals over this domain, we used 5000 simulations to estimate expectations and standard deviations of the Wishart roots. Using weighted linear regression variable selection techniques, symmetric functions of r and c were obtained for approximating the simulated means and standard deviations. Use of these approximating functions will avoid the need for reference to tables for input to computer programs which require these values for tests of significance of sequentially fitted terms in the analyses of balanced linear-bilinear models.     

Further identities for the Wishart distribution with applications in regression

Matrix analogues are given for a known scalar identity which relates certain expectations with respect to the Wishart distribution. (The scalar identity was independently derived by C. Stein and L. Haff.) The matrix analogues are more aptly called “matrix extensions.” They can be derived by using the scalar identity; nevertheless, they are seen (in quite elementary terms) to be more general than the latter. A method of doing multivariate calculations is developed from the identities, and several examples are worked in detail. We compute the first two moments of the regression coefficients and another matrix arising in regression analysis. Also, we give a new result for the matrix analogue of squared multiple correlation: the bias correction of Ezekiel (1930), a result often used in model building, is extended to the case of two or more dependent variables.     

Unbiased Estimator for a Covariance Matrix Under Two-Step Monotone Incomplete Sample

In this article, we consider an inference for a covariance matrix under two-step monotone incomplete sample. The maximum likelihood estimator of the mean vector is unbiased but that of the covariance matrix is biased. We derive an unbiased estimator for the covariance matrix using some fundamental properties of the Wishart matrix. The properties of the estimators are investigated and the accuracies are checked by a numerical simulation.     

New bivariate gamma types with MIMO application

In this paper a bivariate gamma type distribution emanating from the diagonal elements of an inverse Wishart type distribution is developed; which in turn originates from the complex matrix variate elliptical class. From this, a bivariate Weibullised gamma type distribution is also presented, of which the bivariate Nakagami-m type is a special case. The derived results may be applied as decision statistics for a MIMO (multiple input multiple output) system with two transmit antennas. It is proposed that under this elliptical umbrella some performance measures such as the outage probability of MIMO systems can be analyzed in broad generality.     

Equivariant estimators of the covariance matrix

Given a Wishart matrix S [S ∽ W_p(n, Σ)] and an independent multinomial vector X [X ∽ N_p (μ, Σ)], equivariant estimators of Σ are proposed. These estimators dominate the best multiple of S and the Stein-type truncated estimators.     

A Note on the Canonical Correlations From Contingency Tables

Explicit formulae are obtained for the asymptotic variances and covariances of canonical correlations which correspond to non-zero theoretical correlations in (p+ 1) x (q+1) contingency tables, with p≤q. The moments of the roots of a central Wishart matrix distributed as Wp(q; I ) are also given in general, with means, variances and covariances tabulated for p= 2, 3, 4: these may apply to canonical correlations corresponding to zeros.     

Further results on the trace of a noncentral wishart matrix

In an earlier article Mathai (1980) has given compact representations for the moments and cumulants of the trace of a noncentral Wishart matrix. He has also shown that (trA-ntr;∑)/(2ntri∑²)172. is asymptotically standard normal where A is a noncentral wishart matrix with n degrees of freedom and covariance matrix [0, In the present article explicit expressions for the exact density of the trace are given in terms of confluent hypergeometric functions and in terms of zonal polynomials for the general case and as finite sums when the sample size is odd. As a consequence of some of these representations some summation formulae for zonal polynomials are also given     

All Invariant Moments of the Wishart Distribution

Abstract.  In this paper, we compute moments of a Wishart matrix variate U of the form E ( Q ( U )) where Q ( u ) is a polynomial with respect to the entries of the symmetric matrix u , invariant in the sense that it depends only on the eigenvalues of the matrix u . This gives us in particular the expected value of any power of the Wishart matrix U or its inverse U ^{− 1}. For our proofs, we do not rely on traditional combinatorial methods but rather on the interplay between two bases of the space of invariant polynomials in U . This means that all moments can be obtained through the multiplication of three matrices with known entries. Practically, the moments are obtained by computer with an extremely simple Maple program.     

A Characterization of a Matric-Variate Exponential Distribution

Suppose S is a positive definite m x m random matrix and S> Ω denotes the event that S—Ω is positive definite, Ω being a constant positive definite matrix. Under very mild regularity conditions, we show that the constraint Pr(S>Ω₁+Ω₂ | S>Ω₂) = Pr(S>Ω₁) implies that S has a Wishart distribution on m+ 1 degrees of freedom.     

An asymptotic expansion of Wishart distribution when the population eigenvalues are infinitely dispersed

Takemura and Sheena [A. Takemura, Y. Sheena, Distribution of eigenvalues and eigenvectors of Wishart matrix when the population eigenvalues are infinitely dispersed and its application to minimax estimation of covariance matrix, J. Multivariate Anal. 94 (2005) 271–299] derived the asymptotic joint distribution of the eigenvalues and the eigenvectors of a Wishart matrix when the population eigenvalues become infinitely dispersed. They also showed necessary conditions for an estimator of the population covariance matrix to be tail minimax for typical loss functions by calculating the asymptotic risk of the estimator. In this paper, we further examine those distributions and risks by means of an asymptotic expansion. We obtain the asymptotic expansion of the distribution function of relevant elements of the sample eigenvalues and eigenvectors. We also derive the asymptotic expansion of the risk function of a scale and orthogonally equivariant estimator with respect to Stein’s loss. As an application, we prove non-minimaxity of Stein’s and Haff’s estimators, which has been an open problem for a long time.     

On the Distribution of Matrix Quadratic Forms

A characterization of the distribution of the multivariate quadratic form given by X A X′, where X is a p × n normally distributed matrix and A is an n × n symmetric real matrix, is presented. We show that the distribution of the quadratic form is the same as the distribution of a weighted sum of non central Wishart distributed matrices. This is applied to derive the distribution of the sample covariance between the rows of X when the expectation is the same for every column and is estimated with the regular mean.     

Decision-theoretic estimation of generalized variance and generalized precision

Consider a random data matrix X=(X₁,...,X_k):pXk with independent columns [sathik] and an independent p X p Wishart matrix [sathik]. Estimators dominating the best affine equivariant estimators of [sathik] are obtained under four types of loss functions. Improved estimators (Testimators) of generalized variance and generalized precision are also considered under convex entropy loss (CEL).     

Explosion time for some Laplace transforms of the Wishart process

In this article, we focus upon a family of matrix valued stochastic processes and study the problem of determining the smallest time such that their Laplace transforms become infinite. In particular, we concentrate upon the class of Wishart processes, which have proved to be very useful in different applications by their ability in describing non-trivial dependence. Thanks to this remarkable property we are able to explain the behavior of the explosion times for the Laplace transforms of the Wishart process and its time integral in terms of the relative importance of the involved factors and their correlations.     

A solution to the multivariate behrens-fisher problem

Some new algebra on pattern and transition matrices is used to determine the degrees of freedom and the parameter matrix, if the distribution of a linear sum of Wishart matrices is approximated by a single Wishart distribution. This approximation is then used to find a solution to the multivariate Behrens-Fisher problem similar to the Welch (1947) solution in the univariate case.     

Bayesian multivariate normal analysis with a wishart prior

This paper considers the Bayesian analysis of the multivariate normal distribution when its covariance matrix has a Wishart prior density under the assumption of a multivariate quadratic loss function. New flexible marginal posterior distributions of the mean μ and of the covariance matrix Σ are developed and univariate cases with graphical representations are given.     

Improved minimax estimation of a normal precision matrix

Let S_{p × p} have a Wishart distribution with parameter matrix Σ and n degrees of freedom. We consider here the problem of estimating the precision matrix Σ^?1 under the loss functions L₁(σ) tr (σ) - log |σ| and L₂(σ) = tr (σ). James-Stein-type estimators have been derived for an arbitrary p. We also obtain an orthogonal invariant and a diagonal invariant minimax estimator under both loss functions. A Monte-Carlo simulation study indicates that the risk improvement of the orthogonal invariant estimators over the James-Stein type estimators, the Haff (1979) estimator, and the “testimator” given by Sinha and Ghosh (1987) is substantial.     

THE LIMITING FORM OF THE NON-CENTRAL WISHART DISTRIBUTION1,2

Let S (p×p) have a Wishart distribution -with v degrees of freedom and non-centrality matrix θ= [θ_jK] (p×p). Define θ₀= min {| θ_jk |}, let θ₀→∞, and suppose that | θ_jK | = 0(θ_o). Then the limiting form of the standardized non-central distribution, as θ while n? remains fixed, is a multivariate Gaussian distribution. This result in turn is used to obtain known asymptotic properties of multivariate chi-square and Rayleigh distributions under somewhat weaker conditions.