这一页在讲什么
被积函数 在 时趋于无穷(奇点),但曲线下的面积 有限(可积)。最朴素的算法——在 上均匀随机撒点、对 的值求平均——是无偏的:平均来说恰好等于 可是当 时,单个样本 的方差是无穷大:偶尔抽到离 0 很近的点, 的值大得离谱,把平均值拽走。
- 这是第 5 部分「抽样中的应用」的第一页,也是整个第 5 部分要解决的问题。
- 下一页 第 37 页 解释为什么方差这么重要:样本均值的误差大小由 决定,方差无穷时这个「误差尺度」就不存在了。
- 第 38–40 页用第 4 部分的换测度公式(第 35 页)给出解决办法:重要性抽样,把方差变成有限、甚至有界。
先补基础
- 很多积分算不出精确的公式(本页的 就没有初等表达式)。蒙特卡罗方法(Monte Carlo method)的想法是:把积分写成一个期望,再用「随机抽样 + 求平均」去估计它。
- 关键等式:若 服从 上的均匀分布,则 所以「在 上撒很多随机点、算 的平均值」就能近似
- 例:想估计 随机抽到 平均 已经接近 点越多,一般越准。
- 读作「 服从 上的均匀分布」:(也就是 的分布是 上的长度
- 无意识统计学家定律(第 7 页):对非负可测的 期望可以直接对 的分布积分,不必先求 的分布。
- 估计量(estimator):用样本算出来、用来猜一个未知数的量。例如 它依赖于随机样本,所以本身是随机变量。
- 无偏(unbiased):即「重复做很多次试验,估计值的平均恰好是真值」。
- 方差(variance):衡量 围绕均值波动的大小;二阶矩有限时 方差可以是
- 独立同分布(iid): 各自都是 并且互不影响,相当于「把同一个抽样独立地重复做很多次」。正式定义在下一页。
课件原文逐句精读
Example 23 (An integrable singularity) Let , , and .
The uniform-sampling estimator
When the variance is infinite
例 23(可积的奇点) 设
均匀抽样估计量
方差什么时候是无穷
Example 23 (An integrable singularity)|例 23(可积的奇点)
- singularity(奇点):函数在某一点附近「炸掉」的地方。这里 时 而 所以
- integrable(可积):虽然 无界,曲线下的面积仍然有限。理由见下一节:而
- 条件 的两端各管一件事: 让 真的有奇点( 时 是有界的光滑函数); 让奇点可积( 时
- 变量写作 暗示它是「均匀随机数」(uniform)的取值。第 38 页把同一个积分放进测度论的记号里,改用
两个小标题
课件接下来印了两个小标题,内容是板书,下一节补全:
- The uniform-sampling estimator:用均匀随机数构造一个估计 的方法,并说明它无偏;
- When the variance is infinite:计算这个估计量的方差,找出方差为无穷的条件。蓝条标题 An unbiased estimator can have infinite variance(无偏估计量的方差可以是无穷)就是这一部分的结论。
补全课件留白
课件此处留白(课堂板书),以下内容由本站根据标准教材补写,记号与本课件保持一致;如与老师课堂写法不同,以老师为准。
Example 23 的设定:I 是一个有限的正数
- 上界: 时 所以 由积分的单调性和上面的基础补课(
- 下界: 时 所以
- 数值: 没有初等函数的表达式(它是「下不完全伽马函数」用数值积分算得: 时 时 时 确实落在 里。
The uniform-sampling estimator|均匀抽样估计量
When the variance is infinite|什么时候方差是无穷
第 1 步:写出二阶矩。 非负,由无意识统计学家定律,
第 2 步:把 夹住。 时 所以
这个因子不影响「有限还是无穷」,只有 说了算。
第 3 步:用 的判别法。
- :方差有限:例如 时,数值积分得
- :而均值 有限,由下面「推导」一节第(一)条,
第 4 步:对估计量的影响。 当 时,不但单个样本的方差无穷,无论 多大, 都是无穷(「推导」一节第(二)条)。多抽样本也救不回来。
结论:均匀抽样估计量对所有 都无偏,但当 时它的方差无穷。这就是蓝条标题的意思。
为什么「平方」会让奇点变得不可积? 的奇点强度是 平方之后变成 强度翻倍。可积的界线在指数 1: 保证 可积,但只有 才保证 可积。 正好落在「 可积、 不可积」的夹缝里。
在坐标纸上画 时的曲线 横轴从 0 到 1,右端 往左慢慢升高,到 时约为 4.5,再往左就像一面墙一样竖直冲上去。
- 曲线下的面积(有限:墙虽然无限高,但越高的地方越窄,窄得足够快,面积加起来只有约 2.74。
- 均匀撒点:横坐标均匀地落在 上,绝大多数点落在右边平缓的地方,返回 0.4~5 这样的小数;极少数点落在墙脚下,返回成百上千的数。
- 平方之后:把曲线的高度平方,墙变得更高、而宽度不变,窄得「不够快」了,面积变成无穷——这就是二阶矩无穷。估计量的麻烦全部来自那面墙。第 38 页的办法是:让撒点本身也往墙脚下集中,再用权重把多撒的点折算回来。
例子
(没有 因子,
- 有限, 的误差会按正常速度缩小(下一页)。
- 无偏估计 没有问题。
- (对数发散:切掉 后积分是 随 增大没有上界)。
- 所以 已经是方差无穷的情形,条件是 包含等号。
(相当于 没有有限的目标值可估计,「无偏」也无从谈起。这说明例 23 里 这个条件不能去掉。
先看几个具体的样本值(真值
| 抽到的 | |||||
|---|---|---|---|---|---|
再看「 超过 的概率( 在 上递减,所以 等价于 小于某个门槛,这个门槛就是概率):
| 大约每多少次出现一次 | 概率 | 概率 | ||
|---|---|---|---|---|
| 10 | 0.0354 | 28 | 0.354 | 3.54 |
| 100 | 0.00139 | 721 | 0.139 | 13.9 |
| 1000 | 19308 | 0.052 | 51.8 | |
| 10000 | 517949 | 0.019 | 193 |
读表:
- 每乘 10,概率大约除以
- 「概率」越来越小:极大值对均值的贡献越来越小,所以 有限。
- 「概率」越来越大、没有上界:极大值对二阶矩的贡献越来越大,所以
- 实际后果:抽 1000 个点时,很可能碰到一个超过 100 的值,它单独就能把平均值抬高 抽的点越多,越可能碰到更大的值。误差不会像方差有限时那样稳定地缩小。
推导:一步一步来
补上正文用到的两个事实:(一)均值有限而二阶矩无穷,则方差无穷;(二)对非负的独立样本,平均多少个也不能让均方误差变成有限。两者都只用到 这个初等不等式。
(一)若 则
- 记 (有限)。初等不等式 (由 展开得到)用于
- 移项:
- 两边取期望(左边非负,期望有意义):
(二)若 则 其中
- 所有 所以 (和不小于其中一项,两边都非负,平方保持大小关系)。
- 于是 取期望得
- 有限,把(一)用于 得
把(二)用于 就是正文第 4 步: 时,对每个 注意这里连独立性都没用到,只用到了「非负」。
常见疑问与易错点
答:不够。无偏只说「平均起来是对的」,不说「每一次离真值多远」。一个估计量可以一半时候给 0、一半时候给 它也是无偏的,但每次都错一个 衡量「每次大概错多少」要看方差(下一页的均方误差)。方差无穷时,单次运行的结果可能被一个极端样本完全左右。
答:因为平方把奇点放大了: 可积要求 可积要求 例 4 的表格把这件事量化了:极大值对均值的贡献越来越小,对平方的贡献却越来越大。
答:会。强大数定律只要求 就保证 几乎必然(这个定理本讲没有证明,这里只作为背景)。问题在于速度:方差有限时误差大约按 缩小(第 37 页),方差无穷时没有这个保证,收敛又慢又不稳定。所以本部分关心的不是「对不对」,而是「多快、多稳」。
答: 时 而 (例 2)。可积的界线「指数 是严格不等号,指数恰好等于 1 已经发散。
小结与自测
- 例 23:(在 0 处有奇点,但可积:
- 均匀抽样估计量 由无意识统计学家定律和线性性无偏。
- 方差:与 同时有限或无穷,所以
- 方差无穷时,多少个样本求平均,均方误差都是无穷;需要换一种抽样方法(第 38 页)。
1. 时, 的方差有限吗?给出 的一个上界。
点开看答案
有限。由第 2 步的上界,
2. 对 求 并判断 是否有限。
点开看答案
(指数 均值有限,由「推导」(一),方差无穷。
3. 为什么在判断 是否有限时,可以不管 这个因子?
点开看答案
因为在 上 这个因子被两个正的常数夹住。乘上一个夹在两个正常数之间的因子,积分只会变大或变小一个固定的倍数,不会把有限变成无穷,也不会把无穷变成有限。
| 中文 | English | 白话解释 |
|---|---|---|
| 函数 | function / map | f: A → B 给 A 中每个元素指定 B 中唯一一个元素。 |
| 可积 | integrable (L¹) | ∫|f| dμ < ∞。全体可积函数记作 L¹(μ)。 |
| 方差 | variance | Var X = E[(X − E X)²],衡量随机变量围绕均值的波动大小。 |
| 换测度 | change of measure | 用 E_Q[X] = E_P[X · dQ/dP] 在两个概率之间换算期望。 |
| 重要性抽样 | importance sampling | 从另一个更合适的分布抽样,再乘以权重纠正,从而减小估计的方差。 |
| 期望 | expectation | 随机变量对概率测度的积分 E[X] = ∫ X dP,也就是「加权平均」。 |
| 分布 | law / distribution | X 把概率「推」到取值空间上得到的测度:L(X)(B) = P(X ∈ B)。 |
| 无意识统计学家定律 | law of the unconscious statistician (LOTUS) | 求 g(X) 的期望,不必先求出 g(X) 的分布,直接把 g 对 X 的分布 𝓛(X) 积分;X 有密度 f 时就是 E g(X) = ∫ g(x) f(x) dx。 |
| 随机变量 | random variable | 从样本空间到实数的(可测)函数,把结局翻译成数字;随机性来自「实现了哪个 ω」。 |
| 阶矩 | moment | E[|X|ᵖ] 叫 p 阶(绝对)矩;它有限与否反映尾部有多「重」。 |
| 独立同分布 | independent and identically distributed (i.i.d.) | 一列随机变量相互独立,并且分布都相同,例如反复独立地掷同一颗骰子。 |
| 独立 | independence | 任意有限个事件同时发生的概率等于各自概率的乘积。 |
| 勒贝格积分 | Lebesgue integral | 按函数值分块的积分:先定义简单函数的积分,再用下方的简单函数逼近取上确界。 |
| 黎曼积分 | Riemann integral | 大一学的定积分:按 x 轴切竖条求和取极限。 |
| 单调收敛定理 | monotone convergence theorem (MCT) | 非负可测函数单调递增地收敛时,积分也收敛:∫fₙ ↑ ∫f。 |
| 测度 | measure | 给 σ-代数里的每个集合一个 [0, ∞] 中的「大小」,空集为 0,并且可数可加。 |
| 上界 | upper bound | 数集 A 的上界是一个数 M,A 里每个数都 ≤ M。例如 2 和 5 都是 [0,1] 的上界;所有上界中最小的那个叫上确界。 |
| 下界 | lower bound | 数集 A 的下界是一个数 m,A 里每个数都 ≥ m。例如 0 和 −3 都是 (0,1] 的下界;所有下界中最大的那个叫下确界。 |
| 强大数定律 | strong law of large numbers (SLLN) | 独立同分布且 E|X₁| < ∞ 时,样本均值 Sₙ/n 几乎必然收敛到 E X₁。 |
| 几乎必然收敛 | almost sure convergence | 除了概率为 0 的结局以外,Xₙ(ω) 都收敛到 X(ω)。 |