## hgm R

 swMATH ID: 8770 Software Authors: Nakayama, Hiromasa; Nishiyama, Kenta; Noro, Masayuki; Ohara, Katsuyoshi; Sei, Tomonari; Takayama, Nobuki; Takemura, Akimichi Description: Software Packages for Holonomic Gradient Method. The numerical evaluation of the normalizing constant for a given statistical distribution is a fundamental problem in statistics. For example, the normalizing constant of the Gaussian distribution is expressed in terms of a rational expression of a parameter of the distribution named as the standard deviation. However, normalizing constants of many interesting stasistical distributions do not have such closed expressions. The holonomic gradient method, HGM in short, is a general method to evaluate normalizing constant numerically for several parameters in the framework of Zeilberger’s holonomic systems approach. In fact, broad classes of normalizing constants are holonomic functions with respect to parameters. Then, such normalizing constants satisfy holonomic systems of linear partial differential equations. The HGM consists of three steps for a given normalizing constant. (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. The HGM was proposed in 2011 by a group of people inclusing us and has given several new results. For example, the orthant probability is the normalizing constant of the multivariate normal distribution restricted to the first orthant. The HGM can evaluate it in a high accuracy up to the 20 dimensional case when the mean vector is near the origin. In the 20 dimensional case, we numerically solve ordinary differential equation of rank 2^20 =20,148,576. We have developed software packages for the HGM. Packages based on computer algebra systems help us to solve steps (1) and (2). We have implemeted the step (3) for the Fisher-Bingham distribution, the Bingham distribution, the orthant probability, the Fisher distribution on SO(3), some of A-distributions, and the distribution function of the largest root of a Wishart matrix in the language C and/or in the system for statistics R. An implementation for the polyhedral probability is a project in progress. We find an interesting interplay with systems for polytopes in the project. References and current implementations can be found in http://www.math.kobe-u.ac.jp/OpenXM/Math/hgm/ref-hgm.html Homepage: http://www.math.kobe-u.ac.jp/OpenXM/Math/hgm/ref-hgm.html Related Software: Risa/Asir; HolonomicFunctions; SINGULAR; Macaulay2; R; dmod.lib; hgm; SageMath; DLMF; DFVLR-SQP; LBFGS-B; Python; SciPy; D-modules; NAG4M2; Maple; AlgebraicOptimization; crntwin; Julia; ore_algebra Cited in: 26 Documents
all top 5

### Cited by 42 Authors

 8 Sei, Tomonari 8 Takayama, Nobuki 7 Takemura, Akimichi 6 Koyama, Tamio 4 Nakayama, Hiromasa 4 Nishiyama, Kenta 3 Ohara, Katsuyoshi 2 Kume, Alfred 2 Sattelberger, Anna-Laura 1 Betsch, Steffen 1 Breiding, Paul 1 Çelik, Türkü Özlüm 1 Duff, Timothy 1 Ebner, Bruno 1 Görlach, Paul 1 Goto, Yoshiaki 1 Härkönen, Marc 1 Hashiguchi, Hiroki 1 Heaton, Alexander 1 Hibi, Takayuki 1 Hirose, Yoshihiro 1 Jiu, Lin 1 Klar, Bernhard 1 Komaki, Fumiyasu 1 Koutschan, Christoph 1 Kuriki, Satoshi 1 Lehn, Christian 1 Lőrincz, András Cristian 1 Maraj, Aida 1 Marumo, Naoki 1 Matsubara-Heo, Saiei-Jaeyeong 1 Matsumoto, Keiji 1 Nakamoto, Kazuki 1 Noro, Masayuki 1 Oaku, Toshinori 1 Ogawa, Mitsunori 1 Shibata, Hiroki 1 Shimizu, Kunio 1 Takasu, Yuya 1 Venturello, Lorenzo 1 Yano, Keisuke 1 Yürük, Oğuzhan
all top 5

### Cited in 19 Serials

 2 Journal of Multivariate Analysis 2 Advances in Applied Mathematics 2 Japan Journal of Industrial and Applied Mathematics 2 Computational Statistics 2 Statistics and Computing 2 Journal of the Japan Statistical Society. Japanese Issue 1 The Canadian Journal of Statistics 1 Beiträge zur Algebra und Geometrie 1 Advances in Mathematics 1 Funkcialaj Ekvacioj. Serio Internacia 1 Journal of the Mathematical Society of Japan 1 Journal of Pure and Applied Algebra 1 Statistics & Probability Letters 1 Communications in Statistics. Simulation and Computation 1 Kyushu Journal of Mathematics 1 Selecta Mathematica. New Series 1 Mathematics in Computer Science 1 Numerical Algebra, Control and Optimization 1 Information Geometry
all top 5

### Cited in 14 Fields

 18 Statistics (62-XX) 8 Special functions (33-XX) 6 Several complex variables and analytic spaces (32-XX) 5 Commutative algebra (13-XX) 5 Associative rings and algebras (16-XX) 3 Algebraic geometry (14-XX) 3 Operations research, mathematical programming (90-XX) 1 Combinatorics (05-XX) 1 Number theory (11-XX) 1 Linear and multilinear algebra; matrix theory (15-XX) 1 Ordinary differential equations (34-XX) 1 Partial differential equations (35-XX) 1 Probability theory and stochastic processes (60-XX) 1 Numerical analysis (65-XX)