Question: How to calculate q value in rna-seq data?
1
gravatar for --panda--
3.3 years ago by
--panda--30
Evanston
--panda--30 wrote:

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 • 2.6k views
ADD COMMENTlink modified 3.3 years ago by Nicolas Rosewick8.5k • written 3.3 years ago by --panda--30
2
gravatar for Nicolas Rosewick
3.3 years ago by
Belgium, Brussels
Nicolas Rosewick8.5k wrote:

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

ADD COMMENTlink modified 3.3 years ago • written 3.3 years ago by Nicolas Rosewick8.5k

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

ADD REPLYlink written 3.3 years ago by --panda--30

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

ADD REPLYlink modified 3.3 years ago • written 3.3 years ago by Nicolas Rosewick8.5k

Thanks! That's really useful.

ADD REPLYlink written 3.3 years ago by --panda--30
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1816 users visited in the last hour