fdrsample.R
library("fdrtool");
p.adj <- p.adjust(p.values, "fdr");
fdr <- fdrtool(p.values, statistic = "pvalue", plot = TRUE);
fdrs <- cbind(p.adj, fdr$qval, fdr$lfdr);
colnames(fdrs) <- c("adjusted p(BH)", "tail area-based Fdr", "density-based fdr(local fdr)");
write.table(fdrs, "fdrDemo.txt", quote = FALSE, row.names = FALSE, sep = "\t");
参考:
http://cran.r-project.org/web/packages/fdrtool/index.html
http://mametalk.blog32.fc2.com/blog-entry-28.html