References for the Holonomic Gradient Method (HGM) and the Holonomic Gradient Descent Method (HGD)

Papers and Tutorials

  1. J.Hayakawa, A.Takemura, Estimation of exponential-polynomial distribution by holonomic gradient descent arxiv:1403.7852
  2. C.Siriteanu, A.Takemura, S.Kuriki, MIMO Zero-Forcing Detection Performance Evaluation by Holonomic Gradient Method arxiv:1403.3788
  3. T.Koyama, Holonomic Modules Associated with Multivariate Normal Probabilities of Polyhedra, arxiv:1311.6905
  4. T.Hibi, K.Nishiyama, N.Takayama, Pfaffian Systems of A-Hypergeometric Equations I, Bases of Twisted Cohomology Groups, arxiv:1212.6103 (major revision v2 of arxiv:1212.6103)
  5. Intro T.Hibi et al, Groebner Bases : Statistics and Software Systems , Springer, 2013.
  6. Intro Introduction to the Holonomic Gradient Method (movie), 2013. movie at youtube
  7. T.Sei, A.Kume, Calculating the Normalising Constant of the Bingham Distribution on the Sphere using the Holonomic Gradient Method, Statistics and Computing, 2013, DOI
  8. T.Koyama, A.Takemura, Calculation of Orthant Probabilities by the Holonomic Gradient Method, arxiv:1211.6822
  9. T. Koyama, H. Nakayama, K. Nishiyama, N. Takayama, Holonomic Rank of the Fisher-Bingham System of Differential Equations, to appear in Journal of Pure and Applied Algebra
  10. T. Koyama, H. Nakayama, K. Nishiyama, N. Takayama, Holonomic Gradient Descent for the Fisher-Bingham Distribution on the d-dimensional Sphere, Computational Statistics (2013) DOI
  11. Hiroki Hashiguchi, Yasuhide Numata, Nobuki Takayama, Akimichi Takemura, Holonomic gradient method for the distribution function of the largest root of a Wishart matrix, Journal of Multivariate Analysis, 117, (2013) 296-312, DOI
  12. Tomonari Sei, Hiroki Shibata, Akimichi Takemura, Katsuyoshi Ohara, Nobuki Takayama, Properties and applications of Fisher distribution on the rotation group, Journal of Multivariate Analysis, 116 (2013), 440--455, DOI
  13. T.Koyama, A Holonomic Ideal which Annihilates the Fisher-Bingham Integral, Funkcialaj Ekvacioj 56 (2013), 51--61. jstage
  14. Hiromasa Nakayama, Kenta Nishiyama, Masayuki Noro, Katsuyoshi Ohara, Tomonari Sei, Nobuki Takayama, Akimichi Takemura , Holonomic Gradient Descent and its Application to Fisher-Bingham Integral, Advances in Applied Mathematics 47 (2011), 639--658, DOI

Three Steps of HGM

  1. Find a holonomic system satisfied by the normalizing constant. We may use computational or theoretical methods to find it. Groebner basis and related methods are used.
  2. Find an initial value vector for the holonomic system. This is equivalent to evaluating the normalizing constant and its derivatives at a point. This step is usually performed by a series expansion.
  3. Solve the holonomic system numerically. We use several methods in numerical analysis such as the Runge-Kutta method of solving ordinary differential equations and efficient solvers of systems of linear equations.

Software Packages for HGM

Most software packages are experimental and temporary documents are found in "asir-contrib manual" (auto-autogenerated part), or "Experimental Functions in Asir", or "miscellaneous and other documents" of the OpenXM documents or in this folder. The nightly snapshot of the asir-contrib can be found in the Asir-Contrib page below, or look up our cvsweb page.
  1. Command line interfaces are in the folder OpenXM/src/hgm in the OpenXM source tree. See OpenXM distribution page .
  2. hgm package for R for the step 3.
  3. yang (for Pfaffian systems) , nk_restriction (for D-module integrations), tk_jack (for Jack polynomials), ko_fb_pfaffian (Pfaffian system for the Fisher-Bingham system) are for the steps 1 or 2 and in the asir-contrib.
  4. nk_fb_gen_c is a package to generate a C program to perform maximal Likehood estimates for the Fisher-Bingham distribution by HGD (holonomic gradient descent). It is in the asir-contrib.

Programs to try examples of our papers

  1. d-dimensional Fisher-Bingham System 
 $OpenXM: OpenXM/src/hgm/doc/ref-hgm.html,v 1.9 2014/05/15 07:34:05 takayama Exp $