[論文レビュー] An asymptotic formula for the number of non-negative integer matrices with prescribed row and column sums

Alexander Barvinok, J. A. Hartigan|arXiv (Cornell University)|Oct 13, 2009

Advanced Combinatorial Mathematics参考文献 14被引用数 8

ひとこと要約

本稿は、広範な滑らかなマージンのクラスに対して有効な、与えられた行和・列和を持つ $m \times n$ の非負整数行列の個数を計算的に効率的な漸近公式で与える。この公式は、$m$ および $n$ が増加するにつれて相対誤差が消えるように設計されており、エントロピー最大化から得られる典型行列と、エッジワース補正を施したガウス近似を活用している。

ABSTRACT

We count mxn non-negative integer matrices (contingency tables) with prescribed row and column sums (margins). For a wide class of smooth margins we establish a computationally efficient asymptotic formula approximating the number of matrices within a relative error which approaches 0 as m and n grow.

研究の動機と目的

• 広範な条件下で、与えられた行和・列和を持つ非負整数行列（連関表）の個数に対する漸近公式を導出すること。
• 行列の次元が増加するにつれて精度が向上する計算的に効率的な近似手法を確立すること。
• 漸近公式が成り立つ「滑らか」なマージンのクラスを同定し、近似の一様性と安定性を保証すること。
• 漸近的個数とエントロピー最大化による典型行列との関係を明らかにし、それをガウス展開およびエッジワース展開の中心点として用いること。
• 輸送ポリトープの体積近似とスケーリングの議論を用いて、公式の精度を検証すること。

提案手法

• 輸送ポリトープ $P(R,C)$ 上で、厳密に凹関数 $g(X) = \sum_{j,k} \left( (x_{jk}+1)\ln(x_{jk}+1) - x_{jk}\ln x_{jk} \right)$ を最大にする唯一の解 $Z = (\zeta_{jk})$ を典型行列と定義する。
• $\delta$-滑らかなマージンを導入し、次元がバランスしていること（$m, n \geq \delta n, \delta m$）と、典型行列の成分が0および無限大から一様に離れていること（$\delta \tau \leq \zeta_{jk} \leq \tau$）を要件とする。
• スムーズ性の条件のもとで、$\#(R,C)$ の主要項近似として輸送ポリトープの体積 $\operatorname{vol} P(R,C)$ を用い、相対誤差が $O((m+n)^{-1})$ であることを示す。
• マージン $(R,C)$ をスケーリングされたマージン $(\hat{R}, \hat{C})$ に写像するスケーリング変換 $\alpha = (m+n)^9 \tau^{-1}$ を導入し、$\hat{r}_j \approx \alpha r_j$, $\hat{c}_k \approx \alpha c_k$ となるようにする。
• 2次形式 $q(s,t)$ および $\hat{q}(s,t)$ を用いて、漸近公式におけるガウス項を関連付け、$\hat{q} \approx \alpha^2 q$ が相対誤差 $O((m+n)^{-2})$ 以内で成り立つことを示す。
• 6次および4次の同次多項式 $f^2$ および $h$ を用いてエッジワース補正項を組み込み、元のバージョンとスケーリング後のバージョンとの差が制御可能であることを示す。

実験結果

リサーチクエスチョン

• RQ1 $m,n \to \infty$ のとき、計算的に効率的かつ高精度な連関表の個数に対する漸近公式を導出可能か？
• RQ2 行和・列和にどのような条件を課すと、典型行列の成分が0および無限大から一様に離れるようになるか？
• RQ3 エントロピー最大化による典型行列は、固定マージンを持つ整数行列の個数の漸近的カウントとどのように関係するか？
• RQ4 ガウス近似およびエッジワース近似は、典型行列の周辺における連関表の分布をどの程度正確に捉えられるか？
• RQ5 スケーリングパラメータ $\alpha$ を用いて定量化されたマージンの小さな摂動に対しても、漸近公式は安定性を保つか？

主な発見

• マージン $(R,C)$ をもつ $m \times n$ の非負整数行列の個数は、$\delta$-滑らかさの下で、$\operatorname{vol} P(R,C)$ により漸近的に近似可能であり、相対誤差は $O((m+n)^{-1})$ である。
• 漸近公式は、$\alpha = (m+n)^9 \tau^{-1}$ によるマージンのスケーリングに対して不変であり、$\#(R,C) \approx \alpha^{(m-1)(1-n)} \#(\hat{R}, \hat{C})$ と変換される。
• 漸近展開におけるガウス項は、$\exp\left(-\mu/2 + \nu\right) \approx \exp\left(-\hat{\mu}/2 + \hat{\nu}\right)$ と近似可能であり、相対誤差は $O((m+n)^{-1})$ 以内である。
• エッジワース補正項 $f^2$ および $h$ は、元のマージンとスケーリング後のマージンとの間で、相対的な大きさの差が $O((m+n)^{-2})$ 以内である。
• 典型行列の成分 $\zeta_{jk}$ は、$|\hat{\zeta}_{jk} - \alpha \zeta_{jk}| = O(\alpha \zeta_{jk}/(m+n)^2)$ を満たし、スケーリングに対して安定であることが保証される。
• 2次形式 $\psi$ および $\hat{\psi}$ の最小固有値は $\Omega(m+n)$ であり、ガウス近似の有効性を裏付ける。