知识点

大数定律

定义

大数定律刻画样本均值的稳定性:随着样本量增大,样本均值偏离其期望的可能性越来越小。要让“可能性越来越小”有确切含义,先引入依概率收敛。

设 Y1,Y2,…Y_1,Y_2,\ldots 是定义在同一概率空间上的随机变量序列,aa 为常数。若对任意 ε>0\varepsilon>0 都有

lim⁡n→∞P{∣Yn−a∣≥ε}=0,\lim_{n\to\infty}P\{|Y_n-a|\ge\varepsilon\}=0,

则称 YnY_n 依概率收敛于 aa,记作 Yn→PaY_n\xrightarrow{P}a。

它只说明“偏差超过固定误差 ε\varepsilon 的概率”趋于零,不说明某个有限的 nn 之后偏差永远小于 ε\varepsilon,也不要求每个样本结果逐点趋近 aa。

本页记 E(X)E(X) 为数学期望,D(X)=E{[X−E(X)]2}D(X)=E\{[X-E(X)]^2\} 为方差(也写作 Var⁡(X)\operatorname{Var}(X)),Xˉn=1n∑i=1nXi\bar X_n=\dfrac1n\sum_{i=1}^nX_i 为样本均值。期望的定义和性质见数学期望,随机变量的独立性见随机变量的独立性。

于是大数定律就是一类断言:在适当条件下,样本均值与它的期望之差依概率收敛于 00。不同的大数定律差别在于对 XiX_i 的独立性、同分布性和矩条件的要求。

三个大数定律

证明中反复用到切比雪夫不等式:若 D(Y)D(Y) 存在,则对任意 ε>0\varepsilon>0,

P{∣Y−E(Y)∣≥ε}≤D(Y)ε2.P\{|Y-E(Y)|\ge\varepsilon\}\le\frac{D(Y)}{\varepsilon^2}.

切比雪夫大数定律

设 X1,X2,…X_1,X_2,\ldots 相互独立,方差存在且一致有界:存在常数 CC,使 D(Xi)≤CD(X_i)\le C 对所有 ii 成立。则

1n∑i=1nXi−1n∑i=1nE(Xi)→P0.\frac1n\sum_{i=1}^nX_i-\frac1n\sum_{i=1}^nE(X_i)\xrightarrow{P}0.

更一般地,只要

1n2∑i=1nD(Xi)→0,\frac1{n^2}\sum_{i=1}^nD(X_i)\to0,

相同结论也成立。相互独立可以减弱为两两不相关:证明只用到方差的可加性,而两两不相关已足够保证这一点。

证明

记 Xˉn=1n∑i=1nXi\bar X_n=\frac1n\sum_{i=1}^nX_i。由期望的线性性,

E(Xˉn)=1n∑i=1nE(Xi).E(\bar X_n)=\frac1n\sum_{i=1}^nE(X_i).

因为 XiX_i 两两不相关(相互独立时自然成立),方差可以相加:

D(Xˉn)=1n2∑i=1nD(Xi)≤Cn.D(\bar X_n)=\frac1{n^2}\sum_{i=1}^nD(X_i)\le\frac{C}{n}.

对任意 ε>0\varepsilon>0,把切比雪夫不等式用于 Xˉn\bar X_n:

P{∣Xˉn−E(Xˉn)∣≥ε}≤D(Xˉn)ε2≤Cnε2→0.P\{|\bar X_n-E(\bar X_n)|\ge\varepsilon\} \le\frac{D(\bar X_n)}{\varepsilon^2} \le\frac{C}{n\varepsilon^2}\to0.

所以 Xˉn−E(Xˉn)→P0\bar X_n-E(\bar X_n)\xrightarrow{P}0。若只假设 1n2∑i=1nD(Xi)→0\frac1{n^2}\sum_{i=1}^nD(X_i)\to0,最后一步换成

P{∣Xˉn−E(Xˉn)∣≥ε}≤1n2ε2∑i=1nD(Xi)→0,P\{|\bar X_n-E(\bar X_n)|\ge\varepsilon\} \le\frac1{n^2\varepsilon^2}\sum_{i=1}^nD(X_i)\to0,

