DennyQi's Log

04 大数定律

收敛理论

点态收敛

在定义连续随机变量的期望时候,我们是用一列离散的随机变量期望的极限来定义的。一般地,我们也可以定义一列随机变量的极限,这个极限也是一个随机变量,而我们知道随机变量本质上是一个函数,这个极限过程正是数学分析中的函数列的收敛。我们所说的随机变量的极限就是随机变量列的点态收敛。而在概率论中,我们更多时候会用almost surely(a.s.)点态收敛:只要求随机变量在一个测度为1的集合上点态收敛,也即只在一个零测集上不收敛。

关于点态收敛要讨论的一个最重要的问题就是极限和期望的顺序交换问题——是否成立limnE[Xn]=E[limnXn]\lim\limits_{n \to \infty} \mathbb{E}[X_n]=\mathbb{E}[\lim\limits_{n \to \infty} X_n]?(例如在Moment Generating Function一节中我们就默认了这一事实成立而没有加以验证)。

首先我们在([0,1],B([0,1]),PLeb)([0,1],\mathcal{B}([0,1]),P_{\text{Leb}})上有反例Xn=n1[1n,2n]X_n=n \cdot \mathbb{1}_{[\frac{1}{n},\frac{2}{n}]}来说明这一事实并不总是成立。对于任意固定的nnE[Xn]\mathbb{E}[X_n]都为1;而limnXn\lim\limits_{n \to \infty}X_n却a.s.等于0。因此limnE[Xn]=10=E[limnXn]\lim\limits_{n \to \infty} \mathbb{E}[X_n]=1\neq 0=\mathbb{E}[\lim\limits_{n \to \infty} X_n]。那么这一事实在何时成立呢?下面我们给出几个关于充分条件的定理(证明略):

第一个充分条件称为Monotone Convergence Theorem(MCT,单调收敛定理),它指出:如果随机变量列XnX_n非负且递增并收敛到XX(以上条件都只需a.s.成立),那么极限和期望可交换:limnE[Xn]=E[limnXn]=E[X]\lim\limits_{n \to \infty} \mathbb{E}[X_n]=\mathbb{E}[\lim\limits_{n \to \infty} X_n]=\mathbb{E}[X]

第二个充分条件称为Dominated Convergence Theorem(DCT,控制收敛定理),它指出:如果随机变量列XnX_n收敛到XX(a.s.),并且所有的XnX_n都能被一个E[Y]\mathbb{E}[Y]存在的随机变量YYXnY|X_n|\leq Y的方式控制(a.s.),那么极限和期望可交换:limnE[Xn]=E[limnXn]=E[X]\lim\limits_{n \to \infty} \mathbb{E}[X_n]=\mathbb{E}[\lim\limits_{n \to \infty} X_n]=\mathbb{E}[X]。特别地,如果YY取常数函数,那么XnY|X_n|\leq Y恒成立等价于{Xn}\{X_n\}有界,这就得到推论Bounded Convergence Theorem(BCT,有界收敛定理)

依概率收敛,LpL_p收敛,依分布收敛以及强弱关系

就像函数不止点态收敛一种收敛方式一样,点态收敛(a.s.收敛)也不是定义随机变量收敛的唯一方式。一般而言,点态收敛是最强的收敛条件了,但我们很多时候我们需要更弱的收敛条件,因为在许多重要的定理中以强的形式收敛的结论往往是不成立的,只有在更弱时成立。

下面我们依次给出依概率收敛、LpL_p收敛、依分布收敛的定义:

如果ε>0,limnPr[XnX>ε]=0\forall \varepsilon>0,\lim\limits_{n \to \infty}\Pr[|X_n-X|>\varepsilon]=0,就称XnX_n依概率收敛(converge in probability)到XX,记为XnpxX_n \stackrel{p}{\to} x。依概率收敛表示当nn充分大时,XnX_nXX上函数值不同的样本点测度趋向0;

如果limnE[XnXp]=0\lim\limits_{n \to \infty}\mathbb{E}[|X_n-X|^p]=0,就称XnX _n LpL_p收敛到XX,记为XnLpXX_n \stackrel{L_p}{\to} X,表示pp阶矩收敛到同一个值。特别的,当p=1p=1时为L1L_1收敛,它们的收敛到相同的期望;

