統計学第２回 「確率と確率分布」

00

5

10 20 30 40 50 60

10 15 20 25 30 35•http://phi.ypu.jp/statlib/l2-2003.pdf•統計学＝「不確実性を考慮した論理的推論」であった。
•不確実性＝確率的？
•確率とは？
•あらゆる統計的な考え方の基礎

国際看護

栄養 生活環境

学科別出席者数123

学年別出席者数4

どういう現象が確率的？
• サイコロを振ったときの目：振ってみるまでは１から６のどれが出るかはわからない。どの目がでる可能性も等しいから。
• 天気予報：「明日の天気予報は晴れ」といっても「必ず晴れる」とは限らない。「曇ったり雨が降ったりする可能性も少しはあるが，晴れる可能性が高い」ことを意味する。
• 喫煙と肺がんの関係：「タバコを吸うと肺がんになる」とい
うのは，タバコを吸った人と吸わなかった人を比べて，肺がんになった人の割合が吸った人の方で高い，という関係を示す。タバコを吸っても肺がんにならない人もいるし，吸わなくても肺がんになる人もいる。
• 実は自然界のほぼすべての現象は確率的。

統計的事象• こういう「不確かさ」に潜む法則性（長期間繰り返し観察したり，大集団で観察すると見られる）を考える学問を確率論と呼ぶ。
• この種の法則性をもつ現象を，「統計的事象」と呼ぶ。
• 統計的事象の確かさの度合いを示すのに便利なモノサシが「確率」

標本空間• 統計的事象を捉えるには，「どんなことが起こりうるのか」という範囲を定めることが必要。
• 現象は一般に多面的で様々な観察方法がある。以下３
点によって統計的現象を捉えた，記号化された結果の集合のことを「標本空間」と呼ぶ。
– 観察を行う面を特定する– 起こりうる結果の範囲を規定する– その範囲内の各結果に記号を対応させる• 個々の結果の起こりうる可能性を示す数値（これを「確率」という）を考える。一般には「どの結果も同程度に起こる」と考える。各結果に対応付けられた確率は０から１までの数値であり，各確率の値の総和は１にならねばならない。
• サイコロの目では，標本空間は｛１，２，３，４，５，６｝

事象の確率• 問題は，個々の結果の可能性よりも，いくつかの結果が複合された集合（これを「事象」という）の起こる可能性がどのくらいか，ということ。つまり，事象とは，「標本空間の部分集合」である。
• サイコロの例では，「丁（目が偶数）」とか「目が５以上」とか「目が１」とかいうことが事象である。
• ある事象の確率は，その事象に含まれる各結果の生起確率の和である。従って，各結果の生起確率が等しい
場合は，その事象に含まれる結果の場合の数をすべての場合の数で割ると，その事象の確率になる。サイコロの例では，「目が５以上」という事象の確率は，２／６＝０．３３３・・・である。

余事象・和事象・積事象・排反事象• 起こりうるすべての結果の集合を「全事象」という。つまり，全事象は標本空間に等しい。
• 決して起こらない事象を「空事象」といい，空集合 φ で表す。
• 事象 E に対して， E が起こらないという事象を E の「余事象」という。サイコロの例では，「目が偶数」という事象の余事象は「目が奇数」である。 E の余事象を E と書く。
• 事象 E と F の少なくとも一方が起こるという事象を， E とF の「和事象」といい， E∪F で表す。
• 事象 E と F の両方が起こるという事象を， E と F の「積事象」といい， E∩F で表す。
• 事象 E が起これば F は決して起こらないとき， E と F は「排反事象」であるという。 E と F が排反事象なら， E∩F＝ φ である。

事象の確率を表す記号• 事象 E が起こる確率を Pr(E) と書く。
• サイコロで考えると，１回振ったとき「偶数の目が出る」という事象 E が起こる確率 Pr(E) は， {2,4,6} の場合の数３を， {1,2,3,4,5,6} の場合の数６で割った値なので Pr(E)=0.5 。

