你在近五十万个抽屉里找几张关键纸条。逐个检查并不难,难的是:检查次数一多,偶然看起来可疑的纸条也会成批出现。表观基因组研究面对的正是这种局面——研究者想判断吸烟是否先改变某处 DNA 甲基化,再影响肺功能,却必须同时筛查数十万个候选位点。
Haoyu Yang 等人在 Journal of the American Statistical Association(JASA)发表的工作,提出 tail likelihood ratio(TLR,尾部似然比)检验。它把注意力放在统计分布最极端的尾部,试图更准确地挑出显著性很强的中介信号,同时控制海量比较带来的误报。本文依据论文摘要及 JASA 公布的信息;模拟设置、对照方法和具体功效增幅尚未披露。
难点不只是“检验很多次”
中介分析追问的是一条两段式路径:某种暴露是否先影响中间环节,再由中间环节影响结果。放到这项研究里,就是吸烟是否经由某个 CpG 位点——DNA 甲基化测量所对应的位置——影响肺功能。
问题在于,中介效应由两段路径共同构成。只要其中任意一段为零,中介效应就不存在。因此,“没有中介效应”不是一种单一情形,而是由多种情况组成的复合原假设。再把这类检验重复近五十万次,微弱真信号与偶然出现的极端结果就更难分开。
TLR为什么盯住尾部
似然比衡量不同统计模型对已观察数据的支持程度。TLR则针对上述复合原假设,专门构造尾部检验统计量。直观地说,它不只问一个结果是否偏离常态,还更仔细地判断那些最靠近分布边缘、最可能被宣布为“发现”的结果。
这一步很重要,因为大规模筛选最终关心的是 false discovery rate(FDR,错误发现率):在所有被宣布为发现的位点中,预计有多大比例其实是误报。论文称,TLR在原假设比例、信号强度和信号分布变化时仍能保持稳健,并比现有方法取得更高的统计功效——也就是更有机会找出真实信号。不过,这些效果目前都来自作者团队的摘要陈述。
方法还有两个关键安排。第一,作者解析推导了TLR统计量在复合原假设下的分布,而不是完全依靠反复模拟来判断极端结果。第二,甲基化位点往往彼此相关,检验结果并非相互独立。作者用基于 principal component(主成分,把大量相关变化压缩成少数主要方向)的近似,描述检验统计量之间的依赖。
27个位点意味着什么
据JASA公布的信息,作者将TLR用于 Normative Aging Epigenetic Study,在632名男性的484,548个CpG位点中,识别出27个可能介导吸烟对肺功能影响的位点。
这个结果展示了方法的实际用途:面对近五十万个候选对象,它能把后续关注范围压缩到几十个。但“识别出”不等于已经证实因果机制。这27个位点来自一次特定队列分析,材料没有提供外部队列验证。
更值得注意的是理论保证的边界。论文只在各项检验相互独立时,从理论上建立了FDR控制;检验存在依赖时,证据来自模拟。主成分近似用于处理相关性,也不代表所有依赖结构都已有严格保证。
局限与未知
- 摘要没有披露具体模拟参数、比较对象、功效提升幅度和错误率数字,暂时无法判断优势有多大、适用于哪些边界情形。
- 相关检验下的FDR控制由模拟支持,尚非一般性的理论保证。
- 实例仅含632名男性,27个位点没有外部队列验证,对其他性别、年龄和人群的外推仍需谨慎。