Lemma数学ロードマップ

11 確率論 · 第 2 章

大数の法則

目安 9〜13 時間定理など 12演習 7 問

この章の目標

  • 確率収束と概収束の違いを理解し、大数の弱法則・強法則を証明できる
  • 切断・部分列・ボレル–カンテリの補題を組み合わせる強法則の証明技法(エテマディの証明)を身につける
  • コルモゴロフの 0-1 法則を証明し、末尾事象の確率が 0 か 1 であることを使える
  • モンテカルロ法・正規数・ワイエルシュトラスの近似定理・グリヴェンコ–カンテリの定理への応用を説明できる

前提:第1章、06-measure-integration 第4章(測度収束)

2.1 二つの収束

公平なコインを nn 回投げたときの表の回数を SnS_n とすると、経験的に相対頻度 Sn/nS_n/n は 1/21/2 に近づく。確率論の出発点ともいえるこの経験則を定理として述べるには、「確率変数列が収束する」ことの意味を定める必要がある。この章では次の二つを使う(他の収束概念とあわせた比較は第3章で行う)。

定義 2.1(概収束・確率収束)確率変数列 XnX_n と確率変数 XX について

  1. P(lim⁡nXn=X)=1P(\lim_n X_n = X) = 1 のとき、XnX_n は XX に概収束 (almost sure convergence) するといい、Xn→XX_n \to X a.s. と書く。
  2. 任意の ε>0\varepsilon > 0 について P(∣Xn−X∣>ε)→0P(\lvert X_n - X \rvert > \varepsilon) \to 0 のとき、XnX_n は XX に確率収束 (convergence in probability) するといい、Xn→PXX_n \xrightarrow{P} X と書く。

概収束は測度論の a.e. 収束、確率収束は測度収束(06-measure-integration 第4章 定義 4.15)である。全測度が有限なので概収束すれば確率収束するが、逆は成り立たない(第3章)。確率収束は「各時刻 nn で大きくずれる確率が小さい」ことしか言わず、一つの標本路 n↦Xn(ω)n \mapsto X_n(\omega) の挙動については何も言わない。

以下、X1,X2,…X_1, X_2, \dots に対し Sn=X1+⋯+XnS_n = X_1 + \cdots + X_n とおく。

2.2 大数の弱法則

定理 2.2(L2L^2 弱法則)X1,X2,…X_1, X_2, \dots は互いに無相関で、E[Xn]=mE[X_n] = m、Var⁡(Xn)≤C\operatorname{Var}(X_n) \leq C(CC は定数)を満たすとする。このとき Sn/n→mS_n/n \to m は L2L^2 でも確率でも成り立つ。

証明. 無相関なので分散は加法的で、E[(Sn/n−m)2]=Var⁡(Sn)/n2≤C/n→0E[(S_n/n - m)^2] = \operatorname{Var}(S_n)/n^2 \leq C/n \to 0。チェビシェフの不等式より P(∣Sn/n−m∣>ε)≤C/(nε2)→0P(\lvert S_n/n - m \rvert > \varepsilon) \leq C/(n\varepsilon^2) \to 0。□\square

例 2.3 公平なコインでは Var⁡(Xk)=1/4\operatorname{Var}(X_k) = 1/4 なので P(∣Sn/n−1/2∣≥ε)≤1/(4nε2)P(\lvert S_n/n - 1/2 \rvert \geq \varepsilon) \leq 1/(4n\varepsilon^2)。誤差 ε\varepsilon の典型的な大きさは 1/n1/\sqrt{n} の程度である。この n\sqrt{n} のスケールでの揺らぎを精密に述べるのが第4章の中心極限定理である。

分散が存在しなくても、平均さえあれば弱法則は成り立つ。証明の鍵は切断 (truncation) である。

定理 2.4(大数の弱法則, weak law of large numbers)X1,X2,…X_1, X_2, \dots を i.i.d. で E[∣X1∣]<∞E[\lvert X_1 \rvert] < \infty とし、m=E[X1]m = E[X_1] とする。このとき Sn/n→PmS_n/n \xrightarrow{P} m。

証明. nn ごとに Yn,k=Xk1{∣Xk∣≤n}Y_{n,k} = X_k \mathbf{1}_{\lbrace \lvert X_k \rvert \leq n \rbrace}(1≤k≤n1 \leq k \leq n)、Tn=∑k=1nYn,kT_n = \sum_{k=1}^n Y_{n,k}、μn=E[X1;∣X1∣≤n]\mu_n = E[X_1; \lvert X_1 \rvert \leq n] とおく。

(a) P(Sn≠Tn)≤∑k=1nP(∣Xk∣>n)=nP(∣X1∣>n)≤E[∣X1∣;∣X1∣>n]→0P(S_n \neq T_n) \leq \sum_{k=1}^n P(\lvert X_k \rvert > n) = n P(\lvert X_1 \rvert > n) \leq E[\lvert X_1 \rvert; \lvert X_1 \rvert > n] \to 0(優収束定理)。

(b) Yn,1,…,Yn,nY_{n,1}, \dots, Y_{n,n} は i.i.d. なので E[Tn]=nμnE[T_n] = n\mu_n、Var⁡(Tn)≤nE[X12;∣X1∣≤n]\operatorname{Var}(T_n) \leq n E[X_1^2; \lvert X_1 \rvert \leq n]。0<M≤n0 < M \leq n とすると、{∣X1∣≤M}\lbrace \lvert X_1 \rvert \leq M \rbrace では X12≤M2X_1^2 \leq M^2、{M<∣X1∣≤n}\lbrace M < \lvert X_1 \rvert \leq n \rbrace では X12≤n∣X1∣X_1^2 \leq n \lvert X_1 \rvert だから