２回振って「少なくとも１回は偶数の目」の確率は？
•••••0.5+0.5=1.0 ではないのは自明。
偶数の目が１回目に出る事象 E1 と２回目に出る事象 E2 とは排反ではない集合で考えれば， Pr(E1∪E2)= Pr(E 1)+ Pr(E 2) － Pr(E1∩E 2) は自明。この式を「加法法則」と呼ぶ。
「２回とも奇数」の余事象なので 1 － Pr(E1  E2)
事象 E と事象 F が排反なら， Pr(E ∪F)=Pr(E)+Pr(F) という「加法定理」が成立。

１回目が偶数

２回目が偶数(2,1) (2,3) (2,5) (2,2) (2,4) (2,6)(4,1) (4,3) (4,5) (4,2) (4,4) (4,6)(6,1) (6,3) (6,5) (6,2) (6,4) (6,6)(1,2) (3,2) (5,2)(1,4) (3,4) (5,4)(1,6) (3,6) (5,6)(1,1) (1,3) (1,5)(3,1) (3,3) (3,5) ＝ １回目も２回目も奇数(5,1) (5,3) (5,5)

事象の独立性と乗法定理• 事象 E が起こっているときに事象 F が起こる確率を， E が起こったときの F の「条件付き確率」といい， Pr(F|E) と書く。
• Pr(F|E)=Pr(F∩E)/Pr(E) である。
• 事象 E と事象 F が互いに無関係（独立）なら， Fの条件付き確率 Pr(F|E) は， Pr(F) と一致する。
逆にいえば， Pr(F)=Pr(F|E) のときに事象 E と事象 F は互いに独立であるという。独立でないとき「従属である」という。
• 上記２つの式から，事象 E と事象 F が独立なら，Pr(F∩E)=Pr(F)×Pr(E)“乗法定理”が成立。

確率の４つの定義• 操作的接近＝統計的定義：数多く試したときの相対度数の極限。例えば，事象 E が起こる確率 Pr(E) は， N 回試したときにN1 回事象 E が起こるとして， N を無限大にしたときに N1/N が漸近する値である。
• 対称的確率：サイコロの場合，６通りの目の出る確率はどれも等しくなければならず，その和は１でなくてはならないので，例えば１の目が出る確率は 1/6 となる。限定的かつ循環論法。
• 公理的客観確率：標本空間の各要素を ei として， Pr(ei)>=0 かつ Pr(e1)+Pr(e2)+ ・・・ +Pr(eN)=1 かつ Pr(E)=ΣPr(ei) を公理とする。厳密（正確にはもっと長い説明が必要。テキスト参照）。
• 主観確率：ビッグバンとか超新星の爆発といった，観念的にも二度と繰り返すことのできない事象についての「見込み」を扱う。決定理論において重要。

大数の法則（操作的接近の根拠）
n= 100000550101001515020n= 100

1234561245656n= 100000

005005000

1000

10000

150015000

n= 10000

31234561234

確率変数と期待値• 例えば，スロットマシンにコインを入れると，ごくたまに，投入金額の何十倍ものコインが出てくる。
• マシン利用者全員に返ってくる賞金の合計を利用回数で割った値が，１回に期待される賞金額である。これを賭け金で割った値を「賞金還元率」と呼ぶ。すべての賭け事で胴元が儲かるようになっているのは，賞金還元率が 100 ％未満だからである。宝くじでは 40 ％，競馬では75 ％と言われる。
• 一般に，賞金額が x1, x2, x3, ・・・で，その賞金が得られる確率が p1, p2, p3, ・・・のように設定されたスロットマ
シンの期待賞金額 M は， M=x1p1+x2p2+x3p3+ ・・・で与えられる。
• このスロットマシンのようなものを確率変数といい，期待賞金を一般に期待値と呼ぶ。

