ポアソン分布

統計学および確率論で用いられる**ポアソン分布**（）とは、ある事象が一定の時間内に発生する回数を表す離散確率分布である。

数学者シメオン・ドニ・ポアソンが1838年に確率論とともに発表した。

ある離散的な事象について、**ポアソン分布**は所与の時間内での生起回数の確率を示し、指数分布は生起間隔の確率を示す。

## 定義
定数  に対し、 以上の整数を値にとる確率変数  が
:P(X=k)=\frac{\lambda^k e^{-\lambda}}{k!}
を満たすとき、確率変数  は母数  のポアソン分布に従うという。

ここで、 はネイピア数 ()であり、 は  の階乗を表す。また、 は所与の区間内で発生する事象の期待発生回数に等しい。

 は、「所与の時間中に平均で  回発生する事象がちょうど  回（ は非負の整数）発生する確率」に相当する。例えば、事象が平均で10分間に5回発生する場合、10分間の中で事象が発生する回数は、 のポアソン分布モデルを使って求められる。

確率変数 X が母数 \lambda のポアソン分布に従うことを、X\sim\operatorname{Pois}(\lambda) で表すこともある。

## 性質
### 平均・分散
ポアソン分布の平均  および分散  は、 に等しい。
:\begin{align}
\operatorname{E} [X] &=\lambda , \\
\operatorname{V} [X] &=\lambda .
\end{align}

### 最頻値
ポアソン分布の最頻値は、 以下で最大の整数である。

### 積率母関数
平均  のポアソン分布の積率母関数  は、
:M_X (t)=\operatorname{E} \left[ e^{tX} \right] =\sum_{k=0}^\infty e^{tk} P(X=k)=\sum_{k=0}^\infty e^{tk} \frac{\lambda^k e^{-\lambda}}{k!} =e^{\lambda (e^t -1)}
で与えられる。

### モーメント
ポアソン分布の高次モーメントは、 を含むトゥシャール多項式であり、二項係数を持つ。
:\begin{align}
 m_1 &= \operatorname{E} [X] = \lambda , \\
 m_2 &= \operatorname{E} [X^2] = \lambda^2 + \lambda , \\
 m_3 &= \operatorname{E} [X^3] = \lambda^3 + 3\lambda^2 +\lambda , \\
     &\vdots
\end{align}
ポアソン分布の  次のは  である。
:\operatorname{E} [X(X-1) \dotsm (X-n+1)]=\lambda^n .

### キュムラント
ポアソン分布の  次のキュムラント  は全て、平均  と等しい。
:\kappa_n = \left. \frac{\partial^n}{\partial t^n} \log \bigl( M_X(t) \bigr) \right\vert_{t=0} =\lambda.

### 再生性
ポアソン分布は再生性をもつ。すなわち、 と  とが独立な確率変数であり、それぞれパラメータ ,  のポアソン分布に従うとき、 確率変数の和  はパラメータ  のポアソン分布に従う。

### その他
* ポアソン分布は無限分解可能な確率分布である。すなわち、任意の正の整数nに対して、独立同分布な確率変数X_1,\ldots,X_nが存在して、\sum_{i=1}^nX_iの従う分布を母数 \lambda のポアソン分布とできる。具体的には、X_iが母数 \lambda/n のポアソン分布に従うとすれば良い。

* 母数の異なるポアソン分布同士のカルバック・ライブラー情報量はD_\mathrm{KL}(\operatorname{Pois}(\lambda)||\operatorname{Pois}(\lambda_0))=\lambda_0-\lambda+\lambda\log\left(\frac{\lambda}{\lambda_0}\right)で与えられる。

## 近似
 が十分に大きい（たとえば ）ならば、平均 、標準偏差  の正規分布はこのポアソン分布の非常によい近似となる。おおよそ  であれば、適切な連続な分布への修正がなされている場合に限り、正規分布はこのポアソン分布のよい近似となる。例えば  に関して、 が非負の整数ならば、 と置換することができる。

## ポアソン過程
 は、単位時間あたりの事象の平均発生回数などの割合と見なされる場合があり、**到着率**と呼ばれる。このとき、 を時刻  より前に発生した事象の回数とすると、
:P(N_t=k)=\frac{e^{-\lambda t} (\lambda t)^k}{k!}
となる。この式を満たす確率過程を**ポアソン過程**という。さらに、最初の事象が発生するまでの待機時間  は、指数分布による連続確率変数である。この確率分布は、次のように導くことができる。
:P(T>t)=P(N_t =0).
時間を含む場合、すなわち1次元ポアソン過程では、各時間内で事象が発生する回数を確率変数とする離散ポアソン分布と、待機時間を確率変数とする連続アーラン分布の両方を含んでいる。1よりも高い次元のポアソン過程についても同様である。

