This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Access edgeR test statistic

I'd like to know an easy way to access test statistics calculated after using the exactTest() function in edgeR. I'd like to use these with fdrtool to attempt to estimate corrected FDR values due to a skewed p-value distribution. Currently, in the object produced by this function, I can only access: 

> colnames(res$table)
[1] "logFC"  "logCPM" "PValue" "FDR"

Thanks, R

edger
• 1,598 views • link
updated by Gordon Smyth Gordon Smyth • written by rrcutler rrcutler

1 answer

exactTest conducts the exact negative binomial test directly from the data, so the p-values are not computed via a test statistic.

topTags already computes FDR bounds using the Benjamini-Hochberg method. In this context I would recommend Benjamini-Hochberg over fdrtool anyway because it is more conservative and makes fewer assumptions.

• 0 views • link
written by Gordon Smyth Gordon Smyth

Log in to answer this question.

More posts like this