XX的分布函数为F(x)F(x)XnX_n的分布函数为Fn(x)F_n(x)。如果在FF的连续点上始终成立limnFn(x)=F(x)\lim\limits_{n \to \infty}F_n(x)=F(x),就称XnX_n依分布收敛收敛到XX,记为XndxX_n \stackrel{d}{\to} x,表示它们的分布函数收敛到同一个值。

可以证明,r>sr>s时有Lr    LsL_r \implies L_s,也即更高阶的矩收敛可以推出更低阶的。其中最低阶的L1    pL_1 \implies p,这说明LpL_p收敛比依概率收敛更强。同时,a.s.    pa.s. \implies p,几乎处处的点态收敛可以推出依概率收敛。p    dp \implies d,依概率收敛可以推出依分布收敛。可见依分布收敛是最弱的要求。(以上的推出都是不可逆的,构造反例可以说明这一点。)并且我们观察到,a.s.a.s.L1L_1之间的强弱不能直接比较,而这两者正好是a.s.a.s.收敛与期望相等之间的关系——正是我们之前讨论的极限与期望的可交换问题,我们已经知道在特定的充分条件下交换才是成立的。

依分布收敛有以下等价条件:XndX    X_n \stackrel{d}{\to} X \iff对于任意的compactly supported连续函数gg成立limnE[g(Xn)]=E[g(X)]\lim\limits_{n\to\infty}\mathbb{E}[g(X_n)]=\mathbb{E}[g(X)]。其中,compactly supported是指所有使函数值非零的自变量构成的集合是紧集。我们可以证明依分布收敛的Dominated Convergence Theorem:如果XndXX_n \stackrel{d}{\to} X,存在YY使得Xna.s.Y|X_n|\leq_{a.s.} Y恒成立,E[Y]<+\mathbb{E}[Y]<+\infty,则limnE[Xn]=E[X]\lim\limits_{n \to \infty}\mathbb{E}[X_n]=\mathbb{E}[X]

上下极限

虽然有反例说明依概率收敛不能推出a.s.a.s.收敛,但我们可以证明依概率收敛可以推出存在子列几乎处处收敛。为了证明这一点,首先要定义集合列的极限。如果把集合的包含关系看作序关系,那么对于单调的集合列就可以定义极限:对于AiAi+1A_{i} \subseteq A_{i+1},定义limnAn=i1Ai\lim\limits_{n \to \infty} A_n=\bigcup\limits_{i \geq 1}A_i。同理,对于AiAi+1A_{i} \supseteq A_{i+1},定义limnAn=i1Ai\lim\limits_{n \to \infty} A_n=\bigcap\limits_{i \geq 1}A_i。由于是单调的,我们也用上确界或下确界来表示极限。现在,仿照数列的上下极限,定义上极限limsupnAn=limn(supknAk)=n1knAk\lim \sup_n A_n=\lim\limits_{n \to \infty}(\sup\limits_{k \geq n}A_k)=\bigcap\limits_{n \geq 1}\bigcup\limits_{k \geq n}A_k,下极限liminfnAn=limn(infknAk)=n1knAk\lim \inf_n A_n=\lim\limits_{n \to \infty}(\inf\limits_{k \geq n}A_k)=\bigcup\limits_{n \geq 1}\bigcap\limits_{k \geq n}A_k。上下极限也表示一个集合,其中上极限表示所有在{Ai}\{A_i\}中出现次数为无数次的元素构成的集合(如果出现无数次,那么对任意的nn都会落在supknAk\sup\limits_{k \geq n}A_k里,因此最终落在limsupnAn\lim \sup_nA_n中;否则一定存在一个nn使得它不在supknAk\sup\limits_{k \geq n}A_k里,因此最终不在上极限中);下极限表示所有不出现次数为有限次的元素构成的集合。

