第8回 ポアソン分布 確率･統計Ⅰ ここです！ 確率変数と確率分布 確率変数の同時分布、独立性 確率変数の平均 確率変数の分散

Presentation on theme: "第8回 ポアソン分布 確率･統計Ⅰ ここです！ 確率変数と確率分布 確率変数の同時分布、独立性 確率変数の平均 確率変数の分散"— Presentation transcript:

1 第8回 ポアソン分布 確率･統計Ⅰ ここです！ 確率変数と確率分布 確率変数の同時分布、独立性 確率変数の平均 確率変数の分散
確率変数の共分散 ベルヌイ試行、二項分布 二項分布（続き）、幾何分布 ポアソン分布 正規分布 正規分布（続き） 大数の法則、中心極限定理 統計学の基礎1（母集団と標本、確率論との関係） 統計学の基礎2（正規分布を用いた推定・検定） ここです！

2 ポアソン分布 ポアソン分布 二項分布のポアソン近似 ポアソン分布の意味 ポアソン分布の平均・分散

3 ポアソン分布 P( X = x ) が次の式で与えられる確率分布を、パラメータλのポアソン分布 という：
・ポアソン分布は離散分布であるが、x の値は 0, 1, 2, … と理論上は無限個ある。（ただし、大きい x に対する確率は急速に小さくなり、事実上無視できる。そのため、重要なのはいつも x = 0, 1, 2 あたりである。） 問のヒント：eλ のテーラー展開を使う。（分からない者は解析学の本を調べよ。） ( x = 0, 1, 2, … ) 離散型； 値無限個 （問） これが確かに確率分布であることを確かめよ。

4 ポアソン分布 x = 0 1 2 3 4 … （比率はこのまま、全部で 1 になるように、全体に e -λ を掛ける。）
二項分布が、(p+q)n の展開式の各項をバラバラにした分布なのと同様、ポアソン分布は、eλ のテーラー展開の各項をバラバラにした分布になっている。（ただしそれは“比率”の話であって、そのままだと和が 1 でないので、和が 1 になるように定数 eλ で割っておく） ※ eλのテーラー展開の項は、先に行くほど急速に小さくなるので、ポアソン分布の確率も r が（λより）大きくなると急速に小さくなる。（たとえば、λ=1 のとき、r > 1 の確率をすべて合わせても、r = 1 の確率より小さい。） （比率はこのまま、全部で 1 になるように、全体に e -λ を掛ける。）

5 ポアソン分布のグラフ x =λ ポアソン分布 (λ=10 ) （ x = 0, 1, 2, …… ） 最大値
λが 10 で、わりと大きいので、かなり正規分布に近い形をしているが、まだ少し左にゆがんでいる。 ※注: λは整数とは限らないが、x は整数なので、λが整数でない場合の最大値はλの前後どちらかの整数（一番近くとは限らない）でとる。 x =λ 最大値

6 ポアソン分布のグラフとλ λによる変化 x → 1 2 λ= 0.7 λ= 1 λ= 2 λ= 3 λ= 5 λ= 8
1 2 λが小さいときは、x = 0, 1, 2 あたりの確率がほとんどで、あとは急速に減少していく。

7 ポアソン分布 ポアソン分布 二項分布のポアソン近似 ポアソン分布の意味 ポアソン分布の平均・分散

8 二項分布のポアソン近似 X を 二項分布 B(n, p) に従う確率変数とする。 λ=np を一定にして n→∞ のとき、
「ポアソン分布」 証明には、解析学の基本公式のひとつ (1+ x/n)n →ex を使う。 n がある程度大きく、しかも p が小さい場合に、左辺（二項分布） を右辺の式（ポアソン分布）で近似できる。 　　　目安としては、n≧100, p≦0.05

9 二項分布のポアソン近似の様子 p=0.01, n=300 ∴λ=np=3 ポアソン近似の目安は n≧100, p≦0.05

10 二項分布のポアソン近似の様子 二項分布 B(300, 0.01) ∴np=3 （…省略） x → 0 …300
数学的には、横軸（r） は 300 までとらないといけないが、r≧10 の確率は無視できる。

11 二項分布のポアソン近似の様子 二項分布 B (3000, 0.001) ∴np=3 （…省略） x → 0 …3000
数学的には、横軸（r） は 3000 までとらないといけないが、r≧10 の確率は無視できる。

12 二項分布のポアソン近似の様子 二項分布 B (30000, 0.0001) ∴np=3 （…省略） x → 0 …30000
数学的には、横軸（r） は までとらないといけないが、r≧10 の確率は無視できる。

