# 確率分布と統計的な推測：二項分布・正規分布・信頼区間・検定 > 期待値と分散の性質、独立な和、二項分布と幾何分布、正規分布と標準化、二項分布の正規近似、標本平均と大数の法則、区間推定と検定の考え方を、有限の場合は証明し極限定理は大学の定理に結びつける。 > https://rikai.mugen-giken.com/mathematics/hs-probability/distributions-and-inference ## 0. この記事の要点 サイコロを何百回も投げると、結果は毎回違っても、全体のばらつき方にはきまりが見えてきます。この記事は、そのきまりを数で表す道具（平均、ばらつき、釣鐘形の分布）を作り、「この結果は偶然としてありえるか」を判断できるようになることを目指します。以下の箇条書きの用語は、本文で一つずつ説明します。 - [確率変数](/mathematics/probability/random-variables#def-random-variable)の**期待値**と**分散**の基本規則は $E[aX+b]=aE[X]+b$、$V[aX+b]=a^2V[X]$。[**独立な**](/mathematics/probability/random-variables#def-independence)確率変数の和では、期待値だけでなく**分散も足し算**になります。 - この規則から、二項分布 $B(n,p)$ の平均 $np$ と分散 $np(1-p)$ が、「1 回の試行の指示関数（成功なら $1$、失敗なら $0$ をとる変数）の和」として暗算で出ます。幾何分布の平均は $1/p$ で、**無記憶性**をもちます。 - **正規分布**の密度の積分が $1$ になること（Gauss 積分）と、標準化 $Z=(X-\mu)/\sigma$ の意味を確かめます。 - 二項分布の**正規近似**（de Moivre–Laplace）は中心極限定理の特別な場合で、標本平均の分布が $n\to\infty$ で正規分布に近づく理由もこれです。 - 標本平均の期待値は $\mu$、分散は $\sigma^2/n$。Chebyshev の不等式から**大数の法則**（標本が大きいほど平均は母平均に近い）が出ます。 - **区間推定**（信頼区間 $\bar x\pm1.96\,\sigma/\sqrt n$）と**検定**（棄却域・有意水準）の意味と、**やってはいけない解釈**を確認します。 ## 1. 動機 「サイコロを 600 回投げて 6 の目が 80 回出た。このサイコロは歪んでいるか」と聞かれたとき、答えるには 3 つの量が要ります。平均的な回数（$600\cdot\frac16=100$）、ばらつき（標準偏差 $\sqrt{600\cdot\frac16\cdot\frac56}\approx9.1$）、そして「そのずれは偶然としてよくある範囲か」を判断する基準（正規分布の裾の確率）です。この記事はその順に進みます。 連続型の確率変数・正規分布・中心極限定理の厳密な理論は大学の記事（[確率変数](/mathematics/probability/random-variables)、[極限定理](/mathematics/probability/limit-theorems)、[仮説検定と区間推定](/mathematics/mathematical-statistics/hypothesis-testing)）にあります。この記事は、それらを高校の範囲で使うときに、**何が証明されているか**を示します。 ## 2. 期待値・分散の規則 標本空間が有限または可算無限の離散確率空間 $(\Omega,p)$ で考えます（前の記事 [確率](/mathematics/hs-probability/finite-probability) の記号を引き継ぐ）。確率変数 $X$ の期待値は $E[X]=\sum_\omega X(\omega)p(\omega)=\sum_xxP(X=x)$（値ごとの和。、可算無限のときは絶対収束するものだけを考える）です。前の記事の期待値の線形性などの証明は、和の線形性と和の順序の入れ替えだけを使うので、絶対収束する級数については可算無限の場合にもそのまま成り立ちます。 $\mu=E[X]$ とする。$X$ の**分散**を $V[X]=E\bigl[(X-\mu)^2\bigr]$、**標準偏差**を $\sigma(X)=\sqrt{V[X]}$ と定める。 これは、データの分散（）の「各値の確率が $1/n$」の場合で、 の特別な場合でもあります。 $V[X]=E[X^2]-(E[X])^2$ である。 $(X-\mu)^2=X^2-2\mu X+\mu^2$ に期待値の線形性（）を使うと $V[X]=E[X^2]-2\mu E[X]+\mu^2=E[X^2]-\mu^2$。 $a,b$ を実数とすると $E[aX+b]=aE[X]+b$、$V[aX+b]=a^2V[X]$ である。 前半は線形性（、定数の期待値は定数）。後半は $aX+b-E[aX+b]=a(X-\mu)$ より $V[aX+b]=E[a^2(X-\mu)^2]=a^2V[X]$。 確率変数 $X,Y$ が**独立**であるとは、すべての値 $x,y$ について $P(X=x,\ Y=y)=P(X=x)\,P(Y=y)$ が成り立つこと。確率変数 $X_1,\dots,X_n$ が**独立**であるとは、すべての値 $x_1,\dots,x_n$ について $P(X_1=x_1,\dots,X_n=x_n)=\prod_iP(X_i=x_i)$ が成り立つこと。 これは の離散の場合です。 $X,Y$ が独立で、$f,g$ が実数値関数ならば、$f(X)$ と $g(Y)$ も独立である。 $P(f(X)=a,\ g(Y)=b)=\sum_{x:f(x)=a}\sum_{y:g(y)=b}P(X=x,Y=y)=\sum_{x:f(x)=a}P(X=x)\sum_{y:g(y)=b}P(Y=y)=P(f(X)=a)P(g(Y)=b)$。 $X,Y$ が独立で、期待値をもつならば $E[XY]=E[X]\,E[Y]$ である。 $E[XY]=\sum_{x,y}xy\,P(X=x,Y=y)=\sum_{x,y}xy\,P(X=x)P(Y=y)=\bigl(\sum_xxP(X=x)\bigr)\bigl(\sum_yyP(Y=y)\bigr)$（絶対収束する級数の積）。 $X_1,\dots,X_n$ が（どの 2 つも）独立で分散をもつならば、 $$ E\Bigl[\sum_iX_i\Bigr]=\sum_iE[X_i],\qquad V\Bigl[\sum_iX_i\Bigr]=\sum_iV[X_i]. $$ 期待値の式は独立性がなくても成り立つ。 期待値の式は線形性（）。$\mu_i=E[X_i]$、$S=\sum X_i$ とすると $S-E[S]=\sum_i(X_i-\mu_i)$ で $$ V[S]=E\Bigl[\Bigl(\sum_i(X_i-\mu_i)\Bigr)^2\Bigr]=\sum_iE\bigl[(X_i-\mu_i)^2\bigr]+\sum_{i\ne j}E\bigl[(X_i-\mu_i)(X_j-\mu_j)\bigr]. $$ $i\ne j$ では、 により $X_i-\mu_i$ と $X_j-\mu_j$ は独立で、 から $E[(X_i-\mu_i)(X_j-\mu_j)]=E[X_i-\mu_i]\,E[X_j-\mu_j]=0$。よって $V[S]=\sum_iV[X_i]$。 分散の足し算は独立性（正確には共分散が $0$）が要ります。たとえば $Y=X$ なら $V[X+Y]=V[2X]=4V[X]\ne2V[X]$ です。 ## 3. 二項分布と幾何分布 $n$ を正の整数、$0\le p\le1$、$q=1-p$ とする。確率変数 $X$ が $P(X=k)=\binom nkp^kq^{n-k}$（$k=0,\dots,n$）を満たすとき、$X$ は**二項分布** $B(n,p)$ に従うという。 $X\sim B(n,p)$ ならば $E[X]=np$、$V[X]=np(1-p)$ である。 $E[X]$ と $V[X]$ は $X$ の分布だけで決まる（値ごとの和）ので、次の模型で計算してよい。成功確率 $p$ の試行を独立に $n$ 回繰り返す直積確率空間（）で、$X_i$ を第 $i$ 回目が成功なら $1$、失敗なら $0$ とすると、$X=\sum_{i=1}^nX_i$ はちょうど $k$ 回成功する確率が $\binom nkp^kq^{n-k}$（）なので $B(n,p)$ に従う。直積空間の重みの定義から $P(X_1=x_1,\dots,X_n=x_n)=\prod_ip_0(x_i)=\prod_iP(X_i=x_i)$ なので $X_1,\dots,X_n$ は独立。$E[X_i]=p$、$E[X_i^2]=p$（$X_i^2=X_i$）だから より $V[X_i]=p-p^2=pq$。 より $E[X]=np$、$V[X]=npq$。 6 の目が出る回数 $X\sim B(600,\tfrac16)$ は、平均 $100$、分散 $600\cdot\frac16\cdot\frac56=\frac{250}3$、標準偏差 $\approx9.1$。平均から約 $2$ 標準偏差（$\approx18$）離れた $80$ は、「よくある範囲の端」にあたります（§7 の信頼区間・§8 の検定で数値化します）。 $0 $X$ が幾何分布に従うとき（$0m)=q^m$（$m\ge0$）。 3. $E[X]=1/p$。 4. （無記憶性）$m,n\ge0$ について $P(X>m+n\mid X>m)=P(X>n)$。 1. 等比級数（）より $\sum_kq^{k-1}p=p\cdot\frac1{1-q}=1$。 2. $P(X>m)=\sum_{k>m}q^{k-1}p=p\,q^m\cdot\frac1{1-q}=q^m$（同じ等比級数）。 3. $S=\sum_{k\ge1}kq^{k-1}$ は比判定法（、$\frac{(k+1)q^k}{kq^{k-1}}\to q<1$）で絶対収束する。収束する級数の差は項ごとの差なので、$(1-q)S=\sum_{k\ge1}kq^{k-1}-\sum_{k\ge1}kq^k=\sum_{k\ge1}\bigl(kq^{k-1}-(k-1)q^{k-1}\bigr)=\sum_{k\ge1}q^{k-1}=\frac1{1-q}$（右辺の第 2 項を $k\to k-1$ とずらした）。ゆえに $S=\frac1{(1-q)^2}=\frac1{p^2}$、$E[X]=pS=\frac1p$。 4. $\{X>m+n\}\subseteq\{X>m\}$ なので $P(X>m+n\mid X>m)=\dfrac{P(X>m+n)}{P(X>m)}=\dfrac{q^{m+n}}{q^m}=q^n=P(X>n)$。 無記憶性は、「$m$ 回失敗した事実は、次の成功までの待ち時間の分布に影響しない」ということです。サイコロで 6 が出るまでの回数の期待値は $6$ で、ずっと 6 が出ていなくても次の 1 回で 6 が出る確率は $\frac16$ のままです（「そろそろ出る」という錯覚の訂正）。 ## 4. 正規分布 連続型の確率変数 $X$ が**確率密度関数** $f\ge0$ をもつとは、すべての $a\le b$ で $P(a\le X\le b)=\int_a^bf(x)\,dx$ となることです（ の特別な場合）。このとき $E[g(X)]=\int g(x)f(x)\,dx$（[広義積分](/mathematics/calculus/improper-parametric-integrals#def-improper)が絶対収束するとき）と定めます。 実数 $\mu$ と $\sigma>0$ に対し、密度 $$ f(x)=\frac1{\sqrt{2\pi}\,\sigma}\exp\Bigl(-\frac{(x-\mu)^2}{2\sigma^2}\Bigr) $$ をもつ分布を**正規分布** $N(\mu,\sigma^2)$ という。$N(0,1)$ を**標準正規分布**といい、その密度を $\varphi(z)=\dfrac1{\sqrt{2\pi}}e^{-z^2/2}$、$\Phi(x)=\int_{-\infty}^x\varphi(z)\,dz$ と書く。 $\displaystyle\int_{-\infty}^\infty f(x)\,dx=1$、とくに $\Phi(\infty)=\int_{-\infty}^\infty\varphi(z)\,dz=1$ である。 $x=\mu+\sqrt2\,\sigma u$ と置換する（、有限区間 $[-R,R]$ で行って $R\to\infty$）。$dx=\sqrt2\sigma\,du$ で、 $$ \int f(x)\,dx=\frac{\sqrt2\,\sigma}{\sqrt{2\pi}\,\sigma}\int_{-\infty}^\infty e^{-u^2}du=\frac1{\sqrt\pi}\cdot\sqrt\pi=1 $$ （ より $\int_{-\infty}^\infty e^{-u^2}du=\sqrt\pi$）。 $X\sim N(\mu,\sigma^2)$ とすると、$Z=\dfrac{X-\mu}\sigma$ は標準正規分布に従う。すなわち $P(a\le X\le b)=\Phi\Bigl(\dfrac{b-\mu}\sigma\Bigr)-\Phi\Bigl(\dfrac{a-\mu}\sigma\Bigr)$。 $P(a\le X\le b)=\int_a^bf(x)\,dx$ に $x=\mu+\sigma z$ を代入すると（）、$f(\mu+\sigma z)\,\sigma\,dz=\varphi(z)\,dz$ なので $\int_{(a-\mu)/\sigma}^{(b-\mu)/\sigma}\varphi(z)\,dz$。$P(c\le Z\le d)$ はこれと同じ式（$a=\mu+\sigma c$、$b=\mu+\sigma d$）なので $Z$ の密度は $\varphi$。 $X\sim N(\mu,\sigma^2)$ ならば $E[X]=\mu$、$V[X]=\sigma^2$ である。 により $E[g(X)]=\int g(\mu+\sigma z)\varphi(z)\,dz$ と書ける。まず標準正規分布で調べる。$\varphi'(z)=-z\varphi(z)$ である。$\int_0^Rz\varphi(z)\,dz=\dfrac{1-e^{-R^2/2}}{\sqrt{2\pi}}$ は $R\to\infty$ で収束し、$z\varphi(z)$ は奇関数なので $\int_{-\infty}^\infty z\varphi(z)\,dz=0$。次に、部分積分（）で $$ \int_{-R}^Rz^2\varphi(z)\,dz=\bigl[-z\varphi(z)\bigr]_{-R}^{R}+\int_{-R}^R\varphi(z)\,dz=-2R\varphi(R)+\int_{-R}^R\varphi. $$ $e^{R^2/2}\ge R^2/2$ より $R\varphi(R)\le\dfrac2{\sqrt{2\pi}\,R}\to0$、$\int_{-R}^R\varphi\to1$（）なので $\int z^2\varphi=1$。よって標準正規分布は平均 $0$、分散 $1$。$X=\mu+\sigma Z$ なので $E[X]=\mu+\sigma\cdot0=\mu$、$V[X]=E[(\sigma Z)^2]=\sigma^2\cdot1$。 なお、上の証明では、無限区間の積分を区間 $[-R,R]$ での計算と $R\to\infty$ の極限で扱っています。不等式 $e^{R^2/2}\ge R^2/2$ は、$R\varphi(R)$ が $0$ に近づくことを示すために、分母を大きく見積もる道具として使っています。 $\Phi$ は初等関数では書けないので、値は数値積分の表で与えられます。よく使う値は | 区間 | $P(\lvert Z\rvert\le c)$ | |---|---| | $c=1$ | $0.6827$ | | $c=1.96$ | $0.9500$ | | $c=2$ | $0.9545$ | | $c=3$ | $0.9973$ | です。「$\mu\pm\sigma$ に約 $68\%$、$\mu\pm2\sigma$ に約 $95\%$、$\mu\pm3\sigma$ にほぼ全部」が経験則になります。 得点が $N(60,15^2)$ に従うとすると、75 点以上の割合は $P(X\ge75)=P(Z\ge1)=\frac{1-0.6827}{2}\approx0.159$ です。 ## 5. 二項分布の正規近似 中心極限定理（）は、「独立で同分布な確率変数の和を標準化したものの分布関数は、各点で $\Phi$ に収束する」ことを主張します。ここでの「分布関数」は $P(Z_n\le x)$ を $x$ の関数とみたもので、これを区間の確率の形にしておきます。次の補題の証明は、$P(Z_n $X_1,X_2,\dots$ を独立同分布、$\mu=E[X_1]$、$\sigma^2=V[X_1]\in(0,\infty)$ とし、$Z_n=\dfrac{S_n-n\mu}{\sigma\sqrt n}$（$S_n=\sum_{i\le n}X_i$）とおく。任意の $a $F_n(x)=P(Z_n\le x)$ とすると、 により $F_n(x)\to\Phi(x)$（すべての $x$）。$P(a\le Z_n\le b)=F_n(b)-P(Z_n0$ について $F_n(a-\delta)\le P(Z_n $0 成功確率 $p$ の独立な試行の無限列 $\xi_1,\xi_2,\dots$（$P(\xi_i=1)=p$、$P(\xi_i=0)=q$）が存在する（）。$S_n=\sum_{i\le n}\xi_i$ は、 の証明と同じく $B(n,p)$ に従い、$X_n$ と同じ分布をもつ。確率 $P\bigl(a\le(X_n-np)/\sqrt{npq}\le b\bigr)$ は分布だけで決まるので $S_n$ で計算してよい。$E[\xi_i]=p$、$V[\xi_i]=pq>0$ なので が適用でき、$Z_n=(S_n-np)/\sqrt{npq}$ について $P(a\le Z_n\le b)\to\Phi(b)-\Phi(a)$。 $X\sim B(100,\tfrac12)$ は平均 $50$、標準偏差 $5$。$P(45\le X\le55)$ の正確な値は約 $0.7287$ です。そのまま標準化すると $P(-1\le Z\le1)=0.6827$ と少し小さくなります。整数の値を連続な区間に置き換えるときに、両端に $0.5$ ずつ余裕をとる**連続修正**（$44.5\le X\le55.5$、$z=\pm1.1$）を使うと $0.7287$ となって、ほぼ一致します。連続修正は定理ではなく、近似をよくするための経験的な工夫です。 $p=0.3$ の二項分布を標準化し、$n$ を 4 から 400 まで増やしていくと、棒の面積が赤い曲線 $\varphi$ の下の面積に重なっていく様子が見えます。
二項分布 B(n, 0.3) を標準化した棒グラフが、n を 4 から 400 に増やすにつれて標準正規分布の釣鐘型の曲線にぴったり重なっていく動画
## 6. 標本平均と大数の法則 母集団から $n$ 個を**無作為に取り出す**（復元抽出）とは、独立で同じ分布に従う確率変数 $X_1,\dots,X_n$ を考えることです。母集団の平均 $\mu=E[X_1]$、分散 $\sigma^2=V[X_1]$ を**母平均**・**母分散**、$\bar X=\frac1n\sum X_i$ を**標本平均**といいます。 $X_1,\dots,X_n$ が独立で同じ分布をもち、平均 $\mu$、分散 $\sigma^2$ をもつならば、 $$ E[\bar X]=\mu,\qquad V[\bar X]=\frac{\sigma^2}n,\qquad\sigma(\bar X)=\frac\sigma{\sqrt n}. $$ により $E[\sum X_i]=n\mu$、$V[\sum X_i]=n\sigma^2$。 で $a=\frac1n$、$b=0$ とすると $E[\bar X]=\mu$、$V[\bar X]=\frac1{n^2}\cdot n\sigma^2=\frac{\sigma^2}n$。 標本が大きくなると標本平均のばらつきが $1/\sqrt n$ で小さくなる、ということです。「ばらつきを半分にするにはデータが 4 倍要る」の根拠です。 確率変数 $X$ が平均 $\mu$、分散 $V$ をもつとき、任意の $a>0$ について $P(|X-\mu|\ge a)\le\dfrac V{a^2}$ である。 $V=\sum_x(x-\mu)^2P(X=x)\ge\sum_{x:|x-\mu|\ge a}(x-\mu)^2P(X=x)\ge a^2\sum_{x:|x-\mu|\ge a}P(X=x)=a^2P(|X-\mu|\ge a)$。 一般の確率変数の場合は です。 $X_1,X_2,\dots$ が独立で同じ分布をもち、平均 $\mu$、分散 $\sigma^2<\infty$ をもつならば、任意の $\varepsilon>0$ と $n$ について $$ P\bigl(|\bar X_n-\mu|\ge\varepsilon\bigr)\le\frac{\sigma^2}{n\varepsilon^2}\xrightarrow[n\to\infty]{}0. $$ より $\bar X_n$ の平均は $\mu$、分散は $\sigma^2/n$。 に $a=\varepsilon$ を使う。 分散が有限でなくても、平均さえあれば同じ結論が成り立ちます（）。大数の法則は「標本平均は母平均に近づく」ですが、「近づく速さ」は の $1/n$ より精密に、中心極限定理（ の一般形）が与えます。 1 個のさいころの目の平均は $\frac72=3.5$、分散は $\frac{35}{12}$ です。標本平均 $\bar X_n$ の標準偏差は $\sqrt{35/12}/\sqrt n\approx1.71/\sqrt n$ で、$n=10^4$ なら約 $0.017$、$n=10^6$ なら約 $0.0017$ です。次の Python コードで確かめられます。 ```python import random random.seed(2026) def dice_mean(n): return sum(random.randint(1, 6) for _ in range(n)) / n for n in (10, 100, 10_000, 1_000_000): print(n, dice_mean(n)) ``` Python 3.12 で実行すると、`3.8`、`3.91`、`3.4662`、`3.496762` が得られました（乱数の列は実装に依存します）。$n$ が大きくなるほど $3.5$ に近づき、$n=10^4$ で $3.5$ との差が $0.034$（標準偏差の約 $2$ 倍）、$n=10^6$ で $0.003$ 程度です。 ## 7. 区間推定 母集団の平均 $\mu$ は未知で、標本から推定したい。母分散 $\sigma^2$ は既知とします。標本平均 $\bar X_n$ のまわりに $\pm z\,\sigma/\sqrt n$ の幅の区間を作る、というのが高校の「信頼区間」です。 $X_1,X_2,\dots$ を独立同分布、平均 $\mu$（未知）、分散 $\sigma^2\in(0,\infty)$（既知）とし、 $$ I_n=\Bigl[\bar X_n-1.96\frac\sigma{\sqrt n},\ \bar X_n+1.96\frac\sigma{\sqrt n}\Bigr] $$ とする。このとき $n\to\infty$ で $$ P(\mu\in I_n)\longrightarrow\int_{-1.96}^{1.96}\varphi(z)\,dz. $$ この極限の値は、数値積分で $0.9500$（小数第 4 位まで。$0.950004\ldots$）と得られる。 $\mu\in I_n$ は $|\bar X_n-\mu|\le1.96\,\sigma/\sqrt n$ と同値で、$Z_n=\dfrac{S_n-n\mu}{\sigma\sqrt n}=\dfrac{\sqrt n(\bar X_n-\mu)}\sigma$ を使うと $-1.96\le Z_n\le1.96$ と同値。 から $P(-1.96\le Z_n\le1.96)\to\Phi(1.96)-\Phi(-1.96)=\int_{-1.96}^{1.96}\varphi$。定理の最後の数値 $0.9500$ は、この積分の数値積分による値であって、上の論証の一部ではない（$\Phi$ は初等関数で表せないため、値は表や計算機で得る）。 $400$ 人に聞いて $55\%$ が賛成のとき、母比率 $p$ の $95\%$ 信頼区間は、$\sigma=\sqrt{p(1-p)}$ を標本比率 $\hat p=0.55$ で置き換えて $\hat p\pm1.96\sqrt{\hat p(1-\hat p)/400}=0.55\pm0.049$、つまり約 $[0.501,\,0.599]$。区間が $0.5$ をぎりぎり下回らないので、「過半数が賛成」と結論するのは慎重にすべきです。$\sigma$ を推定値で置き換えても被覆確率が $0.95$ に近づくことは、 の考え方（Wald 区間）で保証されます。 $p(1-p)\le\frac14$（$\frac14-p(1-p)=(p-\frac12)^2\ge0$）なので、$p$ が未知でも区間の半幅は $1.96\sqrt{p(1-p)/n}\le\dfrac{1.96\cdot0.5}{\sqrt n}$。半幅を $0.02$ 以下にしたいなら $\sqrt n\ge\dfrac{0.98}{0.02}=49$、$n\ge2401$。誤差を半分にするには $4$ 倍の標本が要ります（ の $\sigma/\sqrt n$）。 信頼区間の定義の一般形（すべての $\theta$ で $P_\theta(\theta\in C(X))\ge1-\alpha$）は にあります。上の定理は被覆確率が極限で $0.95$ に近づくという、漸近的な主張です。 ## 8. 仮説検定の考え方 検定は、「データがある仮説と相容れないほど偶然としては珍しいか」を決める手順です。 **帰無仮説** $H_0$（否定したい仮説）と**対立仮説** $H_1$ を立てる。データから**検定統計量** $T$ を計算し、**棄却域** $R$（データを見る前に決める $T$ の値の集合）に $T$ が入れば $H_0$ を棄却する。**有意水準** $\alpha$（たとえば $0.05$）は、$H_0$ が正しいときに誤って棄却してしまう確率（第 1 種の誤りの確率）の上限で、$P_{H_0}(T\in R)\le\alpha$ となるように $R$ を選ぶ。 一般の定義は にあります。 $00$ について、棄却域 $\{|Z_n|\ge c\}$ の $H_0$ のもとでの確率は $$ P_{H_0}\bigl(|Z_n|\ge c\bigr)\longrightarrow2\bigl(1-\Phi(c)\bigr)\qquad(n\to\infty) $$ である。とくに $c=1.96$ では極限は約 $0.0500$ で、この棄却域は（漸近的に）有意水準 $5\%$ の検定を与える。 と同じく、$X_n$ を二項分布に従う $S_n$ とみなして を使う。任意の $\delta\in(0,c)$ について $$ P(-c+\delta\le Z_n\le c-\delta)\le P(|Z_n| サイコロを $600$ 回投げて 6 の目が $80$ 回出た。$H_0:p=\frac16$ のもとで $E=100$、$\sigma=\sqrt{600\cdot\frac16\cdot\frac56}\approx9.13$。$Z=\dfrac{80-100}{9.13}\approx-2.19$ で $|Z|\ge1.96$ なので、有意水準 $5\%$ で $H_0$ を棄却します（両側の $p$ 値は $2\Phi(-2.19)\approx0.029$）。「6 の目が出にくい」と結論するのではなく、「$p=\frac16$ と矛盾する程度のずれがある」と言うところまでです。 コインを $100$ 回投げて表が $60$ 回。$H_0:p=\frac12$ で $Z=\dfrac{60-50}5=2.0\ge1.96$ なので、正規近似では $5\%$ 水準で棄却されます（$p$ 値 $0.0455$）。ところが連続修正をすると $Z=1.9$（$p$ 値 $0.057$）で棄却されず、二項分布の正確な両側確率も約 $0.0569$ で棄却されません。境界付近の結論は、近似の仕方で変わります。「$5\%$ を超えたかどうか」だけを見ず、$p$ 値の大きさと標本の大きさを一緒に報告するのが誠実です。 次のコードは、「$95\%$ 信頼区間の手順を何度も繰り返すと、約 $95\%$ が真の $p$ を含む」ことを乱数で確かめます（コインの $p=0.5$、$n=100$、標本比率で $\sigma$ を置き換えた Wald 区間）。 ```python import random random.seed(1) trials, n, hit = 20000, 100, 0 for _ in range(trials): x = sum(random.random() < 0.5 for _ in range(n)) p_hat = x / n half = 1.96 * (p_hat * (1 - p_hat) / n) ** 0.5 hit += (p_hat - half <= 0.5 <= p_hat + half) print(hit / trials) ``` Python 3.12 では `0.94145` が出ました。$0.95$ より少し小さいのは、$\sigma$ を推定値で置き換えた Wald 区間の、$n=100$ での近似の誤差です。 ## 9. 大学ではどこへ続くか - 確率変数・期待値・分散・独立の測度論的な定式化と、連続型の確率変数は [確率変数](/mathematics/probability/random-variables)。大数の法則・中心極限定理の証明は [極限定理](/mathematics/probability/limit-theorems)。 - 点推定の良し悪し（不偏性・一致性・有効性）、検定の最適性（Neyman–Pearson の補題）、信頼区間の一般論は [仮説検定と区間推定](/mathematics/mathematical-statistics/hypothesis-testing) で学びます。 - 母分散が未知のときの正確な区間推定（$t$ 分布）や、カイ二乗検定は数理統計学の標準的な話題です。 ## 10. 演習 $X\sim B(20,0.3)$ の平均と分散を求めよ。また $Y=3X-5$ の平均と分散を求めよ。 より $E[X]=6$、$V[X]=20\cdot0.3\cdot0.7=4.2$。 より $E[Y]=3\cdot6-5=13$、$V[Y]=9\cdot4.2=37.8$。 サイコロを繰り返し投げて、初めて 6 の目が出るまでの回数 $X$ について、期待値 $E[X]$ と $P(X>12)$ を求めよ。また、12 回続けて 6 が出なかったとき、あと何回で 6 が出るかの期待値はいくらか。 $X$ は $p=\frac16$ の幾何分布で、 より $E[X]=6$、$P(X>12)=(5/6)^{12}\approx0.112$。無記憶性から、12 回続けて出なくても、あとの回数の分布は最初と同じなので期待値は $6$ 回。 身長が $N(170,6^2)$ に従うとき、$182$ cm 以上の割合を、$P(|Z|\le2)=0.9545$ を用いて求めよ。 より $Z=\dfrac{182-170}6=2$。$P(Z\ge2)=\dfrac{1-0.9545}2\approx0.0228$（標準正規分布の密度は $0$ について対称）。約 $2.3\%$。 標準偏差が $\sigma=8$ とわかっている母集団から $n=100$ 個を取り出して標本平均 $52$ を得た。母平均の $95\%$ 信頼区間を求めよ。母平均がこの区間に入る確率が $95\%$ だと言えるか。 半幅は $1.96\cdot\frac8{\sqrt{100}}=1.568$ で、区間は約 $[50.43,\ 53.57]$（）。母平均は一定の値なので、この区間に入る確率が $95\%$ だとは言えない（§7 の注意）。同じ手順で作る区間のうち、約 $95\%$ が母平均を含む、というのが正しい解釈。 コインを $400$ 回投げて表が $220$ 回出た。$H_0:p=\frac12$ を有意水準 $5\%$ で検定せよ（ の棄却域 $|Z|\ge1.96$ を使う）。また、同じ割合（$55\%$）で標本を $100$ 回に減らすと、結論はどうなるか。 $n=400$: $E=200$、$\sigma=\sqrt{400\cdot\frac14}=10$、$Z=\dfrac{220-200}{10}=2.0\ge1.96$ なので棄却。$n=100$ で表が $55$ 回なら $E=50$、$\sigma=5$、$Z=1.0<1.96$ で棄却されない。割合は同じでも、標本が小さいと偶然のばらつきの範囲に収まる（検出力が足りない）。 ## 参考文献 1. 東京大学教養学部統計学教室（編）『統計学入門』東京大学出版会、1991 — 確率分布・推定・検定の各章。 2. William Feller, *An Introduction to Probability Theory and Its Applications, Volume I*, 3rd ed., Wiley, 1968 — 第 7 章（二項分布の正規近似）。 3. Rick Durrett, *Probability: Theory and Examples*, 5th ed., Cambridge University Press, 2019 — 中心極限定理の章。 4. George Casella, Roger L. Berger, *Statistical Inference*, 2nd ed., Cengage, 2002 — 区間推定・仮説検定の章。