Skip to content

定理Theorem

标量 Freedman 不等式

Scalar Freedman inequality · 可预测方差鞅尾界

用有界鞅差及可预测条件方差控制预算内的全时间越界,证明指数过程并给出适应门控与分层方差预算的实际计算。

形式陈述 ​

一个过程可以根据过去暂停测量或改变下一步的分布。独立求和公式不再适用,但若每次给定过去后的增量均值为零,还能利用其条件方差。Freedman 的任务是:在累计可预测方差不超过预设预算时,控制任何时刻的大偏移。

设 (Sk,Fk)k≥0 是实鞅,S0=0,增量 Di=Si−Si−1 平方可积。设存在确定的 b>0,使 Di≤b 几乎处处。定义

(1)σi2=E[Di2∣Fi−1],Vk=∑i=1kσi2.

因条件均值为零,σi2 也是条件方差。每个 σi2 在本次增量出现前已由过去确定;Vk 因而称为可预测二次变差。它一般是随机过程。

对预先固定的 t>0,v>0,

(2)Pr{∃k≥1:Sk≥t 且 Vk≤v}≤exp⁡{−t22(v+bt/3)}.

单侧上尾只需 Di≤b,不要求 Di≥−b。如果另有 |Di|≤b,对 −S 也能应用,并由并集界得到双侧版本:右边乘2,事件改成 |Sk|≥t。

给定单侧失败概率 0<δ<1,记 ℓ=log⁡(1/δ)。一个方便阈值为

(3)t(v,δ)=2vℓ+2bℓ3.

将其代入式(2),概率至多 δ。若想用较小的反解值,可取

bℓ3+2vℓ+(bℓ/3)2;

式(3)只是便于加总与比较的较松版本,二者不是同一个常数。

直觉

Azuma 的幅度预算让每一轮都按最坏波动付费。Freedman 仍保留单步上界,以排除稀有但任意巨大的正跳跃;与此同时,它用给定过去后的真实二阶矩为常见小波动定价。

式(2)把两个条件放在同一个事件里。它不要求整条路径永远满足 Vk≤v,只控制方差预算仍未超额时是否越过高度 t。一旦预算超过 v,这个固定预算证书便不再保护后来的时刻。

从一个实变量不等式开始 ​

对实数 u,定义

h(u)=∫01(1−s)esuds,

则 eu=1+u+u2h(u),且 h 单调递增。因此 d≤b、λ≥0 时,

eλd≤1+λd+d2b2(eλb−1−λb).

这个论证也处理任意大的负 d,没有在单侧定理的证明里偷偷加入双侧有界。对 0≤z<3,利用 j!≥2⋅3j−2(j≥2)逐项求和,

ez−1−z=∑j≥2zjj!≤z22(1−z/3).

令 0<λ<3/b,定义 ψb(λ)=λ2/[2(1−λb/3)]。给定过去取条件期望,线性项消失,于是

(4)E[eλDi∣Fi−1]≤1+ψb(λ)σi2≤eψb(λ)σi2.

一个预算保护所有时刻 ​

固定这样的 λ,由式(4)逐步相乘可知

Lk(λ)=exp⁡{λSk−ψb(λ)Vk}

是初值1的非负上鞅。非负性和条件期望递推也保证其可积。若某时刻有 Sk≥t,Vk≤v,就有 Lk≥eλt−ψb(λ)v。因此Ville 不等式给

Pr{∃k:Sk≥t,Vk≤v}≤e−λt+ψb(λ)v.

取确定值 λ=t/(v+bt/3);因 v>0,它严格小于 3/b。代回后指数恰为 −t2/[2(v+bt/3)],证明式(2)。这里优化的是事先给定的两个数 t,v;没有根据观察到的路径临时挑 λ。

最后核对方便阈值。令 a=2vℓ、d=2bℓ/3,则 t=a+d,且

t2−2ℓ(v+bt/3)=(a+d)2−a2−d(a+d)=ad≥0,

所以式(3)确实达到要求的指数。

例子与边界

按过去决定是否继续抽样 ​

设 Yi 独立服从 Bernoulli(0.01)。在看到第 i 次结果之前,依据过去选择 Ii∈{0,1},表示是否计入该次观测。令

Di=Ii(Yi−0.01),Nk=∑i≤kIi.

因为 Ii 在给定过去时已经确定,E[Di∣Fi−1]=0,并且

Di≤0.99,Vk=0.0099Nk.

规则可以遇到某种历史信号后暂停,也可以一直测量到指定条件出现;只要计入决定先于新结果,鞅条件就保留。

取 b=0.99,v=9.9,δ=0.01,式(3)给

t=19.8log⁡100+0.66log⁡100≈12.5883583.

所以以至少99%的概率,在所有 Nk≤1000 的时刻,计入成功数减去 0.01Nk 都小于上述精确阈值 t(约12.5883583)。若恰好计入1000项,则成功数至少23会触发该单侧越界;能否看见这种越界不是保证,假设正确时它在整个预算窗口内发生的概率至多1%。

