Question: p-value and FDR
0
gravatar for raj_262920
2.9 years ago by
raj_26292020
Germany
raj_26292020 wrote:

Hi 

How would we explain the diff expression of a gene if it has a p-value and FDR absolutely equal to zero ?

Thank you and Best Regards

rna-seq cufflinks • 1.0k views
ADD COMMENTlink modified 2.9 years ago by fubar1.1k • written 2.9 years ago by raj_26292020
0
gravatar for fubar
2.9 years ago by
fubar1.1k
Australia
fubar1.1k wrote:

tl;dr A tiny p value allows you to confidently reject the null hypothesis of no difference. It doesn't really matter much how tiny. Once the FDR is below 0.05, you reject.

Ignoring the complexities of interpreting such a value, it depends. With (eg) R,  an architecture specific limit (eg 2.225074e-308 - see https://stat.ethz.ch/R-manual/R-devel/library/base/html/zMachine.html) defines the threshold below which you see a zero. It's a sign that the real "true" value is very, very small. 

Assuming you could calculate it correctly using infinite precision arithmetic, the probability that the true value for p from any biological experiment is exactly zero is itself vanishingly small. 

What are the raw counts for the contig in the two groups? The raw data is the key here.

If one group's raw counts for the contig are all consistently very low and the other group's are all consistently very high then you may have good grounds for asserting that the groups have different expression. 

 

 

 

ADD COMMENTlink modified 2.9 years ago • written 2.9 years ago by fubar1.1k
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 16.09
Traffic: 167 users visited in the last hour