我正在测试 DNA 中的突变,我正在使用 Benjamini-Hochberg 方法来修改我正在测试我的 p 值的阈值。该方法基本上是对 p 值进行排序并将它们与定义的新阈值进行比较
,
其中是您的原始阈值,通常为 0.05,是您比较的 p 值的等级,是 p 值的总量(== 测试总量)。
由于我的数据的性质,许多 p 值是相同的。以下有序的 p 值集可能是一个示例:
= 0.01 = 0.03 = 0.03 = 0.03 = 0.09
假设 = 0.05 的阈值,我必须比较我的值是:
= 0.01 = 0.02 = 0.03 = 0.04 = 0.05
因此我接受第二个测试并拒绝第三个和第四个测试,即使它们具有相同的原始 p 值。我试图在文献中寻找解决方案,但没有任何结果。有没有既定的方法呢?如果没有,您首选的临时解决方案是什么?