ひよこ統計

母分散の比の検定

統計検定2級

「2つの母分散が等しい」という主張は、σ12=σ22 \sigma_1^2 = \sigma_2^2 とも σ12σ22=1 \frac{\sigma_1^2}{\sigma_2^2} = 1 とも書けます。母分散の比の検定が 帰無仮説 に置くのは後者の形です。比にしておくと、F分布 に従う量 s12/σ12s22/σ22 \frac{s_1^2/\sigma_1^2}{s_2^2/\sigma_2^2} から未知の母分散が約分で消え、標本の不偏分散だけで計算できる統計量が残るからです。

F=s12s22 ∼ F(ν1, ν2),ν1=n1−1,ν2=n2−1F = \frac{s_1^2}{s_2^2} \ \sim \ F(\nu_1,\ \nu_2), \qquad \nu_1 = n_1 - 1, \quad \nu_2 = n_2 - 1
s12 s_1^2 ・s22 s_2^2 :各群の不偏分散(偏差の2乗和を自由度で割った量)、n1 n_1 ・n2 n_2 :各群の標本の大きさ、ν1 \nu_1 ・ν2 \nu_2 :分子・分母それぞれの自由度、F(ν1, ν2) F(\nu_1,\ \nu_2) :自由度 (ν1, ν2) (\nu_1,\ \nu_2) のF分布

帰無仮説が正しければ F F は 1 の近くに散らばり、母分散が違えば 1 から離れます。「ばらつきが違うか」を問う両側検定では α \alpha を上下へ α2 \frac{\alpha}{2} ずつ振り分け、棄却域を F>Fα/2(ν1,ν2) F > F_{\alpha/2}(\nu_1, \nu_2) または F<F1−α/2(ν1,ν2) F < F_{1-\alpha/2}(\nu_1, \nu_2) とします。ところが 上側確率点 の表にあるのは上側の点だけなので、下側の点は自由度を入れ替えた上側点の逆数 F1−α/2(ν1,ν2)=1Fα/2(ν2,ν1) F_{1-\alpha/2}(\nu_1, \nu_2) = \frac{1}{F_{\alpha/2}(\nu_2, \nu_1)} として読み替えます。分子と分母を取り替えた s22s12 \frac{s_2^2}{s_1^2} は自由度 (ν2,ν1) (\nu_2, \nu_1) のF分布に従います。もとの統計量が小さいことを、逆数が大きいことに置き換えられるわけです。

このページで使う等裾の両側検定では、観測した F F の下側確率と上側確率の小さいほうを2倍し、1を上限として p値 を定めます。F分布は非対称なので、常に上側確率だけを2倍したり、1から等距離の点を反対側に取ったりする方法ではありません。この定義なら、p値法と上下に α/2 \alpha/2 ずつ置く棄却域の判定が一致します。

この検定を使う中心的な場面は、測定器の精度や製品の均一さなど、ばらつきの違い自体が関心事のときです。比がどの範囲にあるかまで知りたい場合は 母分散の比の信頼区間 を作ります。同じ標本で先にこの検定を行い、非有意ならプール型、有意なら Welch と切り替える二段階手続きは、最終的な検定の誤りの確率を名目水準どおりに保てないため推奨されません。平均差の方法は研究計画や事前の仮定に基づいて選びます。この検定には、成り立たなければ結論が意味を失う前提と、計算前に決めておく事柄があります。

同じ区間を走るバスの路線A・Bについて、到着の遅れ(分)を調べました。路線Aは5日、路線Bは9日ぶんで、値は 1, 3, 5, 7, 9(A)と 3, 3, 4, 4, 5, 6, 6, 7, 7(B)です。平均の遅れはどちらも 5 分で差がありません。それでもAは 1 分から 9 分まで散らばっており、日によって当たり外れが大きいように見えます。遅れは正規分布に従い、2つの標本は独立とみなして、この散らばりの違いを偶然で片づけてよいかを有意水準5%の両側検定で確かめます。

偏差の2乗和はAが 16+4+0+4+16=40 16+4+0+4+16 = 40 、Bが 4+4+1+1+0+1+1+4+4=20 4+4+1+1+0+1+1+4+4 = 20 。不偏分散はそれぞれ 40÷4=10 40 \div 4 = 10 と 20÷8=2.5 20 \div 8 = 2.5 で、割る相手は標本の大きさではなく自由度です。

F=s12s22=102.5=4.00F = \frac{s_1^2}{s_2^2} = \frac{10}{2.5} = 4.00

自由度は (5−1, 9−1)=(4, 8) (5-1,\ 9-1) = (4,\ 8) です。F分布表から上側 2.5% 点は F0.025(4,8)=5.05 F_{0.025}(4, 8) = 5.05 、下側 2.5% 点は自由度を入れ替えた F0.025(8,4)=8.98 F_{0.025}(8, 4) = 8.98 の逆数 1÷8.98=0.111 1 \div 8.98 = 0.111 として得ます。上側点 5.05 をそのまま逆数にした 1÷5.05=0.198 1 \div 5.05 = 0.198 を下側点にしてはいけません。自由度を入れ替える手順が抜けています。

有意水準5%
検定の向き棄却域F=4.00 F = 4.00 の位置判定
両側F<0.111 F < 0.111 または F>5.05 F > 5.05 0.111<4.00<5.05 0.111 < 4.00 < 5.05 棄却しない
片側(上側)F>3.84 F > 3.84 4.00>3.84 4.00 > 3.84 棄却する

同じ F=4.00 F = 4.00 でも、問いの向きで結論が分かれました。上側の確率は約 0.0453 なので、片側検定なら 0.0453 < 0.05 からも棄却されると分かります。等裾の両側p値は、小さいほうである上側確率を2倍した約 0.0905 で、5%では棄却しないという表の判定と一致します。そして、結果を見てから片側へ切り替えることはできません。片側にするなら「路線Aのほうが遅れのばらつきが大きいのではないか」という向きを、調べる前に定めておく必要があります。

不偏分散が 4 倍も違うのに両側では棄却されないのは、5 日と 9 日という調査日数の少なさによります。自由度が小さいと F F 自体が大きくばらつくため、境界が 5.05 と遠くに置かれるのです。同じ 4 倍の差でも、各 21 日ずつ調べて自由度が (20,20) (20, 20) になれば上側 2.5% 点は 2.46 まで下がり、両側でも棄却されます。

入門問題

全 4 問
問 1

2つの正規母集団について H0:σ12=σ22 H_0: \sigma_1^2 = \sigma_2^2 を検定します。用いる検定統計量とその帰無分布の組み合わせはどれですか。

問 2

標本の大きさが 6 と 13 の2つの正規母集団から標本をとり、F=s12s22 F = \frac{s_1^2}{s_2^2} (分子は大きさ 6 の群)を計算しました。帰無仮説のもとでこの統計量が従うF分布の自由度はどれですか。

問 3

標本の大きさが 7 と 11 の2群で F=s12s22 F = \frac{s_1^2}{s_2^2} (分子は大きさ 7 の群)を計算し、有意水準5%の両側検定を行います。F0.025(6,10)=4.07 F_{0.025}(6, 10) = 4.07 ・F0.025(10,6)=5.46 F_{0.025}(10, 6) = 5.46 のとき、棄却域の下側の境界はいくつですか。

問 4

同じ設定(自由度 (6, 10) (6,\ 10) 、有意水準5%の両側検定、下側の境界 0.183・上側の境界 4.07)で、統計量が F=3.50 F = 3.50 となりました。結論として正しいものはどれですか。