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.正規分布ふつう
配達時間は、平均 32 min、標準偏差 6 min の正規分布に従います。24.8 から 41.3 min の間にある割合はいくらですか?P を小数で答えなさい。
μ = 32σ = 6x1 = 24.8x2 = 41.3
P =
2.正規分布むずかしい
穴の直径 D は平均 25.04 mm、標準偏差 0.02 mm の正規分布に従います。軸の直径 d は D と独立で、平均 25 mm、標準偏差 0.01 mm の正規分布に従います。無作為に選んだ軸は D − d > 0 のとき穴にはまります。はまる確率 P を求めなさい。z は小数第 2 位に丸めなさい。
穴: μ = 25.04 mm, σ = 0.02 mm軸: μ = 25 mm, σ = 0.01 mm
P =
3.信頼区間ふつう
20 個からなる標本の平均は 118.5、標準偏差は s = 16 です。母集団はほぼ正規分布に従います。 平均の 90% 信頼区間の誤差の幅 E を求めなさい。
n = 20x̄ = 118.5s = 1690%t(19) = 1.729
E =
4.信頼区間むずかしい
母標準偏差は 5 です。平均の 90% 信頼区間の誤差の幅を 1.5 以下にするには、どれだけの標本サイズ n が必要ですか?
σ = 5E = 1.590%
n =
5.仮説検定と信頼区間ふつう
2 つの生産ラインで同じ部品をつくっています。ライン 1:40 個の部品の標本の平均長さは 55.5 mm で、ラインの既知の標準偏差は 8 mm です。ライン 2:60 個の部品の標本の平均は 52.3 mm、既知の標準偏差は 18 mm です。差を x̄₁ − x̄₂(ライン 1 − ライン 2)として、H₀: μ₁ = μ₂ の検定統計量 z を求めなさい。
正規分布:値 x を z スコア z = (x − μ)/σ に変換します。これは平均から標準偏差何個分離れているかを表します。z 表(または Φ(z))は z より左側の面積、つまり P(X < x) を与えます。右側の面積は 1 − Φ(z)、2 つの値の間の面積はそれぞれの左側の面積の差です。
z スコア
z₁ = (24.8 − 32)/6 = −1.2
z₂ = (41.3 − 32)/6 = 1.55
間の面積
P = Φ(1.55) − Φ(−1.2) = 0.9394 − 0.1151 = 0.8244
答え
P = 0.8244
2.正規分布
答え: P = 0.9633
正規変数の和と差: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.04 − 25 = 0.04 mm, σ = 0.02236 mm
標準化する
z = (0 − 0.04)/0.02236 = −1.79
確率
P = 1 − Φ(−1.79) = 0.9633
答え
P = 0.9633
3.信頼区間
答え: E = 6.186
平均の信頼区間: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 = 19, t = 1.729
誤差の幅
E = 1.729·16/√20 = 6.186
答え
E = 6.186
4.信頼区間
答え: n = 31
平均の信頼区間: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 = (1.645·5/1.5)² = 30.07
切り上げる
標本サイズは整数であり、切り捨てると誤差の幅が大きくなりすぎます。
n = 31
答え
n = 31
5.仮説検定と信頼区間
答え: z = 1.209
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 = √(8²/40 + 18²/60) = √7 = 2.646
検定統計量
z = (55.5 − 52.3)/2.646 = 1.209
答え
z = 1.209
|z| < 1.96 なので、有意水準 5%(両側)で H₀ は棄却されません。
6.仮説検定と信頼区間
答え: L = 502.5
平均の検定:検定統計量は、標本平均が H₀ で主張された値から標準誤差何個分離れているかを表します:σ が既知なら z = (x̄ − μ₀)/(σ/√n)、標本から推定するなら t = (x̄ − μ₀)/(s/√n)。信頼区間は x̄ ±(臨界値)·(標準誤差)です。
誤差の幅
E = 1.645·10/√25 = 3.29
区間の端
L = 505.8 − 3.29 = 502.5
答え
L = 502.5
7.線形回帰と相関
答え: r = 0.9947
最小二乗回帰: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 = 3085 − 33·535/6 = 142.5
Sxx = 199 − 33²/6 = 17.5
Syy = 48877 − 535²/6 = 1173
相関係数
r = 142.5/√(17.5·1173) = 0.9947
答え
r = 0.9947
8.線形回帰と相関
答え: x = 24.85
最小二乗回帰: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̄ = 147/6 = 24.5
ȳ = 379/6 = 63.17
平方和と積和
Sxx = Σ(x − x̄)² = 157.5
Sxy = Σ(x − x̄)(y − ȳ) = 376.5
回帰直線
b = 376.5/157.5 = 2.39
a = 63.17 − 2.39·24.5 = 4.6
ŷ = 4.6 + 2.39·x
x について解く
x = (64 − 4.6)/2.39 = 24.85
答え
x = 24.85
9.カイ二乗検定
答え: χ² = 6.061
カイ二乗検定:観測度数 O と H₀ のもとでの期待度数 E を χ² = Σ(O − E)²/E で比べます。適合度検定:各カテゴリーで E = n·p、自由度 df = k − 1(データから推定したパラメータ 1 つごとにさらに 1 を引く)。r × c 表での独立性の検定:E = 行の合計 × 列の合計 / 総計、df = (r − 1)(c − 1)。χ² が選んだ α での臨界値を超えたら H₀ を棄却します。