統計学 II:統計的推測:直前対策パック

公式集、全トピックの解法の要点、10 問の演習問題。完全な解答は巻末にあります。

新しい問題科目に戻る

印刷ウィンドウで「PDF として保存」を選ぶと、コピーを端末に保存できます。

1公式と解法

正規分布

z スコア
z = (x − μ)/σ
x 未満
P(X < x) = Φ(z)
x より大きい
P(X > x) = 1 − Φ(z)
2 つの値の間
Φ(z2) − Φ(z1)
パーセンタイル
x = μ + z·σ

解き方

  1. 値を z スコアに変換します:z = (x − μ)/σ。
  2. z より左側の面積 Φ(z) を表で調べます。
  3. 右側の面積は 1 − Φ(z)、区間の面積は 2 つの左側の面積の差です。
  4. パーセンタイルでは逆にたどります:面積から z を求め、x = μ + z·σ とします。

信頼区間

σ が既知
x̄ ± z·σ/√n
σ が未知
x̄ ± t·s/√n, df = n − 1
z の臨界値
90%: 1.645, 95%: 1.96, 99%: 2.576
標本サイズ
n = (z·σ/E)², rounded up

解き方

  1. z(母集団の σ が既知)と t(標本の s しかなく、自由度 n − 1)のどちらを使うか判断します。
  2. 信頼水準に対応する臨界値を求めます。
  3. 誤差の幅:E = 臨界値 × 標準偏差/√n。
  4. 区間:x̄ ± E。

仮説検定と信頼区間

平均の z 検定
z = (x̄ − μ0)/(σ/√n)
平均の t 検定
t = (x̄ − μ0)/(s/√n), df = n − 1
有意水準 5% の両側検定
|z| > 1.96 なら H0 を棄却
信頼区間
x̄ ± z·σ/√n

解き方

  1. H₀ と H₁ を立て、片側検定か両側検定かを決めます。
  2. 検定統計量を計算します:(x̄ − μ₀)/(σ/√n)、または s と t を使います。
  3. 臨界値と比較するか、p 値を求めます。
  4. 統計量が臨界値を超えていれば(p が α より小さければ)H₀ を棄却し、結論を言葉で述べます。

線形回帰と相関

平方和
Sxx = Σ(x − x̄)² = Σx² − (Σx)²/n, Sxy = Σxy − ΣxΣy/n
傾き
b = Sxy/Sxx
切片
a = ȳ − b·x̄
相関
r = Sxy/√(SxxSyy)
残差
e = y − ŷ

カイ二乗検定

検定統計量
χ² = Σ(O − E)²/E
適合度
E = n·p, df = k − 1 − (estimated parameters)
独立性
E = row total × column total / grand total, df = (r − 1)(c − 1)
臨界値(α = 0.05)
df 1: 3.841, 2: 5.991, 3: 7.815, 4: 9.488
判断
χ² > χ²crit なら H0 を棄却

2演習問題

1.正規分布ふつう
400 g403.8 gT = X₁ + … + X₁₀?

ボルト の質量は互いに独立で、平均 40 g、標準偏差 1.2 g の正規分布に従います。そのうち 10 個の合計質量が 403.8 g を超える確率 P を求めなさい。z は小数第 2 位に丸めなさい。

μ = 40 gσ = 1.2 gn = 10> 403.8 g
P =
2.正規分布むずかしい
μ = 32x = ?5%

配達時間は、平均 32 min、標準偏差 6 min の正規分布に従います。全体の 5% がそれより下にある値 x(第 5 パーセンタイル)を求めなさい。

μ = 32σ = 6P = 5%
x = min
3.信頼区間ふつう
n = 25, s = 10.5LU = ?x̄ = 88E99%

25 個からなる標本の平均は 88、標準偏差は s = 10.5 です。母集団はほぼ正規分布に従います。 平均の 99% 信頼区間の上限を求めなさい。

n = 25x̄ = 88s = 10.599%t(24) = 2.797
上限 =
4.信頼区間むずかしい
n = 11σ = ?Us = 3.295%

機械加工した 11 本の軸の直径の標本標準偏差は 3.2 µm です。直径が正規分布に従うと仮定して、母標準偏差の 95% 信頼区間の下限 σ を求めなさい。df = 10 のとき、χ² 分布の 2.5% は 3.247 より下に、2.5% は 20.483 より上にあります。

n = 11s = 3.2 µmχ² = 3.247, 20.483
σ = µm
5.仮説検定と信頼区間ふつう
−2.2012.201t = ?df = 11

12 人の従業員が研修の前後にタイピングテストを受けます。各人について差 d = 後 − 前(1 分あたりの単語数)を記録します。差の平均は 10.42、差の標準偏差は 9.5 です。H₀: μd = 0 の対応のある検定統計量 t を求めなさい。

n = 12d̄ = 10.42sd = 9.5
t =
6.仮説検定と信頼区間むずかしい
μ₀x̄c

