R 中 fuzzySim 包中的错误发现率 (FDR)

False Discovery Rate (FDR) in fuzzySim package in R

我是 R 的新手,我正在尝试计算我 运行 的一系列 t 检验的错误发现率,我希望 q = 0.05。 FDR function in the fuzzySim package 似乎是一个不错的解决方案,但是当我 运行 fuzzySim 描述中的示例时,我收到了一些错误消息。例如我做...

library(fuzzySim)

df <- data.frame(var = letters[1:5], pval = c(0.02, 0.004, 0.07, 0.03, 0.05))

p_value_info <- FDR(pvalues = df, correction = "fdr", q = 0.05, 
verbose = FALSE, simplif = TRUE)

不幸的是,这是我收到的错误消息...

Error in if (na.loss > 0) message(na.loss, " cases excluded due to missing or non-finite values.") : 
  argument is of length zero

任何关于我做错的想法都将不胜感激,谢谢!

显然在使用这个 fuzzySim 包时,代码中存在一个错误,现在已在可以安装的最新开发版本 (3.5) 中修复

install.packages("fuzzySim", repos="http://R-Forge.R-project.org")

我已经确认它现在可以使用变量名称和 p 值的数据框。

非常感谢此包的作者 A. Márcia Barbosa 更新了代码。 Márcia 告诉我,今年晚些时候将更新包上传到 CRAN。