統計学

ポアソン分布

2021年6月9日


ポアソン分布を学習する前に、相関の強い二項分布についても参照してください。
二項分布

ポアソン分布

ポアソン分布の確率密度関数は以下です。

【ポアソン分布】
\begin{eqnarray}
P(X=k)
&=& \frac{\lambda^{k}e^{-\lambda}}{k!}
\end{eqnarray}

ポアソン分布は滅多に起きない場合の確率密度関数です。
よく題材に上がるのが交通事故の発生確率だったり、さらには隕石が落ちてくる確率とか。
ん?
起きる、起きないを取り扱うということはまさに相反する事象で、ん?二項分布やベルヌーイ分布もそんなようだったような。
そうです!そのため、ポアソン分布は二項分布とはかなり関係のある確率分布となります。
さらに滅多に起きないということは、起きる確率は0にほぼ等しい。起きる確率をpとしたら、pは0にほぼ等しいということになります。
 
確率変数としては、交通事故が起きる回数、隕石が落ちる回数なので、
二項分布と同じように、試行回数を考えるものです。
 
ということなので、二項分布の\(p\)を限りなく0に近いことを仮定することで、
ポアソン分布の確率関数を導出することができます。
滅多に起きないのであれば、試行回数もとにかく多く回数こなさないと1回ですら起きないので、その試行回数を\(n\)とした時、\(n\rightarrow∞\)とします。
そして、起きる確率をpとした時、\(p\rightarrow0\)とします。

\(k\)は二項分布の確率関数をそのまま利用して導出しているので、
二項分布の\(k\)と同じ意味で、試行回数\(n\)回のうち、事象が起きた回数を\(k\)としています。

そこで、ある定数パラメータ\(λ\)を置いて、以下の式を考える。
\begin{eqnarray}
p=\frac{\lambda}{n}
\end{eqnarray}

上記のように式を置くと、
\(n\rightarrow∞\)をとると、\(λ\)は定数なので、右辺は0に近づき、左辺も0に近づくので、仮定の\(p\rightarrow0\)を満たします。
さらに、\(p\rightarrow0\)とすると、右辺も0に近づかないといけないので、\(λ\)が定数なので、\(n\)は\(n\rightarrow∞\)を満たさないといけません。
なので、元々の仮定を満たす式を作りました。

この式を二項分布に代入することで、ポアソン分布を導出することができます。

 

ポアソン分布は二項分布から導出できるので証明してみたいと思います。
\begin{eqnarray}
P(X=k)={}_nC_k p^ k(1-p)^{n-k}
\end{eqnarray}
ここで、\(p= \displaystyle \frac{\lambda}{n}\)と置く。\(\lambda\)はパラメータで定数とします。
そうすると、\(n\rightarrow∞\)とすると、\(\displaystyle \frac{\lambda}{n} \rightarrow 0\)により、\(p\rightarrow0\)になります。
\(p= \displaystyle \frac{\lambda}{n}\)を代入して、

\begin{eqnarray}
P(X=k) &=& {}_nC_k \biggl(\frac{\lambda}{n}\biggl)^ k\biggl(1-\frac{\lambda}{n}\biggl)^{n-k} \\
&=& \frac{n(n-1)(n-2)・・・(n-k-1)}{k!}\biggl(\frac{\lambda}{n}\biggl)^{k} \biggl(1-\frac{\lambda}{n}\biggl)^{n-k} \\
\end{eqnarray}
ここで、分子の各項から\(n\)を抜きます。\(k\)個あるので、\(n \cdot n \cdot \cdot \cdot n = n^{k}\)となり、
\begin{eqnarray}
P(X=k) &=& \frac{n^{k} \cdot 1 \cdot \biggl(1- \displaystyle \frac{1}{n} \biggl) \cdot \biggl(1- \displaystyle \frac{2}{n}\biggl) \cdot \cdot \cdot \biggl(1- \displaystyle \frac{k-1}{n}\biggl)}{k!} \cdot \biggl( \displaystyle \frac{\lambda}{n} \biggl)^{k} \biggl(1- \displaystyle \frac{\lambda}{n} \biggl)^{n-k} \\
&=& \frac{ \biggl(1- \displaystyle \frac{1}{n}\biggl) \cdot \biggl(1- \displaystyle \frac{2}{n}\biggl) \cdot \cdot \cdot \biggl(1- \displaystyle \frac{k-1}{n}\biggl) }{k!} \cdot {\lambda}^{k} \biggl(1- \displaystyle \frac{\lambda}{n} \biggl)^{n-k} \\
&=& \frac{\lambda^{k}}{k!} \cdot \biggl(1- \displaystyle \frac{1}{n} \biggl)\biggl(1- \displaystyle \frac{2}{n}\biggl) \cdot \cdot \cdot \biggl(1- \displaystyle \frac{k-1}{n}\biggl)\biggl(1- \displaystyle \frac{\lambda}{n}\biggl)^{n}\biggl(1- \displaystyle \frac{1}{n}\biggl)^{-k} \\
\end{eqnarray}
 

ここで\(n\rightarrow∞\)をとると、
\( 1- \displaystyle \frac{1}{n} \rightarrow1\)、\( 1- \displaystyle \frac{2}{n} \rightarrow1\)、\( 1- \displaystyle \frac{k-1}{n} \rightarrow1\)、\( \biggl(1- \displaystyle \frac{1}{n}\biggl)^{-k} \rightarrow1\)で、
\( \biggl(1- \displaystyle \frac{\lambda}{n}\biggl)^{n} \)については、\( \biggl(1- \displaystyle \frac{\lambda}{n}\biggl)^{n} \rightarrow e^{-\lambda} \)になるので、
 
\begin{eqnarray}
P(X=k) &=& \frac{\lambda^{k}}{k!} \cdot 1 \cdot 1 \cdot \cdot \cdot \cdot 1 \cdot e^{-\lambda} \cdot 1 \\
&=& \frac{\lambda^{k}e^{-\lambda}}{k!}
\end{eqnarray}

 
ポアソン分布は発生確率がめちゃくちゃ低い場合を想定するので、イメージとしては二項分布での正の確率\(p\)が、
\(p\rightarrow0\)と限りなく0に近い場合を考えるとイメージできれば、簡単です。
 
ポアソン分布も、二項分布と同じで、2つの事柄を考えます。
ただ二項分布と違うところは、滅多に起きない事柄を扱う点です。
滅多に起きないということは、確率が0にほぼ等しいということです。よくある例としては隕石に当たる確率レベルの滅多に起きない事柄です。
 
ポアソン分布のモーメント母関数
期待値、分散

-統計学
-

© 2023 Yosshi Labo. Powered by AFFINGER5