Skip to content

定理Theorem

Sanov 定理

Sanov theorem

以相对熵刻画 IID 经验测度偏离总体分布的指数成本,有限字母表下可由类型计数推导。

形式陈述 ​

样本均值只保留一个数。如果关心整份频率表异常,哪种函数取代 Cramér 定理的均值速率?Sanov 定理的答案是相对熵。

设 X1,X2,… 为Polish 空间 E 上服从 P 的IID 样本,经验测度为 Ln=n−1∑iδXi。在概率测度空间的弱拓扑下,Ln 以速度 n 满足LDP,良速率为

I(Q)=D(Q‖P)={∫log⁡dQdPdQ,Q≪P,+∞,否则.

这里的 D 是采用自然对数的KL 散度,方向是异常经验律 Q 相对于真实抽样律 P。若用比特单位 D2,则 I=D=(ln⁡2)D2,概率指数相应为 e−nD=2−nD2。弱拓扑意味着用有界连续函数的积分检验接近,而不是用总变差检验。

有限字母表 {1,…,r}(r≥1)上,若 P=(p1,…,pr),则 I(q)=∑jqjlog⁡(qj/pj)。任何 qj=0 的项都取零,包括 pj=qj=0 的情形;若 pj=0<qj,该项为无穷,因为原模型不可能产生这个类别。

直觉

一种频率表的概率有两部分:拥有这张表的排列有多少,以及每个排列在原模型下有多大概率。熵计算排列数量,交叉熵计算单个排列的概率,两者相减后恰好留下相对熵。

有限表上,若 nqj 都是非负整数,多项式系数计数相同频率表的有序样本,因此

P(Ln=q)=n!∏j(nqj)!∏jpjnqj.

类型计数给出 (n+1)−re−nD(q‖p)≤P(Ln=q)≤e−nD(q‖p),而类型总数至多 (n+1)r。可以不使用渐近阶乘公式来核验前一个界:在抽样律恰为 q 时,该类型的概率是所有类型中的最大值,所以介于 (n+1)−r 与 1 之间;换回抽样律 p,同类型的每条样本序列恰多出因子 e−nD(q‖p)。最大值性质也可逐项验证:写 kj=nqj,当总计数为 n 时,只需比较 ∏jkjmj/mj!,每个正整数 kj 对应的因子在 mj=kj 处达到最大;kj=0 则只允许 mj=0。这些多项式因子在对数除以 n 后消失。对闭集汇总各类型得上界,对开集选一串逼近目标分布的可实现类型得下界。一般 Polish 空间版本再用有限划分逼近和紧性控制,不能只把有限求和改为积分便算证明完成。

例子与边界

三种结果的真实概率为 p=(1/2,1/4,1/4),研究第一类经验比例不超过 1/4。固定 q1=x 后,剩余质量 1−x 分给第二、三类的最小成本方案保持原比例 1:1,即 q2=q3=(1−x)/2。

代入后,优化降为二元相对熵 D(x‖1/2)。它在 x<1/2 时随 x 增大而减小,所以约束下的最小点是

q∗=(1/4,3/8,3/8),I(q∗)=14log⁡12+34log⁡32.

闭约束与内部 q1<1/4 可逼近同一最小点,因此事件概率的对数除以 n 趋于 −I(q∗)。定理不仅给出成本,也指出异常样本最可能怎样把少掉的第一类质量重新分配。

若 P 连续,有限样本经验测度对 P 奇异,故每个实现的 D(Ln‖P) 都为无穷。这不与 Sanov 矛盾:定理计算的是经验测度落入一个固定弱邻域的概率,不是把每个随机原子测度代入速率后当成其点概率。弱邻域可包含相对熵有限的连续分布。

换成依赖样本时,通常不能保留 D(Q‖P) 作为完整速率。相同的一维边缘不代表相同的频率表概率,转移结构可能改变成本。

推论与应用

在有限字母表上,把 Q 映到 ∑jg(j)qj 连续,收缩原理给出统计量速率 infQ:EQg=xD(Q‖P)。若再条件于稀有频率约束,最小相对熵分布将成为Gibbs 条件化原理的核心对象。

参考资料
关系图谱34 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系