13 二項分布のポアソン近似の様子 （λ=3 のポアソン分布） （…省略） x → 0 ……∞
数学的には、横軸（r） は ∞ までとらないといけないが（どんな大きな r に対する確率も0ではない）、r≧10 の確率は無視できる。

14 二項分布のポアソン近似の方法 X が 二項分布 B(n, p) に従うとき、 nが大きく（n≧100）、 pが小さい（ p≦0.05）ならば
のかわりに を計算 こっちよりはずっと楽 （λ= np ） これもシンドイが

15 二項分布のポアソン近似（例題） X を B(400, 1/400) に従う確率変数とするとき、1 - P(X=0) を求めればよい。
例題： ビジ確率1/400 のパチスロを400回回したとき、少なくとも一回当たる確率を求めよ。 X を B(400, 1/400) に従う確率変数とするとき、1 - P(X=0) を求めればよい。 まじめに計算すると P(X=0) = (399/400)400

16 二項分布のポアソン近似（例題） パラメータλ= 400×1/400 = 1 のポアソン分布として計算すれば 1 - 0.37 = 0.63
例題： ビジ確率1/400 のパチスロを400回回したとき、少なくとも一回当たる確率を求めよ。 パラメータλ= 400×1/400 = 1 のポアソン分布として計算すれば 要するに、eのテーラー展開 e= /2! + 1/3!+ … の第一項×(1/e) が P(X=0) である。 第二項×(1/e) が P(X=1)、第三項×(1/e) が P(X=2)、… だから、P(X=0) = 1/e = 0.37, P(X=1) = 1/e = 0.37, P(X≧2) = 1－ 0.37×2 = 0.26 つまり、確率1/400のパチスロを400回やっても、一回も当たらない確率が37%もあるのである。ちょうど一回当たる確率も同じ37%（そして2回以上当たる確率が残りの26％）。 「平均1回当たる」（二項分布のnpでもポアソン分布のλでも平均はたしかに1）というのは事実であるし、「当たるまでに回す平均回数が400回」（幾何分布の平均は1/pであった）というのも事実であるが、「400回やれば１回ぐらい当たって当然」ではないことに注意！ 「平均」だけですべてを判断してはいけない。分散npq=399/400 がほぼ 1 もあるのだ。分散の理解こそが確率論の鍵…？ = 0.63

17 [4] ある部品が一定期間内に故障を起こさない確率を「精度」と呼ぼう。
［再演習］ 二項分布 [4] ある部品が一定期間内に故障を起こさない確率を「精度」と呼ぼう。 　(1) 精度が0.999の部品1000個のうち、どの部品も故障しない確率を求めよ。 　(2) 精度が0.999の部品10000個のうち、どの部品も故障しない確率を求めよ。 ヒント：(1) 故障する個数 X は、B(1000, 0.001) に従う。(0.999)1000 の計算は大変だが、ポアソン近似すれば、P( X=0 ) = 1/e = 0.37 (2)故障する個数 X は、 B(10000, 0.001) に従う。 ポアソン近似すれば、P( X=0 ) = 1/e10 =

18 ポアソン分布 ポアソン分布 二項分布のポアソン近似 ポアソン分布の意味 ポアソン分布の平均・分散

19 ポアソン分布の意味 空間または時間の１単位あたり、 平均λ回起こる（ことがわかっている） 事象があるとする。
例1：製品１個あたり平均２個の傷が入る工芸品 いきなりこのような例を出されると、どこがベルヌイ試行（二項分布）と関係があるのかわかりにくいだろう。 例2：ぶどうパン１個あたり平均２個の干しぶどうが入る工場 例3：昼間の１時間あたり平均３回の電話がかかる会社

20 ポアソン分布の意味 このとき、特定の１単位に、 実際にそれが 何回起こるかの確率 パラメータλのポアソン分布
例1：製品１個あたり平均２個の傷が入る工芸品 例2：ぶどうパン１個あたり平均２個の干しぶどうが入る工場 例3：昼間の１時間あたり平均３回の電話がかかる会社

21 ポアソン分布の意味 例1：製品１個あたり平均２個の傷が入る工芸品 区画数 N 区画に傷のある確率 p = 2 / N
各区画に傷ができる事象は独立 ∴ 製品1個の傷の個数 X は B(N, p) に従う （N個分のベルヌイ試行とみなせるから） N p = λ = 2

