tecktom.com

Lévy の連続性定理

2026年10月15日25分で読めますAI生成教科書『基礎理論』第13章 中心極限定理
目次
  1. なぜ特性関数を使うか
  2. 分布収束から特性関数の収束へ
  3. Helly の選出定理と緊密性
  4. 連続性定理
  5. 複素数の積の評価
  6. 確認テスト
  7. 木下の層

AI 生成・木下が未勉強

確率論入門 第13章 §13.1 ★★★ 特性関数が各点で収束し、極限が 0 で連続なら、分布が収束する。中心極限定理を特性関数の計算に帰着させる道具

第12章の大数の法則は、独立同分布で期待値 μ\mu をもつ列について、標本平均 Xˉn\bar X_n が期待値に収束することを示した。しかしこの主張は、ずれ Xˉn−μ\bar X_n-\mu が 0 に近づくことを述べるだけで、ずれがどの程度の大きさで、どのような形に分布するかを述べない。分散 σ2\sigma^2 が有限なら、標本平均の分散は σ2/n\sigma^2/n である。したがってずれの大きさはおよそ σ/n\sigma/\sqrt n である。

そこでずれを n\sqrt n 倍に拡大した

Zn=n(Xˉn−μ)σ=Sn−nμσn,Sn=X1+⋯+Xn(13.1)Z_n=\frac{\sqrt n(\bar X_n-\mu)}{\sigma}=\frac{S_n-n\mu}{\sigma\sqrt n},\qquad S_n=X_1+\cdots+X_n\tag{13.1}

を考える。すべての nn で次が成り立つ。

E[Zn]=0,V[Zn]=1.E[Z_n]=0,\qquad V[Z_n]=1.

中心極限定理は、元の XkX_k の分布が何であっても(分散が有限でありさえすれば)、ZnZ_n の分布が標準正規分布に近づくという主張である。この普遍性のために、正規分布は測定誤差・標本平均・多数の小さな寄与の和の近似分布として、確率と統計のあらゆる場面に現れる。

この章で作るものは4つある。

  1. §13.1 では、特性関数の各点収束から分布収束を導く Lévy の連続性定理を証明する。
  2. §13.2 では、この定理と第8章の特性関数の Taylor 展開を組み合わせて中心極限定理を証明する。同分布でない場合の Lindeberg 条件も述べて証明する。
  3. §13.3 では、確率ベクトルの分布収束を1次元の射影へ帰着させる Cramér–Wold の方法を示す。そこから多変量の中心極限定理と Pearson のカイ二乗統計量の極限分布を導く。
  4. §13.4 では、推定量 TnT_n が近似的に正規分布に従うとき、その関数 g(Tn)g(T_n) も近似的に正規分布に従うことを示すデルタ法を扱う。

前の章から借りる結果を、式の形で先に書いておく。第8章 §8.3 から、特性関数 φX(t)=E[eitX]\varphi_X(t)=E[\mathrm{e}^{\mathrm{i} tX}] は次を満たす。

φX(0)=1,∣φX(t)∣≤1,φX+Y=φXφY(X,Y が独立のとき).\varphi_X(0)=1,\qquad\left\vert \varphi_X(t)\right\vert \le1,\qquad\varphi_{X+Y}=\varphi_X\varphi_Y\quad(X,Y\text{ が独立のとき}).

また E[X2]<∞E[X^2]<\infty なら

φX(t)=1+itE[X]−t22E[X2]+o(t2)(t→0)(13.2)\varphi_X(t)=1+\mathrm{i} tE[X]-\frac{t^2}{2}E[X^2]+o(t^2)\qquad(t\to0)\tag{13.2}

が成り立ち、実数 xx について剰余の評価

∣eix−(1+ix−x22)∣≤min⁡(∣x∣36, x2)(13.3)\Bigl|\mathrm{e}^{\mathrm{i} x}-\Bigl(1+\mathrm{i} x-\frac{x^2}{2}\Bigr)\Bigr|\le\min\Bigl(\frac{\left\vert x\right\vert ^3}{6},\ x^2\Bigr)\tag{13.3}

