有没有人知道在BH程序下为FDR调整置信区间的好程序？

Question

我估计了Tajima D值的数据集的置信区间(CI)，使用bootstrapping来导出一个空分布。我想调整这些CI的帐户错误发现率(FDR)使用的程序，因为这允许数据点是非独立的。我已经可以对R包p.adjust中与每个数据点相关联的各个p值执行此操作。然而，我倾向于修改整个数据集的置信区间，而不是每个单独的p值。(这样做可以更快地在图上描述重要的和非重要的值)。

有没有人知道有一个程序能做到这一点(最好是用R或Python编写)？

FDR的详细资料载于: Benjamini和Yekutieli (2005年)。错误发现率-为选定的参数调整多个置信区间。美国统计协会杂志，100 (469)。

Benjamini和Yekutieli (2001年)。依赖条件下多重测试中错误发现率的控制。统计年鉴，29,1165-1188

Answer 1

查看R stats包中的stats(默认情况下已安装和加载)。有各种各样的调整方法，包括一种名为" by“的调整方法，它与你所提到的作者在2001年的引文相同。我想你需要在2005年和2001年的论文中检查这个方法是否相同。

由于p-值可以被“倒置”以生成CDF位置或标准的类似错误的值，因此您应该能够在该函数中使用显示在控制台会话中的现成代码：

p.adjust

..。然后拼凑一个CI来获取你正在按摩的样本统计数据。