The holonomic gradient method for the distribution function of the largest root of a Wishart matrix

Authors:
Hiroki Hashiguchi;Yasuhide Numata;Nobuki Takayama;Akimichi Takemura
Affiliations:
-;-;-;-
Venue:
Journal of Multivariate Analysis
Year:
2013

Citing 5
Cited 2

A holonomic systems approach to special functions identities

Journal of Computational and Applied Mathematics
Regular Article: Algorithms forb-Functions, Restrictions, and Algebraic Local Cohomology Groups ofD-Modules

Advances in Applied Mathematics
Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations

Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations
Distribution of eigenvalues and eigenvectors of Wishart matrix when the population eigenvalues are infinitely dispersed and its application to minimax estimation of covariance matrix

Journal of Multivariate Analysis
Exact distributional computations for Roy's statistic and the largest eigenvalue of a Wishart distribution

Statistics and Computing

Properties and applications of Fisher distribution on the rotation group

Journal of Multivariate Analysis
Approximation of upper percentile points for the second largest latent root in principal component analysis

International Journal of Knowledge Engineering and Soft Data Paradigms

Quantified Score

Hi-index	0.00

Visualization

Abstract

We apply the holonomic gradient method introduced by Nakayama et al. (2011) [23] to the evaluation of the exact distribution function of the largest root of a Wishart matrix, which involves a hypergeometric function "1F"1 of a matrix argument. Numerical evaluation of the hypergeometric function has been one of the longstanding problems in multivariate distribution theory. The holonomic gradient method offers a totally new approach, which is complementary to the infinite series expansion around the origin in terms of zonal polynomials. It allows us to move away from the origin by the use of partial differential equations satisfied by the hypergeometric function. From the numerical viewpoint we show that the method works well up to dimension 10. From the theoretical viewpoint the method offers many challenging problems both to statistics and D-module theory.