が成り立つ(第8章の補題で n=2n=2 とした場合)。第9章 §9.2 の一意性定理から、特性関数が等しい2つの分布は等しい。これは Rd\mathbb{R}^d 値の確率ベクトルについても成り立ち、Nd(μ,Σ)\mathrm N_d(\bm\mu,\Sigma) の特性関数は

φX(t)=E[eit⊤X]=exp⁡(it⊤μ−12t⊤Σt)(13.4)\varphi_{\bm X}(\bm t)=E\bigl[\mathrm{e}^{\mathrm{i}\bm t^\top\bm X}\bigr]=\exp\Bigl(\mathrm{i}\bm t^\top\bm\mu-\frac12\bm t^\top\Sigma\bm t\Bigr)\tag{13.4}

である(第9章 §9.2)。第6章 §6.2 からは、Z∼Nd(μ,Σ)\bm Z\sim\mathrm N_d(\bm\mu,\Sigma) のとき

t⊤Z∼N(t⊤μ, t⊤Σt)\bm t^\top\bm Z\sim\mathrm N(\bm t^\top\bm\mu,\ \bm t^\top\Sigma\bm t)

であることを借りる(分散が 0 のときは定数 t⊤μ\bm t^\top\bm\mu と読む)。収束の言葉は第11章のとおりである。とくに §11.4 の分布収束の言い換え

Xn→dX  ⟺  すべての有界連続関数 g について E[g(Xn)]→E[g(X)](13.5)X_n\xrightarrow{d}X\iff\text{すべての有界連続関数 }g\text{ について }E[g(X_n)]\to E[g(X)]\tag{13.5}

を使い、Rd\mathbb{R}^d 値の確率ベクトルの分布収束は右辺(g:Rd→Rg:\mathbb{R}^d\to\mathbb{R})で定義する。ほかに次の3つも使う。

  • §11.4 の連続写像定理と Slutsky の定理。
  • 分布収束する列は緊密である(定義 13.3)。
  • §11.2 の「定数への分布収束は確率収束を導く」こと。

標準正規分布の密度と分布関数は第6章と同じく ϕ\phi、Φ\Phi と書く。

なぜ特性関数を使うか

独立な確率変数の和 SnS_n の分布は、密度の畳み込み(第7章 §7.3)を n−1n-1 回くり返して得られる。この計算を一般の分布について実行し、n→∞n\to\infty の極限をとることは難しい。一方、特性関数では畳み込みが積 φSn=φX1⋯φXn\varphi_{S_n}=\varphi_{X_1}\cdots\varphi_{X_n} に変わる。積の極限は対数をとれば和の極限になるので、計算が格段に易しい。

残る問題は「特性関数が収束すれば分布も収束する」ことを保証することであり、それがこの節の Lévy の連続性定理である。第9章の一意性定理は「特性関数が等しければ分布が等しい」ことを述べたが、ここで必要なのは「特性関数が近ければ分布が近い」という連続性の主張であり、一意性だけからは出てこない。

分布収束から特性関数の収束へ

易しい向きから示す。

命題 13.1

Xn→dXX_n\xrightarrow{d}X ならば、すべての t∈Rt\in\mathbb{R} で次が成り立つ。

φXn(t)→φX(t).\varphi_{X_n}(t)\to\varphi_X(t).

証明

tt を固定する。x↦cos⁡txx\mapsto\cos tx と x↦sin⁡txx\mapsto\sin tx は有界連続関数である。(13.5) から

E[cos⁡tXn]→E[cos⁡tX],E[sin⁡tXn]→E[sin⁡tX]E[\cos tX_n]\to E[\cos tX],\qquad E[\sin tX_n]\to E[\sin tX]

である。φXn(t)=E[cos⁡tXn]+iE[sin⁡tXn]\varphi_{X_n}(t)=E[\cos tX_n]+\mathrm{i} E[\sin tX_n] なので結論を得る。

逆向きは無条件には成り立たない。図 13.1 の例がそれを示す。

一様分布が広がっていく列

Xn∼U(−n,n)X_n\sim\mathrm U(-n,n) とすると、t≠0t\ne0 で

φXn(t)=12n∫−nneitx dx=sin⁡ntnt\varphi_{X_n}(t)=\frac1{2n}\int_{-n}^{n}\mathrm{e}^{\mathrm{i} tx}\,\mathrm{d} x=\frac{\sin nt}{nt}

