ポアソンの極限定理

ポアソンの極限定理は証明も含めて教科書に載るくらい有名だが,教科書の証明を読んでいてあれ?と思ったので少し考えてみた.

よくある証明

\lambda を適当な正の定数とする. X_n を二項分布 \text{Bn}(n, p) に従う確率変数とし,p = \lambda / nn \rightarrow \infty の極限を考える. このとき,

    \begin{eqnarray*} \mathbb{P}\{X_n = k\}  &=& \binom{n}{k} p^k (1 - p)^{n - k} \\ &=& \frac{n!}{k! (n - k)!} \left( \frac{\lambda}{n} \right)^k \left(1 - \frac{\lambda}{n}\right)^{n - k} \\ &=& \frac{\lambda^k}{k!} \frac{ \left(1 - \frac{\lambda}{n}\right)^n }{ \left(1 - \frac{\lambda}{n}\right)^k } \prod_{i = 0}^{k - 1} \left(1 - \frac{i}{n}\right). \end{eqnarray*}

ここで,

(1)   \begin{eqnarray*} \textstyle \left(1 - \frac{\lambda}{n}\right)^n &\rightarrow& e^{-\lambda} \\ \textstyle \left(1 - \frac{\lambda}{n}\right)^k &\rightarrow& 1 \\ \textstyle \prod_{i = 0}^{k - 1} \left(1 - \frac{i}{n}\right) &\rightarrow& 1   \end{eqnarray*}

であるから

    \[ \mathbb{P} \{X_n = k\} \longrightarrow \frac{\lambda^k}{k!} e^{-\lambda} \]

となる.

気になるところ

(1) では,暗に n \gg k を仮定している. 実際,各点収束をいう場合には,k を止めた後に n を決めればよいので n \gg k とすることができ,上の証明でまったく問題ない.

しかし,k = n として n \rightarrow \infty とすると,(1) は成り立たない. つまり,二項分布の後ろの方の項(少なくとも第 n 項 = 最後の項)は,(1) だけでは収束がいえていないことになる.

実際には,後ろの方の項が収束することは,(1) ではなく,p^k が 0 になることに依る. したがって,真面目に証明するには,k の大小で場合分けして,それぞれが別の理由で収束することを示さないといけない.

ついでに,k = 0 のとき \prod_{i=0}^{k-1} が定義できないという上の証明の小さなほころびも直している.

もう少し真面目な証明

k について場合分けを行う.

まず k = 0 のときは,

    \[ \mathbb{P}\{X = 0\} = (1 - p)^n \rightarrow e^{-\lambda} = \frac{\lambda^0}{0!}e^{-\lambda} \]

であるのでOK.

次に,1 \leq k \leq \sqrt{n} の場合を考える. 1 - \frac{\lambda}{n} < 1 であるから

    \begin{eqnarray*} \left(1 - \frac{\lambda}{n}\right)^k &<& \left(1 - \frac{\lambda}{n}\right)^0 = 1 \\ \left(1 - \frac{\lambda}{n}\right)^k &\geq& \left(1 - \frac{\lambda}{n}\right)^{\sqrt{n}} \rightarrow 1. \end{eqnarray*}

よってはさみうちの原理より

    \[ \left(1 - \frac{\lambda}{n}\right)^k \longrightarrow 1. \]

また,0 \leq i \leq k に対して 1 - \frac{i}{n} \leq 1,かつ

    \[ 1 - \frac{i}{n} \geq 1 - \frac{\sqrt{n}}{n} = 1 - \frac{1}{\sqrt{n}} \rightarrow 1 \]

であるから,はさみうちの原理によって

    \[ 1 - \frac{i}{n} \longrightarrow 1. \]

よって

    \[ \prod_{i = 0}^k \left( 1 - \frac{i}{n} \right) \longrightarrow 1 \]

もいえる.

以上から (1) がいえるので,1 \leq k \leq \sqrt{n} の場合もOK.

最後に,k > \sqrt{n} の場合を考える. このとき,

    \begin{eqnarray*} \textstyle \left(1 - \frac{\lambda}{n}\right)^k &<& \infty \\ \textstyle \prod_{i = 0}^{k - 1} \left(1 - \frac{i}{n}\right) &<& \infty  \end{eqnarray*}

かつ

    \[ \frac{\lambda^k}{k!} \longrightarrow 0 \]

であるから,

    \[ \mathbb{P} \{X = k\} \longrightarrow 0. \]

\frac{\lambda^k}{k!} e^{-\lambda} \rightarrow 0 であるから,この場合も OK となる.

補足1

\left(1 - \frac{\lambda}{n}\right)^{\sqrt{n}} \rightarrow 1 は自明でないかもしれないので,証明しておく.

m = \sqrt{n} とおくと,

    \begin{eqnarray*} \left(1 - \frac{\lambda}{n}\right)^{\sqrt{n}}  &=& \left(1 - \frac{\lambda}{m^2}\right)^{m} \\ &=& \sum_{i = 0}^m (-1)^i \binom{m}{i} \frac{\lambda^i}{m^{2i}} \\ &=& 1 + \frac{1}{m} \sum_{i = 1}^m (-1)^i \binom{m}{i} \frac{\lambda^i}{m^{2i - 1}} \end{eqnarray*}

ここで,a_i^{(m)} = (-1)^i \binom{m}{i} \frac{\lambda^i}{m^{2i - 1}} とすると,

    \[  |a_i^{(m)}|  = \binom{m}{i} \frac{\lambda^i}{m^{2i - 1}}  \leq m^i \cdot \frac{\lambda^i}{m^{2i - 1}}  = \lambda \left(\frac{\lambda}{m}\right)^{i - 1}. \]

よって,m \rightarrow \infty のとき \sum_{i = 1}^m a_i^{(m)} は絶対収束する. したがって,

    \[ \frac{1}{m} \sum_{i = 1}^m (-1)^i \binom{m}{i} \frac{\lambda^i}{m^{2i - 1}} \longrightarrow 0. \]

補足2

今回 kn を同時に飛ばしたが, k > n のとき \mathbb{P}\{X = k\} = 0 として k \rightarrow \infty としてから n \rightarrow \infty すれば, 今回言いたかったことが簡単に言えたんじゃないかと終わってみてから気づいた.

コメントを残す

メールアドレスが公開されることはありません。