Skip to content

定理Theorem

诚实置信带与适应性边界

Honest adaptive confidence bands · Honesty and adaptation · 自适应置信带的不可能性

在固定设计正态回归中构造可复算的诚实带,并用局部尖峰与总变差证明未知光滑性下窄带适应的障碍。

形式陈述 ​

数据看上去很平滑,是否就能把整条曲线的置信带放心收窄?局部平滑区间已经说明偏差和同时覆盖都要校准;本页再问:光滑程度未知时,能否同时保证对较粗糙函数有效、对较光滑函数很窄?

观察固定设计回归

Yi=f(i/n)+σZi,i=1,…,n,Zi独立 N(0,1),σ>0已知.

程序若使用额外随机数,规定其分布与未知函数无关,且 Pf 包含这份共同辅助随机性。乘上同一随机化分布不会增加总变差,因此下面的比较也覆盖这类程序。

这里的“固定设计”只描述观测位置。目标是覆盖整个连续区间 [0,1]上的函数。对 0<r≤1及 L>0,定义

Fr(L)={f∈C[0,1]:‖f‖∞≤L,|f(x)−f(y)|≤L|x−y|r 对所有 x,y}.

L同时限制高度和Hölder半范数,但本页没有把二者的和当作球范数。将置信区间的统一覆盖要求用于整条函数,随机带 Cn(x)=[ℓn(x),un(x)]满足 ℓn≤un,并假定下面的同时事件与宽度都可测。记

Wn=supx∈[0,1][un(x)−ℓn(x)].

若 inff∈FPf{f(x)∈Cn(x)对所有 x}≥1−α,称它在类 F上有限样本诚实。统一渐近版本是 lim infninff∈FPf{f∈Cn}≥1−α。诚实指覆盖承诺统一适用于声明的类,不是描述估计者的品格。

本页证明以下有限样本障碍。设 0<α<1/2,选常数 c>0满足

cσ≤L,κ=1−2α−32c>0.

令 h=n−1/(2r+1)≤1/2。任何在 Fr(L)上诚实的带,都必须在零函数下满足

(1)P0{Wn≥cσn−r/(2r+1)}≥κ,E0Wn≥κcσn−r/(2r+1).

这是一个不带对数因子的下界,已经足以排除更光滑类的快收窄要求;它不是声称给出了尖锐的最优带宽常数或分离半径。

直觉

一条低矮且很窄的峰可以对观测只造成少量总变化。它在某一点的高度却仍超过一条极窄置信带的容纳范围。若数据分布很难区分“完全平坦”和“藏着这条峰”,任何诚实带就不能在平坦函数下总是很窄。

这里的障碍不取决于先用核、小波还是样条估计。它对所有可测的带构造同时成立,包括复杂的调参和随机化。换一个更准确的点估计器,不能自动消除观测本身分不清两种函数的困难。

例子与边界

两条函数怎样迫使带变宽 ​

先给一个适用于任意两条候选函数的步骤。设 f0,f1的观测分布为 P0,P1,d=‖f1−f0‖∞。若带在两条函数下的同时覆盖概率都至少为 1−α,则由两点法的总变差比较,事件 A1={f1∈Cn}满足

P0(A1)≥P1(A1)−TV(P0,P1)≥1−α−TV(P0,P1).

再与 A0={f0∈Cn}交集,在 P0下有概率至少 1−2α−TV。同一个区间若同时包含两个函数在每一点的值,最大宽度至少为 d。所以

(2)P0(Wn≥d)≥1−2α−TV(P0,P1).

这是对宽度的检验型限制。右侧若为负,只给平凡下界;需要选择分布足够接近的函数。

固定网格上的具体局部峰 ​

令 ψ(v)=(1−|v|)+,取 f0=0和

f1(x)=cσhrψ(x−1/2h),h=n−1/(2r+1).

ψ为1-Lipschitz、取值在 [0,1],因此 |ψ(u)−ψ(v)|≤min(|u−v|,1)≤|u−v|r。这给 [f1]r≤cσ≤L,且 ‖f1‖∞=cσhr≤L。两函数确实都在声明的类中。

噪声是已知协方差的多元正态。逐坐标展开正态对数密度比,KL散度在这里一律用自然对数(单位nat),因此为

D(P1‖P0)=12σ2∑i=1nf1(i/n)2.

峰的支撑长为 2h,其中最多有 2nh+1个网格点。利用 nh2r+1=1及 h2r≤1,

D(P1‖P0)≤c22(2nh+1)h2r≤32c2.

Pinsker界给 TV≤D/2≤(3/2)c。函数在 x=1/2处的距离恰为 cσhr,代入式(2)就得到式(1),期望界再由 EWn≥dP(Wn≥d)得到。

即使奇数 n的网格没有 1/2,峰值距离仍然成立,因为带覆盖的是连续域。若目标只覆盖观测点,就必须用最大采样峰值重做距离那一步,不能原样照搬这个高度。

一组可以手算的常数 ​