であり、∣φXn(t)∣≤1n∣t∣→0\left\vert \varphi_{X_n}(t)\right\vert \le\dfrac1{n\left\vert t\right\vert }\to0 である。φXn(0)=1\varphi_{X_n}(0)=1 なので、各点の極限は

φ(t)={1(t=0)0(t≠0)\varphi(t)=\begin{cases}1&(t=0)\\0&(t\ne0)\end{cases}

となる。この φ\varphi は t=0t=0 で連続でないので、どの確率分布の特性関数でもない(第8章 §8.3 で示したとおり、特性関数は一様連続である)。実際、どんな M>0M>0 についても

P(∣Xn∣≤M)=Mn→0P(\left\vert X_n\right\vert \le M)=\frac Mn\to0

である。確率の質量が無限遠へ逃げていくので、XnX_n は分布収束しない。

図

図 13.1 Xn∼U(−n,n)X_n\sim\mathrm U(-n,n) の特性関数 sin⁡(nt)/(nt)\sin(nt)/(nt)。nn が大きくなると山の幅が 1/n1/n の速さで狭まり、t≠0t\ne0 では 0 に近づくが、t=0t=0 では常に 1 である。極限は t=0t=0 で不連続になり、特性関数ではない。

この例が示すように、特性関数が各点収束しても、分布の列が緊密でなければ極限の分布は存在しない。逆向きの証明は、緊密性を特性関数の 0 の近くの振る舞いから引き出すことが要になる。

Helly の選出定理と緊密性

分布関数の列から、収束する部分列を取り出す定理から始める。極限は確率分布の分布関数になるとは限らず、全質量が 1 より小さい「劣分布関数」になることがある。

定理 13.2(Helly の選出定理)

F1,F2,…F_1,F_2,\dots を分布関数の列とする。次を満たす部分列 FnjF_{n_j} と関数 F:R→[0,1]F:\mathbb{R}\to[0,1] が存在する。

  1. (i) FF は非減少で右連続である。
  2. (ii) FF のすべての連続点で、Fnj(x)→F(x)F_{n_j}(x)\to F(x) となる。
証明

有理数を q1,q2,…q_1,q_2,\dots と並べる。Fn(q1)F_n(q_1) は [0,1][0,1] の有界な数列なので、収束する部分列をもつ(Bolzano–Weierstrass の定理)。その部分列の中から Fn(q2)F_n(q_2) が収束する部分列をとる。さらにその中から Fn(q3)F_n(q_3) が収束する部分列をとる、という操作をくり返す。

各回の部分列から、その回数番目の項を取り出して並べる(kk 回目の部分列からは kk 番目の項)。この対角線の列を n1<n2<⋯n_1<n_2<\cdots とする。

qiq_i を1つ固定する。対角線の列の ii 番目から先は、ii 回目にとった部分列の部分列である。したがって Fnj(qi)F_{n_j}(q_i) は j→∞j\to\infty で収束する。

その極限を G(qi)G(q_i) と置く。GG は有理数の上で非減少で、値は [0,1][0,1] にある。

次のように FF を定める。

F(x)=inf⁡{G(q):q∈Q, q>x}.F(x)=\inf\{G(q):q\in\mathbb Q,\ q>x\}.

FF は非減少である。右連続性を示す。xk↓xx_k\downarrow x とすると F(xk)≥F(x)F(x_k)\ge F(x) である。

一方、有理数 q>xq>x を1つとる。十分大きな kk で xk<qx_k<q となる。このとき F(xk)≤G(q)F(x_k)\le G(q) なので、lim⁡kF(xk)≤G(q)\lim_kF(x_k)\le G(q) である。q↓xq\downarrow x について下限をとれば lim⁡kF(xk)≤F(x)\lim_kF(x_k)\le F(x) となる。

xx を FF の連続点とする。有理数 r<x<sr<x<s をとると、次の不等式

Fnj(r)≤Fnj(x)≤Fnj(s)F_{n_j}(r)\le F_{n_j}(x)\le F_{n_j}(s)

から

G(r)≤lim inf⁡jFnj(x)≤lim sup⁡jFnj(x)≤G(s)G(r)\le\liminf_jF_{n_j}(x)\le\limsup_jF_{n_j}(x)\le G(s)

である。ss について下限をとると右端は F(x)F(x) になる。左端については、y<ry<r なら F(y)≤G(r)F(y)\le G(r) である。したがって y<r<xy<r<x を満たすすべての yy で

F(y)≤lim inf⁡jFnj(x)F(y)\le\liminf_jF_{n_j}(x)

である。y↑xy\uparrow x とすると、連続性から次を得る。

F(x)≤lim inf⁡jFnj(x)F(x)\le\liminf_jF_{n_j}(x)

定義 13.3

確率変数の族 {Xn}\{X_n\} が緊密であるとは、次が成り立つことをいう。任意の ε>0\varepsilon >0 に対し、ある M>0M>0 が存在して

P(∣Xn∣>M)<ε(すべての n)P(\left\vert X_n\right\vert >M)<\varepsilon \qquad(\text{すべての }n)

となる。

命題 13.4

{Xn}\{X_n\} が緊密なら、次の2つが成り立つ。

  1. (i) Helly の選出定理の極限 FF は確率分布の分布関数である。つまり F(−∞)=0F(-\infty)=0、F(∞)=1F(\infty)=1 を満たす。
  2. (ii) FF を分布関数にもつ確率変数を XX とする。このとき Xnj→dXX_{n_j}\xrightarrow{d}X となる。
証明

ε>0\varepsilon >0 を与える。緊密性から、sup⁡nP(∣Xn∣>M)<ε\sup_nP(\left\vert X_n\right\vert >M)<\varepsilon となる MM がとれる。FF の不連続点は高々可算個である。したがって x>Mx>M と y<−My<-M を、どちらも連続点にとることができる。このとき

Fnj(x)≥P(∣Xnj∣≤M)>1−ε,Fnj(y)≤P(∣Xnj∣>M)<εF_{n_j}(x)\ge P(\left\vert X_{n_j}\right\vert \le M)>1-\varepsilon ,\qquad F_{n_j}(y)\le P(\left\vert X_{n_j}\right\vert >M)<\varepsilon

である。j→∞j\to\infty として、次を得る。

F(x)≥1−ε,F(y)≤εF(x)\ge1-\varepsilon ,\qquad F(y)\le\varepsilon

これが任意の正の数 ε\varepsilon で成り立つので、次の2式が成り立つ。

F(∞)=1,F(−∞)=0F(\infty)=1,\qquad F(-\infty)=0

第3章 §3.3 から、この性質をもつ非減少で右連続な関数はある確率分布の分布関数である。残りの主張は分布収束の定義そのものである。

次の補題は、特性関数の 0 の近くでの振る舞いから、裾の確率を上から抑える。

補題 13.5(裾の評価)

任意の確率変数 XX と u>0u>0 について

P(∣X∣≥2u)≤1u∫−uu(1−φX(t)) dt(13.6)P\Bigl(\left\vert X\right\vert \ge\frac2u\Bigr)\le\frac1u\int_{-u}^{u}\bigl(1-\varphi_X(t)\bigr)\,\mathrm{d} t\tag{13.6}

が成り立つ(右辺は実数である)。

証明

x≠0x\ne0 について

∫−uu(1−eitx) dt=2u−2sin⁡uxx\int_{-u}^{u}(1-\mathrm{e}^{\mathrm{i} tx})\,\mathrm{d} t=2u-\frac{2\sin ux}{x}

であり、x=0x=0 ではこの値は 0 である。被積分関数は ∣1−eitx∣≤2\left\vert 1-\mathrm{e}^{\mathrm{i} tx}\right\vert \le2 で有界なので、Fubini の定理(第5章)で積分と期待値の順序を入れ替えることができ、

1u∫−uu(1−φX(t)) dt=E[1u∫−uu(1−eitX) dt]=2E[1−sin⁡uXuX]\frac1u\int_{-u}^{u}\bigl(1-\varphi_X(t)\bigr)\,\mathrm{d} t=E\Bigl[\frac1u\int_{-u}^{u}(1-\mathrm{e}^{\mathrm{i} tX})\,\mathrm{d} t\Bigr]=2E\Bigl[1-\frac{\sin uX}{uX}\Bigr]

となる(X=0X=0 では括弧の中を 0 と読む)。∣sin⁡y∣≤∣y∣\left\vert \sin y\right\vert \le\left\vert y\right\vert なので括弧の中は常に 0 以上である。また ∣y∣≥2\left\vert y\right\vert \ge2 では

1−sin⁡yy≥1−1∣y∣≥121-\frac{\sin y}{y}\ge1-\frac1{\left\vert y\right\vert }\ge\frac12

である。したがって

2E[1−sin⁡uXuX]≥2E[(1−sin⁡uXuX)1{∣uX∣≥2}]≥2⋅12P(∣uX∣≥2)=P(∣X∣≥2u)2E\Bigl[1-\frac{\sin uX}{uX}\Bigr]\ge2E\Bigl[\Bigl(1-\frac{\sin uX}{uX}\Bigr)\mathbf 1_{\{\left\vert uX\right\vert \ge2\}}\Bigr]\ge2\cdot\frac12P(\left\vert uX\right\vert \ge2)=P\Bigl(\left\vert X\right\vert \ge\frac2u\Bigr)

である。

連続性定理

証明の最後に使う「部分列の原理」を先に示す。

補題 13.6(部分列の原理)

確率変数の列 XnX_n と確率変数 XX について、次を仮定する。

  • (Xn)(X_n) の任意の部分列が、XX に分布収束するさらに部分列をもつ。

このとき Xn→dXX_n\xrightarrow{d}X である。

証明

gg を有界連続関数とし、次の実数列と実数を考える。

an=E[g(Xn)],a=E[g(X)]a_n=E[g(X_n)],\qquad a=E[g(X)]

an→aa_n\to a でないと仮定する。するとある ε>0\varepsilon >0 と部分列 nkn_k があって、次が成り立つ。

∣ank−a∣≥ε\left\vert a_{n_k}-a\right\vert \ge\varepsilon

仮定から、この部分列のさらに部分列で Xnkj→dXX_{n_{k_j}}\xrightarrow{d}X となる。つまり ankj→aa_{n_{k_j}}\to a であり、これは ∣ankj−a∣≥ε\left\vert a_{n_{k_j}}-a\right\vert \ge\varepsilon と矛盾する。よって an→aa_n\to a であり、(13.5) から結論を得る。

定理 13.7(Lévy の連続性定理)

X1,X2,…X_1,X_2,\dots を確率変数の列とする。

  1. (i) Xn→dXX_n\xrightarrow{d}X ならば、各点で φXn(t)→φX(t)\varphi_{X_n}(t)\to\varphi_X(t) である。

  2. (ii) 次の2つを仮定する。

    • 各点 tt で φXn(t)\varphi_{X_n}(t) が収束する。その極限を φ(t)\varphi(t) と書く。
    • 極限 φ\varphi が t=0t=0 で連続である。

このとき φ\varphi はある確率変数 XX の特性関数である。さらに Xn→dXX_n\xrightarrow{d}X である。

証明

(i) は命題 13.1 である。(ii) を示す。

第1段(緊密性)。φXn(0)=1\varphi_{X_n}(0)=1 なので φ(0)=1\varphi(0)=1 である。ε>0\varepsilon >0 を与える。

φ\varphi は 0 で連続である。したがって、ある正の数 uu がとれて、区間 [−u,u][-u,u] で次が成り立つ。

∣1−φ(t)∣<ε2\left\vert 1-\varphi(t)\right\vert <\frac{\varepsilon }2

このとき

∣1u∫−uu(1−φ(t)) dt∣≤1u⋅2u⋅ε2=ε\Bigl|\frac1u\int_{-u}^{u}\bigl(1-\varphi(t)\bigr)\,\mathrm{d} t\Bigr|\le\frac1u\cdot2u\cdot\frac{\varepsilon }2=\varepsilon

である。

被積分関数は ∣1−φXn(t)∣≤2\left\vert 1-\varphi_{X_n}(t)\right\vert \le2 と有界で、[−u,u][-u,u] の上で各点収束している。有界収束定理(第5章 §5.2)から

1u∫−uu(1−φXn) dt→1u∫−uu(1−φ) dt\frac1u\int_{-u}^{u}(1-\varphi_{X_n})\,\mathrm{d} t\to\frac1u\int_{-u}^{u}(1-\varphi)\,\mathrm{d} t

である。したがって、ある NN から先のすべての番号で左辺は 2ε2\varepsilon より小さくなる。補題 13.5 から、その番号では P(∣Xn∣≥2/u)<2εP(\left\vert X_n\right\vert \ge2/u)<2\varepsilon である。

残りの有限個 X1,…,XN−1X_1,\dots,X_{N-1} については、MM を十分大きくとれば、どれも裾の確率が上の値より小さくなる。さらに M≥2/uM\ge2/u にとれば、次がすべての nn で成り立つ。

P(∣Xn∣>M)<2εP(\left\vert X_n\right\vert >M)<2\varepsilon

第2段(極限の特定)。(Xn)(X_n) の任意の部分列をとる。緊密な族の部分列も緊密である。定理 13.2 と命題 13.4 から、そのさらに部分列 XnjX_{n_j} で、ある確率変数 YY に分布収束するものがある。

(i) から φXnj(t)→φY(t)\varphi_{X_{n_j}}(t)\to\varphi_Y(t) である。仮定から φXnj(t)→φ(t)\varphi_{X_{n_j}}(t)\to\varphi(t) でもあるので、φY=φ\varphi_Y=\varphi である。とくに φ\varphi は特性関数であり、その分布を1つ選んで XX と書く。

一意性定理(第9章 §9.2)から Y=dXY\overset{d}{=}X である。どの部分列から出発しても極限の分布は同じ XX の分布なので、補題 13.6 から Xn→dXX_n\xrightarrow{d}X となる。

注意

  • 定理の (ii) では、極限 φ\varphi が特性関数であることを仮定する必要がない。0 での連続性だけで、緊密性を通じて極限の分布の存在まで保証される。
  • 極限 φ\varphi が初めから特性関数だと分かっているとき(たとえば e−t2/2\mathrm{e}^{-t^2/2})は、連続性の仮定は自動的に満たされる。中心極限定理ではこの形で使う。
  • (i) の収束は、実は tt の有界区間の上で一様である。この節では各点収束だけを使うので、証明は省く。

複素数の積の評価

中心極限定理の証明では、nn 個の複素数の積の極限を扱う。そのための2つの補題を用意する。

補題 13.8

θ≥1\theta\ge1 とする。複素数 z1,…,znz_1,\dots,z_n と w1,…,wnw_1,\dots,w_n の絶対値は、すべてこの値以下であるとする。このとき

∣∏k=1nzk−∏k=1nwk∣≤θn−1∑k=1n∣zk−wk∣(13.7)\Bigl|\prod_{k=1}^nz_k-\prod_{k=1}^nw_k\Bigr|\le\theta^{n-1}\sum_{k=1}^n\left\vert z_k-w_k\right\vert\tag{13.7}

が成り立つ。

証明

nn についての帰納法による。n=1n=1 では等号である。n−1n-1 個で成り立つとすると、

∏k=1nzk−∏k=1nwk=z1(∏k=2nzk−∏k=2nwk)+(z1−w1)∏k=2nwk\prod_{k=1}^nz_k-\prod_{k=1}^nw_k=z_1\Bigl(\prod_{k=2}^nz_k-\prod_{k=2}^nw_k\Bigr)+(z_1-w_1)\prod_{k=2}^nw_k

であり、右辺の絶対値は

θ⋅θn−2∑k=2n∣zk−wk∣+∣z1−w1∣θn−1=θn−1∑k=1n∣zk−wk∣\theta\cdot\theta^{n-2}\sum_{k=2}^n\left\vert z_k-w_k\right\vert +\left\vert z_1-w_1\right\vert \theta^{n-1}=\theta^{n-1}\sum_{k=1}^n\left\vert z_k-w_k\right\vert

以下である。

補題 13.9

複素数の列 cnc_n が cc に収束するならば、次が成り立つ。

(1+cnn)n→ec.\Bigl(1+\frac{c_n}{n}\Bigr)^n\to\mathrm{e}^c.

証明

複素数 ww について

∣ew−1−w∣=∣∑k≥2wkk!∣≤∣w∣2∑j≥0∣w∣j(j+2)!≤∣w∣2e∣w∣\left\vert \mathrm{e}^w-1-w\right\vert =\Bigl|\sum_{k\ge2}\frac{w^k}{k!}\Bigr|\le\left\vert w\right\vert ^2\sum_{j\ge0}\frac{\left\vert w\right\vert ^j}{(j+2)!}\le\left\vert w\right\vert ^2\mathrm{e}^{\left\vert w\right\vert }

である。z=1+cn/nz=1+c_n/n、w=ecn/nw=\mathrm{e}^{c_n/n} と置く。どちらも絶対値は θ=e∣cn∣/n\theta=\mathrm{e}^{\left\vert c_n\right\vert /n} 以下である。実際

∣1+cn/n∣≤1+∣cn∣/n≤e∣cn∣/n\left\vert 1+c_n/n\right\vert \le1+\left\vert c_n\right\vert /n\le\mathrm{e}^{\left\vert c_n\right\vert /n}

である。補題 13.8 で zk=zz_k=z、wk=ww_k=w とすると

∣(1+cnn)n−ecn∣≤e∣cn∣(n−1)/n⋅n⋅∣cn∣2n2e∣cn∣/n≤∣cn∣2e∣cn∣n\Bigl|\Bigl(1+\frac{c_n}{n}\Bigr)^n-\mathrm{e}^{c_n}\Bigr|\le\mathrm{e}^{\left\vert c_n\right\vert (n-1)/n}\cdot n\cdot\frac{\left\vert c_n\right\vert ^2}{n^2}\mathrm{e}^{\left\vert c_n\right\vert /n}\le\frac{\left\vert c_n\right\vert ^2\mathrm{e}^{\left\vert c_n\right\vert }}{n}

である。cnc_n は有界なので右辺は 0 に収束し、ecn→ec\mathrm{e}^{c_n}\to\mathrm{e}^c と合わせて結論を得る。

特性関数による弱法則

第12章 §12.1 の Khinchin の弱法則は、特性関数でも証明できる。X1,X2,…X_1,X_2,\dots を独立同分布とし、次を仮定する。

E[∣X1∣]<∞,E[X1]=μE[\left\vert X_1\right\vert ]<\infty,\qquad E[X_1]=\mu

第8章 §8.3 から、このとき φ=φX1\varphi=\varphi_{X_1} は次のように展開できる。

φ(t)=1+iμt+o(t)(t→0)\varphi(t)=1+\mathrm{i}\mu t+o(t)\qquad(t\to0)

Xˉn\bar X_n の特性関数は

φXˉn(t)=φ(tn)n=(1+cnn)n,cn=n(φ(tn)−1)=iμt+n⋅o(1n)→iμt\varphi_{\bar X_n}(t)=\varphi\Bigl(\frac tn\Bigr)^n=\Bigl(1+\frac{c_n}{n}\Bigr)^n,\qquad c_n=n\Bigl(\varphi\Bigl(\frac tn\Bigr)-1\Bigr)=\mathrm{i}\mu t+n\cdot o\Bigl(\frac1n\Bigr)\to\mathrm{i}\mu t

なので、補題 13.9 から φXˉn(t)→eiμt\varphi_{\bar X_n}(t)\to\mathrm{e}^{\mathrm{i}\mu t} である。これは定数 μ\mu の特性関数なので、定理 13.7 から Xˉn→dμ\bar X_n\xrightarrow{d}\mu である。極限が定数なので Xˉn→Pμ\bar X_n\xrightarrow{P}\mu となる(第11章 §11.2)。

例題 13.1 二項分布の Poisson 極限を特性関数で示す ★★

オリジナルλ>0\lambda>0 を定数とする。n>λn>\lambda について Xn∼Bin(n,λ/n)X_n\sim\mathrm{Bin}(n,\lambda/n) とする。

  1. (1) φXn(t)\varphi_{X_n}(t) を求め、n→∞n\to\infty の極限を計算せよ。
  2. (2) Xn→dPo(λ)X_n\xrightarrow{d}\mathrm{Po}(\lambda) を示せ。また各整数 k≥0k\ge0 で次を導け。

P(Xn=k)→e−λλkk!P(X_n=k)\to\mathrm{e}^{-\lambda}\frac{\lambda^k}{k!} 3. (3) λ=2\lambda=2 とする。n=50n=50 のとき、P(Xn=0)P(X_n=0) を極限の値と比べよ。

解答

(1) p=λ/np=\lambda/n と置く。XnX_n は、成功確率 pp の独立な Bernoulli 確率変数を試行の回数だけ足したものである。1個の特性関数は 1+p(eit−1)1+p(\mathrm{e}^{\mathrm{i} t}-1) である。独立和の特性関数は積なので

φXn(t)=(1+λ(eit−1)n)n\varphi_{X_n}(t)=\Bigl(1+\frac{\lambda(\mathrm{e}^{\mathrm{i} t}-1)}{n}\Bigr)^n

である。cn=c=λ(eit−1)c_n=c=\lambda(\mathrm{e}^{\mathrm{i} t}-1)(nn によらない)として補題 13.9 を使うと

φXn(t)→exp⁡(λ(eit−1))\varphi_{X_n}(t)\to\exp\bigl(\lambda(\mathrm{e}^{\mathrm{i} t}-1)\bigr)

を得る。

(2) 極限は Po(λ)\mathrm{Po}(\lambda) の特性関数(第8章の表)であり、tt について連続である。定理 13.7 (ii) から、Y∼Po(λ)Y\sim\mathrm{Po}(\lambda) について Xn→dYX_n\xrightarrow{d}Y である。YY の分布関数は整数でだけ跳ぶので、半整数 k±12k\pm\frac12 はその連続点である。したがって

P(Xn=k)=FXn(k+12)−FXn(k−12)→FY(k+12)−FY(k−12)=P(Y=k)=e−λλkk!P(X_n=k)=F_{X_n}\Bigl(k+\frac12\Bigr)-F_{X_n}\Bigl(k-\frac12\Bigr)\to F_Y\Bigl(k+\frac12\Bigr)-F_Y\Bigl(k-\frac12\Bigr)=P(Y=k)=\mathrm{e}^{-\lambda}\frac{\lambda^k}{k!}

となる。

(3) 有限の nn での値と極限の値は

P(X50=0)=(1−250)50=0.9650≈0.12989,e−2≈0.13534P(X_{50}=0)=\Bigl(1-\frac2{50}\Bigr)^{50}=0.96^{50}\approx0.12989,\qquad\mathrm{e}^{-2}\approx0.13534

である。差は約 0.0054 で、n=50n=50 でも極限の値に近い。

確認テスト

オリジナル

  1. 13.1-1 Xn∼U(−n,n)X_n\sim\mathrm U(-n,n) について、XnX_n が分布収束しないことを、第11章の「分布収束する列は緊密である」ことを使って示せ。また、定理 13.7 (ii) のどの仮定が満たされないかを答えよ。
  2. 13.1-2 Xn∼Ga(n,n)X_n\sim\mathrm{Ga}(n,n) とする(形もレートも nn)。特性関数の極限を計算し、Xn→P1X_n\xrightarrow{P}1 を示せ。
  3. 13.1-3 0<λ<n0<\lambda<n とする。XnX_n を、成功確率 pn=λ/np_n=\lambda/n の幾何分布に従う確率変数とする。つまり

P(Xn=k)=pn(1−pn)k−1(k=1,2,… )P(X_n=k)=p_n(1-p_n)^{k-1}\qquad(k=1,2,\dots)

である。Xn/n→dEx(λ)X_n/n\xrightarrow{d}\mathrm{Ex}(\lambda) を特性関数で示せ。4. 13.1-4 各番号で XnX_n と YnY_n は独立とする。また Xn→dXX_n\xrightarrow{d}X、Yn→dYY_n\xrightarrow{d}Y とする。独立な確率変数の組 X′,Y′X',Y' で、次を満たすものをとる。

X′=dX,Y′=dYX'\overset d=X,\qquad Y'\overset d=Y

このとき Xn+Yn→dX′+Y′X_n+Y_n\xrightarrow{d}X'+Y' を示せ。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『基礎理論』の記事一覧へ