ひよこ統計

母比率の検定

統計検定2級

袋に「発芽率 80%」と表示された種子を 100 粒まいたところ、芽が出たのは 72 粒でした。表示より 8 ポイント低いこの数字を、たまたま出の悪い袋を引いただけと見るか、表示そのものを疑うか。母比率の検定は、この迷いを「母集団の比率が本当に 80% なら、これだけのずれはどれくらい起こりやすいか」という問いに置き換えます。母比率を特定の値 p0 p_0 とおく 帰無仮説 を立て、ずれの大きさを次の統計量ではかります。

z=p^−p0p0(1−p0)/nz = \frac{\hat{p} - p_0}{\sqrt{p_0(1-p_0)/n}}
p^ \hat{p} :標本比率(該当した個数 ÷ 標本の大きさ)、p0 p_0 :帰無仮説が指定する母比率の値、n n :標本の大きさ

分母は標本比率のばらつきです。母比率の区間推定 では根号の中に標本比率 p^ \hat{p} を入れましたが、検定で入れるのは p0 p_0 のほうです。同じ母比率を扱うのだから分母も同じ、と早合点してはいけません。検定は「帰無仮説が正しいとしたら何が起こるか」を調べるので、ばらつきも仮説が与えた値から計算します。推定のほうには仮説がなく、真の比率の見当が p^ \hat{p} しかない、という違いです。

帰無仮説が正しいとき、z z は近似的に標準正規分布に従います。あとは有意水準を決めて境界と比べるだけで、両側検定なら 5% で ∣z∣>1.96 |z| > 1.96 、1% で ∣z∣>2.576 |z| > 2.576 が棄却域です。p値 で判定しても結論は同じで、帰無分布が左右対称なので両側のp値は片側のちょうど 2 倍になります(両側検定と片側検定)。ただし、この近似がいつでも使えるわけではありません。もとは 0 か 1 しかとらない観測を足し上げた二項分布なので(二項分布の正規近似)、次の目安を満たすだけの標本の大きさが要ります。

二項分布は飛び飛びの値をとるため、境界を半個分ずらす 連続修正 を入れる流儀もあります。統計検定の問題では「連続修正を行わない」と断ってあることが多く、ここでも修正なしの式で通します。

冒頭の種子に戻ります。仮説は H0:p=0.8 H_0: p = 0.8 、H1:p≠0.8 H_1: p \neq 0.8 の両側、有意水準は 5% とします。標本比率は p^=72100=0.72 \hat{p} = \frac{72}{100} = 0.72 で、np0=80 np_0 = 80 ・n(1−p0)=20 n(1-p_0) = 20 はどちらも 5 を超えているので、正規近似の式を使えます。

分母のばらつきは 0.8×0.2100=0.04 \sqrt{\frac{0.8 \times 0.2}{100}} = 0.04 、検定統計量は z=0.72−0.800.04=−2.00 z = \frac{0.72 - 0.80}{0.04} = -2.00 。∣−2.00∣>1.96 |{-2.00}| > 1.96 なので棄却域に入り、表示どおりとは考えにくい、と結論できます。p値でも見ておくと、標準正規分布表から P(Z≤−2.00)=0.0228 P(Z \le -2.00) = 0.0228 なので、両側では 2×0.0228=0.0456 2 \times 0.0228 = 0.0456 。0.05 をわずかに下回る、ぎりぎりの結果です。

根号の中に入れる比率p0=0.80 p_0 = 0.80 (正しい方法)p^=0.72 \hat{p} = 0.72 (誤った方法)
ばらつき0.04000.0449
z z −2.00−1.78
両側のp値0.04560.0750
有意水準 5% の判定棄却する棄却しない

誤った方法の列は、区間推定と同じ感覚で p^ \hat{p} を根号に入れてしまった場合です。0.72 は 0.80 より 0.5 に近く、その分だけ 0.72×0.28=0.2016 0.72 \times 0.28 = 0.2016 は 0.8×0.2=0.16 0.8 \times 0.2 = 0.16 を上回ります。ばらつきは 0.0400 から 0.0449 へ広がり、z z は −1.78 まで縮んで、同じデータなのに棄却できなくなりました。境界に近い問題では、分母の取り違えが結論そのものをひっくり返します。

入門問題

全 4 問
問 1

母比率の検定では、p^−p0 \hat{p} - p_0 を標準誤差で割って検定統計量 z z を求めます。この標準誤差として正しいものはどれですか。

問 2

「表の出る確率は 0.5 である」を帰無仮説として、あるコインを 100 回投げたところ表が 60 回出ました。検定統計量 z z の値はいくつですか。

問 3

H0:p=0.02 H_0: p = 0.02 を大きさ 50 の標本で検定します。正規近似を使ってよいかどうかの判断として最も適切なものはどれですか。

問 4

不良率がこれまでの 5% を上回ったかどうかを H1:p>0.05 H_1: p > 0.05 の片側検定で調べたところ、p値は 0.026 でした。同じデータで H1:p≠0.05 H_1: p \neq 0.05 の両側検定を行った場合のp値はいくつですか。