假设你经营一家大平台,每半年都要汇总成千上万次 A/B 测试,判断“通常一次改版能带来多大影响”。绝大多数实验结果平平,偶尔却冒出一个大得离谱的数字。它可能是真正的突破,也可能只是随机分组失败、指标配置出错。麻烦在于,一种常用估计方法会让这几个数字强烈影响总体判断,继而改变其他实验的效果估计和上线决策。
Manu Burson 和 Ryan Kessler 在这篇 arXiv 预印本中讨论的,正是这个问题。他们研究的对象并非所有 A/B 平台或所有先验,而是 normal-normal empirical Bayes(正态—正态经验贝叶斯)模型中,对真实影响方差所做的最大似然估计。论文提出 bounded influence maximum likelihood estimator,BIMLE:保留每项实验,但限制单项异常结果能够施加的影响。
这项工作的理论结论和 Amazon 应用结果目前都来自作者论文这一项信源,尚无独立交叉验证。
一个实验,为什么会影响其他实验?
经验贝叶斯可以理解成“先从历史中学习常态,再用常态修正单项结果”。普通贝叶斯分析先指定一个先验,也就是看数据前对结果分布的判断;经验贝叶斯则汇总一批历史实验,用它们估计先验参数,然后再回头修正每项实验。
这种修正叫作收缩:如果某项实验噪声很大,它的观测结果会被适度拉向总体平均水平,避免平台把偶然波动误当成真实效果。
论文讨论的正态—正态模型作了两层假设:各项实验的真实效果来自一个正态分布;平台看到的实验结果,又围绕各自的真实效果正态波动。前一个分布的方差很关键。方差大,意味着实验之间本来就可能相差很远,平台应当少收缩;方差小,意味着多数真实效果接近,平台就应更谨慎地对待单项实验的大幅波动。
问题是,真实效果看不见。平台只能从带噪声的观测中估计这个方差。常见做法是 maximum likelihood estimation,MLE(最大似然估计):选择一组参数,让已经看到的数据在模型下尽可能合理。
极端值会把“常态”改写
MLE 的薄弱处来自一个不对称性。论文推导显示,单项实验对先验均值估计的作用随偏差线性增长;对真实影响方差估计的作用却随偏差平方增长。结果越极端,它拉高方差的力量增长得越快,理论上没有上限。
这不只会让报表难看。先验方差决定平台怎样收缩各项实验,因此少数实验可能间接改变大量实验的效果估计,甚至影响上线决策。
最直觉的补救是删掉异常实验。但“异常”并不等于“错误”。一个极端结果既可能来自随机分组失败或指标配置问题,也可能代表真正很大的产品影响。大规模系统很难逐项查清,而且新的数据故障还会不断出现。先定一条阈值再删除,也会把判断变成带有任意性的人工规则;删哪些实验,本身就可能大幅改变结果。
不删除,而是给影响力设上限
BIMLE 采用稳健估计的思路。所谓稳健估计,就是让少量异常数据不至于主宰最终参数;“有界影响”则更直接:无论一项实验多极端,它都不能无限拉动先验。
具体来说,BIMLE 衡量每项结果与当前模型预测有多不一致。符合模型的结果权重为 1;偏离越远,权重越低。它随后用加权后的数据重新估计参数,再根据新参数更新偏离程度和权重,如此反复。论文称,实际计算通常不到 10 次迭代就会收敛。
关键的一步是,权重并非只看常见的检验统计量。相同大小的观测结果,对一个本来变化很大的指标可能很普通,对一个一贯稳定的指标却可能极不寻常。因此,BIMLE 会同时参照采样噪声和正在估计的先验方差。由于方差本身未知,它才需要在“估参数—改权重”之间迭代。
这里的“保留全部样本”需要准确理解:BIMLE 没有删除异常实验,但确实降低了它们的权重。它保留的是信息入口,不是原封不动的影响力。调节常数则控制取舍:数值较低时,方法更重视稳健;数值较高时,它对大效果更敏感。随着所有权重趋近 1,BIMLE 就退回标准 MLE。没有极端观测时,两者一致。
五个替换值,就能让方差膨胀十倍
论文先做了一项“污染实验”:从模型生成 1000 项测试,再把其中少数结果替换成极端值。只替换 5 项,也就是样本的 0.5%,MLE 估出的先验方差便膨胀到原来的十倍以上。
BIMLE 的变化小得多:替换数从 0 增至 5 时,估计方差大约从 1.0 上升到 3.0。它并未完全无视这些结果,而是在吸收一部分信号的同时,阻止五项实验决定整体方差。
论文还比较了 MAP(最大后验估计,即在似然之外加入对参数取值的先验约束)。即使作者给它配置了位置正确、信息较强的超先验,它也只提供了有限保护。原因是超先验只能对参数施加有限的拉力,极端观测对方差估计的作用却没有上限;结果足够极端时,数据仍会压过约束。
Amazon 数据里,极少数降权换来明显稳定
作者又把 MLE、MAP 和 BIMLE 用于 Amazon 的数千项 A/B 测试,考察四个指标,并按半年分别估计。每个指标得到六个时期的结果。
差异最醒目的是 Metric A。MLE 估计的先验标准差从 0.016 变到 0.900,跨期变异系数为 1.86;BIMLE 的范围为 0.016 至 0.026,变异系数降至 0.15。Metric C 中,MLE 的范围是 0.292 至 1.423,变异系数为 0.52;BIMLE 的范围为 0.209 至 0.377,变异系数为 0.23。
综合四个指标,论文称 BIMLE 将先验方差估计的不稳定性降低了 50% 以上。代价看起来很小:平均被降权的观测不到 0.1%,这些观测的平均权重在不同指标上介于 0.37 和 0.50。没有极端结果的时期,BIMLE 又与 MLE 相符。MAP 则几乎复制了 MLE 的结果,未显示出有意义的稳定作用。
这也是论文值得关注之处。它没有另换一套复杂模型,而是保留正态—正态模型容易大规模部署、收缩规则简单的特点,直接处理影响力无上限这个故障点。对生产系统来说,这是一种明确的工程取舍:平台不急着裁定极端结果是真是假,先阻止它在裁定前改写所有实验的“常态”。
局限与未知
- 稳定不等于更接近“真实方差”。极端结果可能包含真实的大效果;任何有界影响方法都会限制这部分信息。BIMLE是在响应极端值与保持稳定之间作取舍。
- 调节常数仍需人为选择,而且它的校准依赖正态分布的尾部概率。若真实效果分布并非正态,这种校准不会完全成立。
- Amazon 结果来自一次内部数据应用。论文没有给出“降低 50% 以上”的统一基线定义、总体样本规模或跨平台复现结果;不到 0.1% 也没有对应的绝对观测数。因此,现有证据不能直接推广到其他平台或其他分布。