ABSTRACT
The hypergeometric function 1F1 of a matrix argument Y is a symmetric entire function in the eigenvalues y1,...,ym of Y. It appears in the distribution function of the largest eigenvalue of a Wishart matrix and its numerical evaluation is important in multivariate distribution theory. Hashiguchi et al. (J. Multivariate Analysis, 2013) proposed an efficient algorithm for evaluating the matrix 1F1 by the holonomic gradient method (HGM). The algorithm is based on the system of partial differential equations (PDEs) satisfied by the matrix 1F1 given by Muirhead (Ann. Math. Statist., 1970) and it cannot be applied to the diagonal cases, i.e. the cases where several yi's are equal because the system of PDEs has singularities on the diagonal region. Hashiguchi et al. derived an ordinary differential equation (ODE) satisfied by 1F1(y,y) in the bivariate case from some relations which are obtained by applying l'opital rule to the system of PDEs for 1F1(y1,y2). In this paper we generalize this approach for computing systems of PDEs satisfied by the matrix 1F1 for various diagonalization patterns. We show that the existence of a system of PDEs for a diagonalized 1F1 is reduced to the non-singularity of the matrices systematically derived from the diagonalization pattern. By checking the non-singularity numerically, we show that there exists a system of PDEs for a diagonalized 1F1 if the size of each diagonal block ≤ 36. We have computed an ODE for 1F1(y,...,y) up to m=22. We made a test implementation of HGM for diagonal cases and we show some numerical results.
- J. C. Butcher. Numerical Methods for Ordinary Differential Equations, 2nd edition. John Wiley and Sons, 2008.Google Scholar
- L. Fousse, G. Hanrot, V. Lefèvre, P. Pélissier, and P. Zimmermann. MPFR: A multiple-precision binary floating-point library with correct rounding. ACM Trans. Math. Softw, 33:00000818, 2007. Google ScholarDigital Library
- H. Hashiguchi and N. Niki. Numerical computation on distributions of the largest and the smallest latent roots of the Wishart matrix. J. Japanese Soc. Comput. Statist., 19(1):45--56, 2006.Google ScholarCross Ref
- H. Hashiguchi, Y. Numata, N. Takayama, and A. Takemura. Holonomic gradient method for the distribution function of the largest root of a Wishart matrix. J. Multivariate Analysis, 117:296--312, 2013. Google ScholarDigital Library
- P. Koev and A. Edelman. The efficient evaluation of the hypergeometric function of a matrix argument. Math. Comp., 75(254):833--846, 2006.Google ScholarCross Ref
- T. Kondo. On a holonomic system of partial differential equations satisfied by the matrix $_1F_1$ (in Japanese). Master's thesis, Kobe university, 2013.Google Scholar
- R. J. Muirhead. System of partial differential equations for hypergeometric functions of matrix argument. Ann. Math. Statist., 41:991--1001, 1970.Google ScholarCross Ref
- R. J. Muirhead. Latent roots and matrix variates: a review of some asymptotic results. Ann. Statist., 6(1):5--33, 1978.Google ScholarCross Ref
- R. J. Muirhead. Aspects of Multivariate Statistical Theory. John Wiley and Sons, 1982.Google Scholar
- H. Nakayama, K. Nishiyama, M. Noro, K. Ohara, T. Sei, N. Takayama, and A. Takemura. Holonomic gradient descent and its application to the Fisher-Bingham integral. Adv. Appl. Math., 47:639--658, 2011. Google ScholarDigital Library
- M. Noro. a Risa/Asir package for computing PDEs for a matrix _1F_1 on diagonal regions. http://www2.rikkyo.ac.jp/web/noro/wishart.Google Scholar
- \relax Risa/Asir developing team. Risa/Asir, a computer algebra system. http://www.math.kobe-u.ac.jp/Asir.Google Scholar
- Y. Sheena and A. Takemura. Inference on eigenvalues of Wishart distribution using asymptotics with respect to the dispersion of population eigenvalues. The Indian J. Statist., 69:717--733, 2007.Google Scholar
Index Terms
- System of Partial Differential Equations for the Hypergeometric Function 1F1 of a Matrix Argument on Diagonal Regions
Recommendations
Sylvester Equations and the numerical solution of partial fractional differential equations
We develop a new matrix-based approach to the numerical solution of partial differential equations (PDE) and apply it to the numerical solution of partial fractional differential equations (PFDE). The proposed method is to discretize a given PFDE as a ...
Hypergeometric Functions and Parabolic Partial Differential Equations
We investigate the summability of formal power series solutions of parabolic linear partial differential equations with exponential coefficients.
Numerical methods for the solution of partial differential equations of fractional order
Anomalous diffusion is a possible mechanism underlying plasma transport in magnetically confined plasmas. To model this transport mechanism, fractional order space derivative operators can be used. Here, the numerical properties of partial differential ...
Comments