ポアソン分布の期待値～最近印象に残った授業～

はじめに

大人のための数学教室和(なごみ)の講師の松中です。最近印象に残った授業シリーズとして不定期に授業の内容を紹介していきたいと思います。今回紹介したいのはポアソン分布の期待値に関してです。

なおこちらの記事の内容は動画でも解説しています。

まずポアソン分に関して簡単に紹介します。例えば、過去の経験から1時間に平均λ回のいたずら電話がかかってくることがわかっている電話があったとします。ある1時間の間にその電話にかかってくるいたずら電話の回数が従う分布がポアソン分布になります。

このように単位時間に平均λ回発生する事象があるとき、ある単位時間にその事象が発生する回数をXとすると、Xはポアソン分布に従います。ポアソン分布に従う確率変数Xが各値をとる確率は次の式で与えられます。

P(X=k)=λke−λk!

具体的にポアソン分布に従う確率変数Xが各値をとる確率を見てみましょう。λを3とすると、各確率は上の表のように計算できます。グラフを書いてみるとたしかにXが3周りの値をとる確率が大きくなっていることが確認できます。λを3より大きくすると、このグラフの山は右の方に動いていきます。

次に期待値について簡単に解説します。確率変数Xの期待値E[X]はXがとる値とその値をとる確率を掛け合わしたものを、全て足すことで計算できます。

中学校や高校で1個のサイコロを投げたときに出る目の期待値を以下のように計算したと思います。

1×16+2×16+⋯+6×16=72

たしかにXがとる値とその確率をかけて足し合わせていますね。

ポアソン分布の期待値も定義通りに以下の式から計算できます。

E[X]=∑k=0∞kλke−λk!

実際に以下のような式変形で計算できます。

途中でk−1=lと変数変換を行っています。また、指数関数の以下のテイラー展開の結果も用いています。

ex=∑n=0∞xnn!

計算の結果、ポアソン分布の期待値はλであることが確認できました。これは元々の確率変数Xの定義から予想できていたことであり、予想通りの結果が得られうれしくなります。

参考までにポアソン分布の分散はX2の期待値を計算することによって得ることができます。結果はλとなり、ポアソン分布では期待値と分散が等しいことがわかります。

さて、このポアソン分布の期待値に関する授業が印象に残っている理由ですが、お客様がお持ちの統計のテキストで、ポアソン分布の期待値の導出にモーメント母関数を用いていたからです。

モーメント母関数は確率変数Xに対して、

MX(t)=E[etX]

で定義される関数のことです。これはtの関数になります。

モーメント母関数MX(t)をtで微分することで、

M′X(t)=ddtE[etX]=E[ddtetX]=E[XetX]

となりますが、ここでt=0とすることで、

E[X]=M′X(0)

であることがわかります。さらにモーメント母関数を微分していくことで、E[Xk]の値も求めることができるのです。

ということで授業の中でお客様とポアソン分布のモーメント母関数を求め、2階まで微分してみました。微分は数Ⅲまでやっている方ならできるかと思います。

その後t=0とすることで、以下のようにポアソン分布の期待値と分散を求めることができました。

たしかに定義に沿って求めた結果と同じく、期待値も分散もλになりました。めでたしめでたし。

今回は最近印象に残った授業としてモーメント母関数を用いたポアソン分布の期待値、分散の求め方を紹介しました。二項分布や正規分布などほかの有名な分布に対してもモーメント母関数を求め、そこから期待値や分散を導出することも可能です。興味を持った方は是非計算してみてください。

モーメント母関数は統計検定準一級や一級で出題されることがあります。このあたりの学習を苦手に感じている方には以下の講座がおすすめです。