tecktom.com

Slutsky の定理と連続写像定理

2026年10月15日28分で読めますAI生成教科書『基礎理論』第11章 確率変数列の収束
目次
  1. Skorokhod の表現
  2. 有界連続関数による言い換え
  3. 連続写像定理
  4. Slutsky の定理
  5. モーメントの収束
  6. 確認テスト
  7. 木下の層

AI 生成・木下が未勉強

確率論入門 第11章 §11.4 ★★★ 分布収束を有界連続関数の期待値の収束で言い換え、連続関数を通しても、確率収束する定数の因子を掛けても、分布収束が保たれることを示す。中心極限定理から使える近似を作るときの道具

第13章の中心極限定理は、独立同分布で平均 μ\mu、分散 σ2\sigma^2 の和について、次を主張する。

Zn=n(Xˉn−μ)σ→dN(0,1)Z_n=\frac{\sqrt n(\bar X_n-\mu)}\sigma\xrightarrow{d}\mathrm N(0,1)

しかし実際のデータでは標準偏差 σ\sigma が分からないことが多い。そこで、それを推定値に置き換えた量の分布が知りたくなる。

また Zn2Z_n^2 や eZn\mathrm{e}^{Z_n} のように、収束する量の関数の分布も知りたくなる。この節の2つの定理は、こうした操作で分布収束が保たれることを保証する。証明の道具として、まず分布収束する列を概収束する列で表す Skorokhod の表現を示す。

Skorokhod の表現

分布収束は別々の確率空間の上の確率変数についても意味を持つが、同じ空間の上で値として近づいている方が扱いやすい。第3章 §3.5 の逆変換法を使うと、分布を変えずに値を近づけることができる。

定理 11.16(Skorokhod の表現(1次元))

Xn→dXX_n\xrightarrow{d}X とする。確率空間 ((0,1),B((0,1)),Leb)((0,1),\mathcal B((0,1)),\mathrm{Leb}) の上に、次の3つを満たす確率変数の列と極限が存在する。

Yn=dXn,Y=dX,Yn→Y a.s.Y_n\overset{d}{=}X_n,\qquad Y\overset{d}{=}X,\qquad Y_n\to Y\ \text{a.s.}

具体的には、分布関数を Fn=FXnF_n=F_{X_n}、F=FXF=F_X と書き、その一般化逆関数を使って次のように置けばよい。

Yn(u)=Fn−1(u),Y(u)=F−1(u)Y_n(u)=F_n^{-1}(u),\qquad Y(u)=F^{-1}(u)

証明

区間 (0,1)(0,1) の上の恒等写像 U(u)=uU(u)=u は、一様分布に従う。逆変換法から、Yn=Fn−1(U)Y_n=F_n^{-1}(U) の分布関数は FnF_n である。同様に Y=F−1(U)Y=F^{-1}(U) の分布関数は FF である。示すべきことは、F−1F^{-1} が連続な点 uu のすべてで、次が成り立つことである。

Fn−1(u)→F−1(u)F_n^{-1}(u)\to F^{-1}(u)

F−1F^{-1} は非減少関数なので、不連続点は高々可算個である(各不連続点で跳び上がる区間に有理数が1つずつとれる)。可算集合の Lebesgue 測度は 0 だから、これで Yn→YY_n\to Y a.s. を結論することができる。以下、F−1F^{-1} が点 u∈(0,1)u\in(0,1) で連続であるとする。その値を xx と置く。

下からの評価:y<xy<x を FF の連続点とする。

  1. F−1(u)≤yF^{-1}(u)\le y ではない。よって (11.3) から、u≤F(y)u\le F(y) ではなく F(y)<uF(y)<u である。
  2. Fn(y)→F(y)F_n(y)\to F(y) なので、十分大きな番号で Fn(y)<uF_n(y)<u となる。再び (11.3) から Fn−1(u)>yF_n^{-1}(u)>y である。
  3. よって lim inf⁡nFn−1(u)≥y\liminf_nF_n^{-1}(u)\ge y である。FF の連続点は稠密なので、yy をこの点にいくらでも近くとることができる。したがって次を得る。

lim inf⁡nFn−1(u)≥x\liminf_nF_n^{-1}(u)\ge x

