[論文レビュー] A full proof of universal inequalities for the distribution function of the binomial law

A. M. Zubkov, Alexander Serov|arXiv (Cornell University)|Jul 16, 2012

Statistical Distribution Estimation and Applications参考文献 4被引用数 4

ひとこと要約

本稿では、標準正規分布の累積分布関数を用いて、二項分布関数 $F_{n,p}(x)$ の普遍的な両側不等式の簡潔で完全な証明を提示する。主な貢献は、アルファースとディンゲスの不等式を、$F_{n,p}(k)$ の上界が $F_{n,p}(k+1)$ の下界と一致する（正確に一致する）形に再定式化したことである。この不等式はすべての $n$、$p$、$k$ に対して有効であり、隣接する整数を用いた正確な分位数推定を可能にする。

ABSTRACT

We present a new form and a short full proof of explicit two-sided estimates for the distribution function F_{n,p}(x) of the binomial law from the paper published by D.Alfers and H.Dinges in 1984. These inequalities are universal (valid for all binomial distribution and all values of argument) and exact (namely, the upper bound for F_{n,p}(k) is the lower bound for F_{n,p}(k+1)). By means of such estimates it is possible to bound any quantile of the binomial law by 2 subsequent integers.

研究の動機と目的

• 普遍的な両側不等式の完全でアクセス可能な証明を提供すること。これらの不等式は以前から知られていたが、長大で不完全な証明のため、検証が困難であった。
• 特に、標準近似が超指数的減衰のために失敗する大偏差領域において、二項分布の尾確率を推定する実用的課題を解決すること。
• 上界 $F_{n,p}(k)$ と下界 $F_{n,p}(k+1)$ が一致するという意味で正確な不等式を確立すること。これにより、連続する整数間の分位数推定がきわめてタイトに保証される。
• これらの不等式が、任意の二項分布の分位数を連続する二つの整数で厳密に挟めることを示し、大偏差問題に対する実用的解決策を提供すること。

提案手法

• 二項分布の累積分布関数の積分表現を、$p$ に関する積分の微分（積分内部微分）を用いて導出し、和をベータ型関数を含む積分に変換する。
• 二項係数に対してスターリングの近似を適用し、正規化係数を $S_n^{k+1}$ として表現する。$S_n^{k+1}$ は対数ガンマ関数の展開を含む補正項である。
• 指数的減衰を支配する関数 $B(z) = \alpha \ln(\alpha/z) + (1-\alpha)\ln((1-\alpha)/(1-z))$ を導入し、$\alpha = (k+1)/n$ とする。
• 被積分関数を標準正規密度 $\varphi$ および累積関数 $\Phi$ を用いて再表現し、二項分布の尾を正規分布に結びつける変換 $a(z)$ を用いる。
• 差関数 $\delta(p) = \mathbf{P}\{X_{n,p} \leq k\} - \Phi(a(p)\sqrt{n})$ の単調性および符号の性質を確立し、すべての $p \in (0,1)$ に対して $\delta(p) < 0$ であることを示す。これにより上界が得られる。
• 対称性と $1-p$ を用いた双対性により、補完分布の上界を用いて下界を導出し、定義域全体にわたるタイトさを保証する。

実験結果

リサーチクエスチョン

• RQ1 アルファースとディンゲスの二項分布関数に関する普遍的両側不等式に対して、簡潔で完全な証明を与えることは可能か？
• RQ2 上界 $F_{n,p}(k)$ と下界 $F_{n,p}(k+1)$ が一致するという正確性が、厳密に確立できるか？
• RQ3 超指数的尾減衰のため標準近似が失敗する大偏差領域において、これらの不等式はどの程度有効に保たれるか？
• RQ4 誤差項 $\delta(p)$ の積分表現を用いて、解析的または数値的に不等式を鋭くできるか？
• RQ5 得られた不等式を用いて、二項分布の任意の分位数を連続する二つの整数で挟むことは可能か？

主な発見

• すべての $n$、$p \in (0,1)$、$k = 0,1,\dots,n-1$ に対して不等式 $C_{n,p}(k) \leq \mathbf{P}\{X_{n,p} \leq k\} \leq C_{n,p}(k+1)$ が成り立ち、等号は $k=0$ または $k=n-1$ のみである。
• 不等式は正確である：$C_{n,p}(k)$ は $F_{n,p}(k)$ の下界であり、$C_{n,p}(k+1)$ は $F_{n,p}(k)$ の上界である。これにより、連続する整数間の遷移においてタイトさが保証される。
• 真の確率と下界との差は、局所確率 $\mathbf{P}\{X_{n,p} = k\}$ よりも厳密に小さい。これにより、不等式の精度が裏付けられる。
• 上界は、すべての $p \in (0,1)$ に対して $\delta(p) < 0$ であることを示すことで得られる。これにより $\mathbf{P}\{X_{n,p} \leq k\} < \Phi(a(p)\sqrt{n}) = C_{n,p}(k+1)$ が導かれる。
• 下界は双対性により得られる：$\mathbf{P}\{X_{n,p} \leq k\} = 1 - \mathbf{P}\{X_{n,1-p} \leq n-k-1\}$ であり、補完分布の上界を適用することで得られる。
• 誤差項 $\delta(p)$ は、被積分関数の成分の単調性を用いて解析的に評価可能であり、特定のパrameter領域において、不等式の鋭さを数値的または解析的に高めることができる。