## 事象
;具体的な例
ポアソン分布は、に関連して発生する。
これは、離散的な自然現象（所与の領域内や所与の時間内において、0回、1回、2回、3回… と発生する現象）に該当するものであり、現象が発生する確率は、時間ないし空間内において一定である。また、時間または空間における発生間隔は指数分布になる。次に、その例を示す。
* 1時間に特定の交差点を通過する車両の台数。
* 1ミリリットルの希釈された水試料中に含まれる特定の細菌の数（細菌数検査における最確法）。
* 単位面積あたりの雨粒の数。
* 1ページの文章を入力するとき、綴りを間違える回数。
* 1日に受け取る電子メールの件数。
* 1時間あたりの電話がかかってくる件数。
* ある一定の時間内の店への来客数。
* 1分間のWebサーバへのアクセス数。
** 例えば、1時間あたりのウィキペディアの最近更新したページの編集数もおおよそポアソン分布。
* 1キロメートルあたりのある通り沿いのレストランの軒数。
* 1ヘクタールあたりのエゾマツの本数。
* 1立方光年あたりの恒星の数。

* 単位時間あたりの放射線の計数値であるカウント毎分やカウント毎秒（半減期による減衰や外部からの放射能などによる変動がないと仮定して）。

;歴史的例
上記の例のほか、歴史的に有名な事例としては、ロシア生まれでドイツで活躍した経済学者、統計学者の  による「プロイセン陸軍で馬に蹴られて死亡した兵士数」の例が知られている。ボルトケヴィッチは著書  において、プロイセン陸軍の14の騎兵連隊の中で、1875年から1894年にかけての20年間で馬に蹴られて死亡する兵士の数について調査しており、1年間当たりに換算した当該事案の発生件数の分布が母数  のポアソン分布によく従うことを示している。

;事象の特徴
上記のように、稀にしか起こらないような現象を大量に観測した結果がポアソン分布に従う例は極めて多く見られる。このようなポアソン分布に従う事象の中で、時間の経過とともに発生する事象の特徴は次のようにまとめられる。
#（希少性）：時間幅  の間に着目している事象がちょうど1回起こる確率が 、2回以上起こる確率が 
#（定常性）：事象の起きる確率は、どの時間帯で同じ
#（独立性）：事象の起きる確率は、それ以前に起こった事象の回数や起こり方には無関係
ここで、 は  に対して高位の無限小を表しており、 のスケールに注目したときに無視できる微小量であることを表す。

## 極限定理
パラメータが  と  である二項分布において、 を一定に保ったまま  を無限大に近づけると、その分布は平均  のポアソン分布に近づく。すなわち、
:\lim_{\lambda=np,~n\to\infty} \binom{n}{k} p^k (1-p)^{n-k} = \frac{\lambda^k e^{-\lambda}}{k!}
が成り立つ。これを**ポアソンの極限定理**という。この定理の名は、数学者シメオン・ドニ・ポアソンが1837年に著書  の中で結果を与えたことに由来する。なお、この中で、二項分布の極限としてポアソン分布が初めて導出されている。

導出の詳細を次に示す。計算には、以下の関係式を用いる。
:\lim_{n\to\infty} \left( 1-\frac{\lambda}{n} \right)^n =e^{-\lambda} .
ここで  とすると、
:\begin{align}
\lim_{n\to\infty} P(X=k) &=\lim_{n\to\infty} \binom{n}{k} p^k (1-p)^{n-k} \\
&=\lim_{n\to\infty} \frac{n!}{(n-k)! \, k!} \left( \frac{\lambda}{n} \right)^k \left( 1 - \frac{\lambda}{n} \right)^{n-k} \\
&=\lim_{n\to\infty} \underbrace{
	\left( \frac{n}{n} \right) 
	\left( \frac{n-1}{n} \right)
	\left( \frac{n-2}{n} \right)
	\dotsm
	\left( \frac{n-k+1}{n} \right)
}
\underbrace{\left( \frac{\lambda^k}{k!} \right)}
\underbrace{\left( 1 - \frac{\lambda}{n} \right)^n}
\underbrace{\left( 1 - \frac{\lambda}{n} \right)^{-k}}.
\end{align}
 を無限大に近づけると、4つの下波括弧のうち、最初の下波括弧の部分は  に近づく。2番目の下波括弧の部分には  が出現しないので、そのままである。3番目の下波括弧の部分は  に近づく。最後の下波括弧の部分は  に近づく。

したがって極限は存在し、
:\frac{\lambda^k e^{-\lambda}}{k!}
となる。

## 少数の法則
**法則**という言葉は、確率分布の同義語として使われることがあり、**法則収束**は分布の収束を意味する。したがってポアソン分布は、滅多に起こり得ない希少な事象の発生数の確率分布であることから、**少数の法則**と呼ばれることがある。

## 脚注

## 関連項目
*指数分布
*
*アーラン分布
*
*再生過程

## 外部リンク
*