H₀: μ = 500 対 H₁: μ < 500 の片側 z 検定を、α = 0.05、標本サイズ 36 で行います。標準偏差は 15 であることがわかっています。実際の平均は 497.5 です。 検定の検出力 1 − β をパーセントで求めなさい。

μ₀ = 500μ = 497.5σ = 15n = 36α = 0.05
1 − β = %
7.線形回帰と相関ふつう
xy125332533738940xy525303540

表は x = 広告費(k€) と y = 売上高(k€) の 5 組の対になった観測値を示しています。 最小二乗直線の決定係数 r² を求めなさい。

n = 5
r² =
8.線形回帰と相関むずかしい
xy140241344447551655xy24640455055x = 1.5

表は x = 広告費(k€) と y = 売上高(k€) の 6 組の対になった観測値を示しています。 最小二乗直線 ŷ = a + b·x を用います。 新しい観測値 (1.5, 36) が得られました。その残差 e = y − ŷ を求めなさい。

n = 6(1.5, 36)
e =
9.カイ二乗検定ふつう
欠陥なし軽微重大合計日勤26193580夜勤642610100合計904545180

分割表は観測度数を示しています。行の変数と列の変数の独立性の検定について、カイ二乗統計量 χ² を計算しなさい。連続性の補正は使いません。

N = 1802 × 3
χ² =
10.カイ二乗検定むずかしい
dfα = 0.05α = 0.0113.8416.63525.9919.21037.81511.34549.48813.277511.07015.086612.59216.812714.06718.475815.50720.090916.91921.6661018.30723.209

2 行 3 列(合計を除く)の分割表について、カイ二乗独立性検定を行います。 表を使って、有意水準 α = 0.01 での臨界値 χ²crit を求めなさい。

α = 0.01
χ²crit =

3溶液

1.正規分布

答え: P = 0.1587

正規変数の和と差:X と Y が独立で正規分布に従うなら、X + Y と X − Y も正規分布に従います。平均は足したり引いたりします:E(X ± Y) = μX ± μY。分散は差の場合でも常に足します:Var(X ± Y) = σX² + σY²、したがって σ = √(σX² + σY²) です。標準偏差を直接足すことはありません。同じ分布に従う独立な n 個の合計は、平均 n·μ、標準偏差 σ·√n です。そのあとはいつもどおり標準化します:z = (x − 平均)/σ とし、Φ(z) を読み取ります。

合計の分布

μT = 400 g, σT = 1.2·√10 = 3.795 g

標準化する

z = (403.8 − 400)/3.795 = 1

確率

P = 1 − Φ(1) = 0.1587

答え

P = 0.1587
2.正規分布

答え: x = 22.13 min

正規分布:値 x を z スコア z = (x − μ)/σ に変換します。これは平均から標準偏差何個分離れているかを表します。z 表(または Φ(z))は z より左側の面積、つまり P(X < x) を与えます。右側の面積は 1 − Φ(z)、2 つの値の間の面積はそれぞれの左側の面積の差です。

パーセンタイルに対応する z

Φ(z) = 0.05 → z = −1.645

x に戻す

x = 32 + −1.645·6 = 22.13 min

答え

x = 22.13 min
3.信頼区間

答え: 上限 = 93.87

平均の信頼区間:x̄ ± 誤差の幅。母集団の σ が既知なら誤差の幅は z·σ/√n、標本標準偏差 s しかなければ自由度 n − 1 の t を使って t·s/√n です。臨界値:z = 1.645(90%)、1.96(95%)、2.576(99%)。誤差の幅を E 以下にするには n ≥ (z·σ/E)² が必要で、端数は切り上げます。

臨界値

df = n − 1 = 24, t = 2.797

誤差の幅

E = 2.797·10.5/√25 = 5.874

区間

88 ± 5.874 → [82.13, 93.87]

答え

上限 = 93.87
4.信頼区間

答え: σ = 2.236 µm

分散の区間:(n − 1)s²/σ² は自由度 n − 1 の χ² 分布に従います。σ² の 95% 区間は、(n − 1)s² を上側の χ² 値(その上に 2.5%)で割った値から、下側の χ² 値(その下に 2.5%)で割った値までです。σ の区間は両方の限界の平方根です。s に関して対称ではありません。

平方和

(n − 1)s² = 10·3.2² = 102.4

σ² の限界

102.4/20.483 = 4.999
102.4/3.247 = 31.54

σ の限界

σ = √4.999 = 2.236 µm

答え

σ = 2.236 µm

σ の 95% 区間:[2.236, 5.616] µm。

5.仮説検定と信頼区間

答え: t = 3.8

Testing a mean: the test statistic measures how many standard errors the sample mean lies from the value claimed in H₀: z = (x̄ − μ₀)/(σ/√n) when σ is known, t = (x̄ − μ₀)/(s/√n) when it is estimated from the sample. A confidence interval is x̄ ± (critical value)·(standard error).

差の平均の標準誤差

SE = 9.5/√12 = 2.742

検定統計量

t = 10.42/2.742 = 3.8

