ポアソン二項分布

ポアソン二項分布

確率質量関数
累積分布関数
母数	[math]\mathbf{p}\in [0,1]^n[/math] - n 回試行のそれぞれに対する成功確率
台	[math] k \in \{0, \ldots, n \}[/math]
確率質量関数	[math]\sum\limits_{A\in F_k} \prod\limits_{i\in A} p_i \prod\limits_{j\in A^c} (1-p_j)[/math]
累積分布関数	[math]\sum\limits_{l=0}^k \sum\limits_{A\in F_l} \prod\limits_{i\in A} p_i \prod\limits_{j\in A^c}{(1-p_j)}[/math]
期待値	[math]\sum\limits_{i=1}^n p_i[/math]
分散	[math] \sigma^2 =\sum\limits_{i=1}^n (1 - {p_i}){p_i}[/math]
歪度	[math]\frac{1}{\sigma^3}\sum\limits_{i=1}^n {\left( 1-2{p_i} \right)\left( 1-{{p}_{i}} \right){{p}_{i}}}[/math]
尖度	[math]\frac{1}{\sigma^4}\sum\limits_{i=1}^n {\left( 1 - 6(1 - p_i){p_i} \right)\left( 1 - p_i \right)p_i}[/math]
モーメント母関数	[math]\prod\limits_{i=1}^n (1-{p_i}+{p_i}{e^{t}})[/math]
特性関数	= [math]\prod\limits_{i=1}^n (1-{p_i}+{p_i}{e^{it}})[/math]

統計学および確率論においてポアソン二項分布 (英: Poisson binomial distribution) とは、独立なベルヌーイ試行の和として定義される離散確率分布である。

別の言い方をすれば、これは成功確率がそれぞれ[math]p_1, p_2, \dots , p_n[/math]であるようなそれぞれ独立なn回の試行を行ったときの成功回数の離散確率分布である。

ここで、すべての成功確率が同じ値[math]p_1 = p_2 = \dots = p_n[/math]とすると、ポアソン二項分布は普通の二項分布になる。すなわち二項分布はポアソン二項分布の特別な場合である。

確率質量関数

n個の確率変数[math]X_i,i \in \{1,2,\ldots,n\}[/math]は、それぞれ独立で成功確率がそれぞれ[math]p_1, p_2, \dots , p_n[/math]であるようなベルヌーイ試行とする。すなわち、

[math]X_i \in \{1,0\} ,\qquad P(X_i=1)=p_i , \qquad P(X_i=0)=1-p_i [/math]

とする。確率変数[math]X= \sum_{i=1}^n X_i[/math]は、このようなn回の試行のうちで成功した回数をあらわす確率変数である。k回成功する確率は次のような和で表現される。

[math]Pr(X=k) = \sum\limits_{A\in {{F}_{k}}}{\prod\limits_{i\in A}{{{p}_{i}}}\prod\limits_{j\in {{A}^{c}}}{(1-{{p}_{j}})}}[/math]

ただし、[math]\ F_k\ [/math]は [math]\{1,2,3, \ldots, n \}[/math]から選べるすべてのk要素部分集合の族である。例えば[math]n= 3[/math]なら、[math]{F_2}=\left\{ \{1,2\},\{1,3\},\{2,3\} \right\}[/math]である。また[math]\ A^c\ [/math]は[math]\ A\ [/math]の補集合。すなわち[math]A^c =\{1,2,3,\dots,n\}\backslash A[/math]である。

これが、定義から直接導かれるポアソン二項分布の確率質量関数である。[math]\ F_k\ [/math]は[math]\ n!/(n-k)!k!\ [/math]要素を含み、この数はnとともに急速に増大するため、試行回数nが小さい場合以外は実際にこの和を計算することは困難である。(例えば[math]\ n=30\ [/math] のとき[math]\ F_{15}\ [/math]は[math]\ 10^{20}\ [/math]もの要素を含む)。 幸いにも、[math]\ Pr(X=k)\ [/math]を計算する非常に効果的な方法がある。一回も成功しない確率が分かれば、n回成功の確率は次のようにして再帰的に計算できる。

[math]\Pr (X=k)=\left\{ \begin{align} & \prod\limits_{i=1}^{n}{(1-{{p}_{i}})}, \qquad k=0 \\ & \frac{1}{k}\sum\limits_{i=1}^{k}{{{(-1)}^{i-1}}\Pr (X=k-i)T(i)}, \qquad k\gt 0 \\ \end{align} \right.[/math]

ただし、[math]T(i)=\sum\limits_{j=1}^{n}{{{\left( \frac{{{p}_{j}}}{1-{{p}_{j}}} \right)}^{i}}}[/math].

他にも離散フーリエ変換を使う次のような計算も可能である。

[math]\Pr (X=k)=\frac{1}{n+1}\sum\limits_{l=0}^{n}{{{C}^{lk}}\prod\limits_{m=1}^{n}{\left( 1+({{C}^{l}}-1){{p}_{m}} \right)}}[/math]

ただし、[math]C=\exp \left( -\frac{2i\pi }{n+1} \right)[/math]である。

さらに他の方法も提案されている。

平均と分散

ポアソン二項分布は独立なベルヌーイ分布に従うn個の確率変数の和だから、その平均と分散は各ベルヌーイ分布における平均および分散の和となる。

[math]\mu = \sum\limits_{i=1}^n p_i[/math]

[math]\sigma^2 =\sum\limits_{i=1}^n (1-p_i) p_i[/math]

レ・カムの定理

次の定理がルーシェン・レ・カム(Lucien le Cam)によって示された。

次のように仮定する

• [math]X_1,\ldots,X_n[/math]はそれぞれベルヌーイ分布に従う独立な確率変数とする。（すなわち0か1の値をとる）ただしそれぞれが同一の分布である必要はない。（発生確率がそれぞれ異なっていてもよい） 各[math] X_i,\ i = 1, 2, 3, \ldots, n[/math]に対して、[math]\ Pr(X_i= 1) = p_i\ [/math]とする。

• [math]\lambda_n = p_1 + \cdots + p_n.\,[/math]

• [math]S_n = X_1 + \cdots + X_n.\,[/math] (すなわち[math]\ S_n\ [/math]はポアソン二項分布に従う。)

このとき、

[math]\sum_{k=0}^\infty \left| \Pr(S_n=k) - {\lambda_n^k e^{-\lambda_n} \over k!} \right| \lt 2 \sum_{i=1}^n p_i^2. [/math]

換言すれば、この和はポアソン分布で近似できる。

各分布がすべて同じ値[math]\ p_i = {\lambda_n}/{n} \ [/math]とすれば、右辺は[math] 2{\lambda_n^2}/{n} [/math]となる。すなわち、この定理は、二項分布の極限がポアソン分布になるというポアソンの極限定理の一般化である。

関連項目

• 二項分布
• ポアソン分布

参考文献

• Steele, J. Michael (1994), “Le Cam's Inequality and Poisson Approximations” (PDF), The American Mathematical Monthly 101 (1): 48-54, doi:10.2307/2325124, http://www-stat.wharton.upenn.edu/~steele/Papers/PDF/LIaPA.pdf . 2010閲覧.