高等概率论逐页精讲

第 3 讲 Borel–Cantelli、收敛与 Radon–Nikodym›第 5 部分 抽样中的应用

PDF 第 37 / 47 页幻灯片 37/49定理

方差与样本均值的误差:均方误差是 σ²/n,误差按 σ/√n 缩小

Variance and the error of a sample average

课件原页 · 点击放大 · ← → 键翻页

这一页在讲什么

人话版一句话版

把 个独立同分布、方差为 的样本求平均,平均值偏离真值的均方误差恰好是 所以误差的典型大小是 样本数乘 4,误差才减半。单个样本的方差 就是一个抽样方案的「单价」。第 36 页 里当 时 这把尺子就失效了。

  • 这是第 5 部分「抽样中的应用」的第二页。上一页 第 36 页 发现:均匀抽样估计量无偏,但方差可能无穷。这一页回答「为什么要在乎方差」:蒙特卡罗方法的精度完全由方差决定。
  • 课件只印了一个小标题 Mean squared error of the sample average(样本均值的均方误差),内容是板书。下面「补全课件留白」一节补全:一个命题(均方误差 它的证明、由它推出的 Chebyshev 概率界,以及方差无穷时会发生什么。
  • 下一页起(第 38 页–第 40 页)的目标由此变得很清楚:换一种抽样方式,让单个样本的方差变成有限,而且尽量小。

先补基础

基础补课均方误差:估计量「平均错多少」
  • 设 是用来估计未知数 的估计量(它依赖随机样本,是随机变量)。它的均方误差(mean squared error,MSE)定义为 误差先平方(去掉正负号,并且大误差罚得更重),再求平均。
  • 均方根误差 把平方开回来,单位和 相同,可以直接读成「典型误差有多大」。
  • 分解:若 则 即「均方误差 = 方差 + 偏差的平方」(证明见「推导」一节第(二)条)。无偏时偏差为 0,均方误差就等于方差。
  • 例:一个永远报告 的「估计量」方差为 0,但均方误差是 1——它稳定地错。另一个一半时候报 一半时候报 偏差为 0,方差和均方误差都是 1。
基础补课方差的三条基本性质

设 这时 因为 (第 32 页 推导一节),所以 是有限的数。

  1. (第 36 页 基础补课)。
  2. 平移不变:因为
  3. 倍数要平方后提出来:因为 平方后多出

例:注意是除以 100,不是除以 10。

基础补课独立随机变量:乘积的期望等于期望的乘积

回忆(第 6 页 补充定义):随机变量 独立,指对任意博雷尔集 本页要用两个事实:

  1. 独立变量的函数仍然独立:若 是博雷尔可测函数,则 与 独立。理由:而 仍是博雷尔集,直接套定义即可。例如 与 独立。
  2. 乘积公式:若 独立并且都可积,则 可积,并且

第 2 条第 1 讲没有讲过。它的证明只用「示性函数 → 简单函数 → 非负函数(单调收敛定理)→ 可积函数」这四步,和 第 32 页–第 33 页 的套路一样,完整写在下面「推导」一节第(一)条。

没有独立性,乘积公式一般不成立:公平硬币 取 则 而 (第 7 页 也讨论过这个例子)。

基础补课Chebyshev 不等式:用方差控制「偏得很远」的概率
  • 广义 Markov 不等式(第 29 页)取 也就是 第 10 页 基础补课里 的情形:对任何随机变量 和
  • 理由:处处有 (在 处左边是 在别处左边是 两边取期望,左边是
  • 用在 上,就是 Chebyshev 不等式:方差小,偏离均值很远的概率就小。
  • 例:则

课件原文逐句精读

课件原文课件印刷文字

Mean squared error of the sample average

中文翻译

样本均值的均方误差

  • 蓝条标题 Variance and the error of a sample average:「方差与样本均值的误差」。它点明本页的主题:样本均值的误差由方差决定。
  • sample average(样本均值):就是把 次抽样的结果加起来除以 第 36 页 的估计量 就是 时的样本均值。
  • mean squared error(均方误差):见基础补课,这里是 其中 是要估计的真值。
  • 课件在这个小标题下面留白。按标题和上下文(上一页刚提出「方差可能无穷」),这里应当写出并证明「均方误差 下一节补全。

补全课件留白

补充说明

课件此处留白(课堂板书),以下内容由本站根据标准教材补写,记号与本课件保持一致;如与老师课堂写法不同,以老师为准。

Mean squared error of the sample average|样本均值的均方误差

命题(样本均值的均方误差)

设 是概率空间 上独立同分布的随机变量,并且 记

则对每个

  1. 无偏:
  2. 均方误差:
  3. Chebyshev 界:对每个 特别地,右边随 趋于 0,所以 依概率收敛(第 14 页 定义 10)。

先看每个条件管什么:

  • :保证 是有限的数,结论才有内容;它也保证 是有限的数。第 36 页 当 时这一条不成立。
  • 同分布:保证每个 的均值都是 方差都是 均值和方差只由分布决定:由无意识统计学家定律(第 7 页), 只依赖于 的分布。
  • 独立:只用在一处——交叉项 (证明里会标出来。
推导证明思路

先把每个样本「中心化」:均值 0、方差 于是 就是 的平均。把 展开成 项: 个平方项各贡献 个交叉项因为独立而期望为 0。合计 再除以 第 3 条是把 Chebyshev 不等式用在 上。

第 1 步:中心化。 令 由同分布,并且

这里把 写成了 个 之和。

第 2 步:无偏。 每个 可积,可积函数的期望是线性的(第 43 页),所以

这就是第 1 条。

第 3 步:把平方展开。 有限和的平方等于所有两两乘积之和:

例如 两个平方项、两个交叉项。

第 4 步:每一项都可积,可以逐项取期望。 由 展开得 所以

有限个可积函数之和的期望等于期望之和(第 43 页),于是

第 5 步:平方项。 每个 一共 项,合计

第 6 步:交叉项(唯一用到独立性的地方)。 设 与 独立,而 分别是它们的函数,所以 与 独立(基础补课第 1 条);两者都可积。由乘积公式,

第 7 步:合起来。 由第 1 步和方差的「倍数平方提出」性质,

又由第 2 步 所以方差的定义 恰好就是 这就是第 2 条。

第 8 步:Chebyshev 界。 固定 把基础补课里的不等式用在 上:

第一个不等号是因为 (概率的单调性);最后一个等号用了第 2 条。 固定时右边随 趋于 0,这正是依概率收敛的定义。这就是第 3 条,它也叫弱大数定律的 版本。

一般形式:独立变量的方差可加

同样的展开对「不同分布」也成立:若 两两独立(相互独立当然更可以),每个 则

证明:令 重复第 3–6 步,交叉项全为 0,平方项之和是 第 6 步只用到了两个变量之间的独立,所以两两独立就够。注意方差可加,标准差不可加:两个独立、标准差都是 1 的变量之和,标准差是 不是 2。

误差的量级:σ/√n

命题第 2 条开平方,得到均方根误差:

它告诉我们典型误差有多大:

样本数
均方根误差

读这张表:

  1. 样本数乘 100,误差只缩小到十分之一。想多得一位准确的小数,要付出 100 倍的计算量。这就是人们常说的「蒙特卡罗方法收敛慢」。
  2. 「几个标准误」的写法:在第 3 条里取 得 叫标准误(standard error)。例如 误差超过 3 个标准误的概率不超过 不超过
  3. 需要多少样本:要让 由第 3 条,只要 即 所需样本数与 成正比。
  4. 两条降低误差的路:(1)加大 代价与 成正比;(2)换一种抽样方式,让单个样本的方差 变小。把 减半,效果等于把样本数乘 4。第 38、39 页走的是第二条路。
补充背景:中心极限定理(本讲不证)

Chebyshev 界对每个 都严格成立,但比较保守。中心极限定理(本讲没有讲,这里只作背景)说: 很大时, 的分布接近标准正态 于是大约 95% 的情形误差不超过 两种说法的共同点是:误差的尺度是 。

方差无穷时:σ/√n 这把尺子失效

回到 第 36 页: 独立同分布、服从

  1. :命题完全适用,误差按 缩小。例如 时 (第 36 页),
  2. :命题的前提不成立。第 36 页「推导」一节第(二)条已经证明:对每个 ,平均多少个样本都一样;Chebyshev 界的右边是 什么也没说。
  3. 估计量还会收敛吗? 会。强大数定律只要求 (本讲不证), 几乎必然。但没有 的速率保证,实际收敛更慢、更不稳定。下面例 3 用模拟展示这一点。

例子

例 1:掷骰子的平均点数

是第 次掷一颗公平骰子的点数,各次独立,

  • 方差:所以
  • :均方误差 均方根误差 Chebyshev:
  • :均方根误差 Chebyshev:

精度提高 10 倍(容忍度从 0.5 缩到 0.05),同样的概率保证需要 100 倍的样本。

例 2:抛硬币——Chebyshev 界是保守的

是独立的公平硬币(正面记 1),(因为

  • 均方根误差
  • Chebyshev:
  • 精确值: 就是「正面次数 或 按二项分布精确计算,概率约为

Chebyshev 给的 0.25 比真实的 0.035 大很多,但它不需要知道分布的具体形状,只用到方差,并且对每个 都成立。

例 3:方差有限与方差无穷,收敛速度对比(模拟)

对 第 36 页 的均匀抽样估计量 每种 独立重复 2000 次,记录每次的误差 取中位数(一半的重复误差比它小,一半比它大)。下面是一次这样的模拟的结果(随机数不同,数字会略有出入):

(方差有限)的中位误差 0.053 0.017 0.0054
(方差无穷)的中位误差 0.61 0.32 0.16
  • 每乘 10,中位误差缩小约 3.1 倍,正是 符合
  • 每乘 10,中位误差只缩小约 2 倍。理论上可以证明(要用本课没有的「稳定分布」理论,这里只作背景),这个例子里典型误差大约按 缩小; 时是 而 与模拟一致。
  • 同一组模拟里, 时约 70% 的重复低估了 (大多数时候抽不到 0 附近的巨大值),而最坏的一次误差超过 100(碰巧抽到了一个离 0 极近的点)。平均来说无偏,是靠极少数巨大的高估撑起来的。
例 4(反例):样本不独立,平均就不降低方差

设 满足 令 (相当于随机数发生器出了故障,把同一个数重复了 次)。

  • 每个 的分布相同,但它们不独立。
  • 所以 不管 多大都不减小。
  • 对照证明:第 6 步的交叉项变成 不再是 0。 个交叉项合计 加上平方项 总共 除以 仍是

所以命题里的「独立」不能去掉。

推导:一步一步来

推导证明思路

补上正文用到的两件事:(一)独立可积变量的乘积公式,走「示性函数 → 简单函数 → 非负函数 → 可积函数」四步;(二)均方误差 = 方差 + 偏差的平方。

(一)若 独立且可积,则 可积,并且

  1. 非负简单函数。 设 是 上的非负简单博雷尔函数(有限和, 是博雷尔集、 也是)。令 因为
    这是 上的非负简单函数。由简单函数积分的线性(第 32 页)和独立性,
  2. 非负的 设 且 取二进逼近 (第 34 页),它是非负简单博雷尔函数,并且对每个 有 令 则 并且 (两个非负递增数列的乘积仍递增,极限等于极限之积)。对三列分别用单调收敛定理(第 36 页),再用第 1 步:
    最后一个等号: 都是有限数的收敛,乘积的极限等于极限的乘积。
  3. 可积的 分别是 的函数,独立且非负,由第 2 步 所以 可积。再把正部、负部(第 43 页) 代入:
    四项都是「 的函数 × 的函数」,非负,并且都 所以可积。由期望的线性和第 2 步,

(二)若 则

  1. 记 (有限的数)。把误差拆成「随机波动」加「系统偏差」:
  2. 平方展开:
  3. 逐项取期望:第一项是 第二项是常数 乘 第三项是常数

所以均方误差 = 方差 + 偏差的平方。对本页的 偏差 均方误差就是方差。

常见疑问与易错点

易错点问:为什么误差是 σ/√n,而不是 σ/n?

答:因为独立误差相加时会部分抵消。 个中心化样本之和 的均方是 (第 7 步前),典型大小是 而不是 有的偏正、有的偏负,互相抵掉一部分。再除以 得到 只有在例 4 那样「完全不抵消」的情形,和才会是 倍大,平均也就一点不改善。

易错点问:命题一定要「相互独立」吗?

答:证明里只用到交叉项 (这叫两两不相关。两两独立就足以保证它(第 6 步只涉及两个变量)。完全相关时(例 4)方差不降;反过来,负相关还能让方差降得更快。例:估计 用一对 和 平均值 恰好是真值,方差为 0。这是「对偶变量」技巧的雏形。

易错点问:均方误差很小,是不是每次都很准?

答:不是。均方误差是平均意义上的;Chebyshev 只说「偏得很远」的概率小,不说它不会发生。例 2 里 时误差超过 0.1 的概率约 3.5%:做 100 次这样的试验,大约有三四次会偏出 0.1。

易错点问:实际计算时不知道 σ²,怎么报告误差?

答:用同一批样本估计它:样本方差 再报告「标准误约 要当心:当真实方差是无穷(第 36 页,时,每次算出的 仍是一个有限的数,看起来很正常,但它会随 增大而不断变大、极不稳定,据此报告的误差是假的。所以要先从理论上确认方差有限,这正是 第 39 页 要做的事。

小结与自测

小结本页小结
  • 命题: 独立同分布、 时,并且均方误差
  • 证明:中心化后展开平方; 个平方项各贡献 交叉项因为独立(乘积公式)期望为 0。
  • Chebyshev:推出 依概率收敛(弱大数定律)。
  • 误差尺度 :样本数乘 100,误差缩小 10 倍; 减半等于样本数乘 4。
  • 第 36 页 中 时 均方误差对每个 都是 这个速率不存在;要换一种抽样方式(第 38 页)。
自测

1. 某抽样方案单个样本的方差 要让均方根误差不超过 至少需要多少个独立样本?

点开看答案

均方根误差是 要求 即

2. 方案 A 单个样本方差 方案 B 单个样本方差 但每个样本的计算时间是 A 的 3 倍。给定同样的总计算时间,哪个方案的均方误差更小?

点开看答案

设总时间能让 A 抽 个样本,则 B 只能抽 个。均方误差分别是 和 B 更小。一般规律:比较「方差 × 每个样本的成本」,A 是 B 是 越小越好。

3. 判断:若 (同一个随机变量复制 份),则

点开看答案

错。此时 命题需要独立(至少两两不相关):这里每个交叉项 证明第 6 步不成立(例 4)。

本页术语中英对照(按出现顺序)
中文English白话解释
独立同分布independent and identically distributed (i.i.d.)一列随机变量相互独立,并且分布都相同,例如反复独立地掷同一颗骰子。
方差varianceVar X = E[(X − E X)²],衡量随机变量围绕均值的波动大小。
随机变量random variable从样本空间到实数的(可测)函数,把结局翻译成数字;随机性来自「实现了哪个 ω」。
平移不变性translation invariance把集合整体平移,测度不变:μ(E + t) = μ(E)。长度就有这个性质。
独立independence任意有限个事件同时发生的概率等于各自概率的乘积。
博雷尔 σ-代数Borel sigma-algebra由实数轴上所有开区间生成的 σ-代数 𝓑(ℝ),包含所有常见的集合。
函数function / mapf: A → B 给 A 中每个元素指定 B 中唯一一个元素。
可测函数measurable function对每个 a,集合 {f > a} 都可测;这样的函数才能积分。
可积integrable (L¹)∫|f| dμ < ∞。全体可积函数记作 L¹(μ)。
单调收敛定理monotone convergence theorem (MCT)非负可测函数单调递增地收敛时,积分也收敛:∫fₙ ↑ ∫f。
示性函数indicator function𝟏_A(x) 在 x ∈ A 时取 1,否则取 0;它的积分就是 A 的测度。
简单函数simple function只取有限个值的可测函数,是若干示性函数的线性组合,图像像台阶。
Markov 不等式Markov's inequalityX ≥ 0 时 P(X ≥ t) ≤ E X / t:期望小,取大值的概率就小。
期望expectation随机变量对概率测度的积分 E[X] = ∫ X dP,也就是「加权平均」。
Chebyshev 不等式Chebyshev's inequalityP(|X − μ| ≥ t) ≤ σ² / t²:方差控制偏离均值的概率。
概率空间probability space三元组 (Ω, 𝓕, P):样本空间、事件的 σ-代数、概率测度。
分布law / distributionX 把概率「推」到取值空间上得到的测度:L(X)(B) = P(X ∈ B)。
依概率收敛convergence in probability对每个 ε > 0,P(|Xₙ − X| > ε) → 0;比几乎必然收敛弱。
无意识统计学家定律law of the unconscious statistician (LOTUS)求 g(X) 的期望,不必先求出 g(X) 的分布,直接把 g 对 X 的分布 𝓛(X) 积分;X 有密度 f 时就是 E g(X) = ∫ g(x) f(x) dx。
弱大数定律weak law of large numbers (WLLN)样本均值依概率收敛到真实均值:对每个 ε > 0,P(|Sₙ/n − μ| > ε) → 0。方差有限时可由 Chebyshev 不等式直接推出。
两两独立pairwise independence任意两个都独立。比相互独立弱:三个事件可以两两独立却不相互独立。
中心极限定理central limit theorem (CLT)独立同分布、方差 σ² 有限时,(Sₙ − nμ)/(σ√n) 依分布收敛到标准正态分布:大量独立小误差之和近似正态。
强大数定律strong law of large numbers (SLLN)独立同分布且 E|X₁| < ∞ 时,样本均值 Sₙ/n 几乎必然收敛到 E X₁。
几乎必然收敛almost sure convergence除了概率为 0 的结局以外,Xₙ(ω) 都收敛到 X(ω)。
正部positive / negative partf⁺ = max(f, 0),f⁻ = max(−f, 0),f = f⁺ − f⁻,|f| = f⁺ + f⁻。