独立性・二項分布・ベルヌーイ分布・大数の法則
樋口さぶろお
龍谷大学理工学部数理情報学科
確率統計☆演習
I L08(2016-11-24 Thu)最終更新: Time-stamp: ”2016-11-24 Thu 06:49 JST hig”
今日の目標
確率変数の独立性が判定できる・利用でき る
塚田確率統計2.4二項分布の母期待値が計算でき る
塚田確率統計4.2 塚田確率統計4.1大数の法則の意味が説明できる
塚田確率統計5.2 http://hig3.net樋口さぶろお (数理情報学科) L08独立性・二項分布・ベルヌーイ分布・大数の法則確率統計☆演習I(2016) 1 / 23
確率不等式・2変数の同時確率分布
L07-Q1
Quiz
解答
:多次元の確率変数の期待値
1 E[X+ 2Y] = 0·(1 + 2·0) +122(2 + 2·0) +121(3 + 2·0) +124 (1 + 2· 2) + 0(2 + 2·2) + 125(3 + 2·2) = 6212.
2 E[1[Y≥1](X, Y)] = 0·0 +122 ·0 +121 ·0 +124 ·1 + 0·1 +125 ·1 = 129.
3
fX(x) =
4/12 (x= 1) 2/12 (x= 2) 6/12 (x= 3) 0 (
他
)fY(y) =
3/12 (y= 0) 9/12 (y= 2) 0 (
他
)4 (1
の別解
)確率不等式・2変数の同時確率分布
L07-Q4 Quiz
解答
:離散型確率変数の独立性
1 fX(2) =fX(4) =fY(2) =fY(4) = 12
であり
, fX(x)·fY(y)̸=fXY(x, y)であり独立でない
.2 E[X] = 12·2 +12·4 = 3,E[Y] = 12·2 +12·4 = 3,E[XY] = 12·4 + 0·8 + 0·8 + 16·12 = 10,E[X+Y] = 3 + 3 = 6,Cov[X, Y] = 10−3×3 = 1.
Cov[X, Y] = E[XY] = E[X]E[Y] = 10−3·3 = 1. Cov[X, Y]̸= 0
からも独立でないことがわかる
.V[X] = V[Y] = 1
より
,ρ[X, Y] = √1 1√1 = 1.
L07-Q6
Quiz
解答
:独立と限らない確率変数の母期待値
1 E[−2X+ 3Y] =−2E[X] + 3E[Y] = 5.
2 V[−2X+ 3Y] = E[(−2X+ 3Y)2]−E[−2X+ 3Y]2=
(−2)2V[X] + 2(−2)(3)Cov[X, Y] + 32V[X] = 20−84 + 99 = 35.
樋口さぶろお (数理情報学科) L08独立性・二項分布・ベルヌーイ分布・大数の法則確率統計☆演習I(2016) 3 / 23
確率不等式・2変数の同時確率分布 独立性
ここまで来たよ
1
確率不等式・
2変数の同時確率分布 独立性
2
独立性・二項分布・ベルヌーイ分布・大数の法則 二項分布
ベルヌーイ分布
独立同分布に従う確率変数の和と大数の法則
確率不等式・2変数の同時確率分布 独立性
独立性
塚田確率統計§2.4,§3.6高校 数学B独立性
確率変数
X, Yが同時分布
fXY(x, y)を持つとする
. X, Yが独立とは
,fXY(x, y) =fX(x)×fY(y)
が成立することをいう
(世の中には
,同値な定義が多数
).独立とは
,X,Yが互いに
「無関係」であること
事象
A, Bが独立
⇐P(Aかつ
B) =P(A)×P(B) 塚田確率統計2.4の特別な 場合
.独立性と母共分散
塚田確率統計定理3.6.3X, Y
が独立なとき
,母共分散
Cov[X, Y] = 0.すぐ後で証明
.母共分散
Cov[X, Y] = 0は
,X, Yが独立であるための
????
条件
.樋口さぶろお (数理情報学科) L08独立性・二項分布・ベルヌーイ分布・大数の法則確率統計☆演習I(2016) 5 / 23
確率不等式・2変数の同時確率分布 独立性
L08-Q1
Quiz(離散型確率変数の独立性)
2
次元の離散型確率変数
(X, Y)を考える
.同時分布
fXY(x, y)は次の表 で与えられる
(現れない
X, Yの確率は
zeroである
).y\x 2 3
3 2/12 1/12
7 A B
X, Y
が独立になるように
,実数
A, Bを定めよう
.確率不等式・2変数の同時確率分布 独立性
X, Y
が独立であるとき
‘だけ
’成立する性質
E[ϕ1(X)×ϕ2(Y)] =E[ϕ1(X)]×E[ϕ2(Y)]
特にCov[X, Y] =(E[XY]−E[X]×E[Y] =)0 V[X+Y] =V[X] + V[Y]
E[XY] =∑
x
∑
y
fXY(x, y)·x·y
=∑
x
∑
y
fX(x)×fY(y)×x×y
=∑
x
fX(x)·x×∑
y
fY(y)·y= E[X]×E[Y] V[X+Y] =E[(X+Y)2]−E[X+Y]2
=E[X2] +2E[XY]+ E[Y2]−(E[X]2+ 2E[X]E[Y]+ E[Y]2)
=V[X] +2Cov[X, Y]+ V[Y]
樋口さぶろお (数理情報学科) L08独立性・二項分布・ベルヌーイ分布・大数の法則確率統計☆演習I(2016) 7 / 23
確率不等式・2変数の同時確率分布 独立性
L08-Q2
Quiz(独立な確率変数の期待値)
独立な確率変数
X, Yを考える
.E[X] = 2,E[Y] = 3,V[X] = 5,V[Y] = 11,
である
.1 E[(−2X+ 3Y)(X+ 5Y)]
を求めよう
.2 V[−2X+ 3Y]
を求めよう
.独立性・二項分布・ベルヌーイ分布・大数の法則 二項分布
ここまで来たよ
1
確率不等式・
2変数の同時確率分布 独立性
2
独立性・二項分布・ベルヌーイ分布・大数の法則 二項分布
ベルヌーイ分布
独立同分布に従う確率変数の和と大数の法則
樋口さぶろお (数理情報学科) L08独立性・二項分布・ベルヌーイ分布・大数の法則確率統計☆演習I(2016) 9 / 23
独立性・二項分布・ベルヌーイ分布・大数の法則 二項分布
二項分布
塚田確率統計§4.2高校 数学B二項分布
離散型確率変数
Xが次の確率分布を持つとき
,Xは二項分布
B(n, p)に 従うという
.P(X=x) =f(x) = {
nCx px(1−p)n−x (x= 0,1,2,3, . . . , n)
0 (
他
)意味
:確率
pで表の出るコインを
n回投げたとき
,x回表が出る確率
.B(40,0.1),B(40,0.5),B(40,0.7),B(4,0.8),B(20,0.8),B(40,0.8)
独立性・二項分布・ベルヌーイ分布・大数の法則 二項分布
二項分布の母平均値と母分散
E[X] =
np
,V[X] =
np(1 − p)
二項定理
高校 数学A(a+b)n=
∑n x=0
nCx(a+b)x
樋口さぶろお (数理情報学科) L08独立性・二項分布・ベルヌーイ分布・大数の法則確率統計☆演習I(2016) 11 / 23
独立性・二項分布・ベルヌーイ分布・大数の法則 ベルヌーイ分布
ここまで来たよ
1
確率不等式・
2変数の同時確率分布 独立性
2
独立性・二項分布・ベルヌーイ分布・大数の法則 二項分布
ベルヌーイ分布
独立同分布に従う確率変数の和と大数の法則
独立性・二項分布・ベルヌーイ分布・大数の法則 ベルヌーイ分布
ベルヌーイ分布
塚田確率統計§4.1ベルヌーイ分布
n= 1
の二項分布
B(1, p)のこと
P(X =x) =f(x) =
1−p (x= 0) p (x= 1) 0 (
他
)意味
:ベルヌーイ試行=(不公平な
)コイン投げ
.表がでる確率
p.ベルヌーイ分布の母平均値と母分散
E[X] =
p
, V[X] =
p(1 − p)
樋口さぶろお (数理情報学科) L08独立性・二項分布・ベルヌーイ分布・大数の法則確率統計☆演習I(2016) 13 / 23
独立性・二項分布・ベルヌーイ分布・大数の法則 ベルヌーイ分布
ベルヌーイ分布と二項分布
X1, X2, . . . , Xn
が独立で
Xi∼B(1, p)のとき
, Un=X1+· · ·+Xnは
Un∼B(n, p).なぜなら
自 分 の 言 葉 で ど う
ぞ
独立性・二項分布・ベルヌーイ分布・大数の法則 ベルヌーイ分布
L08-Q3 塚田確率統計章末問題4.13.1 L08-Q4
Quiz(
二項分布
)確率
p= 23で表のでるいかさまコインがある
. 100回投げる
.1
表が
50回でる確率を求めよう
.2
表がでる回数の母平均値を求めよう
.3
表がでる回数の母分散を求めよう
.樋口さぶろお (数理情報学科) L08独立性・二項分布・ベルヌーイ分布・大数の法則確率統計☆演習I(2016) 15 / 23
独立性・二項分布・ベルヌーイ分布・大数の法則 ベルヌーイ分布
L08-Q5
Quiz(ベルヌーイ分布)
ある宝くじは
,あたりとはずれの
2種類の結果だけがある
.あたりの確率 は
0.05である
.あたりの賞金は
1000円
,はずれの賞金は
0円である
.賞 金を確率変数
Yとする
.1 Y
と
,ベルヌーイ分布
B(1, p)に従う確率変数
Xとの関係を書こう
.2 Y
の母平均値と母分散を求めよう
.独立性・二項分布・ベルヌーイ分布・大数の法則 独立同分布に従う確率変数の和と大数の法則
ここまで来たよ
1
確率不等式・
2変数の同時確率分布 独立性
2
独立性・二項分布・ベルヌーイ分布・大数の法則 二項分布
ベルヌーイ分布
独立同分布に従う確率変数の和と大数の法則
樋口さぶろお (数理情報学科) L08独立性・二項分布・ベルヌーイ分布・大数の法則確率統計☆演習I(2016) 17 / 23
独立性・二項分布・ベルヌーイ分布・大数の法則 独立同分布に従う確率変数の和と大数の法則
復習
チェビシェフの不等式
Chebyshev’s inequality 塚田確率統計3.5X:
離散型または連続型確率変数
µ= E[X]:母平均値
σ2 = V[X]:
母分散
a >0:任意の正の実数 のとき次が成立する
.P(|X−µ| ≥aσ)≤ 1 a2
要するに
自 分 の 言 葉 で ど う ぞ
独立性・二項分布・ベルヌーイ分布・大数の法則 独立同分布に従う確率変数の和と大数の法則
独立同分布の性質
塚田確率統計5.1独立同分布
(i.i.d.)離散型
/連続型確率変数
X1, X2, . . . , Xnが
,たがいに独立で
,すべて同じ 確率分布に従う
(同じ確率関数
f(x))とする
.これを
X1, . . . , Xnは独立同分布に従う
(i.i.d.=independent and identically-distributed)という
.新しい確率変数
: Un=X1+· · ·+Xn母平均値
E[Xi] =µ,母分散
V[Xi] =σ2としたとき
,E[Un] =
∑n i=1
E[Xi] =n×µ.
V[Un] =
∑n i=1
V[Xi] =n×σ2.
Un
の確率密度関数はこん な感じ
?樋口さぶろお (数理情報学科) L08独立性・二項分布・ベルヌーイ分布・大数の法則確率統計☆演習I(2016) 19 / 23
独立性・二項分布・ベルヌーイ分布・大数の法則 独立同分布に従う確率変数の和と大数の法則
新しい確率変数
: Wn= n1Un= 1n(X1+· · ·+Xn) E[Wn] =E[1nUn
]= 1
n×n×µ.
V[Wn] =V[1
nUn
]= (1
n )2
×n×σ2.
Wn
の確率密度関
数はこんな感じ
?独立性・二項分布・ベルヌーイ分布・大数の法則 独立同分布に従う確率変数の和と大数の法則
L08-Q6
Quiz(独立同分布にしたがう変数の和)
確率変数
X1, . . . , X100は
E[Xi] = 3,V[Xi] = 7の独立同分布に従う
.次の確率変数の母平均値と母分散を求めよう
.1
確率変数
A= 1001 (X1+X2+X3+· · ·+X100)2
確率変数
B= 101(X1+X2+X3+· · ·+X100−100·3)3
確率変数
C= 110√
7(X1+X2+X3+· · ·+X100−100·3)
樋口さぶろお (数理情報学科) L08独立性・二項分布・ベルヌーイ分布・大数の法則確率統計☆演習I(2016) 21 / 23
独立性・二項分布・ベルヌーイ分布・大数の法則 独立同分布に従う確率変数の和と大数の法則
大数の
(弱)法則
塚田確率統計5.2Wn= 1 n
∑n i=1
Xi
は
∀ϵ >0 limn→+∞P(|Wn−µX| ≥ϵ) = 0
を満たす
.つまり
,n大で
Wnは母平均値
µX= E[X]に「必ず近い」
(確率収束)証明
µWn =µX, σ2Wn =σX2/n. Wn
に対するチェビシェフの不等式より
, P(|Wn−µX| ≥a× σX√n)≤ 1 a2 a= √σϵX
n
とすると
,n→+∞で
P(|Wn−µX| ≥ϵ)≤ σX2 ϵ2n →0
これが母期待値の直観的意味
.要するに
,自分の言葉でどうぞ
独立性・二項分布・ベルヌーイ分布・大数の法則 独立同分布に従う確率変数の和と大数の法則
連絡
予習問題は
,次々回の授業直前 を締切
(そこまでの最高点を記録
)と
します
.でも
, Trialまでにやったほうが効率いいと思う
.前からそうだけど
,予習問題が満点だと
, Trialの満点の
1/3まで保証 されます
.配布資料は
1-503向かいの引出
,http://hig3.netで再配布
.加減乗除と平方根
(ルート
)の使える電卓持ってきてね
.関数電卓で なくてもいいです
.携帯電話の機能・アプリでもかまいません
.樋口オフィスアワー木
6金昼
(1-502), Mathラウンジ月
-木昼
(1-614)次回は
塚田確率統計4.7 塚田確率統計5.3https://manaba.
ryukoku.ac.jp
樋口さぶろお (数理情報学科) L08独立性・二項分布・ベルヌーイ分布・大数の法則確率統計☆演習I(2016) 23 / 23