结论仍然成立。

独立同分布、方差有限的情形

这是最常用的形式。若 X1,X2,…X_1,X_2,\ldots 独立同分布,E(Xi)=μE(X_i)=\mu,D(Xi)=σ2<∞D(X_i)=\sigma^2<\infty,则

Xˉn→Pμ,\bar X_n\xrightarrow{P}\mu,

并且切比雪夫不等式给出收敛速度的估计:

P{∣Xˉn−μ∣≥ε}≤σ2nε2→0.P\{|\bar X_n-\mu|\ge\varepsilon\}\le\frac{\sigma^2}{n\varepsilon^2}\to0.

同分布使方差自动一致有界,所以这是切比雪夫大数定律的直接推论。

伯努利大数定律

在成功概率为 pp 的 nn 次独立重复试验中,令 NnN_n 为成功次数。则 NnN_n 服从二项分布 B(n,p)B(n,p)(见二项分布),

E(Nnn)=p,D(Nnn)=p(1−p)n≤14n.E\left(\frac{N_n}{n}\right)=p,\qquad D\left(\frac{N_n}{n}\right)=\frac{p(1-p)}{n}\le\frac1{4n}.

由切比雪夫不等式,

P{∣Nnn−p∣≥ε}≤p(1−p)nε2→0,P\left\{\left|\frac{N_n}{n}-p\right|\ge\varepsilon\right\} \le\frac{p(1-p)}{n\varepsilon^2}\to0,

即

Nnn→Pp.\frac{N_n}{n}\xrightarrow{P}p.

频率 Nnn\frac{N_n}{n} 依概率收敛于概率 pp,这正是用大量重复试验的频率估计概率的依据。

辛钦大数定律

若 X1,X2,…X_1,X_2,\ldots 独立同分布,且 E∣Xi∣<∞E|X_i|<\infty,记 μ=E(Xi)\mu=E(X_i),则

Xˉn→Pμ.\bar X_n\xrightarrow{P}\mu.

与有限方差版本相比,辛钦大数定律不要求 D(Xi)D(X_i) 有限,只要求一阶绝对矩存在,因此适用于方差无穷但均值有限的分布。它的证明需要特征函数等工具,这里直接使用结论。

条件的比较

定理条件结论
切比雪夫大数定律相互独立(可减弱为两两不相关),D(Xi)D(X_i) 存在且一致有界;更一般地,1n2∑iD(Xi)→0\frac1{n^2}\sum_iD(X_i)\to01n∑iXi−1n∑iE(Xi)→P0\frac1n\sum_iX_i-\frac1n\sum_iE(X_i)\xrightarrow{P}0
伯努利大数定律nn 次独立重复试验,成功概率为 pp,NnN_n 为成功次数Nnn→Pp\frac{N_n}{n}\xrightarrow{P}p
辛钦大数定律独立同分布,一阶绝对矩有限,均值 μ=E(Xi)\mu=E(X_i)Xˉn→Pμ\bar X_n\xrightarrow{P}\mu

样本量的保守估算

设 X1,X2,…X_1,X_2,\ldots 独立同分布,E(Xi)=μE(X_i)=\mu,D(Xi)=σ2<∞D(X_i)=\sigma^2<\infty。若希望

P{∣Xˉn−μ∣<ε}≥1−α,P\{|\bar X_n-\mu|<\varepsilon\}\ge1-\alpha,

切比雪夫不等式给出一个充分条件。由

P{∣Xˉn−μ∣<ε}≥1−σ2nε2P\{|\bar X_n-\mu|<\varepsilon\}\ge1-\frac{\sigma^2}{n\varepsilon^2}

可知,只要 σ2nε2≤α\frac{\sigma^2}{n\varepsilon^2}\le\alpha,即

n≥σ2αε2,n\ge\frac{\sigma^2}{\alpha\varepsilon^2},

要求就能满足;实际取不小于该值的整数。

