資料置き場 hustat2017 20171201slide

統計学 第９回 正規母集団における標本理論

担当者： 高木 真吾

URL: http://sites.google.com/site/hustat2017/

質問等は， [email protected] ^{までお願いします．}

December 1st, 2017

復習

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

復習

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ 確率変数：実現するまではどの値が出るかわからないが，どの値がどのくらい

の出やすさで実現するかに関するルールは定められているもの

■ 離散型確率変数：個々の取りうる値に確率が付与されている

◆ ベルヌーイ分布に従う確率変数（値１を確率 _p，値０を確率 _{1 − p}）

◆ 二項分布に従う確率変数

■ 連続型確率変数：実現パターンが密度関数として表現されている

◆ ^正規分布

◆ ^一様分布

復習

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ 確率変数の期待値の基本公式： 取りうる値 _× 確率

◆ ベルヌーイ分布に従う確率変数 _X：

■ 平均：E[X] = 1 · p + 0 · (1 − p) = p

■ 分散：V[X] = E[(X − E[X])^{2] = E[X2}] − (E[X])^{2 =} 1² _{· p + 0}² · (1 − p) − p² = p(1 − p)

◆ 正規分布に従う確率変数 _{X (}密度関数は φ(x) = ^{√ 1}

2πσ² exp{−0.5(x − µ)^2/σ2})

■ ^平均：E_{[X] =} ^R∞

−∞ x · φ(x)dx = µ

■ 分散：V[X] = E[(X − E[X])^{2] == R}_−∞^∞ (x − µ)² · φ(x)dx = σ²

復習

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ ^{独立な確率変数} _({Xi}ⁿ_i=1) の和に関する期待値演算（_{Y = β0 +}

Pn

i=1 ^{βi · Xi）}

ただし，_{βi}ⁿ

i=0 は確率変数ではない定数．

◆ ^平均：E_{[Y ] = β0 +} ^Pn_{i=1 βi · E[Xi]}^{，分散：V}_{[Y ] =} ^Pn_{i=1 βi}² _{· V[Xi]}

◆ ^例）「標本平均」という確率変数について考える．_{X =}¯ ¹

n

Pn

i=1 ^Xi

■ ^{「標本平均」とは，}_n ^{個の確率変数を足し，}_n ^{で除す 確率変数}

◆ ^{上の公式で，}_{β0 = 0, β1 = β2} = · · · = β^{n = 1/n とした場合に相当}

■ 平均：E_{[ ¯X] =} ¹

n

Pn

i=1 ^{E[Xi]，分散：V[ ¯X] =} _n¹²

Pn

i=1 ^V[Xi]

◆ さらにすべての確率変数 _{Xi}ⁿ

i=1 が同じ平均，同じ分散を持つなら

（E_{[Xi] = µ, V[Xi] = σ}², for i = 1, 2, . . . , n^）

■ ^平均：E_{[ ¯X] =} ¹

n

Pn

i=1 ^{E[Xi] = µ，分散：V[ ¯X] =} n¹²

Pn

i=1 ^{V[Xi] = σ}

2

n

◆ 無作為抽出によって得られる大きさｎの標本は，ｎ個の確率変数によって 表現され，それぞれの確率変数は独立で同じ分布に従う．そしてその分布 は母集団分布である．

復習

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ 母数：母集団を特徴づけるパラメータ

◆ 母集団平均：母集団全体における中心

◆ 母集団分散：母集団全体における中心からの平均的な乖離の程度

■ 推定量：母数を推定するため，標本から作ることができる確率変数

◆ 「標本平均」という確率変数（_{X =}¯ ¹

n

Pn

i=1 ^Xi）：母集団平均を推定したい

◆ 「標本分散」という確率変数（_S² ₌ ¹

n−1

Pn

i=1^{(Xi − ¯X)2）：母集団分散を}

推定したい

■ 不偏推定量：推定したい母数に対して，偏りがない（その推定量の平均＝起き 方の中心が母数と一致）推定量

◆ ^E_{[ ¯X] = µ}

正規分布関連の確率分布に関する前

提知識

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

必要な確率分布

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

d.f. = degrees-of-freedom = ^自由度

■ ^自由度 _{k (}^{パラメータ} _k) ^{のカイ二乗（}_χ²）分布に従う確率変数の密度関数

f (y; k) = ¹

2 · Γ(k/2)

_y 2

_k/2−1

expⁿ₋^y 2

o, k > 0, y > 0.

■ ^自由度 _{m (}^{パラメータ} _m) ^のｔ分布

f (t; m) = _√ ¹

m · B(m/2, 1/2)



1 + ^t

2

m

_−(m+1)/2

, m > 0, −∞ < t < ∞.

■ _{(n, m) ( (n, m)) :}

正規分布

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

