edgeR: fold change reported by exactTest for zero values of rna-seq
0
0
Entering edit mode
@gordon-smyth
Last seen 1 hour ago
WEHI, Melbourne, Australia

Moderator: This post is in answer to: edgeR: fold change reported by exactTest for zero values of rna-seq

Dear Nick,

See ?exactTest and read the description of the prior.count argument.

For a little more detail, see ?predFC.

Note that division by zero is not undefined in R, it is just infinite (Inf).  If you run exactTest() with prior.count=0, you will get exactly that.

Best wishes
Gordon

 

> Date: Fri, 31 Jan 2014 14:09:23 +0000
> From: Nick N <feralmedic at="gmail.com=">
> To: bioconductor at r-project.org
> Subject: [BioC] edgeR: fold change reported by exactTest for zero values of rna-seq
>
> I have used the exact test in edgeR to compute the log fold changes. Here
> is the snippet:
>
> d <- DGEList(counts=counts, group=samples$Condition)
> d <- calcNormFactors(d)
> d <- estimateCommonDisp(d)
> d <- estimateTagwiseDisp(d)
> de <- exactTest(d)
>
> I've noticed that some genes have zero expression in all samples
> belonging to one of the two conditions. This would make the fold change
> mathematically undefined (division by zero). Yet the FC is reported as
> being ~2^-9. My question is - how does edgeR come up with this value?
> I've checked both the manual and the reference guide but couldn't figure
> out. There are various functions that accept pseudocounts as parameters
> but I have entered none in my snippet. So how does edgeR make up for the
> zero values in this particular case (which seems to be the default usage
> of the exactTest).
>
> Thanks!
>
> Nick

 

edgeR • 1.4k views
ADD COMMENT

Login before adding your answer.

Traffic: 325 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6