This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to calculate q value in rna-seq data?

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): Wiki page

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.

rna-seq

1 answer

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

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

Log in to answer this question.