22 ポアソン分布の意味 例2：パン１個あたり平均２個の干しぶどうが入る工場 × N 個 パンの数 N 干しぶどうの数 2N
干しぶどう１個がこのパンに入る確率 p = 1 / N 例１と異なり、１個のパンを細かく分割すると考えるとなんだかおかしいだろう。 たくさんの小麦粉とたくさんの干しぶどうを混ぜて、パンを無限に生産していくわけだから、パン１個あたりの干しぶどうの平均個数λを一定にして、パンの個数を増やした極限と考えればよい。 各干しぶどうがこのパンに入る事象は独立 ∴ パン1個の干しぶどうの個数 X は B(2N, p) に従う （2N個分のベルヌイ試行とみなせるから） 2N p = λ = 2

23 ポアソン分布の意味 例3：昼間の１時間あたり平均3回の電話がかかる会社 1時間を N等分 1区間に電話のかかる確率 p = 3 / N
t1 t2 t1 t2 t1 t2 1時間を N等分 1区間に電話のかかる確率 p = 3 / N これは例1と同じ考え方の時間版 各区間に電話のかかる事象は独立 ∴ 1時間の電話の回数 X は B(N, p) に従う （N個分のベルヌイ試行とみなせるから） N p = λ = 3

24 ポアソン分布（例） λ=0.7のポアソン分布の値 馬に蹴られて死んだプロシアの兵士 （1875-1894）
これは、ポアソン分布の応用発見にかかわる歴史的な例である（ポーランドの統計家ボルトキーヴィッツが1898年の論文でポアソン分布の適用例として示したもので、ポアソン分布自体の発見からは60年も経っていた）。 pが小さく、nが大きいと考えられる。現代では、交通事故死者数などがこれにあたるだろう。 λ=0.7のポアソン分布の値

25 ポアソン分布（例） 馬に蹴られて死んだプロシアの兵士 （ ）

26 ポアソン分布（例題） X を 1製品あたりの傷の個数とすると、 Xはパラメータλ=2 のポアソン分布に従う。
例題： 1個につき平均2個の傷が普通の工芸品がある。傷が5個以上ある製品は返品を受け付けている。製作した製品の何%が返品されると考えられるか。 X を 1製品あたりの傷の個数とすると、 Xはパラメータλ=2 のポアソン分布に従う。 このように、ポアソン分布の問題では、N や p はぼかされており、「平均個数」のみが与えられる。このため、慣れないと解き方が分かりにくいかもしれない。

27 ポアソン分布（例題） 例題： 1個につき平均2個の傷が普通の工芸品がある。傷が5個以上ある製品は返品を受け付けている。製作した製品の何%が返品されると考えられるか。 約 5% ≒ 0.052

28 ポアソン分布の応用例 細胞内の染色体交替数、バクテリア数など、生物統計への応用 放射線物質の崩壊
在庫管理（いくつ仕入れておけば売り切れの確率をあるレベル以下にできるか） 電話や道路の混雑状況の見積もり→回線数をどれだけ用意すれば、（平常時は）十分やっていけるか etc…

29 ポアソン分布 ポアソン分布 二項分布のポアソン近似 ポアソン分布の意味 ポアソン分布の平均・分散

30 ポアソン分布の平均と分散 ポアソン分布に従うとき、 E(X) =λ V(X) =λ 確率変数 X がパラメータλの
二項分布の平均が np で、λ=np を一定にして n→∞ とした極限がパラメータλのポアソン分布だったから、ポアソン分布の平均がλであることは自然だろう。また二項分布の分散は npq で、qは1に近づくのだから、分散がλ=np になることもうなずける。 問のヒント：二項分布と同様に、階乗の性質を使う方法と、微分を使う方法がある。それぞれ、二項分布の場合と全く同様である。 （問） これらを確かめよ。

31 ポアソン分布のグラフ r =λ ポアソン分布 (λ=10 ) （ r = 0, 1, 2, …… ） 最大値
λが 10 で、わりと大きいので、かなり正規分布に近い形をしているが、まだ少し左にゆがんでいる。 ※注: λは整数とは限らないが、r は整数なので、λが整数でない場合の最大値はλの前後どちらかの整数（一番近くとは限らない）でとる。 r =λ 最大値

32 ポアソン分布のグラフとλ λによる変化 r → 1 2 λ= 0.7 λ= 1 λ= 2 λ= 3 λ= 5 λ= 8
1 2 λが小さいときは、r = 0, 1, 2 あたりの確率がほとんどで、あとは急速に減少していく。

33 メニューに戻る メニューへ