確率と統計 · 大数の法則

大数の法則

同じ試行を独立に繰り返すと、成功の割合は途中で揺れながらも、長期的にはその確率へ収束します。最初の先行分は薄まるのであって、埋め合わされるわけではありません。

01

この図で見る大数の法則

X1,X2,…X_1, X_2, \dots を、有限の平均 μ=E[X]\mu = E[X] を持つ同じ分布からの独立な標本とします(E∣X∣<∞E|X| < \infty で十分です)。最初の n 個の平均を Xˉn\bar X_n とします。大数の法則は Xˉn→μ\bar X_n \to \mu が成り立つことを述べます。弱法則では確率収束、強法則では確率1での収束です。コインなら XiX_i は表で1なので、Xˉn=p^n=Hn/n\bar X_n = \hat p_n = H_n/n は表の割合となり、μ=p\mu = p です。このページでは100,000回の試行を1つ以上の系列で追い、割合が落ち着く理由を探ります。

1つの偏り、2つの尺度
Hn−np=n (p^n−p),SD(p^n)=p(1−p)n,SD(Hn−np)=np(1−p)H_n - np = n\,(\hat p_n - p), \qquad \mathrm{SD}(\hat p_n) = \sqrt{\tfrac{p(1-p)}{n}}, \qquad \mathrm{SD}(H_n - np) = \sqrt{np(1-p)}

上段は p^n\hat p_n、下段は表の回数と期待回数との差 Hn−npH_n - np を示します。公平なコインでは「表の回数 − 裏の回数」の半分です。同じ偏りが上では縮み、下では広がります。各 n の破線の帯は二項分布の厳密な95%範囲です。HnH_n が下側に来る確率、上側に来る確率がそれぞれ2.5%以下となる最も狭い整数区間 [an,bn][a_n, b_n] なので、少なくとも95%を含みます(数値欄には厳密な確率を表示)。割合では [an/n, bn/n][a_n/n,\ b_n/n]、回数の差では [an−np, bn−np][a_n - np,\ b_n - np] となり、同じ事象です。最初に表10回で先行させたコインと、表5回で始めたコインは、それ以降の試行を共有します。表の回数の差は常に厳密に5で、割合の差は 5/n5/n です。先行分は薄まり、埋め合わされるわけではありません。

02

注目するところ

  • 帯は各時点でのものです。固定した各 n では少なくとも95%の系列が中に入りますが、すべての n で中にとどまる系列はずっと少なくなります。ときどき帯の外に出ることは、この法則に反する現象ではありません。
  • 最初の先行分を埋め合わせる仕組みはありません。表が続いた後も、次に表が出る確率は p のままです(そう考えないのがギャンブラーの誤謬です。「連続後の集計」で確かめられます)。割合が落ち着くのは、同じ回数の偏りを、ますます大きくなる n で割るためです。
  • 帯の幅は 1/n1/\sqrt{n} のように縮みます。同じ95%で精度を小数1桁上げるには、約100倍の試行が必要です。この法則には有限の平均が必要です。コーシー分布では n 個の標本の平均も1個の標本より落ち着きません(中心極限定理のページ参照)。
  • ヤコブ・ベルヌーイはコインのような試行の場合を Ars Conjectandi(推測術)(1713年)で証明しました。「大数の法則」という名前はポアソン(1837年)によります。
03

操作方法

  • 赤いキーは固定シードで事前計算した試行を再生します。実行中に押すと停止、再度押すと続行し、最後では次のシードで試行します。図の任意の場所を押すかドラッグするとカーソルが動き、数値欄はその左側の試行だけで計算します。共有はシード、アルゴリズムの版、カーソル位置を保存します。
04

関連項目

  • H16二項分布準備中
  • H19中心極限定理
  • H18標本分布準備中
  • H20信頼区間準備中
  • H05正規分布準備中

さらに読む: ウィキペディア: Law of large numbers (英語); OpenStax, Introductory Statistics 2e, 3.1 Terminology (英語); MacTutor数学史アーカイブ: Jacob Bernoulli (英語).