-4 -2 0 2 4

0.00.10.20.30.4

y

Pr[ -2.0 < Z < 1.0] ^Pr₍^Z< α₎=^0.95

Pr₍Z_{> α)=}0.05

カイ二乗分布に従う確率変数の密度関数

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

0 5 10 15 20

0.00.10.20.30.40.5

カイ二乗分布に従う確率変数の密度関数

cbind(y1, y2, y3)

df₌2 df₌5 df₌10

ｔ分布に従う確率変数の密度関数

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

-4 -2 0 2 4

0.00.10.20.30.4

t分布に従う確率変数の密度関数

x

cbind(y1, y2, y4)

df₌2 df₌5 N₍0_, 1₎

Ｆ分布に従う確率変数の密度関数

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

0 1 2 3 4 5 6

0.00.20.40.60.81.0

F分布に従う確率変数の密度関数

cbind(y1, y2, y3)

df1₌2 , df2₌5 df1₌10 , df2₌20 df1₌20 , df2₌2

正規分布との関係：カイ二乗分布

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ 標準正規分布に従う確率変数 _X：X ∼ N(0, 1)

■ ^{カイ二乗（}_χ^2）分布

◆ _{Y = X}^{2 とおくと，}_Y は自由度１のカイ二乗（_χ²）分布に従う

◆ 独立な標準正規分布に従う確率変数 _X1, . . . , X_n ^を用いて

Y =

n

X

i=1

X_i²

としたとき，_Y は自由度 _n の _χ² 分布に従う

正規分布との関係：ｔ分布

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ ^ｔ分布

◆ _X: ^{標準正規分布／} _{W :} ^自由度 _m ^の _χ^{2 分布}

◆ _X ^と _W が互いに独立であるとき，自由度 _m のｔ分布に従う． X

pW/m ^{∼ t(m)}

◆ ^自由度 _m が大きいとき，標準正規分布に近づく．

正規分布との関係：Ｆ分布

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ ^Ｆ分布

◆ _{V :} ^自由度 _n ^の _χ^{2 分布／} _{W :} ^自由度 _m ^の _χ^{2 分布}

◆ _V ^と _W が互いに独立であるとき，自由度 _{(n, m)} のＦ分布に従う． V /n

W/m ^{∼ F (n, m)}

◆ ^自由度が _{(1, m)} のＦ分布に従う確率変数の平方根は，自由度 _m のｔ分布に

従うということもその構成方法から明らか．

正規母集団における標本理論

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

母集団分布が正規分布

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ ^{母集団分布：正規分布}

■ ^{母数：母集団平均} _µ ^{／ 母集団分散} _σ²

■ ^大きさ _n ^の標本：_{{X1, X2}, . . . , X_n}^．ただし X_{i ∼ N(µ, σ}²) i = 1, 2, . . . , n

■ ^{標本平均：}

X =¯ ¹ n

n

X

i=1

X_i

■ ^{標本分散：}

S² = ¹ n − 1

n

X

i=1

(X_{i − ¯}X)²

標本平均の標本分布

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ ^{標本平均の分布：}_{X = n}^{¯ −1 Pn}_{i=1 Xi}

X ∼ N( µ, σ¯ ^{2/n ) (1)}

■ ここから（標準化によって） X − µ¯

pσ²/n ^{∼ N(0, 1) (2)}

標本平均の問題

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ 天秤による計測：観測には多少の誤差が出る＝誤差を確率変数のように考える

◆ ^{観測誤差は平均} ₀ ^ｇ，分散 _0.1^{（標準偏差 √}_0.1 ^{ｇ）の正規分布}

■ ₁₀ 回の計測を行って，その平均によって「重さ」を確定する．

■ ^{このとき，}「重さ」の取りうる値はどのような分布になるか（標本平均の分布）？ また，真の重さが ₁₀₀ ｇでとき，_{| ¯}X − 100| > 0.3 ^{となる確率は？}

標本平均の問題

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ ^{計測量に関する大きさ} ₁₀ ^の標本：_{{X1, X2}, . . . , X_10}

■ ^{誤差を確率変数} _Ui ^{と表記するとき，}_Ui ∼ N(0, 0.1)

■ ^{計測量の取りうる値を} _Xi ^とすると

X_i = 100 + U_i, i = 1.2. . . . , 10

なので _Xi ∼ N(100, 0.1)^{（真の重さが 100} ｇ，誤差が平均０，分散 _0.1）

■ ^{このとき（}₁^）式より

X =¯ ¹ 10

10

X

i=1

X_i ∼ N(100, 0.1/10)

√

標本平均の問題

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ ^{求めるべき条件式} _{| ¯}X − 100| > 0.3 ^{を Z を含む様に変形する．}