答え

t = 3.8

df = 11 のとき、両側 5% の臨界値は 2.201 です。|t| の方が大きいので、H₀ は棄却されます。

6.仮説検定と信頼区間

答え: 1 − β = 25.95 %

検定の検出力:検出力 = 1 − β は、H₀ が誤っているときに H₀ を棄却する確率です。μ₀ と z_α から境界値 x̄c を求め、真の平均を使って x̄c を超える確率を計算します:z = (x̄c − μ)/(σ/√n)。

臨界標本平均

SE = 15/√36 = 2.5
x̄c = 500 − 1.645·2.5 = 495.9

真の平均のもとでの z

z = (495.9 − 497.5)/2.5 = −0.645

電力

検出力 = P(Z < −0.645) = 25.95%
1 − β = 25.95%

答え

1 − β = 25.95%
7.線形回帰と相関

答え: r² = 0.9446

最小二乗回帰:x̄ と ȳ を平均とし、Sxx = Σ(x − x̄)²、Syy = Σ(y − ȳ)²、Sxy = Σ(x − x̄)(y − ȳ) とすると、直線 ŷ = a + b·x の傾きは b = Sxy/Sxx、切片は a = ȳ − b·x̄ なので、この直線は必ず (x̄, ȳ) を通ります。ピアソンの相関係数 r = Sxy/√(Sxx·Syy) は −1 から 1 の間の値をとります。r² は y の変動のうち直線で説明される割合です。

平均

x̄ = 25/5 = 5
ȳ = 168/5 = 33.6

平方和と積和

Sxx = Σ(x − x̄)² = 40
Sxy = Σ(x − x̄)(y − ȳ) = 72
Syy = Σ(y − ȳ)² = 137.2

決定係数

r² = 72²/(40·137.2) = 0.9446

答え

r² = 0.9446

y の変動の約 94.46% が直線で説明されます。

8.線形回帰と相関

答え: e = −4.162

最小二乗回帰:x̄ と ȳ を平均とし、Sxx = Σ(x − x̄)²、Syy = Σ(y − ȳ)²、Sxy = Σ(x − x̄)(y − ȳ) とすると、直線 ŷ = a + b·x の傾きは b = Sxy/Sxx、切片は a = ȳ − b·x̄ なので、この直線は必ず (x̄, ȳ) を通ります。ピアソンの相関係数 r = Sxy/√(Sxx·Syy) は −1 から 1 の間の値をとります。r² は y の変動のうち直線で説明される割合です。

平均

x̄ = 21/6 = 3.5
ȳ = 278/6 = 46.33

平方和と積和

Sxx = Σ(x − x̄)² = 17.5
Sxy = Σ(x − x̄)(y − ȳ) = 54

回帰直線

b = 54/17.5 = 3.086
a = 46.33 − 3.086·3.5 = 35.53
ŷ = 35.53 + 3.086·x

予測値

ŷ = 35.53 + 3.086·1.5 = 40.16

残差

e = 36 − 40.16 = −4.162

答え

e = −4.162

点は直線より下にあります。

9.カイ二乗検定

答え: χ² = 29.16

カイ二乗検定:観測度数 O と H₀ のもとでの期待度数 E を χ² = Σ(O − E)²/E で比べます。適合度検定:各カテゴリーで E = n·p、自由度 df = k − 1(データから推定したパラメータ 1 つごとにさらに 1 を引く)。r × c 表での独立性の検定:E = 行の合計 × 列の合計 / 総計、df = (r − 1)(c − 1)。χ² が選んだ α での臨界値を超えたら H₀ を棄却します。

期待度数

日勤: 80·90/180 = 40, 80·45/180 = 20, 80·45/180 = 20
夜勤: 100·90/180 = 50, 100·45/180 = 25, 100·45/180 = 25

カイ二乗統計量

χ² = (26 − 40)²/40 + (19 − 20)²/20 + (35 − 20)²/20 + (64 − 50)²/50 + (26 − 25)²/25 + (10 − 25)²/25
χ² = 29.16

答え

χ² = 29.16

df = 2。α = 0.05 での臨界値は 5.991 なので、H₀(独立)は棄却される。

10.カイ二乗検定

答え: χ²crit = 9.21

カイ二乗検定:観測度数 O と H₀ のもとでの期待度数 E を χ² = Σ(O − E)²/E で比べます。適合度検定:各カテゴリーで E = n·p、自由度 df = k − 1(データから推定したパラメータ 1 つごとにさらに 1 を引く)。r × c 表での独立性の検定:E = 行の合計 × 列の合計 / 総計、df = (r − 1)(c − 1)。χ² が選んだ α での臨界値を超えたら H₀ を棄却します。

自由度

df = (2 − 1)·(3 − 1) = 2

表

χ²crit(df = 2, α = 0.01) = 9.210

答え

χ²crit = 9.210

計算した χ² がこの値より大きければ H₀ を棄却します。

ExamProblems · examproblems.com/courses/statistics-2-inference