这一页在讲什么
把 个独立同分布、方差为 的样本求平均,平均值偏离真值的均方误差恰好是 所以误差的典型大小是 样本数乘 4,误差才减半。单个样本的方差 就是一个抽样方案的「单价」。第 36 页 里当 时 这把尺子就失效了。
- 这是第 5 部分「抽样中的应用」的第二页。上一页 第 36 页 发现:均匀抽样估计量无偏,但方差可能无穷。这一页回答「为什么要在乎方差」:蒙特卡罗方法的精度完全由方差决定。
- 课件只印了一个小标题 Mean squared error of the sample average(样本均值的均方误差),内容是板书。下面「补全课件留白」一节补全:一个命题(均方误差 它的证明、由它推出的 Chebyshev 概率界,以及方差无穷时会发生什么。
- 下一页起(第 38 页–第 40 页)的目标由此变得很清楚:换一种抽样方式,让单个样本的方差变成有限,而且尽量小。
先补基础
- 设 是用来估计未知数 的估计量(它依赖随机样本,是随机变量)。它的均方误差(mean squared error,MSE)定义为 误差先平方(去掉正负号,并且大误差罚得更重),再求平均。
- 均方根误差 把平方开回来,单位和 相同,可以直接读成「典型误差有多大」。
- 分解:若 则 即「均方误差 = 方差 + 偏差的平方」(证明见「推导」一节第(二)条)。无偏时偏差为 0,均方误差就等于方差。
- 例:一个永远报告 的「估计量」方差为 0,但均方误差是 1——它稳定地错。另一个一半时候报 一半时候报 偏差为 0,方差和均方误差都是 1。
课件原文逐句精读
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|样本均值的均方误差
先看每个条件管什么:
- :保证 是有限的数,结论才有内容;它也保证 是有限的数。第 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 条开平方,得到均方根误差:
它告诉我们典型误差有多大:
| 样本数 | ||||
|---|---|---|---|---|
| 均方根误差 |
读这张表:
- 样本数乘 100,误差只缩小到十分之一。想多得一位准确的小数,要付出 100 倍的计算量。这就是人们常说的「蒙特卡罗方法收敛慢」。
- 「几个标准误」的写法:在第 3 条里取 得 叫标准误(standard error)。例如 误差超过 3 个标准误的概率不超过 不超过
- 需要多少样本:要让 由第 3 条,只要 即 所需样本数与 成正比。
- 两条降低误差的路:(1)加大 代价与 成正比;(2)换一种抽样方式,让单个样本的方差 变小。把 减半,效果等于把样本数乘 4。第 38、39 页走的是第二条路。
Chebyshev 界对每个 都严格成立,但比较保守。中心极限定理(本讲没有讲,这里只作背景)说: 很大时, 的分布接近标准正态 于是大约 95% 的情形误差不超过 两种说法的共同点是:误差的尺度是 。
方差无穷时:σ/√n 这把尺子失效
回到 第 36 页: 独立同分布、服从
- :命题完全适用,误差按 缩小。例如 时 (第 36 页),
- :命题的前提不成立。第 36 页「推导」一节第(二)条已经证明:对每个 ,平均多少个样本都一样;Chebyshev 界的右边是 什么也没说。
- 估计量还会收敛吗? 会。强大数定律只要求 (本讲不证), 几乎必然。但没有 的速率保证,实际收敛更慢、更不稳定。下面例 3 用模拟展示这一点。
例子
是第 次掷一颗公平骰子的点数,各次独立,
- 方差:所以
- :均方误差 均方根误差 Chebyshev:
- :均方根误差 Chebyshev:
精度提高 10 倍(容忍度从 0.5 缩到 0.05),同样的概率保证需要 100 倍的样本。
是独立的公平硬币(正面记 1),(因为
- 均方根误差
- Chebyshev:
- 精确值: 就是「正面次数 或 按二项分布精确计算,概率约为
Chebyshev 给的 0.25 比真实的 0.035 大很多,但它不需要知道分布的具体形状,只用到方差,并且对每个 都成立。
对 第 36 页 的均匀抽样估计量 每种 独立重复 2000 次,记录每次的误差 取中位数(一半的重复误差比它小,一半比它大)。下面是一次这样的模拟的结果(随机数不同,数字会略有出入):
| (方差有限)的中位误差 | 0.053 | 0.017 | 0.0054 |
| (方差无穷)的中位误差 | 0.61 | 0.32 | 0.16 |
- 每乘 10,中位误差缩小约 3.1 倍,正是 符合
- 每乘 10,中位误差只缩小约 2 倍。理论上可以证明(要用本课没有的「稳定分布」理论,这里只作背景),这个例子里典型误差大约按 缩小; 时是 而 与模拟一致。
- 同一组模拟里, 时约 70% 的重复低估了 (大多数时候抽不到 0 附近的巨大值),而最坏的一次误差超过 100(碰巧抽到了一个离 0 极近的点)。平均来说无偏,是靠极少数巨大的高估撑起来的。
设 满足 令 (相当于随机数发生器出了故障,把同一个数重复了 次)。
- 每个 的分布相同,但它们不独立。
- 所以 不管 多大都不减小。
- 对照证明:第 6 步的交叉项变成 不再是 0。 个交叉项合计 加上平方项 总共 除以 仍是
所以命题里的「独立」不能去掉。
推导:一步一步来
补上正文用到的两件事:(一)独立可积变量的乘积公式,走「示性函数 → 简单函数 → 非负函数 → 可积函数」四步;(二)均方误差 = 方差 + 偏差的平方。
(一)若 独立且可积,则 可积,并且
- 非负简单函数。 设 是 上的非负简单博雷尔函数(有限和, 是博雷尔集、 也是)。令 因为 这是 上的非负简单函数。由简单函数积分的线性(第 32 页)和独立性,
- 非负的 设 且 取二进逼近 (第 34 页),它是非负简单博雷尔函数,并且对每个 有 令 则 并且 (两个非负递增数列的乘积仍递增,极限等于极限之积)。对三列分别用单调收敛定理(第 36 页),再用第 1 步:最后一个等号: 都是有限数的收敛,乘积的极限等于极限的乘积。
- 可积的 分别是 的函数,独立且非负,由第 2 步 所以 可积。再把正部、负部(第 43 页) 代入:四项都是「 的函数 × 的函数」,非负,并且都 所以可积。由期望的线性和第 2 步,
(二)若 则
- 记 (有限的数)。把误差拆成「随机波动」加「系统偏差」:
- 平方展开:
- 逐项取期望:第一项是 第二项是常数 乘 第三项是常数
所以均方误差 = 方差 + 偏差的平方。对本页的 偏差 均方误差就是方差。
常见疑问与易错点
答:因为独立误差相加时会部分抵消。 个中心化样本之和 的均方是 (第 7 步前),典型大小是 而不是 有的偏正、有的偏负,互相抵掉一部分。再除以 得到 只有在例 4 那样「完全不抵消」的情形,和才会是 倍大,平均也就一点不改善。
答:证明里只用到交叉项 (这叫两两不相关。两两独立就足以保证它(第 6 步只涉及两个变量)。完全相关时(例 4)方差不降;反过来,负相关还能让方差降得更快。例:估计 用一对 和 平均值 恰好是真值,方差为 0。这是「对偶变量」技巧的雏形。
答:不是。均方误差是平均意义上的;Chebyshev 只说「偏得很远」的概率小,不说它不会发生。例 2 里 时误差超过 0.1 的概率约 3.5%:做 100 次这样的试验,大约有三四次会偏出 0.1。
小结与自测
1. 某抽样方案单个样本的方差 要让均方根误差不超过 至少需要多少个独立样本?
点开看答案
均方根误差是 要求 即
2. 方案 A 单个样本方差 方案 B 单个样本方差 但每个样本的计算时间是 A 的 3 倍。给定同样的总计算时间,哪个方案的均方误差更小?
点开看答案
设总时间能让 A 抽 个样本,则 B 只能抽 个。均方误差分别是 和 B 更小。一般规律:比较「方差 × 每个样本的成本」,A 是 B 是 越小越好。
3. 判断:若 (同一个随机变量复制 份),则
点开看答案
错。此时 命题需要独立(至少两两不相关):这里每个交叉项 证明第 6 步不成立(例 4)。
| 中文 | English | 白话解释 |
|---|---|---|
| 独立同分布 | independent and identically distributed (i.i.d.) | 一列随机变量相互独立,并且分布都相同,例如反复独立地掷同一颗骰子。 |
| 方差 | variance | Var X = E[(X − E X)²],衡量随机变量围绕均值的波动大小。 |
| 随机变量 | random variable | 从样本空间到实数的(可测)函数,把结局翻译成数字;随机性来自「实现了哪个 ω」。 |
| 平移不变性 | translation invariance | 把集合整体平移,测度不变:μ(E + t) = μ(E)。长度就有这个性质。 |
| 独立 | independence | 任意有限个事件同时发生的概率等于各自概率的乘积。 |
| 博雷尔 σ-代数 | Borel sigma-algebra | 由实数轴上所有开区间生成的 σ-代数 𝓑(ℝ),包含所有常见的集合。 |
| 函数 | function / map | f: 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 inequality | X ≥ 0 时 P(X ≥ t) ≤ E X / t:期望小,取大值的概率就小。 |
| 期望 | expectation | 随机变量对概率测度的积分 E[X] = ∫ X dP,也就是「加权平均」。 |
| Chebyshev 不等式 | Chebyshev's inequality | P(|X − μ| ≥ t) ≤ σ² / t²:方差控制偏离均值的概率。 |
| 概率空间 | probability space | 三元组 (Ω, 𝓕, P):样本空间、事件的 σ-代数、概率测度。 |
| 分布 | law / distribution | X 把概率「推」到取值空间上得到的测度: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 part | f⁺ = max(f, 0),f⁻ = max(−f, 0),f = f⁺ − f⁻,|f| = f⁺ + f⁻。 |