一列事件就是一列集合。我们可以根据定义化简一列事件的上极限的概率:Pr[limsupnAn]=Pr[limnknAk]\Pr[\lim\sup_n A_n]=\Pr[\lim\limits_{n \to \infty} \bigcup \limits_{k \geq n}A_k],根据概率测度的连续性=limnPr[knAk]=limnknPr[Ak]=\lim\limits_{n \to \infty}\Pr[\bigcup\limits_{k \geq n}A_k]=\lim\limits_{n \to \infty}\sum\limits_{k \geq n}\Pr[A_k]。可见,如果k1Pr[Ak]<+\sum\limits_{k\geq 1}\Pr[A_k]<+\infty,那么一定有Pr[limsupnAn]\Pr[\lim\sup_nA_n]。这就是Borel-Cantelli定理,它指出如果一列事件A1,A2,A_1,A_2,\cdots满足n1Pr[An]<+\sum\limits_{n \geq 1}\Pr[A_n]<+\infty,则Pr[limsupnAn]=0\Pr[\lim\sup_n A_n]=0。也即如果所有这些事件发生的概率全部相加是收敛的,那么在这列事件中出现无穷多次的样本点是零测集。它的逆命题不一定成立,然而我们可以验证当AnA_n相互独立时,逆命题成立。此时n1Pr[An]<+    Pr[limsupnAn]=0\sum\limits_{n \geq 1}\Pr[A_n]<+\infty \iff \Pr[\lim\sup_n A_n]=0。我们还可以证明,n1Pr[An]=+    Pr[limsupnAn]=1\sum\limits_{n \geq 1}\Pr[A_n]=+\infty \implies \Pr[\lim\sup_n A_n]=1,这意味着Pr[limsupnAn]\Pr[\lim\sup_nA_n]只能取0或1(我们之后将会用Kolmogorov 0-1 Law这个更高的观点再次看到这个问题),因此也有n1Pr[An]=+    Pr[limsupnAn]=1\sum\limits_{n \geq 1}\Pr[A_n]=+\infty \iff \Pr[\lim\sup_n A_n]=1

根据Borel-Cantelli,我们从一个依概率收敛的随机变量列中挑出一列nmn_m使得Pr[XnmX>1m]<12m\Pr[|X_{n_m}-X|>\dfrac{1}{m}]<\dfrac{1}{2^m},令Am={ωXnm(ω)X(ω)>1m}A_m=\{\omega\mid |X_{n_m}(\omega)-X(\omega)|>\dfrac{1}{m}\},这样就有m1Pr[Am]<m112m<+\sum\limits_{m\geq 1}\Pr[A_m]<\sum\limits_{m \geq 1}\dfrac{1}{2^m}<+\infty,因此Pr[limsupmAm]=0\Pr[\lim\sup_m A_m]=0。在全集中去掉这个零测集以后,我们可以证出点态收敛。因此我们证明了依概率收敛的随机变量列中存在一个a.s.点态收敛的子列。

有了这个定理以后,我们就可以把Dominated Convergence Theorem中的几乎处处收敛放弱到“依概率收敛”。原因是,如果E[Xn]\mathbb{E}[X_n]不收敛到E[X]\mathbb{E}[X],那么由于依概率收敛,它存在子列收敛到aE[X]a \neq \mathbb{E}[X]。而依概率收敛还意味着其任意子序列依概率收敛,因此上面的子序列的子序列必须a.s.a.s.收敛到XX,它的期望必须收敛到E[X]\mathbb{E}[X],矛盾。

大数定律

我们在定义概率空间和随机变量时是从集合和函数出发的,而当我们想要真正理解概率的“意义”时,其实我们已经在使用了大数定律这一事实。硬币正面朝上的概率为1/21/2这句话的意思是,当投掷硬币的次数充分大以至于是一个“大数”时,应当期待有接近一半的次数投掷硬币正面朝上。大数定律描述的就是同一随机事件在被重复足够多次时会收敛到它的期望。

强大数定理与弱大数定理

