[論文レビュー] $A$-Hypergeometric Distributions and Newton Polytopes

Nobuki Takayama, Satoshi Kuriki+1 more|arXiv (Cornell University)|Oct 8, 2015

Advanced Statistical Methods and Models参考文献 10被引用数 3

ひとこと要約

本稿では、ニュートン多面体を用いて、$A$-ハイパージェオメトリック分布のパラメータ空間の商空間とモーメント空間の間の全単射を確立し、ホロノミック勾配法（HGM）を用いたアルゴリズム的逆モーメント写像を可能にする。また、モーメント写像と反復比例スケーリング（IPS）の漸近的同等性を証明し、$A$-ハイパージェオメトリック多項式の新たな近似を提供するとともに、統計学における条件付き最尤推定を可能にする。

ABSTRACT

We give a bijection between a quotient space of the parameters and the space of moments for any $A$-hypergeometric distribution. An algorithmic method to compute the inverse image of the map is proposed utilizing the holonomic gradient method and an asymptotic equivalence of the map and the iterative proportional scaling. The algorithm gives a method to solve a conditional maximum likelihood estimation problem in statistics. Our interplay between the theory of hypergeometric functions and statistics gives some new formulas of $A$-hypergeometric polynomials.

研究の動機と目的

• A-ハイパージェオメトリック分布のパラメータ空間とモーメント空間の間の数学的対応関係を確立すること。
• 条件付き最尤推定に不可欠なモーメント写像の逆写像を計算するアルゴリズム的手法を開発すること。
• モーメント写像と反復比例スケーリング（IPS）の間の漸近的同等性を探索し、統計的推定と代数的幾何学を結びつけること。
• ハイパージェオメトリック関数とニュートン多面体の相乗作用を用いて、$A$-ハイパージェオメトリック多項式の新たな公式を導出すること。

提案手法

• パラメータ$Au = \beta$を固定したポissonモデルにおける条件付き分布として$A$-ハイパージェオメトリック分布を定義する。
• 対数パーティション関数$\psi(\xi) = \log Z(\beta; p(\xi))$を用い、モーメントを$E[U_i] = \partial_{\xi_i} \psi(\xi)$として表現し、$\xi$空間から$\eta$空間へのモーメント写像を確立する。
• Aのトーラス作用を用いてパラメータの商空間を導入し、同一のモーメントをもたらす等価な$\xi$値を特定する。
• モーメント写像の像が$Z(\beta; p)$のニュートン多面体の相対内部に一致することを証明し、これにより商空間との全単射を確立する。
• ホロノミック勾配法（HGM）を用いてモーメント写像を数値的に逆写像し、条件付き最尤推定の計算を可能にする。
• HGM反復と反復比例スケーリング（IPS）の間の漸近的同等性を示し、これによりHGMの初期化および正規化定数$Z(\beta; p)$の近似を実現する。

実験結果

リサーチクエスチョン

• RQ1 モーメント写像$\xi \mapsto E[U]$の像は$\eta$空間でどのように特定され、幾何学的にどのように特徴付けられるか？
• RQ2 A-ハイパージェオメトリックモデルにおける統計的推定のため、モーメント写像の逆写像をどのように効率的に計算できるか？
• RQ3 A-ハイパージェオメトリック分布の文脈において、ホロノミック勾配法と反復比例スケーリング（IPS）の関係は何か？
• RQ4 正規化定数$Z(\beta; p)$はIPS手順を用いて漸近的に近似可能か？その誤差挙動はいかなるものか？

主な発見

• モーメント写像$E[U]$の像は、$A\eta = \beta$, $\eta \geq 0$によって定義されるニュートン多面体の相対内部に一致し、モーメント空間の幾何的特徴付けを確立する。
• パラメータの商空間（Aのトーラス作用による同値類）と$\eta$空間におけるニュートン多面体との間に全単射が確立され、情報幾何学における新たな同型写像が得られる。
• ホロノミック勾配法（HGM）は、A-ハイパージェオメトリック分布の条件付き最尤推定を可能にする効果的なアルゴリズム的アプローチを提供する。
• モーメント写像は反復比例スケーリング（IPS）と漸近的に同等であり、スケーリングされた問題$k\beta$に対して$k \to \infty$で同一の極限に収束する。
• 定理6では、IPS解$m$を用いて$Z(k\beta; p)$の漸近的近似が提供され、$k$が増加するにつれて誤差が減少する。$k=9, 200, 300$について数値的に検証された。
• $k=9$では、HGMによる$\log Z$の近似誤差は1.8052であり、$k=300$では0.3734に減少し、正確な値に収束することが示された。