ひよこ統計

連続修正

統計検定2級

二項分布の正規近似 では、近似して求めた 0.1469 が厳密な 0.1583 と1ポイントほど食い違いました。このずれは偶然ではなく、決まった向きに出ます。とびとびの値をとる分布を、切れ目なくつながった曲線で置き換えたことから来る系統的なずれだからです。

理由は棒グラフを見ると分かります。棒が数えられる大きさのほうが見やすいので、いったん B(20,0.5) B(20, 0.5) に移って「表が 8 回以下」を考えます。X=8 X = 8 の棒は幅 1 をもち、7.5 から 8.5 までの範囲を占めています。曲線で 8 までの面積をとると、この棒の右半分が切り落とされてしまいます。そこで境界を棒の端まで動かす——この手当てが連続修正です。

二項分布の棒グラフに正規分布の曲線を重ね、8以下の棒を濃い色で塗った図。8.0 の位置に破線、8.5 の位置に実線が引かれ、8の棒が7.5から8.5までの幅をもつことが分かる
8 の棒は 7.5 から 8.5 まで。曲線で拾うなら境界は 8.5 に置く

破線の 8.0 で切ると、境界のすぐ内側にある 8 の棒の右半分が抜け落ちます。実線の 8.5 まで伸ばせば棒1本ぶんがまるごと入り、棒の合計と曲線の下の面積がよく合うようになります。どちらへ 0.5 動かすかは、その棒を含めたいかどうかで決まります。

求めたい確率修正後に使う境界
P(X≤a) P(X \le a) a+0.5 a + 0.5 まで
P(X≥a) P(X \ge a) a−0.5 a - 0.5 から
P(X=a) P(X = a) a−0.5 a - 0.5 から a+0.5 a + 0.5 まで

「とりあえず 0.5 を足す」と覚えると必ず間違えます。a a の棒を数えたいなら境界を外側へ広げ、数えたくないなら内側へ狭める、というだけのことです。P(X>a) P(X > a) は a a を含まないので P(X≥a+1) P(X \ge a+1) と同じであり、修正後の境界は a+0.5 a + 0.5 になります。符号を覚えるのではなく、棒を含むかどうかで毎回決めます。

修正によって確率が変わる幅は、n n が大きいほど小さくなります。棒1本ぶんの幅は 1 のまま変わらないのに、分布の広がり np(1−p) \sqrt{np(1-p)} は n n とともに大きくなるので、n n が大きいほど 0.5 のずらしは相対的に小さくなります。n n が十分大きい場面で連続修正が省かれることがあるのは、このためです。

n が小さいほど修正による差が大きい

コインを20回投げて、表が 8 回以下になる確率を求めます。X∼B(20,0.5) X \sim B(20, 0.5) なので平均は 20×0.5=10 20 \times 0.5 = 10 、分散は 20×0.5×0.5=5 20 \times 0.5 \times 0.5 = 5 、標準偏差は 5=2.236 \sqrt{5} = 2.236 です。

まず修正なしで計算します。Z=8−102.236=−0.89 Z = \frac{8 - 10}{2.236} = -0.89 から P(Z≤−0.89)=0.1867 P(Z \le -0.89) = 0.1867 。いっぽう二項分布から直接求めた厳密な値は 0.2517 で、0.065 も小さく出ています。本来の値の4分の1ほどを取りこぼした計算です。

境界を 8.5 に付け替えると Z=8.5−102.236=−0.67 Z = \frac{8.5 - 10}{2.236} = -0.67 、P(Z≤−0.67)=0.2514 P(Z \le -0.67) = 0.2514 。厳密な 0.2517 とほぼ重なりました。0.5 をずらしただけで、誤差が 0.065 から 0.0003 へ縮んでいます。

n=20 n = 20 では標準偏差が 2.236 しかないので、0.5 のずれは標準偏差の2割以上にあたります。小さな n n で修正を省けない理由がここにあります。

近似で残っていたずれを詰める

冒頭に挙げた 0.1469 と 0.1583 の食い違いにも、同じ手当てを当てはめます。10回に1回カウントされる歩数計で1,000歩あるいたときのカウント数 X∼B(1000,0.1) X \sim B(1000, 0.1) について 110 回以上を求めるので、110 の棒を含めるために境界は 109.5 です。Z=109.5−10090=1.00 Z = \frac{109.5 - 100}{\sqrt{90}} = 1.00 から P(Z≥1.00)=0.1587 P(Z \ge 1.00) = 0.1587 。

厳密な値は 0.1583 でしたから、修正なしの 0.1469 に比べてずっと近づきました。n=1000 n = 1000 では標準偏差が 9.49 もあるため、0.5 のずらしが動かす Z Z はわずか 0.05 ほど。それでも確率にすると1ポイントぶんの違いになります。

入門問題

全 4 問
問 1

X∼B(50,0.4) X \sim B(50, 0.4) について P(X≥25) P(X \ge 25) を正規分布で近似します。連続修正を行うとき、正規分布側で使う境界はどれですか。

問 2

二項分布に従う X X について P(X=12) P(X = 12) を正規分布で近似するとき、連続修正後に求める範囲はどれですか。

問 3

X∼B(30,0.5) X \sim B(30, 0.5) について P(X≤12) P(X \le 12) を、連続修正をして正規近似します。標準化した値に最も近いものはどれですか。

問 4

連続修正が必要になる理由として最も適切なものはどれですか。