How to calculate q value in rna-seq data?

Question

[--panda--]

In fact, it is quite a statistical question rather than biological question. After we compare the expression data and get the p value for different genes(iso forms, tss, or any other groups of interest), people sometimes need to control FDR(false discovery rate). Here is the method used in cufflinks to achieve this goal (a wiki page, not from the cufflinks paper):

I understand the steps used here, but I don't know how can the software come up with a q value based on this procedure.

Answer 1

[Nicolas Rosewick]

In R if p is a vector of pvalues then use : q <- p.adjust(p,method="fdr") or use the qvalue package :https://www.bioconductor.org/packages/release/bioc/html/qvalue.html

[--panda--]

Um, is there any description about the math under this transformation?

[Nicolas Rosewick]

In the R vignette you have all the references to the original papers : https://stat.ethz.ch/R-manual/R-devel/library/stats/html/p.adjust.html

And for the q-value package : J. D. Storey. False discovery rates. In Miodrag Lovric, editor, International Encyclopedia of Statistical Science. Springer, 2011. : http://genomine.org/papers/Storey_FDR_2011.pdf

[--panda--]

Thanks! That's really useful.