ポアソン過程―計数分布と到着時間間隔分布―
1/08/2011 熊谷 朝臣
以前、「なぜ降雨イベント間隔とイベント時降雨量は指数分布に従うのか?」なんて文章書きましたが、その中では、ポアソン過程と言っておきながら登場するのは指数分布ばかりでした。その説明が無いのは、なんとも気持ち悪いんではないかな?と思ったので、ずばり、ポアソン分布で表現されるポアソン過程は、表現の仕方次第で指数分布やガンマ分布でも表現されるということを紹介します(ガンマ分布による表現は後編で)。なお、本稿を書くのには、尾崎(1996)による教科書を参考にしました。
ある時刻tまでに起こる事象の数をN(t)として、これを計数過程と呼びます。例えば、ある時間内の、店を訪れる客の数、かかってくる電話の本数、機械の故障の回数・・・など、いくらでも計数過程を想像できそうです。
ポアソン過程は、その計数過程で、これがポアソン分布で表現できることから、その名がつきました。ポアソン過程は、次のような性質を持ちます。と言うか、こういう性質を持つ確率過程をポアソン過程と呼びます。
1.時間ステップ間(tn
− tn-1)での増分(N(tn) −
N(tn-1))は独立であり、tnとtn+h(hは任意)に対する確率変数の分布(N(tn)とN(tn+h))は等しい。
2.微小時間hの間に事象が一回起こる確率は、hに比例する。つまり、
(lは比例定数、o(h)は微小誤差)
3.微小時間hの間に事象が二回以上起こる確率は、ほとんど無い。つまり、
4.N(0) = 0
んんん、言葉による説明の方が、よく意味が分かんないですね。これから数式を追っていった方が分かりやすいと思うので、今、分からなくても無視して進んで下さい。
ずばり、時刻tまでにk回の事象が起きる確率Pk(t)を考えます。すると、iを任意として、上記条件1〜4より(hが微小時間なのを再確認して下さい)、
(1)
となります。余裕ある人には有効な情報ですが、は、だから当然、はのテイラー展開の形になっているのを考えれば、「hが大きいとk = 0となる望みはない。」っていうのと「hが大きいとk = 1は絶対起きる。」ってのが分かりやすいと思います。
で、Pk(t)の漸化式を作ります。k = 0となる場合を分けて考える必要があります。
(2)
(3)
もう、そんなに難しくありません。(2)、(3)式は、hが微小時間であることを考えて、それぞれ次の微分方程式(4)、(5)式になります。ここで、o(h)は無視されました。
(4)
(5)
境界条件としてP0(0) = 1、Pk(0) = 0を利用して、k = 0, 1, 2, 3,・・・と言う風に実験的に解いて行きます。本来なら、数学的帰納法が必要ですが、ま、ここでは無しということで。
すると、こんな式が得られます。
(6)
はい、これは、ずばり、パラメータλtのポアソン分布ですね。もう一度、何を求めたのか確認しましょう。(6)式の表現するものは、「時刻tまでにk回の事象が起きる確率Pk(t)」です。パラメータλtのポアソン分布ということは、平均値と分散がλtになりますから、「時刻tで起きた事象回数の平均はλt回」と言えます。すると、λは「時間当りに起きる事象の平均回数」を表現しているということになります。
さて、もうちょっと頑張りましょう。ポアソン過程に従う計数過程N(t)で、事象が起こる時間間隔をτとすると、その確率分布関数F(τ)は「N(t+τ) −
N(t) = 0とならない」確率(=時間τの間に、必ず何かが起きる確率)で表現されます。すると、(1)式、(6)式を利用して、
(7)
と書けます。すると、τの確率密度関数f(τ)は、
(8)
となります。はい、これは指数分布ですね。くどいですが確認しましょう。ポアソン過程に従う計数過程を考えたとき、その事象の時間間隔τはパラメータをλ(=平均時間間隔)とする指数分布で表現できます。
引用・参考文献
尾崎 (1996) 確率モデル入門.朝倉書店.
□
[HOME] [ワンポイント・レッスン]