DESeq estimateSize factor, if one sample has 0
1
2
Entering edit mode
8.5 years ago
tonja.r ▴ 600

The denominator of this expression can be interpreted as a pseudo-reference sample obtained by taking the geometric mean across samples. Thus, each size factor estimate sˆj is computed as the median of the ratios of the j-th sample's counts to those of the pseudo-reference.

But what if there are zero counts? The geometric mean will appear to be zero. And there is no devision by 0. Who does it handle genes with zero counts in some condition? Or the assumption is that if there are no reads on one condition, then the whole gene should be excluded?

rna-seq • 3.1k views
ADD COMMENT
3
Entering edit mode
8.5 years ago
Martombo ★ 3.1k

if a gene has 0 counts in one sample, it is excluded from this computation. this is not the ideal solution, but usually the number of genes that have some counts in all samples is so high that the estimation of the size factors is good enough.

ADD COMMENT

Login before adding your answer.

Traffic: 1692 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