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.正規分布ふつう
標準正規曲線の下で、z の右側の色つき面積が 0.8997 になるように z をドラッグします。z を小数第 2 位まで求めなさい。
A = 0.8997
z =
2.正規分布むずかしい
穴の直径 D は平均 25.05 mm、標準偏差 0.01 mm の正規分布に従います。軸の直径 d は D と独立で、平均 25 mm、標準偏差 0.02 mm の正規分布に従います。無作為に選んだ軸は D − d > 0 のとき穴にはまります。はまる確率 P を求めなさい。z は小数第 2 位に丸めなさい。
穴: μ = 25.05 mm, σ = 0.01 mm軸: μ = 25 mm, σ = 0.02 mm
P =
3.信頼区間ふつう
12 人の患者の血圧を治療の前後に測定します。差 d = 前 − 後(mmHg)の平均は −4.4、標準偏差は 1.6 です。 差の平均の 90% 信頼区間の誤差の幅 E を求めなさい。
n = 12d̄ = −4.4sd = 1.690%t(11) = 1.796
E = mmHg
4.信頼区間むずかしい
予備調査によると、製品の約 0.7 が不良品です。 割合の 99% 信頼区間の誤差の幅を 0.01 以下にするには、製品をいくつ n 検査する必要がありますか?
面積から z へ戻る:Φ(z) は z の左側の面積です。左側の裾の面積 A なら、表で Φ(z) = A となる z を見つけます。右側の裾の面積 A なら Φ(z) = 1 − A です。−z から z の間の中央の面積 A なら、各裾に (1 − A)/2 が残るので、Φ(z) = (1 + A)/2 です。値 x には x = μ + zσ を使います。
z の左側の面積
Φ(z) = 1 − 0.8997 = 0.1003
z 表を逆に読む
表の中で 0.1003 を見つけ、z を読み取ります。
z = −1.28
答え
z = −1.28
2.正規分布
答え: P = 0.9875
正規変数の和と差:X と Y が独立で正規分布に従うなら、X + Y と X − Y も正規分布に従います。平均は足したり引いたりします:E(X ± Y) = μX ± μY。分散は差の場合でも常に足します:Var(X ± Y) = σX² + σY²、したがって σ = √(σX² + σY²) です。標準偏差を直接足すことはありません。同じ分布に従う独立な n 個の合計は、平均 n·μ、標準偏差 σ·√n です。そのあとはいつもどおり標準化します:z = (x − 平均)/σ とし、Φ(z) を読み取ります。
差の分布
μ = 25.05 − 25 = 0.05 mm, σ = 0.02236 mm
標準化する
z = (0 − 0.05)/0.02236 = −2.24
確率
P = 1 − Φ(−2.24) = 0.9875
答え
P = 0.9875
3.信頼区間
答え: E = 0.8295 mmHg
Interval for a difference: estimate ± critical value × standard error. For two independent means with known σ: SE = √(σ₁²/n₁ + σ₂²/n₂) and z. With equal unknown variances: pooled s_p and t with n₁ + n₂ − 2 degrees of freedom. For paired data, work with the differences: d̄ ± t·s_d/√n.
臨界値
df = 11, t = 1.796
誤差の幅
E = 1.796·1.6/√12 = 0.8295
答え
E = 0.8295 mmHg
4.信頼区間
答え: n = 13936
Confidence interval for a proportion: p̂ ± z·√(p̂(1 − p̂)/n), where p̂ = x/n. To reach a margin E you need n ≥ z²·p(1 − p)/E², rounded up; without a prior estimate use p = 0.5, which gives the largest n.
n について解く
n = 2.576²·0.7·0.3/0.01² = 13935
切り上げる
n = 13936
答え
n = 13936
5.仮説検定と信頼区間
答え: p = 0.04648
Hypothesis test: the test statistic = (estimate − value under H₀)/(standard error under H₀). The p-value is the probability, assuming H₀ is true, of a statistic at least as extreme as the one observed: the area in one tail for a one-sided H₁, twice the tail area for a two-sided H₁. Reject H₀ when p ≤ α.
検定統計量
SE = σ/√n = 16/√64 = 2
z = (503.36 − 500)/2 = 1.68
裾の面積
Φ(1.68) = 0.9535
1 − Φ(1.68) = 0.0465
p 値
p = 0.04648
答え
p = 0.04648
p ≤ 0.05 なので、有意水準 5% で H₀ は棄却されます。
6.仮説検定と信頼区間
答え: L = 96.85
平均の検定:検定統計量は、標本平均が H₀ で主張された値から標準誤差何個分離れているかを表します:σ が既知なら z = (x̄ − μ₀)/(σ/√n)、標本から推定するなら t = (x̄ − μ₀)/(s/√n)。信頼区間は x̄ ±(臨界値)·(標準誤差)です。
誤差の幅
E = 1.96·10/√64 = 2.45
区間の端
L = 99.3 − 2.45 = 96.85
答え
L = 96.85
7.線形回帰と相関
答え: r = 0.6815
最小二乗回帰: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 の変動のうち直線で説明される割合です。
偏差平方和・偏差積和
Sxy = 1014 − 21·285/6 = 16.5
Sxx = 91 − 21²/6 = 17.5
Syy = 13571 − 285²/6 = 33.5
相関係数
r = 16.5/√(17.5·33.5) = 0.6815
答え
r = 0.6815
8.線形回帰と相関
答え: ŷ = 19.92
最小二乗回帰: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
ȳ = 137/6 = 22.83
平方和と積和
Sxx = Σ(x − x̄)² = 17.5
Sxy = Σ(x − x̄)(y − ȳ) = −25.5
回帰直線
b = −25.5/17.5 = −1.457
a = 22.83 − (−1.457)·3.5 = 27.93
ŷ = 27.93 − 1.457·x
予測
ŷ = 27.93 + (−1.457)·5.5 = 19.92
答え
ŷ = 19.92
9.カイ二乗検定
答え: E = 22.82
カイ二乗検定:観測度数 O と H₀ のもとでの期待度数 E を χ² = Σ(O − E)²/E で比べます。適合度検定:各カテゴリーで E = n·p、自由度 df = k − 1(データから推定したパラメータ 1 つごとにさらに 1 を引く)。r × c 表での独立性の検定:E = 行の合計 × 列の合計 / 総計、df = (r − 1)(c − 1)。χ² が選んだ α での臨界値を超えたら H₀ を棄却します。
期待度数
E = 54·41/97 = 22.82
答え
E = 22.82
10.カイ二乗検定
答え: df = 6
カイ二乗検定:観測度数 O と H₀ のもとでの期待度数 E を χ² = Σ(O − E)²/E で比べます。適合度検定:各カテゴリーで E = n·p、自由度 df = k − 1(データから推定したパラメータ 1 つごとにさらに 1 を引く)。r × c 表での独立性の検定:E = 行の合計 × 列の合計 / 総計、df = (r − 1)(c − 1)。χ² が選んだ α での臨界値を超えたら H₀ を棄却します。