Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: How to calculate q value in rna-seq data?
1
gravatar for --panda--
3.7 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.8k views
ADD COMMENT • link •
modified 3.7 years ago by Nicolas Rosewick ♦ 8.7k • written 3.7 years ago by --panda-- • 30
2
gravatar for Nicolas Rosewick
3.7 years ago by
Nicolas Rosewick ♦ 8.7k
Belgium, Brussels
Nicolas Rosewick ♦ 8.7k 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 COMMENT • link modified 3.7 years ago • written 3.7 years ago by Nicolas Rosewick ♦ 8.7k

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

ADD REPLY • link written 3.7 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 REPLY • link modified 3.6 years ago • written 3.6 years ago by Nicolas Rosewick ♦ 8.7k

Thanks! That's really useful.

ADD REPLY • link written 3.6 years ago by --panda-- • 30
Please log in to add an answer.

Similar posts • Search »

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