二項分布のポアソン近似

数理統計

二項分布のポアソン近似

二項分布のポアソン近似

np=\lambdaのもとで、n \to \infty, \quad p \to 0とすると二項分布はポアソン分布に収束する。

証明

積率母関数で証明する

現代数理統計学の基礎では積率母関数で証明する方法が紹介されています。この方法が一番単純なのでまず紹介します。

n,pが変数なので、p=\frac{\lambda}{n}とし変数をnだけにします。

    \[(pe^t + 1 -p )^n = (\frac{\lambda}{n}e^t+1 -\frac{\lambda}{n})^n = (1+\frac{\lambda(e^t-1)}{n})^n\]

n \to \inftyにして考えると、eの定義から下のように収束します。

    \[\lim_{n \to \infty}(1+\frac{\lambda(e^t-1)}{n})^n = e^{\lambda(e^t-1)}\]

これはポアソン分布の積率母関数と一致するので、ポアソン分布に収束します。

確率関数で証明する

積率母関数で考える方法が一番簡単ではありますが、確率関数で考える方法もあります。データ解析のための数理統計入門ではこの方法が紹介されています。

二項分布の確率関数を下のように変形します。

    \[p(k)=\frac{n!}{k!(n-k)!} (\frac{\lambda}{n})^k(1-\frac{\lambda}{n})^{n-k}=\frac{\lambda^k}{k!}(1-\frac{\lambda}{n})^n \frac{n!}{(n-k)!n^k}(1-\frac{\lambda}{n})^{-k}\]

n \to \inftyのとき、それぞれ4つがどうなるかというと、

\frac{\lambda^k}{k!}nが含まれていないので\frac{\lambda^k}{k!}はそのままになります。

(1-\frac{\lambda}{n})^neの定義からもわかる通り、e^{-\lambda} に収束します。

\frac{n!}{(n-k)!n^k}は下のように変形できます。

    $$\begin{align} \frac{n!}{(n-k)!n^k} &= \frac{n \cdot (n-1) \cdot (n-2) \cdots (n-k+1)}{n \cdot n \cdot n \cdots n}\\ &= \frac{n}{n} \frac{n-1}{n} \frac{n-2}{n} \cdots \frac{n-k+1}{n} \\ &= 1 (1-\frac{1}{n})( 1-\frac{2}{n}) \cdots (1-\frac{k-1}{n})\end{align}$$

それぞれ1に収束するので、積も1に収束する。

(1-\frac{\lambda}{n})^{-k}\frac{\lambda}{n}は0に収束するので、1の-k乗と考えることができるので、1です。

よって確率関数は

    \[\lim_{n \to \infty} p(k) = \lim_{n \to \infty} \frac{\lambda^k}{k!}(1-\frac{\lambda}{n})^n \frac{n!}{(n-k)!n^k}(1-\frac{\lambda}{n})^{-k} = \frac{\lambda^k}{k!} e^{-\lambda}\]

このように収束し、この確率関数はポアソン分布と一致します。

タイトルとURLをコピーしました