| ¯X − 100| > 0.3 ⇔ 10 · | ¯X − 100| > 3 ⇔ |Z| > 3

■ 標準正規分布の数表より

Pr[Z > 3.00] ≈ 0.00135, Pr[|Z| > 3.00] ≈ 0.0027

標本平均に関する問題２

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ 天秤による計測：観測誤差は平均 ₀ ｇ，分散 _0.1 ｇ（標準偏差

√0.1 ^{ｇ）の正規}

分布

■ _X^{¯ の誤差が} _0.1 ^{ｇ以下となる確率が} ₉₀ ％以上となるようにしたい

■ 何回くらい計測すればよいか？

■ ^{計測量に関する大きさ} _n ^の標本：_{{X1, X2}, . . . , X_n}

■ ^{このとき（}₂^{）式を参照して，}

Z ≡ ^{X − 100¯} p0.1/n ⁼

√10n · ( ¯X − 100) ∼ N(0, 1)

Pr[| ¯X − 100| < 0.1] ≥ 0.90 ⁿ

標本平均に関する問題２

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

| ¯X − 100| < 0.1 ⇔ ^√10n · | ¯X − 100| < 0.1^√10n ⇔ |Z| < 0.1^√10n

■ 標準正規分布の数表より

Pr[|Z| < 1.65] = 0.90

なので臨界点を _1.65 よりも大きくすると題意を満たす．つまり 0.1^√10n > 1.65 −→ n ≥ ^16.5

2

10 ^{= 27.225}

つまり _{n = 28} 以上に設定すればよい．

Pr[|Z| < c] = 0.90

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

-4 -2 0 2 4

0.00.10.20.30.4

Pr[ -1.65 < Z < 1.65 ]

y1

µ =^{0 ,} σ²=¹

Pr[Z>1.65] = 0.050 Pr[Z<-1.65] = 0.050

標本分散の標本分布

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ ^{標本分散：}_S² _{= (n − 1)}^{−1 Pn}_{i=1(Xi − ¯X)}²

■ ^{標本分散は自由度} _{n − 1} ^{のカイ二乗（}_χ^{2）分布に従う}

(n − 1) · ^S

2

σ^{2 =}

n

X

i=1

 X_{i − ¯}X σ

²

∼ χ²(n − 1) ⁽³⁾

c.f.

n

X

i=1

 X_{i − µ} σ

2

=

n

X

i=1

Z_i², Z_{i ∼ N(0, 1)}

■ ^自由度が _n ^{ではなく，}_{n − 1} ^{であることに注意．}

■ ^データは _n 個あって，本来これらがすべてどのような値でもとりうるはずだが， X¯ による縛りが１つあるので，_X¯ が与えられた下で「自由に」動けるのは

n − 1 ^個．

■ 証明は厄介なので省略．

標本分散に関する問題

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ 母集団分布は正規分布／母集団平均 ₄，母集団分散 ₁₅．

■ ^大きさ ₁₀ ^の標本：_{{X1, X2}, . . . , X_n}^（X_{i ∼ N(4, 15)}^）

■ ^{標本分散について} _S² _{> a} ^{となる確率が} _0.05 ^{となるような} _a ^{はいくつか？}

■ ^（₃）式より，以下の関係が成り立つ (n − 1) · ^S

2

σ^{2 = 9 ·} S²

15 ^{∼ χ}

2₍₉₎

■ ^また

S² > a ⇔ 9 · ^S

2

> 9 · ^a

標本分散に関する問題

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ カイ二乗分布表の自由度９の欄の５％点をみると Pr[V > 16.9190] = 0.05 ^なの でなので臨界点を _16.9190 とすると題意を満たす．つまり

9 · ^a

15 = 16.9190 −→ a = 16.9190 · 15

9 ^{≈ 28.19833}

Pr[V > c] = 0.05

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

0 5 10 15 20

0.000.020.040.060.080.10

自由度９のカイ二乗分布: Pr[ V > 16.9190 ]

y1

df₌9

Pr[V>16.9] = 0.050

標本平均・標本分散に関して

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

X − µ¯

pσ²/n ^{∼ N(0, 1),}

X − µ¯

pS²/n ^{∼ t(n − 1)}

■ ^（₂^）式は _µ ^にも _σ^{2 にも依存．}

◆ _µ についてだけ調べたいときには _σ² が厄介．

◆ ^何とか _σ² に依存しないようにできないか（_S² で置換）？

◆ _X^{¯ と} _S² は独立になることを利用すると次の定理が成立する．

標本平均・標本分散に関して

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

X − µ¯

pσ²/n ^{∼ N(0, 1),}

X − µ¯

pS²/n ^{∼ t(n − 1)}

■ ^{重要な結果１：}_{( ¯X − µ)/pσ}²_/n^{：標準正規分布}

