「你看到的效应量是被望远镜口径切出来的」——Malmquist 偏倚与元分析里的赢家诅咒
【锚点】这一篇只回答一个问题:为什么把「达到 p<0.05 的研究」汇总起来,平均效应量会系统性地大于真实值——以及 1922 年天文学里一个修正公式,为什么比 trim-and-fill 更诚实地承认「选择函数不可观测」。
一、望远镜口径如何制造虚假的「亮星」
【显式】Malmquist 偏倚(Malmquist 1922):在一个有流量极限的巡天里,越远的星系你只能看到越亮的。于是「被探测到的样本」的平均光度随距离系统性偏高,偏差随距离增大。修正需要知道选择函数(flux limit)。
【隐式】把「探测」换成「显著」:元分析里你只收录 p<0.05 的研究。这不是在精度(SE)上随机抽样,而是在 (效应量, SE) 的联合空间里沿一条对角线 d_hat/SE = 1.96 切割。小样本研究必须效应大才显著。于是被收录的研究偏向大效应——这就是【赢家诅咒】。
二、一个能算的数字
【细微】设真实 d=0.2,单研究 SE=0.1。显著(双侧)要求 |d_hat|>0.196。 - P(显著) ≈ P(d_hat>0.196) ≈ 0.516(负向显著几乎为 0) - 在显著且正向的研究里,E[d_hat | 显著] ≈ 0.2 + 0.1·φ(-0.04)/0.516 ≈ 0.277
【隐式】也就是说,只看显著研究,平均效应被高估约 39%。复现研究若用同样样本量、但不再筛选显著性,平均会回落到 0.2 附近——这不是「原研究造假」,而是【选择效应回吐】(winner’s curse 的回归)。
三、为什么 trim-and-fill 修不准
【细微】trim-and-fill 假设漏斗图的不对称来自「对精度(SE)的选择」,并据此把缺失研究镜像补上。但真实选择是对 p 值(=效应/SE)的选择——在 (效应, SE) 平面上是一条斜线,不是一条水平线。镜像补出来的「缺失研究」形状根本不对,所以它可能补多、可能补少。
【极难发现】所有选择修正都要求知道选择函数。Malmquist 知道自己的(flux limit 可测);元分析作者通常不知道(抽屉里的阴性研究不可观测)。trim-and-fill 假装能从对称性里「估」出选择函数——而这恰恰是它假设成立时才成立的东西。这是循环。
四、可执行的一点
【流程节点链路】真实效应 d=0.2 → 抽样 d_hat~N(0.2,0.1²) → 显著性筛选 |d_hat|>1.96·SE → 条件分布 E[d_hat|显著]≈0.28 → 汇总报告虚高效应 → 复现回归 0.2(被误读为「复现失败」)
【结论】不要把「汇总显著研究」当成总体估计;它是对一个被切窄的参考类的估计,而切法(选择函数)通常没被写进方法学部分。
结尾问题:如果选择函数不可观测,我们该继续报告「修正后的点估计」,还是改报「在给定选择函数假设下的效应量区间」?换句话说——发表偏倚修正,到底是在修正数据,还是在修正我们的无知?







0 评论