version 1.8, 2014/04/03 07:34:31 |
version 1.10, 2014/05/16 11:30:31 |
Line 88 Advances in Applied Mathematics 47 (2011), 639--658, |
|
Line 88 Advances in Applied Mathematics 47 (2011), 639--658, |
|
|
|
<h2> Three Steps of HGM </h2> |
<h2> Three Steps of HGM </h2> |
<ol> |
<ol> |
<li> Find a holonomic system satisfied by the normalizing constant. |
<li> Finding a holonomic system satisfied by the normalizing constant. |
We may use computational or theoretical methods to find it. |
We may use computational or theoretical methods to find it. |
Groebner basis and related methods are used. |
Groebner basis and related methods are used. |
<li> Find an initial value vector for the holonomic system. |
<li> Finding an initial value vector for the holonomic system. |
This is equivalent to evaluating the normalizing constant and its derivatives |
This is equivalent to evaluating the normalizing constant and its derivatives |
at a point. |
at a point. |
This step is usually performed by a series expansion. |
This step is usually performed by a series expansion. |
<li> Solve the holonomic system numerically. We use several methods |
<li> Solving the holonomic system numerically. We use several methods |
in numerical analysis such as the Runge-Kutta method of solving |
in numerical analysis such as the Runge-Kutta method of solving |
ordinary differential equations and efficient solvers of systems of linear |
ordinary differential equations and efficient solvers of systems of linear |
equations. |
equations. |
|
|
<a href="http://www.math.kobe-u.ac.jp/OpenXM/Current/doc/index-doc.html"> |
<a href="http://www.math.kobe-u.ac.jp/OpenXM/Current/doc/index-doc.html"> |
OpenXM documents</a> |
OpenXM documents</a> |
or in <a href="./"> this folder</a>. |
or in <a href="./"> this folder</a>. |
The nightly snapshot of the asir-contrib can be found in the Asir-Contrib page below, |
The nightly snapshot of the asir-contrib can be found in the asir page below, |
or look up our <a href="http://www.math.sci.kobe-u.ac.jp/cgi/cvsweb.cgi/"> |
or look up our <a href="http://www.math.sci.kobe-u.ac.jp/cgi/cvsweb.cgi/"> |
cvsweb page</a>. |
cvsweb page</a>. |
<ol> |
<ol> |
<li> <a href="http://www.math.kobe-u.ac.jp/OpenXM/Math/hgm/"> hgm package for R </a> for the step 3. |
<li> Command line interfaces are in the folder OpenXM/src/hgm |
<li> yang (for Pfaffian systems) , nk_restriction (for D-module integrations), |
in the OpenXM source tree. See <a href="http://www.math.kobe-u.ac.jp/OpenXM"> |
tk_jack (for Jack polynomials), ko_fb_pfaffian (Pfaffian system for the Fisher-Bingham system) |
OpenXM distribution page </a>. |
are for the steps 1 or 2 and in the |
<li> <a href="http://www.math.kobe-u.ac.jp/OpenXM/Math/hgm/"> hgm package for R </a> (hgm_*tar.gz, hgm-manual.pdf) for the step 3. |
<a href="http://www.math.kobe-u.ac.jp/Asir"> asir-contrib</a>. |
<li> The following packages are |
<li> nk_fb_gen_c is a package to generate a C program to perform |
for the computer algebra system |
|
<a href="http://www.math.kobe-u.ac.jp/Asir"> Risa/Asir</a>. |
|
They are in the asir-contrib collection. |
|
<ul> |
|
<li> yang.rr (for Pfaffian systems) , |
|
nk_restriction.rr (for D-module integrations), |
|
tk_jack.rr (for Jack polynomials), |
|
ko_fb_pfaffian.rr (Pfaffian system for the Fisher-Bingham system), |
|
are for the steps 1 or 2. |
|
<li> nk_fb_gen_c.rr is a package to generate a C program to perform |
maximal Likehood estimates for the Fisher-Bingham distribution by HGD (holonomic gradient descent). |
maximal Likehood estimates for the Fisher-Bingham distribution by HGD (holonomic gradient descent). |
It is in the |
<li> ot_hgm_ahg.rr (HGM for A-distributions, very experimental). |
<a href="http://www.math.kobe-u.ac.jp/Asir"> asir-contrib</a>. |
</ul> |
</ol> |
</ol> |
|
|
<h2> Programs to try examples of our papers </h2> |
<h2> Programs to try examples of our papers </h2> |
|
|
<li> <a href="http://www.math.kobe-u.ac.jp/OpenXM/Math/Fisher-Bingham-2"> d-dimensional Fisher-Bingham System </a> |
<li> <a href="http://www.math.kobe-u.ac.jp/OpenXM/Math/Fisher-Bingham-2"> d-dimensional Fisher-Bingham System </a> |
</ol> |
</ol> |
|
|
<pre> $OpenXM: OpenXM/src/hgm/doc/ref-hgm.html,v 1.7 2014/03/31 07:23:09 takayama Exp $ </pre> |
<pre> $OpenXM: OpenXM/src/hgm/doc/ref-hgm.html,v 1.9 2014/05/15 07:34:05 takayama Exp $ </pre> |
</body> |
</body> |
</html> |
</html> |