The distribution of the ratio of independent central wishart determinants

Ratios of independent central Wishart determinants are useful statistics in multivariate analyses, particularly in the study of multivariate linear models. A method based on the inversion of characteristic functions is outlined for deriving new experessions for the probability distribution functions of the logarithms of these statistics. Accurate tables of the percentiles of these distributions have been obtained covering many bivariate and trivariate cases which have been computed by approximating these expression.     

An approximation for the hotelling-lawley trace in the noncentral case

A procedure for estimating power in conjunction with the Hotelling-Lawley trace is developed. By approximating a non-central Wishart distribution with a central Wishart, and using McKeon's (1974) F-type approximation, a relatively simple procedure for obtaining power estimates is obtained. The accuracy of the approximation is investigated by comparing the approximate results with those for a wide range of conditions given in Olson's (1973) extensive Monte Carlo study. Siotani's (1971) asymptotic expansion is used to provide further comparative assessments. It is demonstrated that the approximation is of sufficient accuracy to be used in practical applications.     

Wishart and Pseudo-Wishart distributions under elliptical laws and related distributions in the shape theory context

In this paper, we determine the density of a singular elliptically contoured matrix. From this, the study of Wishart and Pseudo-Wishart distributions, whether central or noncentral, whether singular or nonsingular, is extended to the case of elliptical models. Some elated distributions are studied in the context of shape theory. Particular attention is paid to singular size-and-shape and size-and-shape cone densities.     

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.     

Moments of generalized quadratic forms and distributions of certain multivariate test statistics

Summary Moments and distributions of quadratic forms or quadratic expressions in normal variables are available in literature. Such quadratic expressions are shown to be equivalent to a linear function of independent central or noncentral chi-square variables. Some results on linear functions of generalized quadratic forms are also available in literature. Here we consider an arbitrary linear function of matrix-variate gamma variables. Moments of the determinant of such a linear function are evaluated when the matrix-variate gammas are independently distributed. By using these results, arbitrary non-null moments as well as the non-null distribution of the likelihood ratio criterion for testing the hypothesis of equality of covariance matrices in independent multivariate normal populations are derived. As a related result, the distribution of a linear function of independent matrix-variate gamma random variables, which includes linear functions of independent Wishart matrices, is also obtained. Some properties of generalized special functions of several matrix arguments are used in deriving these results.     

Saddlepoint Approximations to the Density and the Distribution Functions of Linear Combinations of Ratios of Partial Sums

A class of ratios of partial sums, including Normal, Weibull, Gamma, and Exponential distributions, is considered. The distribution of a linear combination of ratios of partial sums from this class is characterized by the distribution of a linear combination of Dirichlet components. This article presents two saddlepoint approaches to calculate the density and the distribution function for such a class of linear combinations. A simulation study is conducted to assess the performance of the saddlepoint methods and shows the great accuracy of the approximations over the usual asymptotic approximation. Applications of the presented approximations in statistical inferences are discussed.     

Some decompositions of spherical distributions

Certain rotationally symmetric or spherical multivariate distributions can be factored into independent components that produce in a natural manner various uniform, chi, triangular chi, root F, triangular root F, t, disguised t, triangular root beta, Wishart and disguised Wishart distributions.     

The noncentral bivariate chisquared distribution and extensions

The distribution of certain correlated noncentral chisquared variates P, Q, is termed the noncentral bivariate chisquared distribution. Moment generating functions of the distributions of (P, Q), (P+Q) and other quadratic forms have been obtained. A relationship to the linear case of the noncentral Wishart distribution is indicated. Convolution properties and applications are presented.     

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.     

On the real representation of quaternion random variables

For the first time, the matrix-variate quaternion normal and quaternion Wishart distributions are derived from first principles, that is, from their real counterparts, exposing the relations between their respective densities and characteristic functions. Applications of this theory in hypothesis testing are presented, and the density function of Wilks’ statistic is derived for quaternion Wishart matrices.     

