高等概率论逐页精讲

作业 2 收敛方式、大数定律、中心极限定理与鞅›大数定律与中心极限定理

第 6 题 / 共 12 题解答原页 第 6–7 页证明题

Glivenko–Cantelli 定理:经验分布函数一致收敛到真实分布

The Glivenko–Cantelli theorem

解答原页第 6 页(本题所在部分)
解答原页第 7 页(本题所在部分)
解答原页第 6–7 页(本题部分)· 点击图片看整页

这道题在问什么

人话版一句话版

从分布 中独立抽 个样本,统计「不超过 的样本所占的比例」,得到经验分布函数 强大数定律说:对每个固定的 几乎必然。Glivenko–Cantelli 定理更强:所有 一起, 几乎必然,也就是整条阶梯曲线一致地贴近真实的分布函数。它被称为「统计学基本定理」:样本足够多时,数据能还原出整个分布。证明的关键是把无穷多个 化成有限个分位点,再用分布函数的单调性夹逼。

  • 考点:强大数定律(讲义之外,作为已知结论)、分布函数的单调性与右连续性、左极限 分位点。
  • 出处:Durrett《Probability: Theory and Examples》定理 2.4.9,原解答几乎逐字照搬了它的证明。
  • 和讲义的关系:分布函数(第 2 讲 第 8 页)、分位数与逆变换(第 2 讲 第 14 页)、几乎必然收敛(第 3 讲 第 9 页)。

先补基础

基础补课分布函数的三条性质与左极限

满足(第 2 讲 第 8 页):单调不减;右连续;

左极限 在 处有跳跃时,跳跃的高度是 例:公平骰子在 处

基础补课经验分布函数