■ ^{重要な結果２：}_{(n − 1)S}²_/σ^{2 ：自由度} _{n − 1} ^の _χ^{2 分布}

■ ^{したがって}

T = ^{( ¯X − µ)/pσ}

2_/n

p{(n − 1)S^2/σ2}/(n − 1) ⁼

X − µ¯

pS²/n ^{∼ t(n − 1) (4)}

■ _{2 σ}² _S²

標本平均・標本分散の問題

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ 母集団分布：正規分布／母集団平均 ₃，母集団分散 _σ²

■ ^大きさ ₁₅ ^の標本：_{{X1, X2}, . . . , X_15}

■ 標本平均と標本分散を用いて，_{( ¯X − 3)/}

√S² > a ^{となる確率が} 0.01 ^となる a ^は いくらか？

■ ^{問題から大きさ} ₁₅ ^の標本：_{{X1, X2}, . . . , X_15} ^{において，}X_{i ∼ N(3, σ}²)

■ このとき標本平均と標本分散は X =¯ ¹

15

15

X

i=1

X_{i ∼ N(3, σ}²/15), S² = ¹ 15 − 1

15

X

i=1

(X_{i − ¯}X)²

であり，（₄）式より

標本平均・標本分散の問題

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

T ≡ ^{X − 3¯}

pS²/15 ^{∼ t(14)}

■ ^ところで

X − 3¯

√S^{2 > a ⇔}

X − 3¯ pS²/15 ^>

a

p1/15 ^{= a}

√15

つまり

Pr

 _¯

X − 3

√S^{2 > a}



= Pr

"

X − 3¯

pS²/15 ^{> a}

√15

#

= 0.01

Pr[T > c] = 0.05

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

-4 -2 0 2 4

0.00.10.20.30.4

自由度14のｔ分布: Pr[ T > 2.624 ]

x

y1

df₌14

Pr[T>2.624] = 0.010

正規分布とｔ分布の乖離

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ ｎ回の計測で得られるであろう結果が _{{X1, X2}, . . . , X_n}

■ ₆ 回の計測で標本平均と母集団平均の乖離が，標本標準偏差の _1.0 倍となる確率 はいくらか

■ ₆ 回の計測で標本平均と母集団平均の乖離が，（母集団）標準偏差の _1.0 倍となる 確率はいくらか

◆ ^{母集団分散} _σ^{2 は既知とする}

正規分布とｔ分布の乖離

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ ^前提：

T ≡ ^{X − µ¯}

pS²/n ∼ t(n − 1), Z ≡ ^{X − µ¯}

pσ²/n ^{∼ N(0, 1)}

■ ６回の計測で標本平均と母集団平均の乖離が，標本標準偏差の _1.0 倍となる確率 Pr[| ¯X − µ| ≥ ^{√S2] = Pr}

"

| ¯X − µ| pS²/n ^≥

√S² pS²/n

#

= Pr_{|T | ≥} ^√n^

n = 6 なのでｔ分布表の自由度５の欄から

√6 ≈ 2.45 ^{となる確率は５％∼ １０}

％の間．（正確には _5.78 ％くらい）

正規分布とｔ分布の乖離

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ ^前提：

T ≡ ^{X − µ¯}

pS²/n ∼ t(n − 1), Z ≡ ^{X − µ¯}

pσ²/n ^{∼ N(0, 1)}

■ 母集団分散が分かっている _σ² として，₆ 回の計測で標本平均と母集団平均の乖 離が，（母集団）標準偏差の _1.0 倍となる確率はいくらか

Pr[| ¯X − µ| ≥ ^{√σ2] = Pr}

"

| ¯X − µ| pσ²/n ^≥

√σ² pσ²/n

#

= Pr_{|Z| ≥} ^√n^

n = 6^（√_{6 ≈ 2.45}）のとき正規分布表から求める確率は約 _0.014

（= 0.007 × 2^）．

本日の演習問題

復習 正規分布関連の確率分布に関する前提知識 正規母集団における標本理論

■ ^{母集団平均が} _µ 分（母集団分散・標準偏差については未知）の正規分布であっ たとするならば、「調査による平均視聴時間」と「母集団における平均視聴時 間」が標本標準偏差の１倍以上も乖離してしまうという結果が出る確率が _0.05 以下で収まる¹のは，何人について調査を実施する場合か？

■ ^{標本標準偏差の} _0.5 倍以上も乖離してしまうという結果が出る確率が _0.05 以下 となるのは，何人について調査を実施する場合か？

■ 上の結果は母集団分散が既知と考えた場合に比べてどの程度異なっているか？

1

調査結果が，ここで想定している調査誤差の許容範囲（標本標準偏差の１倍）を超えて しまう確率が５％以内であると述べている