Var⁡(Tn/n)≤1nE[X12;∣X1∣≤n]≤M2n+E[∣X1∣;∣X1∣>M]\operatorname{Var}(T_n/n) \leq \frac{1}{n} E[X_1^2; \lvert X_1 \rvert \leq n] \leq \frac{M^2}{n} + E[\lvert X_1 \rvert; \lvert X_1 \rvert > M]

n→∞n \to \infty、次に M→∞M \to \infty として Var⁡(Tn/n)→0\operatorname{Var}(T_n/n) \to 0 を得る。

(c) 優収束定理より μn→m\mu_n \to m。ε>0\varepsilon > 0 に対し nn が十分大きく ∣μn−m∣<ε/2\lvert \mu_n - m \rvert < \varepsilon/2 となれば

P(∣Sn/n−m∣>ε)≤P(Sn≠Tn)+P(∣Tn/n−μn∣>ε/2)≤P(Sn≠Tn)+4Var⁡(Tn/n)ε2→0□P(\lvert S_n/n - m \rvert > \varepsilon) \leq P(S_n \neq T_n) + P(\lvert T_n/n - \mu_n \rvert > \varepsilon/2) \leq P(S_n \neq T_n) + \frac{4\operatorname{Var}(T_n/n)}{\varepsilon^2} \to 0 \qquad \square

注意 2.5 証明の骨格は (a) nP(∣X1∣>n)→0n P(\lvert X_1 \rvert > n) \to 0 と (b) 切断した和の分散の評価であり、可積分性は (b) の評価と μn→m\mu_n \to m に使った。(b) を問題 2.7 のように評価し直すと、可積分性を仮定しなくても xP(∣X1∣>x)→0x P(\lvert X_1 \rvert > x) \to 0 だけから Sn/n−μn→P0S_n/n - \mu_n \xrightarrow{P} 0 が従う。この条件は可積分性より真に弱く、弱法則は成り立つが強法則は成り立たない例がある。

2.3 大数の強法則

弱法則は「nn を一つ固定すれば、Sn/nS_n/n が mm から離れている確率は小さい」と言うだけである。コイン投げを延々と続けたとき、一つの試行列の上で相対頻度が 1/21/2 に収束するか、という問いに答えるのが強法則である。まずボレル–カンテリの補題が直接使える場合を扱う。

定理 2.6(4 次モーメントによる強法則)X1,X2,…X_1, X_2, \dots は独立で、E[Xn]=mE[X_n] = m、E[(Xn−m)4]≤KE[(X_n - m)^4] \leq K(KK は定数)を満たすとする。このとき Sn/n→mS_n/n \to m a.s.。

証明. Xn−mX_n - m を改めて XnX_n とおき m=0m = 0 とする。Sn4=(∑iXi)4S_n^4 = \left(\sum_i X_i\right)^4 を展開すると、独立性と E[Xi]=0E[X_i] = 0 より、ある添字がちょうど 1 回だけ現れる項の期待値は 0 になる。残るのは Xi4X_i^4 と Xi2Xj2X_i^2 X_j^2(i≠ji \neq j)の項だけで

E[Sn4]=∑i=1nE[Xi4]+6∑i<jE[Xi2]E[Xj2]≤nK+3n(n−1)K≤3Kn2E[S_n^4] = \sum_{i=1}^n E[X_i^4] + 6\sum_{i < j} E[X_i^2] E[X_j^2] \leq nK + 3n(n-1)K \leq 3Kn^2

である(E[Xi2]2≤E[Xi4]≤KE[X_i^2]^2 \leq E[X_i^4] \leq K を使った)。マルコフの不等式より P(∣Sn/n∣>ε)≤E[Sn4]/(n4ε4)≤3K/(ε4n2)P(\lvert S_n/n \rvert > \varepsilon) \leq E[S_n^4]/(n^4 \varepsilon^4) \leq 3K/(\varepsilon^4 n^2) で、これは nn について総和有限である。ボレル–カンテリの第 1 補題より、a.s. で有限個の nn を除き ∣Sn/n∣≤ε\lvert S_n/n \rvert \leq \varepsilon。ε=1/l\varepsilon = 1/l(l∈Nl \in \mathbb{N})について確率 1 の事象の共通部分をとれば Sn/n→0S_n/n \to 0 a.s.。□\square

コイン投げのように有界な確率変数ならこれで十分である。一般の場合は、次の定理が最良である(逆が注意 2.8)。ここでは独立性を「対ごとの独立性」に弱めたエテマディ(N. Etemadi, 1981)の証明を与える。

定理 2.7(大数の強法則, strong law of large numbers)X1,X2,…X_1, X_2, \dots は対ごとに独立で同分布、E[∣X1∣]<∞E[\lvert X_1 \rvert] < \infty とし、m=E[X1]m = E[X_1] とする。このとき Sn/n→mS_n/n \to m a.s.。

証明. Xn=Xn+−Xn−X_n = X_n^{+} - X_n^{-} と分けると、(Xn+)n(X_n^{+})_n、(Xn−)n(X_n^{-})_n もそれぞれ対ごとに独立で同分布だから、Xn≥0X_n \geq 0 としてよい。

(1) 切断。Yk=Xk1{Xk≤k}Y_k = X_k \mathbf{1}_{\lbrace X_k \leq k \rbrace}、Tn=∑k=1nYkT_n = \sum_{k=1}^n Y_k とおく。命題 1.14 より ∑kP(Xk≠Yk)=∑kP(X1>k)≤E[X1]<∞\sum_k P(X_k \neq Y_k) = \sum_k P(X_1 > k) \leq E[X_1] < \infty なので、第 1 補題より a.s. で有限個の kk を除き Xk=YkX_k = Y_k となり、(Sn−Tn)/n→0(S_n - T_n)/n \to 0 a.s.。よって Tn/n→mT_n/n \to m a.s. を示せばよい。