X1,,Xn,X_1,\cdots,X_n,\cdots是相互独立且同分布(independent and identically distributed, i.i.d.)的随机变量,大数定理要描述i[n]Xin\dfrac{\sum_{i \in [n]}X_i}{n}(记为Snn\dfrac{S_n}{n})以何种方式收敛到E[Xi]\mathbb{E}[X_i](记为μ\mu)。我们已经知道随机变量的收敛是有许多不同强弱的种类的。Snnpμ\dfrac{S_n}{n} \stackrel{p}{\to} \mu这一事实称为弱大数定理(Weak Law of Large Numbers, WLLN),Snna.s.μ\dfrac{S_n}{n} \stackrel{a.s.}{\to} \mu这一事实称为强大数定理(Strong Law of Large Numbers, SLLN)。

我们首先在附加上二阶矩有限(E[Xi2]σ2\mathbb{E}[X_i^2] \leq \sigma^2)的前提下证明弱大数定理,这只需用Markov不等式说明Pr[Snn>ε]=Pr[Snn2>ε2]E[(Snn)2]ε2σ2ε2n\Pr[\left|\dfrac{S_n}{n}\right|>\varepsilon]=\Pr[\left|\dfrac{S_n}{n}\right|^2>\varepsilon^2]\leq \dfrac{\mathbb{E}[\left(\frac{S_n}{n}\right)^2]}{\varepsilon^2}\leq\dfrac{\sigma^2}{\varepsilon^2n},因此Snn\dfrac{S_n}{n}依概率收敛。在同样的前提下,为了证明强大数定理,我们也想用Markov不等式,结合n=1Pr[Snn>ε]<+\sum\limits_{n=1}^{\infty}\Pr[\left|\dfrac{S_n}{n}\right|>\varepsilon]<+\infty用Borel-Cantelli说明a.s.点态收敛,此时我们发现仅规定二阶矩有限是不够的,为此我们附加四阶矩有限的条件,用相同的方法得到证明。

现在我们要去掉二阶矩有限的条件,证明真正的弱大数定理。此时我们不再能直接运用Markov不等式了,因为二阶矩可能是无界的。这里我们要用到称为truncation(截断)的证明思路:我们把随机变量拆分成>M>MM\leq M两种情形,于是Pr[Snnμ>ε]Pr[Sn,Mnμ>ε]+Pr[Sn,>M0]\Pr[\left|\dfrac{S_n}{n}-\mu\right|>\varepsilon]\leq \Pr[\left|\dfrac{S_{n,\leq M}}{n}-\mu\right|>\varepsilon]+\Pr[S_{n,>M}\neq 0]。取M=nM=n,前者我们把随机变量的取值控制在了有限范围内,后者在nn \to \infty时显然趋向0,于是我们发现我们能够证明这两个概率都趋向0,这样就证明了弱大数定理。

我们暂时还不能给出强大数定理的证明。

Kolmogorov 0-1 Law

从更一般的观点来看大数定律,它其实指出了当nn趋向无穷时,Pr[Snnμ>ε]\Pr[\left|\dfrac{S_n}{n}-\mu\right|>\varepsilon]总为0(弱大数定理),Pr[Snn=μ]\Pr[\dfrac{S_n}{n}=\mu]总为1(强大数定理)。在Borel-Cantelli中,我们也看到了Pr[limsupnAn]\Pr[\lim\sup_n A_n]总是只能取0或者1。。事实上这是一个更为普遍的规律,我们能够证明一列相互独立事件的极限事件(tail event)发生的概率总是0或1的。这就是Kolmogorov 0-1 Law。

