我需要一些帮助来找到一种方法来计算 R 中的 Welch p 值的 q 值，使用置换 FDR 方法，就像 SAM 技术一样。
假设我的数据如下
x = matrix(rnorm(1000*3),nrow=1000)
y = matrix(rnorm(1000*3),nrow=1000)
我可以获得每个“行 x-vs-y”的韦尔奇 p 值列表，如下所示
pvalues = sapply(1:1000, function(n) t.test(x[n,], y[n,])$p.val)
然后如何使用置换 fdr 调整 p 值？据我所见， p.adjust不能这样做，samr只允许 t-statistic 或 wilcoxon ，而且我找不到任何可以帮助我的包。有什么建议么？