取 r=1/2,L=1,σ=1,α=0.05,c=0.2,n=256,则 h=1/16、峰高为 0.05。共有31个严格位于支撑内部的非零采样值。按对称三角形平方求和,可得

D(P1‖P0)=0.013359375,D/2≈0.0817294.

所以式(2)至少给 P0(Wn≥0.05)≥0.8182706,期望宽度至少约 0.0409135。使用不计算实际网格的统一上界也可得到较松的 0.7267949概率下界;两者都是下界,不是实际覆盖概率。

推论与应用

为什么更光滑的自适应带受到阻碍 ​

固定 r<s≤1。本页定义下 Fs(L)⊆Fr(L),因为 |x−y|≤1时 |x−y|s≤|x−y|r。若希望同一程序在大类上诚实,同时在更光滑类上满足

(3)supf∈Fs(L)EfWn≤K(log⁡nn)s/(2s+1),

其中 K不随 n变化,那么零函数也必须满足该上界。但 s/(2s+1)>r/(2r+1),式(3)右侧与式(1)下界的比值趋零,产生矛盾。零函数已经最光滑,仍无法在较粗糙类的统一覆盖承诺下获得所要求的宽度。

统一渐近诚实也有相同障碍:对任意小 ε>0,足够大的 n对整个类的覆盖都至少为 1−α−ε,在式(2)中多扣 2ε即可。若只知道每个固定函数的点态渐近覆盖,不能把它用于随 n改变的峰 f1;这正是统一量词发挥作用的位置。

光滑界已知时,一条诚实带确实能算出来 ​

负结果并不排除诚实带本身。设 r,L已知、n=Mq,把 [0,1]分成 M个宽 1/M的箱,第 j箱包含观测 i=(j−1)q+1,…,jq。函数评价时第一箱包含0,其余按右闭端点分配。每箱响应均值记为 Y―j。

任一箱中目标 f(x)与该箱观测真均值之差至多 LM−r,噪声标准差为 σ/q。因此在该箱输出

(4)Cn(x)=Y―j ±[LM−r+z1−α/(2M)σq].

对 M个正态箱均值用并集界,以至少 1−α的概率同时控制它们的噪声;在同一事件上,已知Hölder界控制每个箱内的全部连续位置,故式(4)在整个 Fr(L)上诚实。这里不需要在无穷多个位置分别累加错误率。

例如 n=64,M=4,q=16,r=1,L=σ=1,α=0.05,半宽为 1/4+z0.99375/4≈0.8744。它可能宽,但保证明确。平衡 M−r与 Mlog⁡M/n会提示 M≍(n/log⁡n)1/(2r+1);这是在已知 r,L下选择分箱尺度,不能把数据估出的 r直接代入就继承统一覆盖。

可以改变什么要求 ​

若预先排除与较光滑类统计上难区分的粗糙函数,或加入可检测的自相似性条件,某些模型允许自适应诚实带。排除范围与可检测尺度必须由对应定理指定,不能只从图像上判断曲线“没有坏峰”。Massé–Meiniel研究固定设计回归中的分离条件,其尖锐必要性定理另要求 r>1/2;Hoffmann–Nickl的相关结果则在密度估计模型中。本页给出的两点证明具有自己的 0<r≤1范围,不替代这些更精确结论。

改变损失或覆盖对象也可能改变结论。L2置信球容许窄而高的局部变化在积分下贡献较小,其适应性范围不能由上确界带的负结果直接判定。共形预测覆盖未来响应,通常是输入分布下的边际事件,也不是本页整条均值函数的同时覆盖。问题对象变了,需要重新说明承诺。

自测。 在较粗糙类上保持95%诚实,却只在若干典型平滑函数上展示窄带,是否证明了式(3)?没有。式(3)对整个平滑类取上确界,且零函数也受式(1)约束;有限演示既不检查统一覆盖,也不检查统一宽度。

自测二。 式(4)若把 L改成由同一批数据估出的较小值,还能引用原证明吗?不能。原证明需要对所有允许真函数成立的确定偏差上界;一个可能低估的拟合曲率或连续模不是这个上界。可以另构造具有联合错误预算的有效上界,但需要新的保证。

参考资料
  • Pierre-Yves Massé and William Meiniel,Adaptive Confidence Bands in the Nonparametric Fixed Design Regression Model,arXiv:1207.3975v3,Definitions 1–3、Theorem 2及§2.2:固定网格、已知正态方差、诚实性与分离条件;必要性部分要求 r>1/2,α<1/2。本文式(1)是独立证明的较弱无对数障碍。
  • Marc Hoffmann and Richard Nickl,On Adaptive Inference and Confidence Bands,Annals of Statistics 39(5),2011,pp.2383–2409,§§1–2:IID密度估计中的诚实性、适应性及可区分子类。本文不将其密度模型定理直接当作回归定理。
  • Adam D. Bull and Richard Nickl,Adaptive Confidence Sets in L2,Probability Theory and Related Fields 156,2013:说明积分损失下的置信集合具有不同的适应性边界;本页未陈述该文全部条件。
关系图谱9 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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