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.正規分布ふつう
X は平均 250 g、標準偏差 4 g の正規分布に従い、Y は X と独立で、平均 500 g、標準偏差 9 g の正規分布に従います。X + Y が 742.1 g を超える確率 P を求めなさい。z は小数第 2 位に丸めなさい。
μX = 250, σX = 4μY = 500, σY = 9X + Y > 742.1 g
P =
2.正規分布むずかしい
穴の直径 D は平均 25.02 mm、標準偏差 0.01 mm の正規分布に従います。軸の直径 d は D と独立で、平均 25 mm、標準偏差 0.01 mm の正規分布に従います。無作為に選んだ軸は D − d > 0 のとき穴にはまります。はまる確率 P を求めなさい。z は小数第 2 位に丸めなさい。
穴: μ = 25.02 mm, σ = 0.01 mm軸: μ = 25 mm, σ = 0.01 mm
正規変数の和と差:X と Y が独立で正規分布に従うなら、X + Y と X − Y も正規分布に従います。平均は足したり引いたりします:E(X ± Y) = μX ± μY。分散は差の場合でも常に足します:Var(X ± Y) = σX² + σY²、したがって σ = √(σX² + σY²) です。標準偏差を直接足すことはありません。同じ分布に従う独立な n 個の合計は、平均 n·μ、標準偏差 σ·√n です。そのあとはいつもどおり標準化します:z = (x − 平均)/σ とし、Φ(z) を読み取ります。
和の分布
μ = 750 g, σ = √(4² + 9²) = 9.849 g
標準化する
z = (742.1 − 750)/9.849 = −0.8
上側確率
P = 1 − Φ(−0.8) = 1 − 0.2119 = 0.7881
答え
P = 0.7881
2.正規分布
答え: P = 0.9207
正規変数の和と差: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.02 − 25 = 0.02 mm, σ = 0.01414 mm
標準化する
z = (0 − 0.02)/0.01414 = −1.41
確率
P = 1 − Φ(−1.41) = 0.9207
答え
P = 0.9207
3.信頼区間
答え: 下限 = −1.486 h
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.
点推定値
x̄₁ − x̄₂ = 104 − 102 = 2
標準誤差
SE = √(7²/100 + 20²/100) = 2.119
誤差の幅
E = 1.645·2.119 = 3.486
区間
2 ± 3.486 → [−1.486, 5.486]
答え
下限 = −1.486 h
区間は 0 を含むので、データは有意な差を示していません。
4.信頼区間
答え: n = 612
平均の信頼区間:x̄ ± 誤差の幅。母集団の σ が既知なら誤差の幅は z·σ/√n、標本標準偏差 s しかなければ自由度 n − 1 の t を使って t·s/√n です。臨界値:z = 1.645(90%)、1.96(95%)、2.576(99%)。誤差の幅を E 以下にするには n ≥ (z·σ/E)² が必要で、端数は切り上げます。
n について解く
n = (2.576·24/2.5)² = 611.6
切り上げる
標本サイズは整数であり、切り捨てると誤差の幅が大きくなりすぎます。
n = 612
答え
n = 612
5.仮説検定と信頼区間
答え: t = 1.68
平均の検定:検定統計量は、標本平均が H₀ で主張された値から標準誤差何個分離れているかを表します:σ が既知なら z = (x̄ − μ₀)/(σ/√n)、標本から推定するなら t = (x̄ − μ₀)/(s/√n)。信頼区間は x̄ ±(臨界値)·(標準誤差)です。
標準誤差
SE = 15/√81 = 15/9 = 1.667
検定統計量
t = (102.8 − 100)/1.667 = 1.68
答え
t = 1.68
t0.025; df = 80 = 1.99
|t| < 1.99 なので、有意水準 5%(両側)で H₀ は棄却されません。
6.仮説検定と信頼区間
答え: n = 174
平均の検定:検定統計量は、標本平均が H₀ で主張された値から標準誤差何個分離れているかを表します:σ が既知なら z = (x̄ − μ₀)/(σ/√n)、標本から推定するなら t = (x̄ − μ₀)/(s/√n)。信頼区間は x̄ ±(臨界値)·(標準誤差)です。
n について解く
n = (1.645·20/2.5)² = 173.2
切り上げる
n = 174
答え
n = 174
7.線形回帰と相関
答え: r = −0.9452
最小二乗回帰: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 = 484 − 21·147/6 = −30.5
Sxx = 91 − 21²/6 = 17.5
Syy = 3661 − 147²/6 = 59.5
相関係数
r = −30.5/√(17.5·59.5) = −0.9452
答え
r = −0.9452
8.線形回帰と相関
答え: e = 5.35
最小二乗回帰: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̄ = 30/5 = 6
ȳ = 75/5 = 15
平方和と積和
Sxx = Σ(x − x̄)² = 40
Sxy = Σ(x − x̄)(y − ȳ) = 98
回帰直線
b = 98/40 = 2.45
a = 15 − 2.45·6 = 0.3
ŷ = 0.3 + 2.45·x
予測値
ŷ = 0.3 + 2.45·3 = 7.65
残差
e = 13 − 7.65 = 5.35
答え
e = 5.35
点は直線より上にあります。
9.カイ二乗検定
答え: E = 12.76
カイ二乗検定:観測度数 O と H₀ のもとでの期待度数 E を χ² = Σ(O − E)²/E で比べます。適合度検定:各カテゴリーで E = n·p、自由度 df = k − 1(データから推定したパラメータ 1 つごとにさらに 1 を引く)。r × c 表での独立性の検定:E = 行の合計 × 列の合計 / 総計、df = (r − 1)(c − 1)。χ² が選んだ α での臨界値を超えたら H₀ を棄却します。
期待度数
E = 37·40/116 = 12.76
答え
E = 12.76
10.カイ二乗検定
答え: χ²crit = 5.991
カイ二乗検定:観測度数 O と H₀ のもとでの期待度数 E を χ² = Σ(O − E)²/E で比べます。適合度検定:各カテゴリーで E = n·p、自由度 df = k − 1(データから推定したパラメータ 1 つごとにさらに 1 を引く)。r × c 表での独立性の検定:E = 行の合計 × 列の合計 / 総計、df = (r − 1)(c − 1)。χ² が選んだ α での臨界値を超えたら H₀ を棄却します。