Linear generalizations of various matric-t distributions

The concept of a matric-t variate is extended to cases where the positive (definite) part of the variate, which is usually Wishart distributed independently of the normal part, is a linear sum of positive (definite) variates with positive coefficients. These distributions and their quadratic forms are of importance i.a, for the exact solution to the multi¬variate Behrens-Fisher problem. A few useful identities con¬cerning the invariant polynomials with matrix arguments are derived     

p-Value for two-Sample Linear Rank Tests Using Permutation Simulations and the Normal Approximation Method

Test statistics from the class of two-sample linear rank tests are commonly used to compare a treatment group with a control group. Two independent random samples of sizes m and n are drawn from two populations. As a result, N = m + n observations in total are obtained. The aim is to test the null hypothesis of identical distributions. The alternative hypothesis is that the populations are of the same form but with a different measure of central tendency. This article examines mid p-values from the null permutation distributions of tests based on the class of two-sample linear rank statistics. The results obtained indicate that normal approximation-based computations are very close to the permutation simulations, and they provide p-values that are close to the exact mid p-values for all practical purposes.     

A Hybrid Approximation Bayesian Test of Variance Components for Longitudinal Data

The test of variance components of possibly correlated random effects in generalized linear mixed models (GLMMs) can be used to examine if there exists heterogeneous effects. The Bayesian test with Bayes factors offers a flexible method. In this article, we focus on the performance of Bayesian tests under three reference priors and a conjugate prior: an approximate uniform shrinkage prior, modified approximate Jeffreys' prior, half-normal unit information prior and Wishart prior. To compute Bayes factors, we propose a hybrid approximation approach combining a simulated version of Laplace's method and importance sampling techniques to test the variance components in GLMMs.     

Regression via order statistics and their concomitants

We consider a five-dimensional normal distribution and derive the exact joint distribution one variable, linear combinations of order statistics from two other variables, and linear combinations of the corresponding concomitants of these order statistics. We show that this joint distribution is a mixture of trivariate unified skew-normal distributions. This mixture representation enables us to predict one variable based on linear combinations of order statistics from two other variables and linear combinations of the corresponding concomitants. We finally illustrate the usefulness of these results by using a real data.     

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.     

Generalized Inverse Gaussian Distributions and their Wishart Connections

The matrix generalized inverse Gaussian distribution (MGIG) is shown to arise as a conditional distribution of components of a Wishart distributio n. In the special scalar case, the characterization refers to members of the class of generalized inverse Gaussian distributions (GIGs) and includes the inverse Gaussian distribution among others     

A probabilistic approach to Wallis’ formula

We present an alternative proof of Wallis’ formula from the probabilistic point of view. Based on the classical central limit theorem, some discrete distributions with additive property, such as binomial, negative binomial, Poisson and multinomial distributions, are considered to derive π/2.     

On characterizing the normal and poisson distributions

Multivariate normal, correlated multivariate Poisson and multiple Poisson distributions are characterized, in the class of exponential-type distributions, by the properties of the linear combinations of the variables, the properties of their cumulants and the recurance relation between the cumulants.     

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.     

Multivariate stochastic volatility with Bayesian dynamic linear models

This paper develops a Bayesian procedure for estimation and forecasting of the volatility of multivariate time series. The foundation of this work is the matrix-variate dynamic linear model, for the volatility of which we adopt a multiplicative stochastic evolution, using Wishart and singular multivariate beta distributions. A diagonal matrix of discount factors is employed in order to discount the variances element by element and therefore allowing a flexible and pragmatic variance modelling approach. Diagnostic tests and sequential model monitoring are discussed in some detail. The proposed estimation theory is applied to a four-dimensional time series, comprising spot prices of aluminium, copper, lead and zinc of the London metal exchange. The empirical findings suggest that the proposed Bayesian procedure can be effectively applied to financial data, overcoming many of the disadvantages of existing volatility models.