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