形式陈述
样本均值只保留一个数。如果关心整份频率表异常,哪种函数取代 Cramér 定理的均值速率?Sanov 定理的答案是相对熵。
设 为Polish 空间公理库Polish 空间Polish space可分且可完全度量化的拓扑空间;用换度量、可数乘积和子空间例子区分拓扑性质与指定度量的完备性。 上服从 的IID 样本公理库独立同分布样本IID sample · Independent and identically distributed sample以乘积分布描述来自同一总体的独立重复观测。,经验测度公理库经验测度Empirical measure · Empirical distribution把有限观测的重复频数编码为原子概率测度。为 。在概率测度空间的弱拓扑公理库测度的弱收敛Weak convergence of measures · Weak convergence of probability measures以所有有界连续测试函数的积分收敛来定义有限 Borel 测度的拓扑收敛。下, 以速度 满足LDP公理库大偏差原理Large deviation principle · LDP用开集下界和闭集上界刻画概率的指数衰减,区分速率函数、拓扑与精确概率。,良速率为
这里的 是采用自然对数的KL 散度公理库KL 散度Kullback–Leibler divergence · Relative entropy同一可测空间上分布 P 相对于 Q 的对数 Radon–Nikodym 导数在 P 下的积分。,方向是异常经验律 相对于真实抽样律 。若用比特单位 ,则 ,概率指数相应为 。弱拓扑意味着用有界连续函数的积分检验接近,而不是用总变差检验。
有限字母表 ()上,若 ,则 。任何 的项都取零,包括 的情形;若 ,该项为无穷,因为原模型不可能产生这个类别。
直觉
一种频率表的概率有两部分:拥有这张表的排列有多少,以及每个排列在原模型下有多大概率。熵计算排列数量,交叉熵计算单个排列的概率,两者相减后恰好留下相对熵。
有限表上,若 都是非负整数,多项式系数公理库多项式系数Multinomial coefficient把 n 个可区分对象分入若干有标号组时的计数系数 n!/(n1!⋯nk!)。计数相同频率表的有序样本,因此
类型计数给出 ,而类型总数至多 。可以不使用渐近阶乘公式来核验前一个界:在抽样律恰为 时,该类型的概率是所有类型中的最大值,所以介于 与 之间;换回抽样律 ,同类型的每条样本序列恰多出因子 。最大值性质也可逐项验证:写 ,当总计数为 时,只需比较 ,每个正整数 对应的因子在 处达到最大; 则只允许 。这些多项式因子在对数除以 后消失。对闭集汇总各类型得上界,对开集选一串逼近目标分布的可实现类型得下界。一般 Polish 空间版本再用有限划分逼近和紧性控制,不能只把有限求和改为积分便算证明完成。
例子与边界
三种结果的真实概率为 ,研究第一类经验比例不超过 。固定 后,剩余质量 分给第二、三类的最小成本方案保持原比例 ,即 。
代入后,优化降为二元相对熵 。它在 时随 增大而减小,所以约束下的最小点是
闭约束与内部 可逼近同一最小点,因此事件概率的对数除以 趋于 。定理不仅给出成本,也指出异常样本最可能怎样把少掉的第一类质量重新分配。
若 连续,有限样本经验测度对 奇异,故每个实现的 都为无穷。这不与 Sanov 矛盾:定理计算的是经验测度落入一个固定弱邻域的概率,不是把每个随机原子测度代入速率后当成其点概率。弱邻域可包含相对熵有限的连续分布。
换成依赖样本时,通常不能保留 作为完整速率。相同的一维边缘不代表相同的频率表概率,转移结构可能改变成本。
推论与应用
在有限字母表上,把 映到 连续,收缩原理公理库大偏差收缩原理Contraction principle沿连续映射转移大偏差原理,以所有原像中最小的速率作为输出成本。给出统计量速率 。若再条件于稀有频率约束,最小相对熵分布将成为Gibbs 条件化原理公理库Gibbs 条件原理Gibbs conditioning principle在稀有经验约束下,以最小相对熵分布描述固定少量样本的渐近条件规律。的核心对象。
参考资料