(2) 分散の和の評価。x>0x > 0 について ∑k∈N,k≥xk−2≤2/x\sum_{k \in \mathbb{N}, k \geq x} k^{-2} \leq 2/x である(x≥1x \geq 1 なら k0=⌈x⌉k_0 = \lceil x \rceil として ∑k≥k0k−2≤k0−2+∫k0∞t−2 dt≤2/k0\sum_{k \geq k_0} k^{-2} \leq k_0^{-2} + \int_{k_0}^\infty t^{-2}\ dt \leq 2/k_0、0<x<10 < x < 1 なら ∑kk−2<2<2/x\sum_k k^{-2} < 2 < 2/x)。トネリの定理より

∑k=1∞E[Yk2]k2=E[X12∑k≥X11k2]≤E[X12⋅2X1;X1>0]=2E[X1]<∞\sum_{k=1}^\infty \frac{E[Y_k^2]}{k^2} = E\left[X_1^2 \sum_{k \geq X_1} \frac{1}{k^2}\right] \leq E\left[X_1^2 \cdot \frac{2}{X_1}; X_1 > 0\right] = 2E[X_1] < \infty

(3) 幾何級数的な部分列。α>1\alpha > 1 を固定し k(n)=⌊αn⌋k(n) = \lfloor \alpha^n \rfloor とおく。YkY_k は対ごとに独立なので無相関で、Var⁡(TN)=∑j≤NVar⁡(Yj)\operatorname{Var}(T_N) = \sum_{j \leq N} \operatorname{Var}(Y_j)。y≥1y \geq 1 で ⌊y⌋≥y/2\lfloor y \rfloor \geq y/2 だから ∑n ⁣:k(n)≥jk(n)−2≤∑n ⁣:αn≥j4α−2n≤4(1−α−2)j2\sum_{n \colon k(n) \geq j} k(n)^{-2} \leq \sum_{n \colon \alpha^n \geq j} 4\alpha^{-2n} \leq \frac{4}{(1 - \alpha^{-2}) j^2}。チェビシェフの不等式と和の順序交換により

∑n=1∞P(∣Tk(n)−E[Tk(n)]∣>εk(n))≤1ε2∑n1k(n)2∑j≤k(n)Var⁡(Yj)≤4ε2(1−α−2)∑jVar⁡(Yj)j2<∞\sum_{n=1}^\infty P\left(\lvert T_{k(n)} - E[T_{k(n)}] \rvert > \varepsilon k(n)\right) \leq \frac{1}{\varepsilon^2} \sum_{n} \frac{1}{k(n)^2} \sum_{j \leq k(n)} \operatorname{Var}(Y_j) \leq \frac{4}{\varepsilon^2 (1 - \alpha^{-2})} \sum_{j} \frac{\operatorname{Var}(Y_j)}{j^2} < \infty

第 1 補題と ε=1/l\varepsilon = 1/l の共通部分により (Tk(n)−E[Tk(n)])/k(n)→0(T_{k(n)} - E[T_{k(n)}])/k(n) \to 0 a.s.。単調収束定理より E[Yk]=E[X1;X1≤k]↑mE[Y_k] = E[X_1; X_1 \leq k] \uparrow m なので、チェザロ平均も E[Tn]/n→mE[T_n]/n \to m。よって Tk(n)/k(n)→mT_{k(n)}/k(n) \to m a.s.。

(4) 補間。Yk≥0Y_k \geq 0 より TlT_l は ll について単調非減少だから、k(n)≤l≤k(n+1)k(n) \leq l \leq k(n+1) のとき

k(n)k(n+1)⋅Tk(n)k(n)≤Tll≤Tk(n+1)k(n+1)⋅k(n+1)k(n)\frac{k(n)}{k(n+1)} \cdot \frac{T_{k(n)}}{k(n)} \leq \frac{T_l}{l} \leq \frac{T_{k(n+1)}}{k(n+1)} \cdot \frac{k(n+1)}{k(n)}

k(n+1)/k(n)→αk(n+1)/k(n) \to \alpha なので、a.s. で m/α≤lim inf⁡lTl/l≤lim sup⁡lTl/l≤αmm/\alpha \leq \liminf_l T_l/l \leq \limsup_l T_l/l \leq \alpha m。α=1+1/j\alpha = 1 + 1/j(j∈Nj \in \mathbb{N})について共通部分をとれば Tl/l→mT_l/l \to m a.s.。□\square

注意 2.8(強法則の逆)XnX_n が i.i.d. で E[∣X1∣]=∞E[\lvert X_1 \rvert] = \infty ならば lim sup⁡n∣Sn∣/n=∞\limsup_n \lvert S_n \rvert / n = \infty a.s. である。実際、問題 1.7 より lim sup⁡n∣Xn∣/n=∞\limsup_n \lvert X_n \rvert/n = \infty a.s. であり、∣Xn∣≤∣Sn∣+∣Sn−1∣\lvert X_n \rvert \leq \lvert S_n \rvert + \lvert S_{n-1} \rvert だから ∣Sn∣/n\lvert S_n \rvert / n が有界なら ∣Xn∣/n\lvert X_n \rvert / n も有界になってしまう。したがって i.i.d. 列について、Sn/nS_n/n が a.s. で有限の極限をもつことと E[∣X1∣]<∞E[\lvert X_1 \rvert] < \infty は同値である。たとえばコーシー分布の i.i.d. 列の標本平均は収束しない。

2.4 コルモゴロフの 0-1 法則

コイン投げで「相対頻度が収束する」という事象は、最初の有限回の結果をどう変えても起こるかどうかが変わらない。このような事象の確率は 0 か 1 に限られる。

定義 2.9(末尾 σ-加法族, tail σ-algebra)確率変数列 X1,X2,…X_1, X_2, \dots に対し Tn=σ(Xn+1,Xn+2,… )\mathcal{T}_n = \sigma(X_{n+1}, X_{n+2}, \dots)、T=⋂nTn\mathcal{T} = \bigcap_{n} \mathcal{T}_n とおく。T\mathcal{T} を末尾 σ-加法族、その元を末尾事象 (tail event) という。

