多重检验——FWER,FDR,FCR

A/B实验设计——如何避免多重检验错误中,介绍了什么是多重检验,并且介绍了Bonferroni correctionHolm–Bonferroni method。本篇将更深入讨论应对多重检验的策略。

符号定义

  • m:总检验假设数
  • m0:零假设正确的数量,我们无法得知
  • m - m0:备择假设正确的数量
  • V:假阳性结论数量
  • S:真阳性数量
  • T:假阴性数量
  • U:真阴性数量
  • R = V + S:拒绝零假设数量

在m个假设检验中,m0个零假设为真,R是观察到的显著情况的随机变量,S、T、U、V都是不可观测的随机变量。

FWER(Family-wise error rate)

代表一系列检验中,至少存在一个假阳性的概率:
FWER = P_r(V >= 1)

控制过程

无论检验间是否独立的,\alpha \leq m * \alpha_{sub}都成立。
利用这个不等式,可以通过Bonferroni correctionHolm–Bonferroni method来对FWER进行控制。

缺点

由于FWER限制过于严格,会导致power相对比较低,容易错失正确的决策机会。
例如当两个比较是完全相关,多次比较并不会增加假阳性水平,但是矫正后却增加了假阴性。

FDR(False discovery rate)

发现中假阳性错误的比例期望:
Q = V / R,
FDR = Q_E = E[Q] = E[V/R|R>0]*P(R > 0)

控制过程

Benjamini–Hochberg procedure(BH step-up procedure)
  1. 将多重比较的P值排序,找到满足P_{(k)} \leq \frac{k}{m} {\alpha} 的最大的k
  2. 拒绝1 ~ k对应的原假设。

检验间独立或者正相关情况下,HB过程控制结果满足:
E(Q) \leq \frac{m_0}{m}\alpha \leq \alpha

如何理解?


  1. 设共有M个假设,M_0个零假设为真,它们的P值为均匀分布,显著水平为h,则期望的假阳性数量为h * M_0;
  2. 红线的斜率为\alpha / M,红线下方最大的P值对应的序号为L
  3. 拒绝零假设中,期望的假阳性数为h * M_0 = M_0\frac{\alpha * L}{M},因此:
    FDR \leq \frac{M_0\frac{\alpha * L}{M}}{L} = \alpha *{ M_0} / M \leq \alpha

HB过程在每次比较独立或者正相关时是有效的

Benjamini–Yekutieli procedure

此过程在任意情况下,都能控制假阳性。方式为在BH过程中,引入参数c,找到最大k满足P_{(k)} \le \frac{k}{m * {c(m)}}\alpha

  • 如果检验间独立或者正相关,c(m) = 1
  • 其他情况,c(m) = \sum _{i=1}^{m}\frac{1}{i}

缺点

相对于FWER,有较高的假阳性率。

FCR(False coverage rate)

错误覆盖率。FWER、FDR用于控制P值,FCR用于控制置信区间。

控制过程

BH过程对应的置信区间修正
  1. 将多重比较的P值排序,找到满足P_{(k)} \leq \frac{k}{m} {\alpha} 的最大的k
  2. 拒绝1 ~ k对应的原假设;
  3. 为每个比较中的参数,构建 1 - \frac{k}{m} {\alpha} 水平的置信区间。

结语

我认为需要根据不同情况,选择适合自己的标准与方法。
目的明确,决策指标明确,测试干预是否有效,也许FWER来进行控制比较合适;探索性实验,指标很多的情况下,可能FDR会更适合。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 214,128评论 6 493
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,316评论 3 388
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 159,737评论 0 349
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,283评论 1 287
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,384评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,458评论 1 292
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,467评论 3 412
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,251评论 0 269
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,688评论 1 306
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,980评论 2 328
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,155评论 1 342
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,818评论 4 337
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,492评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,142评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,382评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,020评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,044评论 2 352