Computing and approximating multivariate chi-square probabilities

We consider computational methods for evaluating and approximating multivariate chi-square probabilities in cases where the pertaining correlation matrix or blocks thereof have a low factorial representation. To this end, techniques from matrix factorization and probability theory are applied. We outline a variety of statistical applications of multivariate chi-square distributions and provide a system of MATLAB programs implementing the proposed algorithms. Computer simulations demonstrate the accuracy and the computational efficiency of our methods in comparison with Monte Carlo approximations, and a real data example from statistical genetics illustrates their usage in practice.     

A spatial Bayesian semiparametric mixture model for positive definite matrices with applications in diffusion tensor imaging

Studies on diffusion tensor imaging (DTI) quantify the diffusion of water molecules in a brain voxel using an estimated 3 × 3 symmetric positive definite (p.d.) diffusion tensor matrix. Due to the challenges associated with modelling matrix‐variate responses, the voxel‐level DTI data are usually summarized by univariate quantities, such as fractional anisotropy. This approach leads to evident loss of information. Furthermore, DTI analyses often ignore the spatial association among neighbouring voxels, leading to imprecise estimates. Although the spatial modelling literature is rich, modelling spatially dependent p.d. matrices is challenging. To mitigate these issues, we propose a matrix‐variate Bayesian semiparametric mixture model, where the p.d. matrices are distributed as a mixture of inverse Wishart distributions, with the spatial dependence captured by a Markov model for the mixture component labels. Related Bayesian computing is facilitated by conjugacy results and use of the double Metropolis–Hastings algorithm. Our simulation study shows that the proposed method is more powerful than competing non‐spatial methods. We also apply our method to investigate the effect of cocaine use on brain microstructure. By extending spatial statistics to matrix‐variate data, we contribute to providing a novel and computationally tractable inferential tool for DTI analysis.     

Spatial prediction and temporal backcasting for environmental fields having monotone data patterns

The authors develop a methodology for predicting unobserved values in a conditionally lognormal random spatial field like those commonly encountered in environmental risk analysis. These unobserved values are of two types. The first come from spatial locations where the field has never been monitored, the second, from currently monitored sites which have been only recently installed. Thus the monitoring data exhibit a monotone pattern, resembling a staircase whose highest step comes from the oldest monitoring sites. The authors propose a hierarchical Bayesian approach using the lognormal sampling distribution, in conjunction with a conjugate generalized Wishart distribution. This prior distribution allows different degrees of freedom to be fitted for individual steps, taking into account the differential amounts of information available from sites at the different steps in the staircase. The resulting hierarchical model is a predictive distribution for the unobserved values of the field. The method is demonstrated by application to the ambient ozone field for the southwestern region of British Columbia.     

Moments of the trace of a noncentral wishart matrix

Several articles have appeared on the moments of the trace of a noncentral Wishart matrix. Partial results are given in these papers but the representations of these partial results are often too cumbersome for any practical use. In this paper general results are given in compact form which are readily usable in practical problems.     

On the Dirichlet distributions on symmetric matrices

In this paper, we introduce a generalization of the Dirichlet distribution on symmetric matrices which represents the multivariate version of the Connor and Mosimann generalized real Dirichlet distribution. We establish some properties concerning this generalized distribution. We also extend to the matrix Dirichlet distribution a remarkable characterization established in the real case by Darroch and Ratcliff.     

Polynomial estimation of eigenvalues

In estimating the eigenvalues of the covariance matrix of a multivariate normal population, the usual estimates are the eigenvalues of the sample covariance matrix. It is well known that these estimates are biased. This paper investigates obtaining improved eigenvalue estimates through improved estimates of the characteristic polynomial, which is a function of the sample eigenvalues. A numerical study investigates the improvements evaluated under both a square error and an entropy loss function.     

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).     

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.