分散• 期待賞金が同じでも，値動きの幅が小さいと一喜一憂の程度が小さく，逆に幅が大きいと滅多に当たらないが当たったときの喜びは大きくなる。
つまり，ギャンブル性は，値動きの幅と，チャンスの大きさに依存している。
• 各賞金がどれくらい期待賞金から隔たりがあり，それを獲得できる可能性がどれくらいあるのかを見積もれば，ギャンブル性が表せる。
• V= （マシンのギャンブル性） =Σ （期待値からの隔たり） × （可能性）という値が定義できる。このV を「分散」と呼ぶ。なお，各賞金額 x と期待値
M の隔たりは，普通，差の平方値 D=(xM)2 で表す。

確率変数と確率分布• 一般に，とりうる値の集合 x=(x1,x2,x3, ・・・ ) と，それぞれの値が実現する確率 p=(p1,p2,p3, ・・・) が与えられていて，事象として x のうちどれか１つの値のみ実現するとき， (x,p) という１セットを，「確率変数」と呼んで， X で表す。
• 期待値は E(X)=μ=Σxipi• 分散は V(X)=σ2=Σ(xi-μ)2pi• 分散の平方根 σ を標準偏差と呼ぶ。
• 横軸に x の各々の値を示す位置に， p の各々の可能性を示す高さの棒を立ててみれば，これが確率変数の「確率分布」ということになる。

ベルヌーイ試行と２項分布• １回の実験で S か F かのどちらかが起こり，しかもそれらが起こる可能性が， Pr(S)=p ， Pr(F)=1p=q で何回実験しても変わらないとき，これを「ベルヌーイ試行」という。
• ベルヌーイ試行を n 回行って， S がちょうど k回起こる確率は， Pr(X=k)=nCkpkqn-k• nCk は２項係数と呼ばれる。このような確率変数 X は，「２項分布に従う」といい， B(n,p) と表す。 E(X)=np ， V(X)=npq である。

二項分布のシミュレーション
n= 100050

100200300

100 150 200 250 300400

n= 4012340123568

n= 50002050406010080150120200

n= 2040 12 3 45 6 78 911357911 13 15 17 19▼ 正二十面体（各面には１から 20 までの数字が割り振られている）サ
イコロを n 回（ n=4, 10,20, 50 ）投げたときの，１から４までの目が出る回数を１試行と考えれば，これはベルヌーイ試行である。１回投げたときに１から４までの目が出る確率は 0.2 であるとして（＝母比率を 0.2 とする），左
図は試行 1000 セットの度数分布

二項分布の理論値
各 n についての理論的な確率分布は， Pr(X=k)=nCk0.2k0.8n-k より下図のようになる。
n= 10

0.00.00

0.1

0.10

0.2

0.20

0.3

0.40.30

n= 40.04

0.080.12

n= 500.00

0.00

0.05

0.10

0.15

0.20

n= 20

正規分布
• n が非常に大きい場合には，２項分布 B(n,p) の確率
Pr(X=np+d) という値が，1/√(2πnpq) ・ exp(-d2/(2npq)) という値で近似できる。
• 一般にこの極限である，Pr(X=x)= 1/√(2πσ2) ・ exp(-(x-μ)2/(2σ2))という形をもつ確率分布を正規分布と呼び，N(μ, σ2) と書く。
• z=(x-μ)/σ と置けば，Pr(Z=z)= 1/√(2π) ・ exp(-z2/2)となる。これを標準正規分布と呼び， N(0,1) と書く。
• 統計学でよく使われる確率分布であるカイ二乗分布とか t 分布とか F 分布は，正規分布から導かれる。

練習問題
８頭で出走する競馬のレースがあり，「どの馬が勝つチャンスも等しい」と仮定した場合，ある特定の馬が勝つと予想して当たる確率は 1/8 となるが，２回のレースの少なくともどちらか一方に当たる確率はいくらか？

解答例• １回目に当たる事象と２回目に当たる事象は排反ではないから， 1/8+1/81/64=15/64• 「１回目も２回目も当たらない」ことはない，という確率なので， 1/8)·(11/8)=15/64