统计学的下界证明有一个几乎统一的套路:先把估计问题换成一个检验问题,再证明这个检验问题很难。
设 X1,…,Xn∼N(θ,1),要估计 θ。样本均值给出风险 1/n。这个 1/n 是不是最优的?也就是说,是否存在某个估计量,能让最坏情况下的风险比 1/n 更小?Le Cam 的两点方法给出的回答是:不存在。它的论证只用了两个候选参数 θ0=0 与 θ1=Δ——如果连这两个点都分不清,那就更不可能把 θ 估准。
本文介绍这套语言里最常用的三个工具:Le Cam 的两点方法、Fano 不等式,以及 Assouad 的超立方体引理。三者的共同点是:先选一个有限的「困难子问题」,把估计下界转化为检验下界,再用总变差、KL 散度或互信息来度量「有多难分」。
一个估计量 θ^ 的优劣,通常用它的风险衡量:在真实参数为 θ 时,θ^ 的平均损失。但真实参数是未知的,所以我们不能只关心某一个 θ 下的风险。极小极大(minimax)的立场是:先让大自然挑一个最不利的 θ,再让统计学家挑估计量,看统计学家能保证什么。
用博弈的话说,极小极大风险是
θ^infθ∈Θsup EPθ[ℓ(θ,θ^)].这个量回答的是:在所有可能的参数中,最坏情况下的风险最小能压到多少。 一个上界对应一个具体的估计量;一个下界则对应一句否定:无论你多聪明,都不可能比这个数更好。
下界为什么难?因为你要同时对所有估计量说话。逐一检查是不可能的。所以下界证明总是走同一条路:把参数空间缩小到一个有限的、精心挑选的子集,使得在这个子集上「区分参数」本身就已经很难;如果区分都做不到,估计就更做不到。
这就是「归约到检验」。接下来的三种方法,区别只在于子集有多大、用什么距离度量「难分」。
设 {Pθ:θ∈Θ} 是样本空间 X 上的一族分布,参数空间 Θ 可以是任意集合。我们观测 n 个独立同分布的样本
Xn=(X1,…,Xn)∼Pθ⊗n,一个估计量是任意可测映射 θ^:Xn→Θ。损失函数 ℓ:Θ×Θ→[0,∞) 度量估计值 θ^ 与真值 θ 的差距,常用的有:
- 平方损失 ℓ(θ,θ^)=∥θ−θ^∥22;
- 绝对损失 ℓ(θ,θ^)=∥θ−θ^∥1;
- 0-1 损失 ℓ(θ,θ^)=1[θ=θ^]。
估计量 θ^ 在 θ 处的风险是
Rθ(θ^)=EPθ⊗n[ℓ(θ,θ^(Xn))].注意这里的期望是对数据取的,θ 固定。极小极大风险则再对 θ 取最坏情况:
Mn(Θ,ℓ)=θ^infθ∈ΘsupRθ(θ^).极小极大风险有一个贝叶斯下界。对任意 Θ 上的先验 π,记贝叶斯风险为
Bn(π)=θ^inf∫Rθ(θ^)π(dθ).由于「积分 ≤ 最大值」,
Bn(π)≤Mn(Θ,ℓ).所以,任何一个先验下的贝叶斯风险,都是极小极大风险的下界。这个观察非常有用:下界证明往往先选一个先验(通常是有限集上的均匀分布),再算贝叶斯风险。Le Cam 与 Fano 本质上都是在做这件事,只是把贝叶斯风险进一步转化成了检验误差。
对任意有限子集 Θ0={θ1,…,θM}⊆Θ,
Mn(Θ,ℓ)≥Mn(Θ0,ℓ).所以我们只需要在一个有限问题上证明下界。有限问题可以看成一场 M 元检验:自然以某个先验选出 θK,统计学家看到数据后要猜出 K。如果连「猜出是哪一个」都做不到,那么「把 θ 估到很小的损失」当然也做不到。
把这句话写成不等式,就是 Le Cam 与 Fano 的工具。
先看最简单的情形:只有两个候选参数。
定理 3.1(Le Cam 两点方法) 设 θ0,θ1∈Θ,ℓ(θ0,θ1)≥2δ,且 ℓ 满足三角不等式。则对任意估计量 θ^,
i=0,1max EPi⊗n[ℓ(θi,θ^)]≥2δ(1−P0⊗n−P1⊗nTV).
证明只用一个测试。令
ψ=1[ℓ(θ0,θ^)≥δ].在事件 {ψ=0} 上,ℓ(θ0,θ^)<δ,由三角不等式
ℓ(θ1,θ^)≥ℓ(θ0,θ1)−ℓ(θ0,θ^)>2δ−δ=δ.于是
EP0ℓ(θ0,θ^)≥δP0(ψ=1),EP1ℓ(θ1,θ^)≥δP1(ψ=0).两式相加:
EP0ℓ(θ0,θ^)+EP1ℓ(θ1,θ^)≥δ[1+P0(ψ=1)−P1(ψ=1)]≥δ(1−∥P0−P1∥TV).最后,最大值不小于平均值,得到 δ(1−TV)/2。证毕。
如果取 0-1 损失、δ=1/2,就得到最常见的检验形式:
θ^infi=0,1max Pi⊗n(θ^=θi)≥21(1−P0⊗n−P1⊗nTV).这个式子的意思是:两个分布的总变差越小,任何检验的错误率就越接近 1/2。 当 P0=P1 时,错误率至少是 1/2,这当然——数据里根本没有信息。
定理 3.1 要求 ℓ 满足三角不等式,平方损失并不满足。对平方损失,有一个常数稍有不同的版本。设 Δ=∥θ0−θ1∥22,则
i=0,1maxEPi⊗n∥θi−θ^∥22≥4Δ(1−∥P0⊗n−P1⊗n∥TV).证明用贝叶斯视角。取均匀先验,平方损失下的贝叶斯估计是后验均值
θ^B(Xn)=π(Xn)θ0+(1−π(Xn))θ1,π(Xn)=Pr(θ=θ0∣Xn).其风险为 ΔE[π(Xn)(1−π(Xn))]。由 π(1−π)≥21min(π,1−π) 与均匀先验下贝叶斯检验误差的下界 E[min(π,1−π)]≥21(1−TV),风险至少是 4Δ(1−TV)。极小极大风险不小于贝叶斯风险,得证。
总变差通常不好直接算,尤其对乘积分布。Pinsker 不等式把它换成 KL 散度:
∥P−Q∥TV≤21DKL(P∥Q).而 KL 对乘积分布是可加的:
DKL(P0⊗n∥P1⊗n)=nDKL(P0∥P1).代入定理 3.1:
i=0,1maxEPi⊗nℓ(θi,θ^)≥2δ(1−2nDKL(P0∥P1)).这个式子把下界拆成了两个因子:损失间隔 δ,和信息量 nDKL。要让下界为正,必须 2nDKL<1;换句话说,如果两个分布在这个样本量下已经能被可靠地区分,两点方法就给不出非平凡的下界了。
这也解释了为什么下界里的样本量总以「n× 单样本信息」的形式出现:独立样本的信息是线性累加的。
设 Xi∼Ber(p),p∈{1/2,1/2+Δ},损失为平方损失。两个分布的单样本 KL 是
DKL(Ber(21)Ber(21+Δ))=−21log(1−4Δ2)≈2Δ2(Δ→0).由 Pinsker,n 个样本的总变差至多 Δn。取 Δ=1/(2n),则 Δn=1/2,平方损失的版本给出
imaxE[(pi−p^)2]≥4Δ2(1−21)=32n1.所以 Bernoulli 均值的极小极大速率是 1/n。这个结论不稀奇,但推导方式值得记住:我们只用了两个点,就把「1/n 不可能被突破」这件事证明完了。
两点方法简单、常数清楚,但它只能看见一个方向。考虑 d 维高斯均值:
X1,…,Xn∼N(θ,σ2Id),θ∈Rd.如果只取 θ0=0、θ1=ae1,那么 KL 是 na2/(2σ2),Le Cam 给出的下界是 a2/4 的量级,取 a≍σ/n 后只有 1/n。但真实的极小极大风险是 σ2d/n——差了整整一个 d。
原因很直接:两点方法只探测了一个坐标方向,而 d 维问题的困难来自 d 个方向同时存在的不确定性。要捕捉这一点,需要让假设集合随着 d 指数增长。这就是 Fano 不等式和 Assouad 引理的用武之地。
定理 4.1(Fano 不等式) 设 K 在 {1,…,M} 上均匀分布,X 的分布依赖于 K,K^=K^(X) 是任意估计。记 Pe=Pr(K^=K),则
Pe≥1−logMI(K;X)+log2,其中 I(K;X) 是互信息。
证明分两步。第一步是条件熵的界:
H(K∣K^)≤h(Pe)+Pelog(M−1)≤log2+PelogM,其中 h(Pe)≤log2 是二元熵的最大值。第二步是数据处理不等式:
H(K∣K^)≥H(K∣X)=H(K)−I(K;X)=logM−I(K;X).两式合并,整理即得。
这个不等式有一个很直观的读法:如果 X 几乎不包含关于 K 的信息(I≈0),那么任何猜测的错误率至少是 1−log2/logM,当 M 很大时接近 1,也就是接近「瞎猜」。 反过来,要让错误率小,互信息必须达到 logM 的量级——你需要至少 logM 奈特的信息才能区分 M 个等可能的选择。
把 Fano 用在极小极大下界上,需要一个把「损失」翻译成「猜错」的分离条件。
定理 4.2(Fano 方法,度量形式) 设 d 是 Θ 上的度量,θ1,…,θM∈Θ 满足
d(θi,θj)≥2δ(i=j).令 K∼Unif{1,…,M},X∼PθK。则对任意估计量 θ^,
M1i=1∑MEPi[d(θi,θ^)]≥δ(1−logMI(K;X)+log2).特别地,若 I(K;X)≤β,则
imaxEPi[d(θi,θ^)]≥δ(1−logMβ+log2).
证明的关键是一个解码论证。定义最近点解码
K^=argimind(θ^,θi).若 d(θK,θ^)<δ,则对任意 j=K,
d(θ^,θj)≥d(θK,θj)−d(θK,θ^)>2δ−δ=δ>d(θK,θ^),所以 K 是唯一的最小值点,K^=K。因此
{d(θK,θ^)<δ}⊆{K^=K},从而 Pe≤Pr(d(θK,θ^)≥δ)。于是
E[d(θK,θ^)]≥δPr(d≥δ)≥δPe≥δ(1−logMI+log2).左边正是 M 个风险的平均,平均不超过最大,定理得证。
对平方损失 ℓ=d2,同样的解码论证把 d 换成 ℓ。若平方间隔为 Δ,即 d(θi,θj)≥Δ,则取 δ=Δ/2,得到
M1i∑E[ℓ(θi,θ^)]≥4Δ(1−logMI(K;X)+log2).这里用了 {d<δ} 蕴含解码正确,因此 E[d2]≥δ2Pr(d≥δ)≥δ2Pe。平方损失的常数是 Δ/4,与 Le Cam 的平方损失版本一致。本文后面的例子只关心速率,常数不影响结论。
剩下的任务是把互信息控制住。由定义,
I(K;X)=M1i=1∑MDKL(Pi∥Pˉ),Pˉ=M1i=1∑MPi.再用 KL 对第二个变量的凸性,
I(K;X)≤M21i,j=1∑MDKL(Pi∥Pj)≤i,jmaxDKL(Pi∥Pj).如果每个 Pi 都是 n 个独立样本的分布,那么成对 KL 是单样本 KL 的 n 倍。这样,Fano 下界就变成了「logM 与 n×KL 的竞争」。
定理 4.2 常被称为广义 Fano 不等式。它的使用可以总结成一句口诀:
只要 logM ≫ n⋅i,jmaxDKL(Pi∥Pj), 就有 Mn≳δ. 这里 M 是「用 δ 分开的假设个数」,也就是损失空间里的打包数(packing number)。于是下界问题变成了一个几何问题:在参数空间里能放多少个半径为 δ 的小球,同时又让它们在 KL 意义下足够接近?
Yang 与 Barron 的版本直接用 I(K;X)=M1∑iDKL(Pi∥Pˉ),而不是把它放缩成最大成对 KL。这在局部问题里更精确,常常能去掉多余的 log 因子。其思想是:不要要求所有假设两两都接近,只要它们整体上离它们的「重心」接近就够了。
Fano 还有一个局部版本:不要求 θi 覆盖整个参数空间,只要求它们落在一个半径 r 的小球里。非参数速率 n−2β/(2β+d) 正是这样来的——局部的打包数随 r 增长,而 KL 随 r 增长,两者平衡给出最优的 r。
Fano 处理的是「M 个任意假设」,Assouad 处理的是一种特殊但极有用的结构:超立方体。
定理 5.1(Assouad 引理,平方损失) 设 Θ={0,1}m,Pθ 是观测分布,损失为平方欧氏距离。对每个坐标 j,令
Qj,0=2m−11θ:θj=0∑Pθ,Qj,1=2m−11θ:θj=1∑Pθ.则
θ^infθ∈ΘmaxEPθ[∥θ^−θ∥22]≥4m(1−jmax∥Qj,0−Qj,1∥TV).
证明用贝叶斯视角。取 Θ 上的均匀先验。平方损失下的贝叶斯估计是后验均值 θ^B(X)=E[θ∣X],其贝叶斯风险为
j=1∑mE[πj(X)(1−πj(X))],πj(X)=Pr(θj=1∣X).固定 j。边际后验 πj 正是 Qj,0 与 Qj,1 在均匀先验下的后验概率。对任意 π∈[0,1],
π(1−π)≥21min(π,1−π),而均匀先验下的贝叶斯检验误差满足
E[min(πj,1−πj)]≥21(1−∥Qj,0−Qj,1∥TV),这是定理 3.1 的检验形式。于是每个坐标贡献至少 41(1−TVj),求和得
Bn≥4m(1−jmax∥Qj,0−Qj,1∥TV).极小极大风险不小于贝叶斯风险,定理得证。
Assouad 的威力在于它把 m 维问题拆成了 m 个一维问题。只要每个坐标上的两个混合分布足够接近,m 个坐标就会贡献 m 倍的下界——这正是高维速率里 d 的来源。
对 Θ={0,ε}m,令 θ=εθ′,其中 θ′∈{0,1}m。平方损失满足
∥θ^−θ∥22=ε2θ^/ε−θ′22,所以定理 5.1 直接给出
θ^infθ∈{0,ε}mmaxE∥θ^−θ∥22≥4ε2m(1−jmax∥Qj,0−Qj,1∥TV).如果损失换成 Hamming 距离,同样的逐坐标论证给出常数 1/2:把后验均值换成阈值检验即可。
模型:
X1,…,Xn∼N(θ,σ2Id),θ∈Rd,损失为 ∥θ^−θ∥22。
上界。 样本均值 Xˉ 的风险是
E∥Xˉ−θ∥22=nσ2d.下界。 把参数限制到超立方体 Θ0={0,ε}d。对坐标 j,混合分布 Qj,0 与 Qj,1 只在这一维的均值上相差 ε。由 KL 的联合凸性,
DKL(Qj,0∥Qj,1)≤2σ2nε2.Pinsker 给出
∥Qj,0−Qj,1∥TV≤2σεn.取 ε=σ/n,则 TV 至多 1/2。Assouad:
θ∈Θ0maxE∥θ^−θ∥22≥4ε2d(1−21)=8nσ2d.上下界只差常数 8,速率 σ2d/n 被夹住了。
这个例子也说明了 Fano 与 Assouad 的分工:Fano 需要构造一个 δ-分离的打包,Assouad 则直接利用乘积结构,把打包问题化成了逐坐标的 TV 计算。
模型:
Yi=f(Xi)+ξi,Xi∼Unif([0,1]d),ξi∼N(0,σ2),f 属于 β-Hölder 球 Σ(β,L),损失为 ∥f^−f∥L22。
上界。 核估计或局部多项式估计可以达到
n−2β+d2β(见 Tsybakov 第 1 章)。下界用 Fano。
取 h∈(0,1),把 [0,1]d 分成 m=h−d 个边长为 h 的小方块。取一个光滑的鼓包函数 φ,支撑在单位方块内,令
φj(x)=φ(hx−xj),其中 xj 是第 j 个方块的中心。对任意子集 S⊆{1,…,m},定义
fS=aj∈S∑φj,a=chβ.当 c 足够小时,fS∈Σ(β,L):鼓包的高度是 hβ、宽度是 h,所以 β 阶导数的量级是 hβ/hβ=1。
两个函数 fS,fS′ 的 L2 距离只取决于它们不同的那些鼓包。若 ∣S△S′∣=k,则
∥fS−fS′∥L22=a2hdk∥φ∥L22=c2h2β+dk∥φ∥L22.现在取一个长度为 m、相对 Hamming 距离至少 1/4 的二元码。这样的码有 M=2c0m 个码字(c0>0 是常数),任意两个码字至少在 m/4 个位置不同。于是
∥fS−fS′∥L22≥c1h2β,即平方间隔 Δ=c1h2β。
另一方面,两个回归分布的 KL 是
DKL(PfS⊗n∥PfS′⊗n)=2σ2n∥fS−fS′∥L22≤σ2nc2h2β.Fano 要求 logM≫nh2β,即
h−d≳nh2β⟺h≲n−2β+d1.取 h≍n−1/(2β+d),下界就是
Mn≳h2β≍n−2β+d2β,与上界同阶。证毕。
这个构造值得反复看:鼓包的宽度 h 控制两件事——损失间隔(h2β)和 KL(nh2β),而鼓包的个数 h−d 控制假设的个数。 三者平衡,就得到非参数速率。维度 d 出现在指数上,正是因为它决定了「同样宽度下能放多少个鼓包」。
| 问题 | 参数类 | 极小极大速率 |
|---|
| d 维高斯均值 | ∥θ∥2≤R | σ2d/n |
| 非参数回归 | β-Hölder,L | n−2β/(2β+d) |
| 密度估计 | β-Hölder,L | n−2β/(2β+d) |
| s-稀疏均值 | ∥θ∥0≤s≪d | nσ2slogsd |
| 一维单调回归 | 单调、有界 | n−2/3 |
这张表里的每一个速率,都可以用本文的三种方法之一推出来。稀疏均值的下界通常用 Fano 加一个随机支撑的打包;单调回归则用 Le Cam 加一个局部的「台阶」构造。
| 方法 | 假设个数 | 距离 | 优势 | 局限 |
|---|
| Le Cam | 2 | TV / KL | 简单,常数清楚 | 只探测一个方向,高维问题会松 |
| Fano | M | 互信息 / KL | 通用,适合高维与非参数 | 需要构造打包,常数较松 |
| Assouad | 2m | 逐坐标 TV | 乘积结构下最锋利 | 只适用于超立方体结构 |
三者的关系可以这样理解:
- Le Cam 是 Fano 的两点版本。 当只有两个假设时,Fano 退化为两点检验;Le Cam 用 TV 给出了更干净的常数。
- Fano 是 Le Cam 的多假设推广。 它把「两个点分不清」推广为「M 个点平均分不清」。
- Assouad 是 Fano 在乘积结构上的精细化。 它不直接对 2m 个假设做 Fano,而是逐坐标做 Le Cam,再把 m 个下界加起来。
实际使用时,选择哪种方法取决于参数空间的结构:
- 一维、局部、两样本问题:Le Cam 通常足够,而且常数最干净。
- 高维、非参数、需要 logM 的地方:Fano 是标准工具。
- 乘积结构、ℓ2 或 ℓ∞ 速率:Assouad 最方便。
回到开头的问题:为什么 1/n 不能再小?因为两个足够接近的分布,在 n 个样本下仍然难以区分;而任何估计器都必须先区分它们。这个论证只需要两个点,所以叫两点方法。
但当困难来自许多方向时,两个点就不够了。Fano 把「两个点」换成「M 个点」,用互信息度量整体的不可区分性;Assouad 则把 m 维问题拆成 m 个一维问题,逐坐标地累加下界。
如果只记三句话:
- 下界证明的第一步永远是归约:把估计问题换成一个有限检验问题。
- Le Cam 用两个点,Fano 用 M 个点,Assouad 用 2m 个点;工具在变,思想都是「构造一个在损失上离得远、在信息上离得近的假设集合」。
- 信息论下界告诉你统计上不可能;它不告诉你计算上不可能。深度学习时代真正难的,往往是后一个问题。