Skip to content

算法Algorithm

全共形预测

Full conformal prediction · Transductive conformal prediction · 全保形预测

把候选响应加入数据后对称地重新拟合,以全部分数的秩检验候选值,在可交换性下构造有限样本预测集合。

形式陈述 ​

全共形预测逐个询问:“若新响应等于 y,它在包含自己的这份数据中是否显得异常?”设 n≥1,N=n+1,观测对 Zi=(Xi,Yi) 与未来对 ZN 可交换。取确定、可测且不依赖样本排列的拟合算法 A,以及对各位置相同的实值评分规则 s(z,f);大分数表示不相容。算法 A 输出的预测器可以拟合得很差,覆盖结论不要求模型正确。

给定输入 x、候选响应 y∈Y 和错误率 0<α<1,构造增广数据及分数:

Dx,y=(Z1,…,Zn,(x,y)),f^x,y=A(Dx,y),Six,y=s(Zi,f^x,y)(i≤n),SNx,y=s((x,y),f^x,y).

每换一个候选值,都按这个定义重新拟合,并重新计算历史点和候选点的分数。定义上尾计数及预测集合

py=1N∑i=1N1{Six,y≥SNx,y},C(x)={y∈Y:py>α}.

计数包含候选点自身,也包含与它并列的分数。在上述条件及相关事件可测的前提下,

P{YN∈C(XN)}≥1−α.

概率对历史数据与未来观测共同取平均;它不是给定某份历史数据或每个固定 x 的条件保证。若真实响应处的 N 个分数几乎必然无并列,覆盖恰为 ⌈N(1−α)⌉/N;有并列时保留覆盖下界。

直觉

分割共形预测先冻结模型,让校准点和测试点都站在训练过程之外。全共形采用另一种对称性:把候选测试点也放进训练过程,让全部点共同参与拟合,再用同一把尺子评分。历史点可以重复使用,关键是不能只让历史点享受训练内优势。

只有当候选值恰好是真实的 YN 时,增广数据才恢复成真正的可交换样本。此时交换两个点的位置不会改变拟合结果,只会交换相应分数的位置。覆盖证明只需考察这个真实候选值,无需假设所有虚构候选数据也可交换。

候选响应改变拟合和全部比较分数

图中最后一项始终是候选点的分数。上行有两个分数达到候选分数,下行只有候选点自己达到;是否并列直接决定端点的去留。

例子与边界

均值重拟合:完整反解出 [−1,3] ​

忽略输入,只用截距预测所有响应。历史响应为 0,1,2,取 α=1/4。加入候选 y 后,四个点共同拟合的均值是 μy=(3+y)/4,使用绝对残差得到

(S1y,S2y,S3y,S4y)=14(|y+3|,|1−y|,|5−y|,3|y−1|).

因为 py>1/4,候选自身之外至少还须有一个历史分数不小于它。令 t=y−1,最大的历史分数为

max(S1y,S2y,S3y)=14max{|t+4|,|t|,|4−t|}=4+|t|4.

最后一个等式可按 t≥0 与 t<0 分别检查。因此接受条件恰为

3|t|4≤4+|t|4⟺|t|≤2⟺−1≤y≤3,C=[−1,3].

在右端点 y=3,μ3=1.5,四个分数是 (1.5,0.5,0.5,1.5),故 p3=2/4>1/4,端点保留。在 y=4,μ4=1.75,分数变成 (1.75,0.75,0.25,2.25),故 p4=1/4,候选被拒绝。左端点同理由 |t|=2 保留。

这个计算说明给定这份数据时输出哪个集合。算法在可交换采样下具有至少 75% 的边际覆盖;本例含并列,不能据此宣称覆盖恰为 75%,更不能把覆盖概率解释为固定区间 [−1,3] 对任意响应分布都成立。

对称性保住有效性,却未必保住信息量 ​

设输入和响应各自独立且连续,训练算法记住每个训练输入的标签,在未见输入上输出零。直接用训练残差校准,会因残差全为零而在新输入上输出 {0};连续响应落入它的概率为零。全共形若对每个候选重新执行同样的记忆拟合,则候选也被记住,全部残差都是零,py=1,输出整个 R。这里修复了对称性,但没有得到有用的预测精度。

更一般地,计数至少包含自身,故 py≥1/N。若 α<1/N,所有候选都被接受。反过来,一般评分规则下集合未必是区间,也未必能解析求出端点。对连续响应只检查有限网格并返回网格点,不能自动继承原来整个响应空间上的覆盖结论。

推论与应用

含并列分数的秩证明 ​

在真实候选 (x,y)=(XN,YN) 处,记分数为 S1,…,SN。输入可交换且拟合对排列不变,所以分数向量可交换。给每个分数附上随机数 Ui;这些随机数相互独立、均服从 Uniform(0,1),且整个随机数向量独立于原始数据及分数向量。按 (Si,Ui) 从大到小排序;它们仍可交换且几乎必然没有并列。记第 N 个点的降序秩为 RN,则 RN 在 {1,…,N} 上均匀。

令 MN=#{i:Si≥SN}。随机打破并列只会把第 N 个点放在其并列组的某个位置,故逐点有 RN≤MN。于是对任意 u∈[0,1],

P{pYN≤u}=P{MN≤Nu}≤P{RN≤Nu}=⌊Nu⌋N≤u.

真实候选的计数因此具有p 值的超均匀性质;取 u=α,再对拒绝事件取补集即得覆盖。若没有并列,则 RN=MN,覆盖为 1−⌊Nα⌋/N=⌈N(1−α)⌉/N。辅助随机数只用于证明,上述算法本身并不需要随机打破并列。

与候选相关的分位数形式 ​

令 k=⌈N(1−α)⌉=N−⌊Nα⌋,并令 qy 为 S1x,y,…,Snx,y,+∞ 的第 k 个顺序统计量。则

py>α⟺SNx,y≤qy.

当 k≤n 时,右侧等价于至少 n−k+1=⌊Nα⌋ 个历史分数不小于候选分数,恰好就是左侧计数条件;k=N 时两侧对所有候选成立。这也解释了为何使用非严格分数比较,却使用严格的 py>α 接受规则。与分割法不同,qy 本身随候选改变,不能先用原模型算一次历史残差,再把这个阈值冒充全共形阈值。

有限标签空间可以枚举全部候选。若有 m 个标签,一次对 N 个点拟合的成本为 TA(N),逐点评分成本为 cs,直接计数的总成本为 O(m[TA(N)+Ncs]),不必为了求 py 对分数排序。连续回归的难点则常在反解整个集合;均值例子的绝对值结构使这一步可以精确完成。

自测 ​

保持历史响应 0,1,2 和同一拟合规则,把错误率改为 α=0.2,集合如何变化?此时 N=4,所有 py≥1/4>0.2,所以输出 R。这不是均值估计失灵,而是四个可交换位置提供的离散秩尚不足以排除任何候选。

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

拖动节点调整位置。

显示关系

显示:依赖

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