Skip to content

原则Principle

统计与数值误差预算

Statistical and numerical error budget

将样本量、局部曲率与求解残差联动,判定近似算法是否保留目标估计量的一阶推断。

形式陈述 ​

一项推断程序有两个不同输出:理论上定义的估计量 θ^n,和计算机真正返回的 θ~n。设确定尺度 an→∞ 下

an(θ^n−θ0)⇒Z.

这里的 Z 可以是渐近正态极限,也可为其他分布。若计算误差满足

an‖θ~n−θ^n‖→P0,

则 an(θ~n−θ0) 依分布收敛到同一个 Z。证明是在统计误差与计算误差的恒等分解中应用Slutsky 定理。本页研究如何从实际可检查的数值证书推出这个条件,而非把它当成求解器默认保证。

从平均准则差到一阶可忽略性 ​

设 Mn 为按样本量平均的随机准则。在概率趋于一的事件上,假定精确解 θ^n 和计算输出 θ~n 位于同一凸集合,且在该集合有已验证的二次增长

Mn(θ)−Mn(θ^n)≥μn2‖θ−θ^n‖2,μn>0.

若可认证 0≤Mn(θ~n)−Mn(θ^n)≤rn,则

an‖θ~n−θ^n‖≤an2rn/μn.

因此 an2rn/μn→P0 是充分条件。固定非退化曲率、an=n 时,它成为平均准则差 rn=oP(n−1)。若程序优化未平均的和 Ln=nMn,同一要求写成总准则差 oP(1);不能在两个定标之间沿用同一数字容差。

从梯度或方程残差到证书 ​

若 Mn 可微、μn-强凸,且精确极小点位于内点,则强单调性与 Cauchy–Schwarz 给出

‖θ~n−θ^n‖≤‖∇Mn(θ~n)‖/μn.

故 an‖∇Mn(θ~n)‖/μn→P0 足以保持极限。约束极小点的普通梯度不必为零,需要适当的投影梯度或 KKT 误差界,不能直接套用这个内点版本。

对一般向量方程,应另证局部逆界 ‖θ~n−θ^n‖≤Cn‖Ψn(θ~n)‖,再检验 anCn‖Ψn(θ~n)‖→P0。可逆 Jacobian 在解处只给局部信息;如果算法未进入该邻域,这不是全局证书。一般残差与误差估计解释这种数值转换,本页增加的是它与统计尺度的联动。

直觉

“计算误差比统计误差小”必须比较同一种单位。准则差、梯度残差、参数误差和标准误各用不同标尺;曲率或逆 Jacobian 负责转换。固定一个很小的数字容差,可能在小样本过度计算,在大样本又不足以保留区间的中心。

这项预算也不能改变统计目标。如果准则估计的是错配模型的伪真参数,数值精确只意味着更准确地求出那个参数,不意味着模型偏差消失。

例子与边界

一致性的容差不足以保护正态中心 ​

取 Xi∼N(μ,1) 独立同分布,Mn(θ)=12n∑i(Xi−θ)2。精确解为 X¯,曲率恒为1。若返回 θ~n=X¯+n−1/4,则准则差恰为 rn=1/(2n)→0,所以计算输出仍一致。

但其标准化误差为 n(X¯−μ)+n1/4,偏移发散,围绕该输出使用通常标准误 1/n 的固定水平区间覆盖趋于零。若偏移改为 c/n,准则差恰为 c2/(2n),标准化极限为 N(c,1)。这说明 OP(n−1) 差也不够保证不变中心;若只凭这个 gap 上界保护同一中心,可采用小 o 这一充分条件;否则需要进一步分析并把可能的偏移纳入分布。

一次可审计的停止选择 ​

仍取曲率1、标准误 1/n。希望计算误差最多占一个标准误的 τn,可要求

rn≤τn22n或‖∇Mn(θ~n)‖≤τnn.

n=10000、τn=0.01 时,参数误差预算为 10−4,平均准则差上限为 5×10−9,平均梯度残差上限为 10−4。总准则的 gap 上限则是 5×10−5。实际程序还需证明 gap 上界,例如用可行对偶下界;只把当前目标减去上一轮目标不能当成 rn。

固定的 τ=0.01 给“至多百分之一标准误”的数值控制,并不字面蕴含渐近可忽略。若要沿整个样本序列保持完全相同的一阶极限,需安排 τn→0。有限精度下若无法进一步下降,应报告当前误差预算和舍入平台,而非无限迭代。

弱识别时必须重新计算尺度 ​

令第 n 个实验观察 Yi=κnθ0+εi,εi∼N(0,1) IID,κn=n−1/4 已知。平均二次准则的精确解为 θ^n=Y¯/κn,曲率为 μn=κn2=n−1/2,统计标准误为 1/(nκn)=n−1/4。常规根号样本量方差稳定的假设在此失效,但有效尺度 an=nκn=n1/4 仍给精确标准正态。

返回 θ~n=θ^n+n−1/4 时,平均梯度残差为 κn2n−1/4=n−3/4=o(n−1/2)。它看似满足常规固定曲率下的残差要求,却恰好产生一个统计标准误的偏移。正确条件是

‖∇Mn(θ~n)‖=o(κn/n),

因为 an/μn=n/κn。当 κn=0 时,数据分布完全不依赖 θ,连识别都不存在;把容差继续收紧无法解决。

数据复用还会改变选择后的分布 ​

取正偶数 n,把 Xi IID N(μ,1) 的观测平均分成两半,得到两个独立半样本均值 T1,T2。每一个预先固定的方法都满足 n/2(Tj−μ)∼N(0,1)。若看过两者后挑大的那个报告为 T=max(T1,T2),则

P{n/2(T−μ)≤z}=Φ(z)2,

不是标准正态。例如用标准正态上侧5%临界值 z0.95 检验,真实拒绝概率为 1−0.952=9.75%。这个直接计算定位了错误校准;此特例的对称双侧区间可能碰巧保持覆盖,不能由此恢复完整的标准正态枢轴。优化器即使完全精确,选择步骤仍改变了推断对象的分布。把同一数据用于选择阈值、模型或估计规则后,固定配置下的点态定理不能直接套给选出的配置;需要隔离选择、联合分析或选择后的校准。

推论与应用

评价算法停止状态时,应同时记录:平均还是求和准则、目标参数或对比、统计抽样单位、可用误差证书、曲率或逆界、已达到的数值阈值,以及最大迭代数或非有限值退出。求解成本由算法另行报告;本页的阈值不预言达到它需要多少步。

近似 ERM把准则差加入总体风险界;这里则要求计算误差不改变参数的一阶分布。风险保证足够的容差可能不足以做参数区间,反之继续求解到机器精度也不修复泛化、识别或数据选择问题。

自测与答案 ​

  1. 若 an=n、μn≥1/2,认证平均 gap 为 n−3/2,是否足够?答案:足够,nrn/μn≤2n−1/2→0。
  2. 若只认证当前点和上轮点的目标差为 10−12,能否声称 gap 已达 10−12?答案:不能。缓慢变化、过小步长或舍入停滞都可能发生在远离最优处,尚缺最优值下界或有效残差误差界。
参考资料
  • Whitney K. Newey 与 Daniel McFadden,Large Sample Estimation and Hypothesis Testing,1994,§2.1(印刷2121–2122页)近似极值的一致性,§3.1 局部渐近正态;本页的容差换算由二次增长不等式直接推导。
  • A. W. van der Vaart,Asymptotic Statistics,1998,Ch.5:近似极值与近似方程估计。本页不把一般收敛容差等同于保留一阶推断的容差。
关系图谱13 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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