对成功频率,D(Nnn)=p(1−p)nD\left(\frac{N_n}{n}\right)=\frac{p(1-p)}{n} 且 p(1−p)≤14p(1-p)\le\frac14,所以未知 pp 时可用

n≥14αε2,n\ge\frac1{4\alpha\varepsilon^2},

再向上取整。这是保守保证:切比雪夫不等式不依赖具体分布,代价是通常比用正态近似算出的样本量大。

例子

用切比雪夫不等式估算样本量

设总体方差 σ2=4\sigma^2=4,要求样本均值与总体均值之差小于 0.50.5 的概率至少为 0.950.95。取 ε=0.5\varepsilon=0.5、α=0.05\alpha=0.05,则

n≥40.05×0.52=320.n\ge\frac{4}{0.05\times0.5^2}=320.

所以至少需要 320320 个样本。这里的 320320 是充分条件给出的保守值,实际所需样本量可能更小。

未知成功概率时的样本量

要求频率与成功概率 pp 的偏差小于 0.010.01 的概率至少为 0.950.95,但不知道 pp。利用 p(1−p)≤14p(1-p)\le\frac14,取 ε=0.01\varepsilon=0.01、α=0.05\alpha=0.05,得

n≥14×0.05×0.012=50000.n\ge\frac{1}{4\times0.05\times0.01^2}=50000.

即使不知道 pp,也能给出保证,代价是样本量较大。

切比雪夫界有多保守

设 p=12p=\frac12,做 n=10000n=10000 次独立重复试验,取 ε=0.01\varepsilon=0.01。切比雪夫不等式给出

P{∣Nnn−p∣≥0.01}≤0.2510000×0.0001=0.25.P\left\{\left|\frac{N_n}{n}-p\right|\ge0.01\right\} \le\frac{0.25}{10000\times0.0001}=0.25.

而 Nnn\frac{N_n}{n} 的标准差为 0.0050.005。记 Φ\Phi 为标准正态分布函数(见正态分布),用正态近似估计同一概率约为 2(1−Φ(2))≈0.04552(1-\Phi(2))\approx0.0455。两者相差一个数量级:切比雪夫不等式只提供宽松的上界,不适合用来精确计算概率。

方差无穷时仍可使用辛钦大数定律

设 XiX_i 独立同分布,服从尾指数 1<α<21<\alpha<2 的帕累托分布。此时 E∣Xi∣<∞E|X_i|<\infty 但 D(Xi)=∞D(X_i)=\infty。切比雪夫大数定律不适用,因为方差不存在;辛钦大数定律仍给出 Xˉn→Pμ\bar X_n\xrightarrow{P}\mu。这说明“方差有限”只是切比雪夫版本的条件,不是大数定律本身必需的。

交互探索

用伪随机数模拟成功概率为 1/21/2 的独立重复试验,调整次数可查看同一序列的不同前缀,刷新可观察另一条模拟路径。频率会波动;模拟帮助理解现象,不能代替定理的证明。

试验次数增加时,频率怎样变化?

通过可刷新的模拟,观察成功频率随试验次数的变化。

正在准备演示…

题集

练习

  1. 某总体方差为 99。用切比雪夫不等式估计,要使样本均值与总体均值之差小于 0.30.3 的概率至少为 0.950.95,样本量至少为多少?

    参考解答

    由 n≥σ2αε2=90.05×0.32=2000n\ge\dfrac{\sigma^2}{\alpha\varepsilon^2}=\dfrac{9}{0.05\times0.3^2}=2000。

  2. 下列条件中,哪些能保证大数定律的结论成立?

    题 conditions · options-conditions
    参考解答

    前三项分别对应切比雪夫大数定律、辛钦大数定律和切比雪夫大数定律的独立同分布情形。第四项不满足切比雪夫大数定律的方差条件,不能由此保证结论成立。

  3. 关于伯努利大数定律,下列哪一项正确?

    选择答案
    参考解答

    依概率收敛只控制偏差超过 ε\varepsilon 的概率,不保证单次试验序列的偏差,也不保证单调或最终相等;定理对任意 0<p<10<p<1 都成立。