我们首先要定义什么是极限事件。为此,我们要定义关于随机变量的σ\sigma-algebra。对于随机变量XX,定义σ(X)\sigma(X)为能使得XX可测的最小σ\sigma-algebra。在定义随机变量时,我们已经要求它在所有Borel Set下的原像落在事件集里,那么我们直接取出所有这些原像X1(B(R))X^{-1}(\mathcal{B}(\R)),可以证明这本身就是一个σ\sigma-algebra,因此直接有σ(X)=X1(B(R))\sigma(X)=X^{-1}(\mathcal{B}(\R))。我们可以这样理解“最小可测”,我们知道2Ω2^\Omega总是一个使得XX可测的事件集,但有时XX的特性使得它并不会用到全部这些子集,例如当XX仅仅只是骰子是奇数还是偶数时,我们便无需关心{1,2},{1,3,4,5}\{1,2\},\{1,3,4,5\}这样的集合,而只需关心{1,3,5},{2,4,6}\{1,3,5\},\{2,4,6\}这两个集合。换言之,使得不同的XX可测需要的其实是不同大小的σ\sigma-algebra,这和XX本身包含的“信息”有关。如果我们关心骰子的具体取值,那么我们需要一个相对庞大的σ\sigma-algebra;而如果只关心骰子的奇偶,则只需要一个较小的σ\sigma-algebra。而一旦知道了具体取值,我们就一定知道了奇偶,因此我们说前者包含了后者的信息。σ(X)\sigma(X)刻画了XX包含的信息。如果σ(Y)σ(X)\sigma(Y)\subseteq \sigma(X),说明可以用σ(X)\sigma(X)来测YY,也就说明XX中包含比YY更多的信息。对于多个随机变量,我们定义σ(X1,X2)=σ(σ(X1)σ(X2))\sigma(X_1,X_2)=\sigma(\sigma(X_1)\cup \sigma(X_2)),也就是使得X1,X2X_1,X_2都可测的最小σ\sigma-algebra。定义σ\sigma-algebra F,G\mathcal{F},\mathcal{G}独立当且仅当AF,BG\forall A \in \mathcal{F},B\in\mathcal{G}都有AABB独立,容易根据定义证明XY    σ(X)σ(Y)X \bot Y\iff \sigma(X)\bot \sigma(Y)

现在我们定义极限事件。对一列相互独立的随机变量 X1,X2,X_1,X_2,\cdots ,定义Fn=σ(X1,X2,...,Xn)\mathcal{F} _n=\sigma(X_1,X_2,...,X_n)F=σ(X1,X2,...)\mathcal{F}_{\infty}=\sigma(X_1,X_2,...)。容易验证F=σ(n1Fn)\mathcal{F}_\infty=\sigma(\bigcup\limits_{n\ge 1}\mathcal{F}_n) 。定义 Fn=σ(Xn+1,Xn+2,...)\mathcal{F}_n^*=\sigma(X_{n+1},X_{n+2},...)F=n0Fn\mathcal{F} _{\infty}^*=\bigcap\limits_{n\ge 0}\mathcal{F}_n^* ,其中F\mathcal{F} _{\infty}^*被称为tail algebra。任何F\mathcal{F}^*_\infty中的事件就称为极限事件。极限事件与任意有限的XnX_n中的信息无关,只与极限过程中的随机变量的信息有关。

Kolmogorov 0-1 Law指出,AF,P(A)=0\forall A\in \mathcal{F} _{\infty}^*,P(A)=011。Pf:“P(A)=0P(A)=011”可以转化为AAA\bot A,因为AAA\bot A的定义恰好是P(AA)=P(A)P(A)P(A\cap A)=P(A)\cdot P(A),也即P(A)=P(A)2P(A)=P(A)^2,解得P(A)=01P(A)=0或1。我们把所有与AA独立的事件收集进集合H\mathcal{H},那么只需证AHA \in \mathcal{H}。显然任意有限的Fn\mathcal{F}_n都与Fn\mathcal{F}^*_n独立(一个描述前nn项的信息,一个描述nn以后的信息),而AFnA \in \mathcal{F}_n^*,因此对任意的nn总有FnHF_n\in \mathcal{H}。也即n1FnH\bigcup\limits_{n \geq 1}\mathcal{F}_n\subseteq \mathcal{H}。于是可以证明(从略)FH\mathcal{F}_\infty \in \mathcal{H}。而AFA \in \mathcal{F}_{\infty},因此AHA\in \mathcal{H},证毕。

大数定理中的Snn\dfrac{S_n}{n}收敛就是极限事件,因为数列的收敛与任意有限项都无关。因此它要么一概率收敛,要么一概率不收敛(E[Xi]\mathbb{E}[X_i]不收敛);上极限与任意有限项无关,它也是一个极限事件,因此Pr[limsupnAn]\Pr[\lim\sup_n A_n]只能取0或1。