MathLabs

確率と統計

推定

標本データを使って母集団の未知の値を、点推定や区間推定によって推し量る手法。

直観一握りから全体へ

大きな鍋のスープをひとさじ味見して、鍋全体の塩加減を判断する場面を想像してほしい。鍋のスープを全部飲む必要はなく、代表的なひとさじ(標本)だけで鍋全体(母集団)の塩加減を推測できる。推定とは、そのひとさじを信頼できる推測に変え、その推測がどれくらいずれうるかまで示す数学である。

標本サイズが大きくなるにつれて平均の周りに狭まっていく標本分布を示す釣鐘型の曲線。
標準正規標本分布 N(0,1)\mathcal{N}(0,1) と中央の信頼区間 [−z,+z][-z, +z](緑色領域):pp を動かして信頼水準の幅を変えてみよう。

中高点推定:ただ一つの最良の推測

定義: 点推定量

点推定量とは、標本データから計算される単一の数値であり、母集団の未知のパラメータに対する最良の推測として用いられる。母平均 μ\mu を推定する場合、自然な点推定量は標本平均 Xˉ\bar X、すなわち観測値 X1,…,XnX_1,\dots,X_n の平均である。

Xˉ=1n∑i=1nXi\bar X = \frac{1}{n}\sum_{i=1}^{n} X_i

ここで nn は標本サイズ、X1,…,XnX_1,\dots,X_n は個々の標本観測値、Xˉ\bar X はそれらの平均である。標本が異なれば Xˉ\bar X の値も異なるため、それ自体が独自の分布を持つ確率変数であり、それを平均の標本分布と呼ぶ。この分布の広がりは nn が大きくなるほど狭まっていく。

Xˉ−zα/2σn≤μ≤Xˉ+zα/2σn\bar X - z_{\alpha/2}\frac{\sigma}{\sqrt{n}} \le \mu \le \bar X + z_{\alpha/2}\frac{\sigma}{\sqrt{n}}

これは母標準偏差 σ\sigma が既知のときの μ\mu に対する 100(1−α)%100(1-\alpha)\% 信頼区間である:zα/2z_{\alpha/2} は標準正規分布の各裾に確率 α/2\alpha/2 を残す臨界値であり、σn\frac{\sigma}{\sqrt n} は Xˉ\bar X の標準誤差、すなわち Xˉ\bar X と μ\mu の典型的な隔たりである。

よく使われる信頼水準とその zz 臨界値
信頼水準zα/2z_{\alpha/2}
90%zα/2=1.645z_{\alpha/2}=1.645
95%zα/2=1.960z_{\alpha/2}=1.960
99%zα/2=2.576z_{\alpha/2}=2.576

大学信頼区間が機能する理由

X1,…,XnX_1,\dots,X_n が独立同分布で、すべての ii について E[Xi]=μE[X_i] = \mu であるとき、標本平均 Xˉ=1n∑i=1nXi\bar X = \frac{1}{n}\sum_{i=1}^{n} X_i は E[Xˉ]=μE[\bar X] = \mu を満たす。

なぜ正しいのか?

平均を取ることで系統的な偏りが取り除かれる:多数の仮想的な標本にわたって平均すると、Xˉ\bar X は真の母平均 μ\mu にちょうど一致し、系統的に高すぎたり低すぎたりすることはない。

証明

定義より Xˉ=1n∑i=1nXi\bar X = \frac{1}{n}\sum_{i=1}^{n} X_i である。期待値は線形作用素なので、和と定数係数 1n\frac{1}{n} の上に分配できる:E[Xˉ]=1n∑i=1nE[Xi]E[\bar X] = \frac{1}{n}\sum_{i=1}^{n} E[X_i]。

各 XiX_i は同じ母集団から得られるため、nn 個の項それぞれについて E[Xi]=μE[X_i] = \mu となり、和は 1n∑i=1nμ=1n(nμ)=μ\frac{1}{n}\sum_{i=1}^{n} \mu = \frac{1}{n}(n\mu) = \mu に帰着する。

したがって、任意の標本サイズ nn に対して E[Xˉ]=μE[\bar X] = \mu となる:標本が小さくても大きくても標本平均は不偏である。これは標本分布の中心についての主張にすぎない -- その広がり、すなわち Var⁡(Xˉ)=σ2/n\operatorname{Var}(\bar X) = \sigma^2/n は nn が大きくなるにつれて依然として狭まっていき、それが以下の信頼区間で使われる。

X1,…,XnX_1,\dots,X_n が独立同分布で、平均 μ\mu、有限分散 σ2\sigma^2 を持つとき、nn が大きければ信頼区間 Xˉ−zα/2σn≤μ≤Xˉ+zα/2σn\bar X - z_{\alpha/2}\frac{\sigma}{\sqrt n} \le \mu \le \bar X + z_{\alpha/2}\frac{\sigma}{\sqrt n} は確率およそ 1−α1-\alpha で μ\mu を含む。

なぜ正しいのか?

中心極限定理によれば、nn が大きくなると標準化された標本平均は標準正規変数のように振る舞うため、正規分布の分位点 zα/2z_{\alpha/2} を使って、元の母集団の形に関わらず既知の長期的成功率で μ\mu を挟み込むことができる。

証明