若忽略小方差,固定1000次的 Hoeffding 单侧计数阈值为 1000log⁡100/2≈47.9853。这个数用于说明幅度预算可能更粗;它本身只有固定时刻含义,不能不经证明就作为此处任意时刻的竞争边界。

四条路径重算适应过程 ​

把上例改为公平硬币,最多看三次,并在首次正面后不再计入。于是 I1=1,随后仅当此前全为反面时 Ii=1。四种终局如下:

首次正面 概率 计入数 终局 S 终局 V
第1次 1/2 1 1/2 1/4
第2次 1/4 2 0 1/2
第3次 1/8 3 −1/2 3/4
三次都没有 1/8 3 −3/2 3/4

加权终局均值为0,虽然后来的增量显然不是独立的。事件“某时刻 S≥0.4 且 V≤0.25”仅在第一条路径发生,真实概率为 1/2。Freedman 取 b=0.5 给上界

exp⁡{−0.422(0.25+0.5⋅0.4/3)}≈0.776754,

确实覆盖真实概率。公式在这种小预算下较松,不应把上界当成精确尾概率。

看完当前结果才选择,不是可预测 ​

若改成 Ii=1{Yi=1},仍写 Di=Ii(Yi−p),则 Di=(1−p)Yi,条件均值为 p(1−p)>0。只收录成功记录引入了漂移,第一步的鞅条件已经失败。把所得计数代回方差公式不能修复这种选择偏差。

此外,Vk=∑iE[Di2∣Fi−1] 不是观测平方和 ∑iDi2,也不是围绕样本平均的样本方差。它可能含未知参数;实际报告前要有模型下可计算的上界,或在候选参数下构造合法检验。经验 Bernstein用校准后的样本方差,承担的是不同任务。

推论与应用

何时真的可以适应观察到的方差 ​

式(2)没有直接证明 Sk<t(Vk,δ) 对所有 k 同时成立。要让阈值随随机方差变化,可以提前列出预算并分配错误概率。

固定 v0>0,设 vj=2jv0,j=0,1,…,并设

δj=δ(j+1)(j+2),∑j=0∞δj=δ.

对每个 j 使用式(2)–(3),再对这些可数事件取并集。令 j(V) 为满足 V≤vj 的最小非负整数,则以至少 1−δ 的概率,对所有 k,

(5)Sk<t(vj(Vk),δj(Vk)).

实际 Vk 可以随机;现在有效性来自整张预算表共同覆盖,不是未经校准地把随机值代入一个固定预算定理。例如 b=1,v0=10,δ=0.05,观测到 Vk=13 时属于 j=1,用 v1=20,δ1=1/120,阈值为17.0300017。不能仍使用 v=13 和完整的0.05预算冒充同一证明。

若所有增量独立,并采用它们的自然滤过(或另证每个 Di 与整个 Fi−1 独立),则 σi2=Var(Di) 为确定数。在固定时刻 n,取 v=∑i≤nVar(Di),式(2)包含该时刻事件并恢复标量 Bernstein型界。独立情形的原证明更短;这里新增的是条件方差、适应性和预算内全时间控制。

本页不要求条件对称。若没有幅度界、但有条件对称结构,可考察混合边界的自归一化支路;若既无幅度、也无其他尾部结构,仅有限条件方差不能凭空给指数尾。

自测与答案 ​

  1. 在稀有故障例中,把计入读数变成原值的5倍,b,V,t 各如何变化?若要双侧99%保证,对数项怎样调整?
  2. Vk≤v 只在某条已观察路径上成立,能否使用式(2)?这是否意味着任意随机 v=Vk 都可直接代入?

答案:第一题 b 乘5,V 乘25,阈值乘5。双侧还要有绝对增量界;原中心化读数可用 b=0.99 的双侧界,因为负幅度0.01较小;总预算0.01时每侧用0.005,对数变为 log⁡200。第二题式(2)本来就控制偏差与预算的联合事件,因此不要求所有路径满足预算;但数字 v 必须预先固定,随机选择预算需要式(5)这样的共同校准。

参考资料
  • Joel A. Tropp,Freedman’s inequality for matrix martingales,arXiv:1101.3039v1,2011-01-16;§1.2 Theorem 1.1(PDF第2页)明确给出标量、增量单侧上界和存在时刻的联合事件。本文只用这一标量入口,完整重证指数上鞅,不依赖矩阵 Lieb 定理。
  • David A. Freedman,On Tail Probabilities for Martingales,Annals of Probability 3(1), 1975,100–118,Theorem 1.6 为原始来源;本文常数形式依据上面已核验的 Tropp 陈述。
关系图谱9 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

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

上位 / 更一般

暂未标注直接上位概念。

下位 / 直接特例

类型化关系