Please enable JavaScript to view this site.

对较低的P值进行解读并非易事

假设您正在筛选药物,以验证它们是否能降低血压。 您采用通常的P<0.05阈值作为统计学显著的标准。根据预期观察到的数据离散程度以及您所关注的最小变化量,您已为每次实验确定了样本量,以确保在P值小于0.05的情况下,具有80%的功效来检测您所寻找的差异。

如果确实得到了小于0.05的P值,那么该药物真正有效的概率是多少?

答案是:这取决于实验的具体背景。让我们先考虑这样一种情况:根据研究背景,您估计该药物真正有效的概率为10%。当您进行1000次实验时,会发生什么? 根据您10%的估计值,下方的两列总和分别为100和900。由于功效为80%,您预计80%的真正有效的药物在实验中会产生小于0.05的P值,因此左上角的单元格数值为80。 由于您将统计显著性的阈值设定为0.05,预计5%的无效药物在实验中会产生小于0.05的P值,因此右上角单元格的数值为45。  

 

药物确实有效

药物确实无效

总计

P<0.05,“显著”

80

45

125

P>0.05,“无统计学意义”

20

855

875

总计

100

900

1000

总而言之,您预计会有125项实验得出“统计学显著”的结果,而其中只有80项实验中该药物确实有效。 其余 45 项实验虽然得出了“统计学显著”的结果,但实际上是假阳性或错误发现。错误发现率(简称 FDR)为 45/125,即 36%。不是 5%,而是 36%。这也被称为假阳性率(FPR)。

下表摘自《生物统计学精要》第12章,展示了该情景及其他三种情景的FDR。


先验概率

P<0.05时的FDR

0.045 < P < 0.050 时的 FDR

在临床试验中比较治疗前随机分组的组别

0%

100%

100%

测试一种可能有效的药物

10%

36%

78%

测试一种成功概率为50:50的药物

50%

6%

27%

阳性对照组

100%

0%

0

 

上表中的每一行对应一种不同的情景,这些情景由不同的先验概率(即在收集数据之前)定义,即存在真实效应的概率。中间一列显示了如上所述计算出的预期FDR(也称为FPR)。 该列回答了以下问题:“如果P值小于0.05,那么实际上并不存在效应、结果仅仅是随机抽样所致的情况有多大概率?”请注意,该答案并非5%。FDR与α(用于定义统计学显著性的阈值P值)截然不同。

最右侧一列数据由模拟计算得出,基于 Colquhoun(1,2) 的研究,提出了一个略有不同的问题: “如果 P 值仅略低于 0.05(介于 0.045 和 0.050 之间),那么实际不存在效应、结果仅仅是随机抽样所致的概率是多少?”这些数值要高得多。 请关注第三行,其中先验概率为50%。在此情况下,如果P值仅略低于0.05,则该效应归因于随机性的概率为27%。 注意:是27%,不是5%!而在更具探索性的情境中,如果您认为先验概率为10%,那么对于刚好低于0.05的P值,错误发现率高达78%。这种情况下,一个(按常规定义的)统计学显著的结果几乎毫无意义。  

归根结底:您不能脱离具体情境来解释统计显著性(或P值)。您的解释依赖于实验的具体背景。错误发现率可能远高于α值(通常为5%)。解释结果需要常识、直觉和判断力。

参考文献

1.Colquhoun, D. (2014). 关于错误发现率与P值误解的研究。《皇家学会开放科学》,1(3),140216–140216。 http://doi.org/10.1098/rsos.140216

2.Colquhoun, D (2019). 《假阳性风险:关于如何处理p值的建议》。《美国统计学家》,第73卷,增刊1。