中心極限定理により、Xˉ−μσ/n\frac{\bar X - \mu}{\sigma/\sqrt n} は nn→∞\to\infty のとき分布収束で N(0,1)N(0,1) に近づく。よって nn が大きいとき、−zα/2≤Xˉ−μσ/n≤zα/2-z_{\alpha/2} \le \frac{\bar X - \mu}{\sigma/\sqrt n} \le z_{\alpha/2} が確率およそ 1−α1-\alpha で成り立つ。ここで zα/2z_{\alpha/2} は N(0,1)N(0,1) の各裾に確率 α/2\alpha/2 を残す値である。

不等式の三辺すべてに σn\frac{\sigma}{\sqrt n}>0>0 を掛ける(不等号の向きは変わらない):−zα/2σn≤Xˉ−μ≤zα/2σn-z_{\alpha/2}\frac{\sigma}{\sqrt n} \le \bar X - \mu \le z_{\alpha/2}\frac{\sigma}{\sqrt n}。

Xˉ\bar X をすべての辺から引き、−1-1 を掛ける(不等号の向きが反転する):Xˉ−zα/2σn≤μ≤Xˉ+zα/2σn\bar X - z_{\alpha/2}\frac{\sigma}{\sqrt n} \le \mu \le \bar X + z_{\alpha/2}\frac{\sigma}{\sqrt n}。これはまさに信頼区間の式である:元の確率の主張が確率およそ 1−α1-\alpha で成り立っていたので、並べ替えたこの区間も同じ確率で μ\mu を含む。

誤差の限界 E=zα/2σnE = z_{\alpha/2}\frac{\sigma}{\sqrt n} は 1/n1/n ではなく 1/n1/\sqrt n のようにしか縮まらない:誤差の限界を半分にするには標本サイズ nn を4倍にしなければならず、三分の一にするには観測数を9倍にする必要がある。精度は高くつく:精度の桁を一つ増やすごとに、その前よりもはるかに多くのデータが必要になる。

大学実世界での応用と具体例

信頼区間は、不完全な標本から判断を下さなければならないあらゆる場面に現れる:工場ラインの品質管理、選挙前の世論調査、医学における投与量の試験、繰り返し測定による望遠鏡の較正などである。どの場合も、同じ公式が標本平均を、自身の不確実性に正直な範囲へと変える。

例: ボルトの平均直径の信頼区間

ある品質管理エンジニアが n=64n=64 本のボルトを標本として抽出し、標本平均直径 xˉ=12.02\bar x = 12.02 mm を得た。母標準偏差は σ=0.16\sigma = 0.16 mm と既知である。z0.025=1.96z_{0.025}=1.96 を用いて、真の平均直径 μ\mu に対する95%信頼区間を構成せよ。

解答

標準誤差は 0.1664=0.02\frac{0.16}{\sqrt{64}} = 0.02 mm である。

誤差の限界は 1.96×0.02=0.03921.96 \times 0.02 = 0.0392 mm である。

95%信頼区間は (12.02−0.0392, 12.02+0.0392)=(11.981, 12.059)(12.02 - 0.0392,\ 12.02 + 0.0392) = (11.981,\ 12.059) mm である:真のボルト平均直径がこの範囲に入ると95%の確信を持って言える。

例: ある政策を支持する有権者の割合の推定

ある世論調査員が無作為に選んだ n=400n=400 人の有権者を調査し、ある政策に賛成する標本比率 p^=0.53\hat p = 0.53 を得た。標本比率を0/1回答の標本平均とみなし、z0.025=1.96z_{0.025}=1.96 を用いて、賛成する真の母集団比率 μ\mu に対する95%信頼区間を構成せよ。

解答

単一の0/1回答の推定標準偏差は p^(1−p^)=0.53×0.47≈0.499\sqrt{\hat p(1-\hat p)} = \sqrt{0.53\times 0.47} \approx 0.499 であるため、標本比率の標準誤差は 0.499400≈0.025\frac{0.499}{\sqrt{400}} \approx 0.025 である。

誤差の限界は 1.96×0.025≈0.0491.96 \times 0.025 \approx 0.049 である。

95%信頼区間は (0.53−0.049, 0.53+0.049)=(0.481, 0.579)(0.53 - 0.049,\ 0.53 + 0.049) = (0.481,\ 0.579) である:この調査では真の支持率が約48.1%から57.9%の間にあると95%の確信を持って言えるが、政策が過半数の支持を得ているとは確実に言えないほど広い -- より明確な結論には、より大きな標本が必要である。

誤差の限界の公式 E=zα/2σnE = z_{\alpha/2}\frac{\sigma}{\sqrt n} を用いると、zα/2=1.96z_{\alpha/2}=1.96、σ=8\sigma=8、n=64n=64 のときの誤差の限界はいくらか。

ある1つの標本から、平均気温に対する95%信頼区間が計算された。「95%の信頼」の正しい解釈は何か。

すべての標本サイズ nn について成り立つ性質 E[Xˉ]=μE[\bar X] = \mu は、推定量 Xˉ\bar X のどの特徴を表すか。

標本サイズ nn が100から400に増え、σ\sigma と zα/2z_{\alpha/2} が変わらないとき、誤差の限界は何倍になるか。

参考文献

  1. NIST/SEMATECH (2013). Confidence Limits for the Mean
  2. Diez, D.; Cetinkaya-Rundel, M.; Barr, C. (2019). OpenIntro Statistics