上からの評価:u<u′<1u<u'<1 をとる。y>F−1(u′)y>F^{-1}(u') を FF の連続点とする。

  1. (11.3) から u′≤F(y)u'\le F(y) である。よって F(y)>uF(y)>u である。
  2. 十分大きな番号で Fn(y)>uF_n(y)>u となる。(11.3) から Fn−1(u)≤yF_n^{-1}(u)\le y である。
  3. よって lim sup⁡nFn−1(u)≤y\limsup_nF_n^{-1}(u)\le y である。y↓F−1(u′)y\downarrow F^{-1}(u') として次を得る。

lim sup⁡nFn−1(u)≤F−1(u′)\limsup_nF_n^{-1}(u)\le F^{-1}(u') 4. F−1F^{-1} はこの点で連続なので、u′↓uu'\downarrow u として次となる。

lim sup⁡nFn−1(u)≤F−1(u)=x\limsup_nF_n^{-1}(u)\le F^{-1}(u)=x

例題 11.1 の Yn=n(1−Mn)Y_n=n(1-M_n) の表現

分布関数は次のとおりである。

Fn(x)=1−(1−xn)n(0≤x≤n)F_n(x)=1-\Bigl(1-\frac xn\Bigr)^n\qquad(0\le x\le n)

これは連続で狭義に増加するので、一般化逆関数は普通の逆関数である。u=Fn(x)u=F_n(x) を xx について解いて

Fn−1(u)=n(1−(1−u)1/n)F_n^{-1}(u)=n\bigl(1-(1-u)^{1/n}\bigr)

を得る。展開

(1−u)1/n=exp⁡(1nlog⁡(1−u))=1+1nlog⁡(1−u)+O(n−2)(1-u)^{1/n}=\exp\Bigl(\frac1n\log(1-u)\Bigr)=1+\frac1n\log(1-u)+O(n^{-2})

から、各 u∈(0,1)u\in(0,1) で Fn−1(u)→−log⁡(1−u)F_n^{-1}(u)\to-\log(1-u) である。これは Ex(1)\mathrm{Ex}(1) の分布関数 1−e−x1-\mathrm{e}^{-x} の逆関数である。たとえば u=12u=\frac12 では次のようになる。

F10−1(12)=10(1−2−1/10)≈0.6697,極限は log⁡2≈0.6931F_{10}^{-1}\Bigl(\frac12\Bigr)=10(1-2^{-1/10})\approx0.6697,\qquad\text{極限は}\ \log2\approx0.6931

有界連続関数による言い換え

定理 11.17(分布収束の言い換え)

Xn→dXX_n\xrightarrow{d}X であることと、すべての有界連続関数 g:R→Rg:\mathbb{R}\to\mathbb{R} について

E[g(Xn)]→E[g(X)](11.17)E[g(X_n)]\to E[g(X)]\tag{11.17}

が成り立つことは同値である。

証明

⇒\Rightarrow:定理 11.16 の表現 Yn,YY_n,Y をとる。期待値 E[g(Xn)]E[g(X_n)] は分布だけで決まる(第5章 §5.3 の変数変換公式)。よって

E[g(Xn)]=E[g(Yn)],E[g(X)]=E[g(Y)]E[g(X_n)]=E[g(Y_n)],\qquad E[g(X)]=E[g(Y)]

である。gg は連続なので、g(Yn)→g(Y)g(Y_n)\to g(Y) a.s. である。さらに ∣g∣≤C\left\vert g\right\vert \le C なので、有界収束定理(第5章 §5.2)から E[g(Yn)]→E[g(Y)]E[g(Y_n)]\to E[g(Y)] である。

⇐\Leftarrow:点 xx を極限の分布関数の連続点とする。正の数 ε\varepsilon を1つとる。関数 gg を次のように定める。

  • (−∞,x](-\infty,x] で 1 とする。
  • [x+ε,∞)[x+\varepsilon ,\infty) で 0 とする。
  • そのあいだでは直線的に減少させる。

すると gg は有界連続で

1(−∞,x]≤g≤1(−∞,x+ε]\mathbf 1_{(-\infty,x]}\le g\le\mathbf 1_{(-\infty,x+\varepsilon ]}

を満たす。左の不等式から FXn(x)≤E[g(Xn)]F_{X_n}(x)\le E[g(X_n)] である。右の不等式から E[g(X)]≤FX(x+ε)E[g(X)]\le F_X(x+\varepsilon ) である。よって

lim sup⁡nFXn(x)≤lim⁡nE[g(Xn)]=E[g(X)]≤FX(x+ε)\limsup_nF_{X_n}(x)\le\lim_nE[g(X_n)]=E[g(X)]\le F_X(x+\varepsilon )

である。次に、gg を ε\varepsilon だけ左にずらした関数を使う。これは (−∞,x−ε](-\infty,x-\varepsilon ] で 1、[x,∞)[x,\infty) で 0 となる関数である。同様に次を得る。

lim inf⁡nFXn(x)≥FX(x−ε)\liminf_nF_{X_n}(x)\ge F_X(x-\varepsilon )

ε↓0\varepsilon \downarrow0 として、この点での連続性から FXn(x)→FX(x)F_{X_n}(x)\to F_X(x) となる。

注意

  • 多次元への拡張。(11.17) は分布関数を使わないので、Rd\mathbb{R}^d 値の確率ベクトルにもそのまま意味を持つ。Rd\mathbb{R}^d 値の列については、すべての有界連続関数 g:Rd→Rg:\mathbb{R}^d\to\mathbb{R} で次が成り立つことを、分布収束の定義とする。

E[g(Xn)]→E[g(X)]E[g(\bm X_n)]\to E[g(\bm X)]

第13章 §13.3 の多変量の中心極限定理は、この意味の分布収束を主張する。

  • なめらかな関数に限ってよい。定理の「⇐\Leftarrow」の証明が示すように、gg はなめらかな関数に限っても同値性は保たれる。たとえば何回でも微分可能で、導関数がすべて有界なものでよい。証明の関数の傾きの部分を、滑らかな曲線に置き換えればよいからである。
  • 有界性は外せない。有界でない連続関数 g(x)=xg(x)=x については、(11.17) は一般に成り立たない(確認テスト 11.4-5)。モーメントの収束には §11.3 の一様可積分性が要る(命題 11.22)。

連続写像定理

定理 11.18(連続写像定理)

h:R→Rh:\mathbb{R}\to\mathbb{R} を Borel 可測関数とし、それが連続でない点の集合を DhD_h とする(これが Borel 集合になることは示せる。ここでは認める)。さらに P(X∈Dh)=0P(X\in D_h)=0 とする。このとき次が成り立つ。

  1. (i) 概収束は保たれる。

Xn→X a.s. ⟹ h(Xn)→h(X) a.s.X_n\to X\ \text{a.s.}\ \Longrightarrow\ h(X_n)\to h(X)\ \text{a.s.} 2. (ii) 確率収束は保たれる。

Xn→PX ⟹ h(Xn)→Ph(X)X_n\xrightarrow{P}X\ \Longrightarrow\ h(X_n)\xrightarrow{P}h(X) 3. (iii) 分布収束は保たれる。

Xn→dX ⟹ h(Xn)→dh(X)X_n\xrightarrow{d}X\ \Longrightarrow\ h(X_n)\xrightarrow{d}h(X)

証明

(i):事象 {Xn→X}∩{X∉Dh}\{X_n\to X\}\cap\{X\notin D_h\} は確率 1 である。その上では hh が X(ω)X(\omega) で連続なので、次が成り立つ。

h(Xn(ω))→h(X(ω))h(X_n(\omega))\to h(X(\omega))

(ii):系 11.9 を使う。任意の部分列 (Xnk)(X_{n_k}) について、確率収束から、概収束するさらに部分列 XnkjX_{n_{k_j}} がとれる。

(i) から h(Xnkj)→h(X)h(X_{n_{k_j}})\to h(X) a.s. である。こうして h(Xn)h(X_n) の任意の部分列が、概収束する部分列を持つ。系 11.9 から h(Xn)→Ph(X)h(X_n)\xrightarrow{P}h(X) である。

(iii):定理 11.16 の表現 YnY_n、YY をとる。

  1. 極限は同分布なので、P(Y∈Dh)=P(X∈Dh)=0P(Y\in D_h)=P(X\in D_h)=0 である。
  2. (i) から h(Yn)→h(Y)h(Y_n)\to h(Y) a.s. である。定理 11.7 (i)(iv) から、次が成り立つ。

h(Yn)→dh(Y)h(Y_n)\xrightarrow{d}h(Y) 3. 同分布の確率変数を同じ Borel 関数に通すと、同分布になる。よって

h(Yn)=dh(Xn),h(Y)=dh(X)h(Y_n)\overset{d}{=}h(X_n),\qquad h(Y)\overset{d}{=}h(X)

であり、h(Xn)→dh(X)h(X_n)\xrightarrow{d}h(X) を得る。

連続写像定理の使い方

  • 2乗。Zn→dN(0,1)Z_n\xrightarrow{d}\mathrm N(0,1) とし、Z∼N(0,1)Z\sim\mathrm N(0,1) とする。関数 h(x)=x2h(x)=x^2 は連続なので、次が成り立つ。

Zn2→dZ2Z_n^2\xrightarrow{d}Z^2

Z2Z^2 の分布は、第7章 §7.1 で求めた自由度 1 のカイ2乗分布、すなわち Ga(12,12)\mathrm{Ga}(\frac12,\frac12) である。

  • 逆数。0 でない点では h(x)=1/xh(x)=1/x とし、h(0)=0h(0)=0 と定める。この関数の不連続点は原点だけである。よって Xn→dXX_n\xrightarrow{d}X で P(X=0)=0P(X=0)=0 なら、次が成り立つ。

1Xn→d1X\frac1{X_n}\xrightarrow{d}\frac1Xたとえば X∼N(0,1)X\sim\mathrm N(0,1) ならこの条件が成り立つ。

落とし穴

極限が hh の不連続点に確率をもつと、結論は成り立たない。次の設定で確かめる。

Xn=1n,X=0,h=1(0,∞)X_n=\frac1n,\qquad X=0,\qquad h=\mathbf 1_{(0,\infty)}

列はすべての意味で極限に収束する。ところが、すべての番号で h(Xn)=1h(X_n)=1 であり、h(X)=0h(X)=0 に収束しない。ここでは不連続点の集合が Dh={0}D_h=\{0\} で、P(X∈Dh)=1P(X\in D_h)=1 である。

Slutsky の定理

まず、分布収束する列は「確率的に有界」であることを示す。

補題 11.19(分布収束する列の緊密性)

Xn→dXX_n\xrightarrow{d}X とする。任意の ε>0\varepsilon >0 について、ある M>0M>0 が存在して、すべての番号で次が成り立つ。

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

証明

極限の分布関数 F=FXF=F_X は、左の端で 0、右の端で 1 に近づく。また連続点は稠密である。よって、2点 ±M0\pm M_0 がともに連続点で、次を満たす M0M_0 がとれる。

F(−M0)<ε4,1−F(M0)<ε4F(-M_0)<\frac{\varepsilon }4,\qquad 1-F(M_0)<\frac{\varepsilon }4

FXn(±M0)→F(±M0)F_{X_n}(\pm M_0)\to F(\pm M_0) なので、ある NN から先の番号で次が成り立つ。

FXn(−M0)<ε2,1−FXn(M0)<ε2F_{X_n}(-M_0)<\frac{\varepsilon }2,\qquad 1-F_{X_n}(M_0)<\frac{\varepsilon }2

よって

P(∣Xn∣>M0)≤FXn(−M0)+1−FXn(M0)<εP(\left\vert X_n\right\vert >M_0)\le F_{X_n}(-M_0)+1-F_{X_n}(M_0)<\varepsilon

である。残りの有限個 X1,…,XN−1X_1,\dots,X_{N-1} を考える。各項は実数値なので、M→∞M\to\infty で P(∣Xn∣>M)→0P(\left\vert X_n\right\vert >M)\to0 である。よって M≥M0M\ge M_0 を大きくとれば、すべての番号で P(∣Xn∣>M)<εP(\left\vert X_n\right\vert >M)<\varepsilon となる。

補題 11.20

  1. (i) XnX_n が補題 11.19 の結論を満たすとする。Zn→P0Z_n\xrightarrow{P}0 ならば、XnZn→P0X_nZ_n\xrightarrow{P}0 である。
  2. (ii) Xn→dXX_n\xrightarrow{d}X とする。Zn→P0Z_n\xrightarrow{P}0 ならば、Xn+Zn→dXX_n+Z_n\xrightarrow{d}X である。
証明

(i):ε,η>0\varepsilon ,\eta>0 を与える。sup⁡nP(∣Xn∣>M)<η\sup_nP(\left\vert X_n\right\vert >M)<\eta となる MM をとる。∣XnZn∣>ε\left\vert X_nZ_n\right\vert >\varepsilon ならば、次のどちらかが成り立つ。

∣Xn∣>Mまたは∣Zn∣>εM\left\vert X_n\right\vert >M\quad\text{または}\quad\left\vert Z_n\right\vert >\frac{\varepsilon } M

よって

P(∣XnZn∣>ε)≤η+P(∣Zn∣>ε/M)P(\left\vert X_nZ_n\right\vert >\varepsilon )\le\eta+P(\left\vert Z_n\right\vert >\varepsilon /M)

である。番号を大きくすると第2項は 0 に近づく。よって lim sup⁡nP(∣XnZn∣>ε)≤η\limsup_nP(\left\vert X_nZ_n\right\vert >\varepsilon )\le\eta であり、η\eta は任意である。

(ii):定理 11.7 (iv) の証明と同じ議論を使う。点 xx を FXF_X の連続点とする。さらに、2点 x±εx\pm\varepsilon も連続点となる ε>0\varepsilon >0 をとる。

不連続点は高々可算なので、そうした ε\varepsilon はいくらでも小さくとれる。次の2つが成り立つ。

  • Xn+Zn≤xX_n+Z_n\le x ならば、次のどちらかである。

Xn≤x+εまたは∣Zn∣>εX_n\le x+\varepsilon \quad\text{または}\quad\left\vert Z_n\right\vert >\varepsilon

  • Xn≤x−εX_n\le x-\varepsilon ならば、次のどちらかである。

Xn+Zn≤xまたは∣Zn∣>εX_n+Z_n\le x\quad\text{または}\quad\left\vert Z_n\right\vert >\varepsilon

よって

FXn(x−ε)−P(∣Zn∣>ε)≤P(Xn+Zn≤x)≤FXn(x+ε)+P(∣Zn∣>ε)F_{X_n}(x-\varepsilon )-P(\left\vert Z_n\right\vert >\varepsilon )\le P(X_n+Z_n\le x)\le F_{X_n}(x+\varepsilon )+P(\left\vert Z_n\right\vert >\varepsilon )

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

FX(x−ε)≤lim inf⁡nP(Xn+Zn≤x)≤lim sup⁡nP(Xn+Zn≤x)≤FX(x+ε)F_X(x-\varepsilon )\le\liminf_nP(X_n+Z_n\le x)\le\limsup_nP(X_n+Z_n\le x)\le F_X(x+\varepsilon )

あとは ε↓0\varepsilon \downarrow0 とすればよい。

定理 11.21(Slutsky の定理)

Xn→dXX_n\xrightarrow{d}X とする。定数 cc について Yn→PcY_n\xrightarrow{P}c とする(2つの列は同じ確率空間の上にあるとする)。このとき

  1. (i) Xn+Yn→dX+cX_n+Y_n\xrightarrow{d}X+c、
  2. (ii) YnXn→dcXY_nX_n\xrightarrow{d}cX、
  3. (iii) c≠0c\ne0 なら Xn/Yn→dX/cX_n/Y_n\xrightarrow{d}X/c。ただし Yn=0Y_n=0 となる標本点では、比を 0 などと任意に定める。
証明

(i):分布関数は FXn+c(x)=FXn(x−c)F_{X_n+c}(x)=F_{X_n}(x-c) を満たす。ずらした分布関数の連続点 xx は、もとの分布関数の連続点 x−cx-c にちょうど対応する。よって Xn+c→dX+cX_n+c\xrightarrow{d}X+c である。さらに次のように分ける。

Xn+Yn=(Xn+c)+(Yn−c),Yn−c→P0X_n+Y_n=(X_n+c)+(Y_n-c),\qquad Y_n-c\xrightarrow{P}0

補題 11.20 (ii) から結論を得る。

(ii):YnXn=cXn+(Yn−c)XnY_nX_n=cX_n+(Y_n-c)X_n と分ける。h(x)=cxh(x)=cx は連続なので、連続写像定理 (iii) から cXn→dcXcX_n\xrightarrow{d}cX である。補題 11.19 から、XnX_n は確率的に有界である。Yn−c→P0Y_n-c\xrightarrow{P}0 なので、補題 11.20 (i) から次が成り立つ。

(Yn−c)Xn→P0(Y_n-c)X_n\xrightarrow{P}0

補題 11.20 (ii) から YnXn→dcXY_nX_n\xrightarrow{d}cX となる。

(iii):0 でない点では h(y)=1/yh(y)=1/y とし、h(0)=0h(0)=0 と定める。この関数は点 c≠0c\ne0 で連続である。連続写像定理 (ii) を定数への確率収束に使って、1/Yn→P1/c1/Y_n\xrightarrow{P}1/c を得る。(ii) を 1/Yn1/Y_n に使えば、Xn/Yn→dX/cX_n/Y_n\xrightarrow{d}X/c となる。

落とし穴

YnY_n の極限が定数でなければ、Slutsky の定理は成り立たない。確認テスト 11.2-1 の列で確かめる。標準正規分布に従う確率変数を1つとり、それ自身の列と、符号を反転した列を並べた例である。

どちらの列も XX に分布収束する。しかし和 Xn+Yn=0X_n+Y_n=0 は 2X2X に分布収束しない。2つの列の和の分布は、それぞれの分布だけでは決まらず、同時分布によって決まるからである。極限が定数なら同時分布の問題が起きないことが、定理の仮定の意味である。

モーメントの収束

分布収束から期待値の収束を引き出すには、§11.3 の一様可積分性を Skorokhod の表現と組み合わせる。

命題 11.22(分布収束とモーメントの収束)

Xn→dXX_n\xrightarrow{d}X とし、ある r>0r>0 について次を仮定する。

sup⁡nE[∣Xn∣r]<∞\sup_nE[\left\vert X_n\right\vert ^r]<\infty

このとき 0<p<r0<p<r について E[∣X∣p]<∞E[\left\vert X\right\vert ^p]<\infty で、

E[∣Xn∣p]→E[∣X∣p]E[\left\vert X_n\right\vert ^p]\to E[\left\vert X\right\vert ^p]

である。さらに k<rk<r を満たす正の整数について、その次数のモーメントも収束する。

証明

定理 11.16 の表現 Yn→YY_n\to Y a.s. をとり、Wn=∣Yn∣pW_n=\left\vert Y_n\right\vert ^p と置く。

  1. 一様可積分である。Wn→∣Y∣pW_n\to\left\vert Y\right\vert ^p a.s. である。また E[Wnr/p]=E[∣Xn∣r]E[W_n^{r/p}]=E[\left\vert X_n\right\vert ^r] は番号について有界である。r/p>1r/p>1 なので、命題 11.14 (iii) から {Wn}\{W_n\} は一様可積分である。
  2. 前半。Vitali の定理(概収束は確率収束を導く)から、E∣Y∣p<∞E\left\vert Y\right\vert ^p<\infty かつ E[Wn]→E[∣Y∣p]E[W_n]\to E[\left\vert Y\right\vert ^p] である。同分布から E[Wn]=E[∣Xn∣p]E[W_n]=E[\left\vert X_n\right\vert ^p]、E[∣Y∣p]=E[∣X∣p]E[\left\vert Y\right\vert ^p]=E[\left\vert X\right\vert ^p] なので、前半を得る。
  3. 後半。Ynk→YkY_n^k\to Y^k a.s. である。∣Ynk∣=∣Yn∣k\left\vert Y_n^k\right\vert =\left\vert Y_n\right\vert ^k は同じ理由で一様可積分である。よって Vitali の定理から E[Ynk]→E[Yk]E[Y_n^k]\to E[Y^k] となる。

例題 11.1 の Yn=n(1−Mn)Y_n=n(1-M_n) のモーメント

最大値の1からのずれの密度は、区間 (0,1) の上で次のとおりである。

n(1−t)n−1n(1-t)^{n-1}

よって

E[(1−Mn)k]=n∫01tk(1−t)n−1 dt=k! n!(n+k)!E[(1-M_n)^k]=n\int_0^1t^k(1-t)^{n-1}\,\mathrm{d} t=\frac{k!\,n!}{(n+k)!}

である。したがって

E[Yn]=nn+1→1,E[Yn2]=2n2(n+1)(n+2)→2E[Y_n]=\frac{n}{n+1}\to1,\qquad E[Y_n^2]=\frac{2n^2}{(n+1)(n+2)}\to2

である。これは極限 Ex(1)\mathrm{Ex}(1) の1次と2次のモーメント(1 と 2)に一致する。命題の仮定は、P(Yn>x)=(1−x/n)n≤e−xP(Y_n>x)=(1-x/n)^n\le\mathrm{e}^{-x} と裾の公式(第5章 §5.3)から確かめられる。すべての r>0r>0 で次が成り立つからである。

E[Ynr]=∫0∞rxr−1P(Yn>x) dx≤∫0∞rxr−1e−x dx=Γ(r+1)E[Y_n^r]=\int_0^\infty rx^{r-1}P(Y_n>x)\,\mathrm{d} x\le\int_0^\infty rx^{r-1}\mathrm{e}^{-x}\,\mathrm{d} x=\Gamma(r+1)

例題 11.4 最大値の対数を Slutsky の定理で扱う ★★

オリジナル例題 11.1 の記号を使う(Mn=max⁡(U1,…,Un)M_n=\max(U_1,\dots,U_n)、Yn=n(1−Mn)→dEx(1)Y_n=n(1-M_n)\xrightarrow{d}\mathrm{Ex}(1))。

  1. (1) Wn=−nlog⁡MnW_n=-n\log M_n の分布をすべての nn について正確に求めよ。
  2. (2) (1) を使わずに、Slutsky の定理と連続写像定理から Wn→dEx(1)W_n\xrightarrow{d}\mathrm{Ex}(1) を導け。
  3. (3) Vn=n(1Mn−1)V_n=n\bigl(\frac{1}{M_n}-1\bigr) について Vn→dEx(1)V_n\xrightarrow{d}\mathrm{Ex}(1) を示せ。
解答

(1) x≥0x\ge0 について

P(Wn>x)=P(log⁡Mn<−x/n)=P(Mn<e−x/n)=(e−x/n)n=e−xP(W_n>x)=P\bigl(\log M_n<-x/n\bigr)=P\bigl(M_n<\mathrm{e}^{-x/n}\bigr)=\bigl(\mathrm{e}^{-x/n}\bigr)^n=\mathrm{e}^{-x}

である。よって WnW_n はすべての番号で正確に Ex(1)\mathrm{Ex}(1) に従う。これは、−log⁡Ui∼Ex(1)-\log U_i\sim\mathrm{Ex}(1) の最小値が Ex(n)\mathrm{Ex}(n) に従うことの言い換えでもある。

(2) 次の関数を置く。

h(m)=−log⁡m1−m(0<m<1),h(1)=1h(m)=\frac{-\log m}{1-m}\quad(0<m<1),\qquad h(1)=1

m→1m\to1 では次の展開が成り立つ。

−log⁡m=(1−m)+12(1−m)2+⋯-\log m=(1-m)+\frac12(1-m)^2+\cdots

よって h(m)→1h(m)\to1 であり、関数 hh は点 1 で連続である。Mn∈(0,1)M_n\in(0,1) のとき、次のように書ける(MnM_n が 0 または 1 となる確率は 0)。

Wn=n(1−Mn)⋅h(Mn)=Yn h(Mn)W_n=n(1-M_n)\cdot h(M_n)=Y_n\,h(M_n)

例題 11.1 (1) と定理 11.7 (i) から、Mn→P1M_n\xrightarrow{P}1 である。連続写像定理 (ii) から h(Mn)→Ph(1)=1h(M_n)\xrightarrow{P}h(1)=1 である。Slutsky の定理 (ii) を、分布収束する列 YnY_n と定数 1 に確率収束する列 h(Mn)h(M_n) に使う。すると次となる。

Wn→d1⋅Ex(1)W_n\xrightarrow{d}1\cdot\mathrm{Ex}(1)

これは (1) の正確な結果と一致する。

(3) 次のように書ける。

Vn=n 1−MnMn=YnMnV_n=n\,\frac{1-M_n}{M_n}=\frac{Y_n}{M_n}

分母は Mn→P1≠0M_n\xrightarrow{P}1\ne0 を満たす。よって Slutsky の定理 (iii) から、次となる。

Vn→dEx(1)/1=Ex(1)V_n\xrightarrow{d}\mathrm{Ex}(1)/1=\mathrm{Ex}(1)

例題 11.5 比率の推定と分散の置き換え ★★★

オリジナルSn∼Bin(n,p)S_n\sim\mathrm{Bin}(n,p)(0<p<10<p<1)とする。標本比率を p^n=Sn/n\hat p_n=S_n/n と置く。第13章 §13.2 で示す中心極限定理の特別な場合

Zn=Sn−npnp(1−p)→dN(0,1)Z_n=\frac{S_n-np}{\sqrt{np(1-p)}}\xrightarrow{d}\mathrm N(0,1)

を認める。

  1. (1) Markov の不等式(p=2p=2 の場合、Chebyshev の不等式)を使って p^n→Pp\hat p_n\xrightarrow{P}p を示せ。
  2. (2) 次のように置く。

Tn={p^n−pp^n(1−p^n)/n(p^n∈(0,1))Zn(それ以外)T_n=\begin{cases}\dfrac{\hat p_n-p}{\sqrt{\hat p_n(1-\hat p_n)/n}}&(\hat p_n\in(0,1))\\[2ex] Z_n&(\text{それ以外})\end{cases}

Tn→dN(0,1)T_n\xrightarrow{d}\mathrm N(0,1) を示せ。3. (3) n=400n=400 で、p^n=0.3\hat p_n=0.3 が観測されたとする。近似 P(∣Tn∣≤1.96)≈0.95P(\left\vert T_n\right\vert \le1.96)\approx0.95 から、pp の範囲を求めよ。

解答

(1) 平均と分散は E[p^n]=pE[\hat p_n]=p、V[p^n]=p(1−p)/nV[\hat p_n]=p(1-p)/n である。(11.4) を Y=p^n−pY=\hat p_n-p、指数 2 として使うと

P(∣p^n−p∣≥ε)≤E[(p^n−p)2]ε2=p(1−p)nε2→0P(\left\vert \hat p_n-p\right\vert \ge\varepsilon )\le\frac{E[(\hat p_n-p)^2]}{\varepsilon ^2}=\frac{p(1-p)}{n\varepsilon ^2}\to0

である。

(2) 次の関数を置く。

h(q)=p(1−p)q(1−q)(q∈(0,1)),h(q)=1(q∈{0,1})h(q)=\sqrt{\frac{p(1-p)}{q(1-q)}}\quad(q\in(0,1)),\qquad h(q)=1\quad(q\in\{0,1\})

この関数は点 pp で連続で、h(p)=1h(p)=1 である。p^n∈(0,1)\hat p_n\in(0,1) のとき

Tn=p^n−pp(1−p)/n⋅p(1−p)p^n(1−p^n)=Zn h(p^n)T_n=\frac{\hat p_n-p}{\sqrt{p(1-p)/n}}\cdot\sqrt{\frac{p(1-p)}{\hat p_n(1-\hat p_n)}}=Z_n\,h(\hat p_n)

である。p^n∈{0,1}\hat p_n\in\{0,1\} のときも Tn=Zn=Znh(p^n)T_n=Z_n=Z_nh(\hat p_n) である。よって常に Tn=Zn h(p^n)T_n=Z_n\,h(\hat p_n) である。

(1) と連続写像定理 (ii) から h(p^n)→P1h(\hat p_n)\xrightarrow{P}1 である。Slutsky の定理 (ii) から、次を得る。

Tn→d1⋅Z=N(0,1)T_n\xrightarrow{d}1\cdot Z=\mathrm N(0,1)

(3) 条件 ∣Tn∣≤1.96\left\vert T_n\right\vert \le1.96 は、次と同じである。

∣p^n−p∣≤1.96p^n(1−p^n)/n\left\vert \hat p_n-p\right\vert \le1.96\sqrt{\hat p_n(1-\hat p_n)/n}

0.3×0.7/400=0.000525≈0.02291\sqrt{0.3\times0.7/400}=\sqrt{0.000525}\approx0.02291 なので、右辺は約 0.0449 である。よって

0.2551≤p≤0.34490.2551\le p\le0.3449

が近似的に確率 0.95 で成り立つ範囲である。分母の p(1−p)p(1-p) を観測値 p^n(1−p^n)\hat p_n(1-\hat p_n) で置き換えても、極限の分布は変わらない。これを保証しているのが、(2) の Slutsky の定理である。

確認テスト

オリジナル

  1. 11.4-1 Xn→dN(0,1)X_n\xrightarrow{d}\mathrm N(0,1) とする。P(Xn2≤1)P(X_n^2\le1) の極限を求めよ。数値 Φ(1)≈0.8413\Phi(1)\approx0.8413 を使ってよい。
  2. 11.4-2 Xn∼Po(n)X_n\sim\mathrm{Po}(n) とし、第13章の結果 (Xn−n)/n→dN(0,1)(X_n-n)/\sqrt n\xrightarrow{d}\mathrm N(0,1) を認める。次のように置く。

Rn={(Xn−n)/Xn(Xn>0)0(Xn=0)R_n=\begin{cases}(X_n-n)/\sqrt{X_n}&(X_n>0)\\0&(X_n=0)\end{cases}

Rn→dN(0,1)R_n\xrightarrow{d}\mathrm N(0,1) を示せ。3. 11.4-3 定理 11.16 の逆関数 Fn−1(u)F_n^{-1}(u) を、次の2つの列について求めよ。極限の逆関数に収束しない点があるかどうかを調べよ。(a) Xn∼U(0,1+1n)X_n\sim\mathrm U(0,1+\frac1n) と極限 U(0,1)\mathrm U(0,1)。(b) n≥3n\ge3 で次の分布と、極限 P(X=1)=P(X=0)=12P(X=1)=P(X=0)=\frac12。

P(Xn=1)=12+1n,P(Xn=0)=12−1nP(X_n=1)=\frac12+\frac1n,\qquad P(X_n=0)=\frac12-\frac1n 4. 11.4-4 同じ確率空間の上で Xn→PXX_n\xrightarrow{P}X、Yn→PYY_n\xrightarrow{P}Y とする。次の2つを示せ。

Xn+Yn→PX+Y,XnYn→PXYX_n+Y_n\xrightarrow{P}X+Y,\qquad X_nY_n\xrightarrow{P}XY 5. 11.4-5 列の分布を次で定める。

P(Xn=n)=1n,P(Xn=0)=1−1nP(X_n=n)=\frac1n,\qquad P(X_n=0)=1-\frac1n

Xn→d0X_n\xrightarrow{d}0 であるが E[Xn]↛0E[X_n]\not\to0 であることを確かめよ。命題 11.22 の仮定のどこが成り立たないかを述べよ。

木下の勉強ステータス

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

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

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