样本中不超过 的比例。它是一个阶梯函数,在每个样本点处向上跳 例:样本 (则 时 时 时 时 同样定义

基础补课强大数定律(这里作为已知结论使用)

i.i.d.,则 几乎必然。(Durrett 定理 2.4.1。)本题只用到 是示性函数(取 0、1)的情形,此时它也叫 Borel 强大数定律,可以用 Borel–Cantelli 引理较直接地证明。

基础补课可数个概率为 1 的事件,交仍然概率为 1

若 (则 因为 本题要让可数多个点上的收敛同时成立,靠的就是这一条。

题目与原解答

课件原文题目
  1. We define the empirical distribution functions as the following. Let be i.i.d. with distribution and let

= the observed frequency of values that are , hence the name given above.

Prove the Glivenko-Cantelli theorem. As ,

(Hence, converges uniformly to as .)

中文翻译
  1. 定义经验分布函数如下:设 i.i.d.,分布函数为 令 即观测值中不超过 的频率,因此得名。

证明 Glivenko–Cantelli 定理: 时 几乎必然(于是 一致收敛到

  • 题目中求和写成 上限 和下标 都不一致,应为

原解答的思路(左栏,即 Durrett 定理 2.4.9 的证明):

  1. 固定 由强大数定律得 几乎必然;
  2. 取分位点 在这有限个点上同时让误差小于
  3. 对两个分位点之间的 用 的单调性夹逼,得到

证明是正确的,但跳过了几处需要说明的细节(见「原解答逐行对照」)。

思路:怎么想到的

  • 难点:「对每个 收敛」和「对所有 一致收敛」差别很大,因为有不可数多个 不能简单地把所有 的零测例外集并起来(不可数个零测集的并可以不是零测集)。
  • 办法:只看有限个精心选择的点。在两个相邻的点之间, 和 都单调,它们的值被夹在端点的值之间;只要相邻端点处 的增长不超过 中间的误差就被端点的误差控制住。
  • 跳跃的麻烦: 若有跳跃,一个跳跃可能比 还高,相邻分位点之间的增长就控制不住。解决办法:在分位点处同时考虑 和左极限 跳跃发生在这两者之间,不影响夹逼。

完整解答

第 1 步(逐点收敛,包括左极限):固定 (i.i.d.,均值 由强大数定律, 几乎必然。同理 的均值是 所以 几乎必然。

第 2 步(分位点):固定正整数 对 令

(第 分位点,第 2 讲 第 14 页),并令 由于 这些 都是有限实数,并且

分位点的两个性质(

  • 对 由下确界的定义 取极限即得。( 时右边是 1,自然成立。)
  • 由右连续性,下确界处取得 ( 时左边按约定为

两式相减:

第 3 步(可数个点同时收敛):由第 1 步,对每个 和每个 和 都几乎必然成立。所有 只有可数多对,由「可数个概率为 1 的事件之交仍概率为 1」,存在一个概率为 1 的事件 在 上它们同时成立。下面固定 对每个 只有有限个点,所以存在 使 时对所有

( 时两边都是 0 或 1,自动成立。)

第 4 步(夹逼):设 任取实数

  • 若 恰好是某个 :由第 3 步,
  • 否则 对某个 由单调性和

(第一行第一个不等号:所以 最后一个: 单调。第二行:第一个不等号用 和 单调;最后一个用 所以

所以对所有 即 (

第 5 步(结论): 任意,所以在 上 对所有 成立,即 定理得证。

原解答逐行对照

  1. 「Fix and let ... the strong law of large numbers implies that a.s.」:意思是「固定 令 由强大数定律, 几乎必然收敛到 求和上限 应为 正确。
  2. 「In general, if is a sequence of nondecreasing functions that converges pointwise to a bounded and continuous limit then . However, the distribution function may have jumps」:意思是「一般地,若一列单调不减的函数 逐点收敛到一个有界连续的极限 则 可是分布函数 可能有跳跃」。这句话说明了为什么要处理左极限: 连续时,只看 就够了(这个结论叫 Pólya 定理);有跳跃时还要看
  3. 「Again, fix and let ... a.s.」:下标混用(意思正确。
  4. 「For let . The pointwise convergence ... imply that we can pick ...」:意思是「对 令 (在这些点上的)逐点收敛……说明可以选出 这里需要一个公共的概率为 1 的事件,使所有 和 的收敛同时成立(第 3 步)。原解答没有说明。
  5. 「If ... using ... :两处遗漏:① 不等式 没有证明(它依赖分位点的两个性质:下确界给出 右连续性给出 ② 只处理了开区间, 恰好等于分位点的情形没有说(它由第 3 步的逐点估计直接覆盖)。
  6. 两行夹逼不等式与结论 :正确。
易错点原解答的问题:「同时成立」需要一个公共的概率 1 事件

对每个固定的点,收敛都是几乎必然的,但例外的零测集依赖于这个点。证明里要在无穷多个 每个 下的 个点上同时使用收敛,必须先说明:这些点总共只有可数个,所以把可数个例外集并起来仍是零测集。原解答直接说「pick 省略了这一步。这一步正是「把不可数个 化成可数个点」这个思路的意义所在。

易错点原解答的问题:关键不等式 (∗) 没有证明

是整个夹逼的基础,也是处理跳跃的关键:跳跃的高度可能远大于 但它发生在 和 之间,不在两个相邻分位点之间。原解答直接使用,没有给出理由(见第 2 步)。

例子与数值验证

例 1:一致误差随样本量的变化

从 上的均匀分布抽样,计算 的平均值(Python 模拟,固定随机种子;四列分别重复 4 万、2 万、4000、800 次,误差都在最后一位数字的 ±1 以内):

10 100 1000 10000
的平均值(模拟)

趋于 0,而且大约按 的速度:样本量乘 100,误差缩小到十分之一。

例 2:有跳跃的分布——掷骰子

掷 60 次公平骰子(固定随机种子),得到 在 处的值依次是 0.183、0.350、0.483、0.633、0.833、1,真实值 最大误差 骰子的分布函数每处跳 比证明中 较大时的 高得多;证明之所以仍然有效,是因为跳跃夹在 与 之间,两者都被单独控制了。

例 3:逐点收敛不等于一致收敛(非概率的例子)

函数列 对每个 可 对所有 成立。这说明「每个 都收敛」推不出一致收敛。Glivenko–Cantelli 之所以成立,是因为 和 都单调有界,单调性把有限个点的信息「传播」到了所有

拓展延伸

DKW 不等式:收敛速度的定量版本

Dvoretzky–Kiefer–Wolfowitz 不等式(Massart 给出了最优常数):对任意

例如 时,右边约为 这比 Glivenko–Cantelli 更强:它给出了任意有限 下的概率上界,而且与 是什么分布无关。由它和 Borel–Cantelli 引理也能推出 Glivenko–Cantelli 定理。

Kolmogorov–Smirnov 检验

更精细地, 依分布收敛到一个与 无关的极限分布(Kolmogorov 分布, 连续时),例 1 里的 就是这个极限分布的均值 统计学中用它检验「数据是否来自某个指定的分布」:若 太大,就拒绝这个分布。

经验过程与 bootstrap

把 看成一个随机函数,它收敛到「布朗桥」(Donsker 定理)。这是经验过程理论的起点。统计中的 bootstrap 方法,正是把未知的 换成经验分布 去重复抽样;Glivenko–Cantelli 定理保证了 是 的好替身。

常见疑问与易错点

易错点问:为什么不能对每个 x 用强大数定律,然后取 sup?

答:对每个 例外的零测集 不同; 有不可数多个, 可能不是零测集,甚至可能不是事件。所以只能在可数多个点上使用强大数定律,再用单调性推广到所有

易错点问:F 连续时证明能简化吗?

答:能。连续时 只需在 处用 并取 使 这时就是原解答提到的 Pólya 定理。

易错点问:一致收敛有什么实际意义?

答:它保证用 估计任何区间的概率都同时准确:如果只有逐点收敛,就不能保证所有区间同时可靠。

小结与自测

小结本题小结
  • 逐点:强大数定律给出 几乎必然。
  • 分位点 满足 (下确界 + 右连续)。
  • 可数个点同时收敛(一个公共的概率 1 事件)+ 单调性夹逼 ⇒
  • 原解答:下标笔误;没有说明公共的零测集、没有证明 没有处理分位点本身。
  • 拓展:误差约为 DKW 不等式
自测

1. 样本为 (写出 和

点开看答案

( 都 (只有

2. 对 上的均匀分布,取 写出分位点

点开看答案

所以

3. 用 DKW 不等式估计:要使 大约需要多少样本?

点开看答案

需要 即

本页术语中英对照(按出现顺序)
中文English白话解释
分布law / distributionX 把概率「推」到取值空间上得到的测度:L(X)(B) = P(X ∈ B)。
独立independence任意有限个事件同时发生的概率等于各自概率的乘积。
经验分布函数empirical distribution functionFₙ(x) = 样本 X₁, …, Xₙ 中不超过 x 的个数 ÷ n,即样本里 ≤ x 的比例;它是真实分布函数 F 的估计。
强大数定律strong law of large numbers (SLLN)独立同分布且 E|X₁| < ∞ 时,样本均值 Sₙ/n 几乎必然收敛到 E X₁。
Glivenko–Cantelli 定理Glivenko–Cantelli theorem独立同分布样本的经验分布函数几乎必然一致收敛到真实分布函数:supₓ |Fₙ(x) − F(x)| → 0。
几乎必然收敛almost sure convergence除了概率为 0 的结局以外,Xₙ(ω) 都收敛到 X(ω)。
累积分布函数cumulative distribution function (CDF)F(x) = P(X ≤ x),单调不减、右连续,完全决定分布。
函数function / mapf: A → B 给 A 中每个元素指定 B 中唯一一个元素。
左极限left limit自变量从左边(比 x 小的一侧)靠近 x 时,函数值趋向的数,记作 f(x−)。分布函数在 x 处的跳跃 F(x) − F(x−) 就是 P(X = x)。
分位数quantileQ(u) = inf{x : F(x) ≥ u},是 CDF 的「广义反函数」。
单调递增nondecreasing后一项不小于前一项;单调有界的数列一定有极限。
样本点sample point一个完整的结局 ω;每次试验恰好实现一个样本点。
独立同分布independent and identically distributed (i.i.d.)一列随机变量相互独立,并且分布都相同,例如反复独立地掷同一颗骰子。
Borel–Cantelli 引理Borel–Cantelli lemma概率之和有限 ⇒ 只发生有限次;独立且概率之和无穷 ⇒ 发生无穷多次。
示性函数indicator function𝟏_A(x) 在 x ∈ A 时取 1,否则取 0;它的积分就是 A 的测度。
可数countable元素能排成第 1 个、第 2 个、第 3 个……的一列(有限也算),例如自然数、有理数。
一致收敛uniform convergence所有点上的最大差距 supₓ |fₙ(x) − f(x)| → 0。比逐点收敛强:不允许某些点收敛得特别慢。
不可数uncountable无论怎么排都排不成一列,例如区间 [0,1] 里的全体实数。
零测集null set测度为 0 的集合,例如单点、有理数全体、康托尔集。
逐点收敛pointwise convergence对定义域里的每一个点 x,数列 fₙ(x) 都收敛到 f(x);这样得到的函数 f 叫逐点极限。
下确界infimum (inf)最大的下界。例如 (0,1] 没有最小值,但下确界是 0。
事件event样本空间的(允许谈概率的)子集;实现的 ω 落在里面,就说事件发生。
区间interval数轴上的一段,[a,b] 包含端点,(a,b) 不包含端点。
上界upper bound数集 A 的上界是一个数 M,A 里每个数都 ≤ M。例如 2 和 5 都是 [0,1] 的上界;所有上界中最小的那个叫上确界。
依分布收敛convergence in distribution在极限分布函数 F 的每个连续点 x 上,P(Xₙ ≤ x) → F(x)。它只比较分布,不要求 Xₙ 本身靠近 X。