例 2.10 {∑nXn が収束}\lbrace \sum_n X_n \text{ が収束} \rbrace、{Sn/n が収束}\lbrace S_n/n \text{ が収束} \rbrace は末尾事象であり、lim sup⁡nSn/n\limsup_n S_n/n は T\mathcal{T}-可測である(どの mm についても、(Sn−Sm)/n(S_n - S_m)/n の挙動で決まり、それは Tm\mathcal{T}_m-可測)。一方 {Sn>0 i.o.}\lbrace S_n > 0 \text{ i.o.} \rbrace は、X1X_1 の値を変えると変わりうるので一般には末尾事象ではない。

定理 2.11(コルモゴロフの 0-1 法則, Kolmogorov's zero-one law)X1,X2,…X_1, X_2, \dots が独立ならば、任意の A∈TA \in \mathcal{T} について P(A)=0P(A) = 0 または 11 である。

証明. An=σ(X1,…,Xn)\mathcal{A}_n = \sigma(X_1, \dots, X_n) とおく。グループ化(系 1.24)により An\mathcal{A}_n と Tn\mathcal{T}_n は独立で、T⊂Tn\mathcal{T} \subset \mathcal{T}_n だから T\mathcal{T} は各 An\mathcal{A}_n と独立である。An\mathcal{A}_n は増大列なので ⋃nAn\bigcup_n \mathcal{A}_n は π-系であり、定理 1.23 より T\mathcal{T} は σ(⋃nAn)=σ(X1,X2,… )\sigma\left(\bigcup_n \mathcal{A}_n\right) = \sigma(X_1, X_2, \dots) と独立である。T\mathcal{T} はこれに含まれるので、A∈TA \in \mathcal{T} は自分自身と独立で P(A)=P(A∩A)=P(A)2P(A) = P(A \cap A) = P(A)^2。□\square

系 2.12 XnX_n が独立ならば、T\mathcal{T}-可測な [−∞,∞][-\infty, \infty] 値確率変数 YY は a.s. で定数である。

証明. 各 xx について P(Y≤x)∈{0,1}P(Y \leq x) \in \lbrace 0, 1 \rbrace なので、c=inf⁡{x∈R∣P(Y≤x)=1}∈[−∞,∞]c = \inf\lbrace x \in \mathbb{R} \mid P(Y \leq x) = 1 \rbrace \in [-\infty, \infty](inf⁡∅=∞\inf \emptyset = \infty)とおけば Y=cY = c a.s.。□\square

したがって独立な確率変数の級数は、確率 1 で収束するか確率 1 で発散するかのどちらかである。どちらになるかを判定する道具が次節の定理である。

2.5 コルモゴロフの不等式と級数の収束

定理 2.13(コルモゴロフの不等式, Kolmogorov's inequality)X1,…,XnX_1, \dots, X_n は独立で E[Xk]=0E[X_k] = 0、E[Xk2]<∞E[X_k^2] < \infty とする。λ>0\lambda > 0 について

P(max⁡1≤k≤n∣Sk∣≥λ)≤Var⁡(Sn)λ2P\left(\max_{1 \leq k \leq n} \lvert S_k \rvert \geq \lambda\right) \leq \frac{\operatorname{Var}(S_n)}{\lambda^2}

チェビシェフの不等式は ∣Sn∣\lvert S_n \rvert だけを評価するが、この不等式は途中の最大値まで同じ上界で評価する。

証明. Ak={∣Sk∣≥λ, ∣Sj∣<λ (j<k)}A_k = \lbrace \lvert S_k \rvert \geq \lambda,\ \lvert S_j \rvert < \lambda\ (j < k) \rbrace は互いに交わらず、和集合が左辺の事象である。Sk1AkS_k \mathbf{1}_{A_k} は σ(X1,…,Xk)\sigma(X_1, \dots, X_k)-可測、Sn−SkS_n - S_k はそれと独立で平均 0 なので E[Sk1Ak(Sn−Sk)]=0E[S_k \mathbf{1}_{A_k}(S_n - S_k)] = 0。よって

E[Sn2]≥∑k=1nE[Sn2;Ak]≥∑k=1nE[Sk2+2Sk(Sn−Sk);Ak]=∑k=1nE[Sk2;Ak]≥λ2∑k=1nP(Ak)□E[S_n^2] \geq \sum_{k=1}^n E[S_n^2; A_k] \geq \sum_{k=1}^n E[S_k^2 + 2S_k(S_n - S_k); A_k] = \sum_{k=1}^n E[S_k^2; A_k] \geq \lambda^2 \sum_{k=1}^n P(A_k) \qquad \square

定理 2.14 X1,X2,…X_1, X_2, \dots は独立で E[Xn]=0E[X_n] = 0、∑nVar⁡(Xn)<∞\sum_n \operatorname{Var}(X_n) < \infty を満たすとする。このとき ∑nXn\sum_n X_n は a.s. で収束する。

証明. m<nm < n とし、Xm+1,…,XnX_{m+1}, \dots, X_n に定理 2.13 を適用して n→∞n \to \infty とすれば

P(sup⁡k>m∣Sk−Sm∣>ε)≤1ε2∑k>mVar⁡(Xk)P\left(\sup_{k > m} \lvert S_k - S_m \rvert > \varepsilon\right) \leq \frac{1}{\varepsilon^2}\sum_{k > m} \operatorname{Var}(X_k)

Wm=sup⁡j,k>m∣Sj−Sk∣W_m = \sup_{j, k > m} \lvert S_j - S_k \rvert は mm について単調非増加で、Wm≤2sup⁡k>m∣Sk−Sm∣W_m \leq 2\sup_{k > m} \lvert S_k - S_m \rvert。極限 W=lim⁡mWmW = \lim_m W_m について P(W>2ε)≤P(Wm>2ε)≤ε−2∑k>mVar⁡(Xk)→0P(W > 2\varepsilon) \leq P(W_m > 2\varepsilon) \leq \varepsilon^{-2}\sum_{k > m} \operatorname{Var}(X_k) \to 0 なので W=0W = 0 a.s.、すなわち (Sn)(S_n) は a.s. でコーシー列である。□\square

例 2.15(ランダムな符号の調和級数)εn\varepsilon_n を P(εn=±1)=1/2P(\varepsilon_n = \pm 1) = 1/2 の i.i.d. とすると、∑n1/n2<∞\sum_n 1/n^2 < \infty なので ∑nεn/n\sum_n \varepsilon_n / n は a.s. で収束する。調和級数は発散し交代調和級数は収束するが、符号を公平なコインで決めると確率 1 で収束する。

定理 2.16(コルモゴロフの 3 級数定理, three-series theorem)X1,X2,…X_1, X_2, \dots を独立とし、A>0A > 0、Yn=Xn1{∣Xn∣≤A}Y_n = X_n \mathbf{1}_{\lbrace \lvert X_n \rvert \leq A \rbrace} とおく。∑nXn\sum_n X_n が a.s. で収束するための必要十分条件は、次の三つの級数がすべて収束することである。

(i) ∑nP(∣Xn∣>A),(ii) ∑nE[Yn],(iii) ∑nVar⁡(Yn)\text{(i) } \sum_n P(\lvert X_n \rvert > A), \qquad \text{(ii) } \sum_n E[Y_n], \qquad \text{(iii) } \sum_n \operatorname{Var}(Y_n)

証明の概略. 十分性:(i) と第 1 補題より a.s. で有限個を除き Xn=YnX_n = Y_n。(iii) と定理 2.14 より ∑n(Yn−E[Yn])\sum_n (Y_n - E[Y_n]) は a.s. で収束し、(ii) と合わせて ∑nYn\sum_n Y_n、したがって ∑nXn\sum_n X_n が a.s. で収束する(ここまでは完全な証明である)。必要性:∑nXn\sum_n X_n が収束すれば Xn→0X_n \to 0 なので、第 2 補題の対偶から (i) が従う。(iii) は、独立なコピー Yn′Y_n' との差 Yn−Yn′Y_n - Y_n'(対称化)を考え、有界な項の和が収束するなら分散の和も有限であること(コルモゴロフの不等式の逆向きの評価、または中心極限定理による)を示して得る。(ii) は (iii) と定理 2.14 から従う。詳細は Durrett の教科書を参照。□\square

2.6 応用

モンテカルロ法

f ⁣:[0,1]d→Rf\colon [0,1]^d \to \mathbb{R} が可積分なら、[0,1]d[0,1]^d 上の一様分布に従う i.i.d. U1,U2,…U_1, U_2, \dots について、強法則より

1n∑k=1nf(Uk)→∫[0,1]df(x) dxa.s.\frac{1}{n}\sum_{k=1}^n f(U_k) \to \int_{[0,1]^d} f(x)\ dx \quad \text{a.s.}

例 2.17(π\pi の推定)d=2d = 2、f=4⋅1{x12+x22≤1}f = 4 \cdot \mathbf{1}_{\lbrace x_1^2 + x_2^2 \leq 1 \rbrace} とすれば積分は π\pi である。Var⁡(f(U))=16⋅π4(1−π4)=4π−π2≈2.70\operatorname{Var}(f(U)) = 16 \cdot \frac{\pi}{4}\left(1 - \frac{\pi}{4}\right) = 4\pi - \pi^2 \approx 2.70 なので、チェビシェフの不等式より n=106n = 10^6 点で P(∣π^n−π∣≥0.01)≤2.70/(106⋅10−4)=0.027P(\lvert \hat{\pi}_n - \pi \rvert \geq 0.01) \leq 2.70/(10^6 \cdot 10^{-4}) = 0.027。誤差は Var⁡(f(U))/n\sqrt{\operatorname{Var}(f(U))/n} の程度で、次元 dd によらない。格子点を使う数値積分では、1 軸あたり N1/dN^{1/d} 点しか置けず高次元で精度が急速に落ちるのと対照的である。

ボレルの正規数定理

x∈[0,1)x \in [0, 1) の bb 進展開(b≥2b \geq 2)で、長さ kk のすべての数字列 ww が極限頻度 b−kb^{-k} で現れるとき、xx は bb 進正規数 (normal number) であるという。

定理 2.18(ボレルの正規数定理, Borel 1909)ルベーグ測度についてほとんどすべての x∈[0,1)x \in [0,1) は、すべての b≥2b \geq 2 について bb 進正規数である。

証明. bb を固定する。例 1.32 と同様に、bb 進展開の各桁 d1,d2,…d_1, d_2, \dots はルベーグ測度のもとで {0,…,b−1}\lbrace 0, \dots, b-1 \rbrace 上の一様分布に従う i.i.d. である。長さ kk の数字列 ww を固定し、ηj=1{(dj,…,dj+k−1)=w}\eta_j = \mathbf{1}_{\lbrace (d_j, \dots, d_{j+k-1}) = w \rbrace} とおく。ηj\eta_j たちは重なりのために独立ではないが、剰余 r∈{1,…,k}r \in \lbrace 1, \dots, k \rbrace ごとに ηr,ηr+k,ηr+2k,…\eta_{r}, \eta_{r+k}, \eta_{r+2k}, \dots は重ならないブロックで決まるので Be⁡(b−k)\operatorname{Be}(b^{-k}) の i.i.d. である。定理 2.6 より各 rr について 1N∑i=0N−1ηr+ik→b−k\frac{1}{N}\sum_{i=0}^{N-1} \eta_{r + ik} \to b^{-k} a.s. であり、j≤nj \leq n を剰余で kk 組に分けて和をとれば 1n∑j=1nηj→b−k\frac{1}{n}\sum_{j=1}^n \eta_j \to b^{-k} a.s.。bb、kk、ww は可算個なので、例外集合の和も零集合である。□\square

ほとんどすべての数が正規数であるにもかかわらず、2\sqrt{2}、π\pi、ee が 10 進正規数かどうかは知られていない。確率論的な存在証明は、具体例を与えずに「ほとんどすべて」を示せるという特徴をもつ。

ワイエルシュトラスの近似定理

定理 2.19(バーンスタイン多項式による近似)f∈C([0,1])f \in C([0,1]) に対し

Bnf(p)=∑k=0nf(kn)(nk)pk(1−p)n−kB_n f(p) = \sum_{k=0}^n f\left(\frac{k}{n}\right)\binom{n}{k} p^k (1-p)^{n-k}

とおくと、Bnf→fB_n f \to f は [0,1][0,1] 上一様収束する。特に連続関数は多項式で一様近似できる。

証明. Sn∼B(n,p)S_n \sim B(n, p) とすると Bnf(p)=E[f(Sn/n)]B_n f(p) = E[f(S_n/n)] である(コインの表の出る確率 pp を相対頻度 Sn/nS_n/n で推定して ff に代入する)。ff は一様連続なので、ε>0\varepsilon > 0 に対し ∣x−y∣<δ\lvert x - y \rvert < \delta なら ∣f(x)−f(y)∣<ε\lvert f(x) - f(y) \rvert < \varepsilon となる δ>0\delta > 0 がある。M=max⁡∣f∣M = \max \lvert f \rvert とすると、チェビシェフの不等式と p(1−p)≤1/4p(1-p) \leq 1/4 より

∣Bnf(p)−f(p)∣≤E[∣f(Sn/n)−f(p)∣]≤ε+2MP(∣Snn−p∣≥δ)≤ε+2Mp(1−p)nδ2≤ε+M2nδ2\lvert B_n f(p) - f(p) \rvert \leq E\left[\lvert f(S_n/n) - f(p) \rvert\right] \leq \varepsilon + 2M P\left(\left\lvert \frac{S_n}{n} - p \right\rvert \geq \delta\right) \leq \varepsilon + \frac{2M p(1-p)}{n\delta^2} \leq \varepsilon + \frac{M}{2n\delta^2}

右辺は pp によらないので lim sup⁡nmax⁡p∣Bnf(p)−f(p)∣≤ε\limsup_n \max_p \lvert B_n f(p) - f(p) \rvert \leq \varepsilon。□\square

解析的な証明は 01-calculus 第6章 も参照。確率論的な証明は近似多項式を具体的に与える点に特徴がある。

2.7 経験分布関数とグリヴェンコ–カンテリの定理

未知の分布 FF に従う i.i.d. な観測値 X1,…,XnX_1, \dots, X_n から FF を推定したい。自然な推定量は経験分布関数 (empirical distribution function)

Fn(x)=1n∑k=1n1{Xk≤x}F_n(x) = \frac{1}{n}\sum_{k=1}^n \mathbf{1}_{\lbrace X_k \leq x \rbrace}

である。各 xx を固定すれば、強法則より Fn(x)→F(x)F_n(x) \to F(x) a.s.。しかし例外集合は xx に依存し、xx は非可算個あるので、これだけでは「関数として」の収束は言えない。

定理 2.20(グリヴェンコ–カンテリの定理, Glivenko–Cantelli theorem)X1,X2,…X_1, X_2, \dots を分布関数 FF をもつ i.i.d. とすると

sup⁡x∈R∣Fn(x)−F(x)∣→0a.s.\sup_{x \in \mathbb{R}} \lvert F_n(x) - F(x) \rvert \to 0 \quad \text{a.s.}

証明. k∈Nk \in \mathbb{N} を固定し、x0=−∞x_0 = -\infty、xk=∞x_k = \infty、1≤j≤k−11 \leq j \leq k-1 について xj=inf⁡{y∣F(y)≥j/k}x_j = \inf\lbrace y \mid F(y) \geq j/k \rbrace とおく。右連続性から F(xj)≥j/kF(x_j) \geq j/k、下限の定義から F(xj−)≤j/kF(x_j -) \leq j/k である(F(−∞)=0F(-\infty) = 0、F(∞−)=1F(\infty -) = 1 と約束する)。よって各 jj で F(xj−)−F(xj−1)≤1/kF(x_j -) - F(x_{j-1}) \leq 1/k。強法則を 1{Xi≤xj}\mathbf{1}_{\lbrace X_i \leq x_j \rbrace} と 1{Xi<xj}\mathbf{1}_{\lbrace X_i < x_j \rbrace} に適用すると、a.s. で

εn:=max⁡j(∣Fn(xj)−F(xj)∣+∣Fn(xj−)−F(xj−)∣)→0\varepsilon_n := \max_{j} \left(\lvert F_n(x_j) - F(x_j) \rvert + \lvert F_n(x_j -) - F(x_j -) \rvert\right) \to 0

x∈[xj−1,xj)x \in [x_{j-1}, x_j) とすると、単調性から

Fn(x)≤Fn(xj−)≤F(xj−)+εn≤F(xj−1)+1k+εn≤F(x)+1k+εnF_n(x) \leq F_n(x_j -) \leq F(x_j -) + \varepsilon_n \leq F(x_{j-1}) + \frac{1}{k} + \varepsilon_n \leq F(x) + \frac{1}{k} + \varepsilon_n

同様に Fn(x)≥Fn(xj−1)≥F(xj−1)−εn≥F(xj−)−1k−εn≥F(x)−1k−εnF_n(x) \geq F_n(x_{j-1}) \geq F(x_{j-1}) - \varepsilon_n \geq F(x_j -) - \frac{1}{k} - \varepsilon_n \geq F(x) - \frac{1}{k} - \varepsilon_n。よって lim sup⁡nsup⁡x∣Fn(x)−F(x)∣≤1/k\limsup_n \sup_x \lvert F_n(x) - F(x) \rvert \leq 1/k a.s.。k∈Nk \in \mathbb{N} について共通部分をとればよい。□\square

この定理は「データを集めれば真の分布がわかる」ことの数学的な保証であり、統計学の基本定理とも呼ばれる。誤差 nsup⁡x∣Fn−F∣\sqrt{n}\sup_x \lvert F_n - F \rvert の極限分布(コルモゴロフ–スミルノフ検定の基礎)は第7章のブラウン運動(ブラウン橋)と関係する。

まとめ

  • 概収束(a.e. 収束)は標本路ごとの収束、確率収束(測度収束)は各時刻での分布の集中であり、前者の方が強い。
  • L2L^2 弱法則はチェビシェフの不等式から、一般の弱法則は切断から従う。
  • 4 次モーメントがあればボレル–カンテリの補題で強法則が直ちに示せる。一般の強法則(E[∣X1∣]<∞E[\lvert X_1 \rvert] < \infty)は、切断・幾何級数的な部分列・単調性による補間で示す(エテマディ)。i.i.d. では可積分性は強法則の必要十分条件である。
  • 独立列の末尾事象の確率は 0 か 1(コルモゴロフの 0-1 法則)。
  • コルモゴロフの不等式から、分散の和が有限な独立級数は a.s. で収束する。収束の判定は 3 級数定理で完全に与えられる。
  • 応用:モンテカルロ法(誤差は次元によらず n−1/2n^{-1/2} 程度)、ほとんどすべての実数は正規数、バーンスタイン多項式による一様近似、経験分布関数の一様収束。

演習問題

問題 2.1 ★ 公平なコインを nn 回投げる。チェビシェフの不等式を用いて、表の相対頻度と 1/21/2 との差が 0.010.01 未満である確率が 0.950.95 以上となることを保証するには、nn をいくつ以上にすればよいか。

解答

P(∣Sn/n−1/2∣≥0.01)≤1/4n⋅10−4=2500nP(\lvert S_n/n - 1/2 \rvert \geq 0.01) \leq \frac{1/4}{n \cdot 10^{-4}} = \frac{2500}{n} なので、2500/n≤0.052500/n \leq 0.05、すなわち n≥50000n \geq 50000 ならよい。(中心極限定理を使えば約 96009600 回で足りることがわかる。チェビシェフの不等式は安全側の粗い評価である。)

問題 2.2 ★ X1,X2,…X_1, X_2, \dots を i.i.d. で P(X1=±1)=1/2P(X_1 = \pm 1) = 1/2 とする。lim sup⁡nSn/n\limsup_n S_n/\sqrt{n} は a.s. で定数(±∞\pm\infty も許す)であることを示せ。

解答

mm を固定すると Sm/n→0S_m/\sqrt{n} \to 0 なので lim sup⁡nSn/n=lim sup⁡n(Sn−Sm)/n\limsup_n S_n/\sqrt{n} = \limsup_n (S_n - S_m)/\sqrt{n} であり、右辺は σ(Xm+1,Xm+2,… )\sigma(X_{m+1}, X_{m+2}, \dots)-可測である。mm は任意なので T\mathcal{T}-可測で、系 2.12 より a.s. で定数である。(実は +∞+\infty である。第4章 問題 4.3 を参照。)

問題 2.3 ★★ XnX_n を i.i.d. で E[X1+]=∞E[X_1^{+}] = \infty、E[X1−]<∞E[X_1^{-}] < \infty とする。Sn/n→∞S_n/n \to \infty a.s. を示せ。

解答

M>0M > 0 について Xn∧MX_n \wedge M は i.i.d. で可積分(∣Xn∧M∣≤M+Xn−\lvert X_n \wedge M \rvert \leq M + X_n^{-})だから、強法則より 1n∑k≤n(Xk∧M)→E[X1∧M]\frac{1}{n}\sum_{k \leq n} (X_k \wedge M) \to E[X_1 \wedge M] a.s.。Sn≥∑k≤n(Xk∧M)S_n \geq \sum_{k \leq n}(X_k \wedge M) なので lim inf⁡nSn/n≥E[X1∧M]\liminf_n S_n/n \geq E[X_1 \wedge M] a.s.。E[X1∧M]=E[X1+∧M]−E[X1−]→∞E[X_1 \wedge M] = E[X_1^{+} \wedge M] - E[X_1^{-}] \to \infty(M→∞M \to \infty、単調収束定理)なので、M∈NM \in \mathbb{N} について共通部分をとればよい。

問題 2.4 ★★ (1)(クロネッカーの補題)実数列 ana_n について ∑nan/n\sum_n a_n/n が収束するならば 1n∑k=1nak→0\frac{1}{n}\sum_{k=1}^n a_k \to 0 であることを示せ。(2) XnX_n が独立で E[Xn]=0E[X_n] = 0、∑nVar⁡(Xn)/n2<∞\sum_n \operatorname{Var}(X_n)/n^2 < \infty ならば Sn/n→0S_n/n \to 0 a.s. であることを示せ。

解答

(1) b0=0b_0 = 0、bn=∑k≤nak/k→bb_n = \sum_{k \leq n} a_k/k \to b とおく。ak=k(bk−bk−1)a_k = k(b_k - b_{k-1}) だから、アーベルの総和法により

1n∑k=1nak=1n(nbn−∑k=0n−1bk)=bn−1n∑k=0n−1bk→b−b=0\frac{1}{n}\sum_{k=1}^n a_k = \frac{1}{n}\left(n b_n - \sum_{k=0}^{n-1} b_k\right) = b_n - \frac{1}{n}\sum_{k=0}^{n-1} b_k \to b - b = 0

(チェザロ平均は元の数列と同じ極限をもつ)。(2) Xn/nX_n/n に定理 2.14 を適用すると ∑nXn/n\sum_n X_n/n は a.s. で収束し、(1) より Sn/n→0S_n/n \to 0 a.s.。

問題 2.5 ★★ εn\varepsilon_n を P(εn=±1)=1/2P(\varepsilon_n = \pm 1) = 1/2 の i.i.d. とする。∑nεnn−α\sum_n \varepsilon_n n^{-\alpha} が a.s. で収束するための必要十分条件は α>1/2\alpha > 1/2 であることを示せ。

解答

α>1/2\alpha > 1/2 なら ∑nVar⁡(εnn−α)=∑nn−2α<∞\sum_n \operatorname{Var}(\varepsilon_n n^{-\alpha}) = \sum_n n^{-2\alpha} < \infty なので定理 2.14 より収束する。α≤1/2\alpha \leq 1/2 のとき 3 級数定理を A=2A = 2 で適用すると、Yn=Xn=εnn−αY_n = X_n = \varepsilon_n n^{-\alpha} で (i) (ii) は成り立つが (iii) ∑nn−2α=∞\sum_n n^{-2\alpha} = \infty が成り立たない。よって a.s. 収束はせず、0-1 法則により確率 1 で発散する。

問題 2.6 ★★(再生定理)ξ1,ξ2,…\xi_1, \xi_2, \dots を正の値をとる i.i.d. で E[ξ1]=μ<∞E[\xi_1] = \mu < \infty とし、Tn=ξ1+⋯+ξnT_n = \xi_1 + \cdots + \xi_n(電球を nn 個使い切る時刻)、N(t)=max⁡{n≥0∣Tn≤t}N(t) = \max\lbrace n \geq 0 \mid T_n \leq t \rbrace(T0=0T_0 = 0)とする。N(t)/t→1/μN(t)/t \to 1/\mu a.s.(t→∞t \to \infty)を示せ。

解答

強法則より Tn/n→μT_n/n \to \mu a.s.。この事象の上では Tn→∞T_n \to \infty なので N(t)<∞N(t) < \infty、また各 Tn<∞T_n < \infty より N(t)→∞N(t) \to \infty(t→∞t \to \infty)。N(t)≥1N(t) \geq 1 のとき TN(t)≤t<TN(t)+1T_{N(t)} \leq t < T_{N(t)+1} だから

TN(t)N(t)≤tN(t)<TN(t)+1N(t)+1⋅N(t)+1N(t)\frac{T_{N(t)}}{N(t)} \leq \frac{t}{N(t)} < \frac{T_{N(t)+1}}{N(t)+1} \cdot \frac{N(t)+1}{N(t)}

両端は μ\mu に収束するので t/N(t)→μt/N(t) \to \mu。

問題 2.7 ★★★ (1) XnX_n が i.i.d. で xP(∣X1∣>x)→0x P(\lvert X_1 \rvert > x) \to 0(x→∞x \to \infty)ならば、μn=E[X1;∣X1∣≤n]\mu_n = E[X_1; \lvert X_1 \rvert \leq n] について Sn/n−μn→P0S_n/n - \mu_n \xrightarrow{P} 0 を示せ。(2) X1X_1 が対称分布で x≥ex \geq e について P(∣X1∣>x)=e/(xlog⁡x)P(\lvert X_1 \rvert > x) = e/(x \log x) を満たすとき、Sn/n→P0S_n/n \xrightarrow{P} 0 だが Sn/nS_n/n は a.s. 収束しないことを示せ。

解答

(1) 定理 2.4 の証明と同じ記号を使う。(a) P(Sn≠Tn)≤nP(∣X1∣>n)→0P(S_n \neq T_n) \leq nP(\lvert X_1 \rvert > n) \to 0。(b) 命題 1.14 の証明と同様に E[X12;∣X1∣≤n]=∫0∞2yP(y<∣X1∣≤n) dy≤∫0n2yP(∣X1∣>y) dyE[X_1^2; \lvert X_1 \rvert \leq n] = \int_0^\infty 2y P(y < \lvert X_1 \rvert \leq n)\ dy \leq \int_0^n 2y P(\lvert X_1 \rvert > y)\ dy。g(y)=2yP(∣X1∣>y)g(y) = 2yP(\lvert X_1 \rvert > y) は有界で y→∞y \to \infty で 0 に収束するので、1n∫0ng(y) dy→0\frac{1}{n}\int_0^n g(y)\ dy \to 0。よって Var⁡(Tn/n)≤1nE[X12;∣X1∣≤n]→0\operatorname{Var}(T_n/n) \leq \frac{1}{n}E[X_1^2; \lvert X_1 \rvert \leq n] \to 0。(a)(b) とチェビシェフの不等式から P(∣Sn/n−μn∣>ε)≤P(Sn≠Tn)+Var⁡(Tn/n)/ε2→0P(\lvert S_n/n - \mu_n \rvert > \varepsilon) \leq P(S_n \neq T_n) + \operatorname{Var}(T_n/n)/\varepsilon^2 \to 0。

(2) xP(∣X1∣>x)=e/log⁡x→0xP(\lvert X_1 \rvert > x) = e/\log x \to 0 で、対称性より μn=0\mu_n = 0 だから (1) より Sn/n→P0S_n/n \xrightarrow{P} 0。一方 E[∣X1∣]=∫0∞P(∣X1∣>x) dx≥∫e∞e dxxlog⁡x=∞E[\lvert X_1 \rvert] = \int_0^\infty P(\lvert X_1 \rvert > x)\ dx \geq \int_e^\infty \frac{e\ dx}{x \log x} = \infty なので、注意 2.8 より lim sup⁡n∣Sn∣/n=∞\limsup_n \lvert S_n \rvert/n = \infty a.s.。

この章を読み終えたら

「読了」にすると学習記録とロードマップに反映されます。演習の自己採点もお忘れなく。

この章の誤りを報告GitHub で見る