Question: Statistical question (Deseq, Cuffdis) when one condition is zero?
1
gravatar for manekineko
3.7 years ago by
manekineko130
Bulgaria
manekineko130 wrote:

Hi,
I always am asking myself and can't answear then you have differential expression and one of the conditions are zero? Most of the time the software get division by zero and discard these results (or if the value is close to zero you get "inf" of something), is this correct?

I mean if on one condition the gene is 0 (not expressed) and in other 1000, you will get it discard it, but in general you can have switching on a gene, or probably, your sequencing, if you have it a bit more deep you will get 10 and 10000, and than will not be discarded, and get a nice log2, upregulation?

So my question is what is the right thing you do in such cases, and what is a biologically relevant/right (not statistically)?

ADD COMMENTlink modified 3.7 years ago by dariober10.0k • written 3.7 years ago by manekineko130

why would u discard it when you see Inf to -Inf ?

ADD REPLYlink modified 3.7 years ago • written 3.7 years ago by geek_y9.4k

Ok the Inf probably was bad example, but I'm more interested in the main case when one is zero?

ADD REPLYlink modified 3.7 years ago • written 3.7 years ago by manekineko130
1

Cuffdiff seems to discard things for no reason. I don't recall DESeq2 or edgeR doing that.

ADD REPLYlink written 3.7 years ago by Devon Ryan89k
2
gravatar for dariober
3.7 years ago by
dariober10.0k
WCIP | Glasgow | UK
dariober10.0k wrote:

A popular strategy to cope with zeros is to add a small number to all counts so that you avoid division by zero and at the same time you don't bias the results (e.g. 1000:0 is reasonably equivalent to 1001:1). Having said that, this is an issue that bugs me sometime when interpreting fold change ratios since small numbers can have a large effect which is not consistent with the biological interpretation. For example, if you add 1 to all your counts you could get log2(1001/1)= 9.97; if instead you add 0.1(biologically the same, I would argue) you get log2(1000.1/0.1)= 13.29, which is a big difference.

ADD COMMENTlink written 3.7 years ago by dariober10.0k
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: 979 users visited in the last hour