DESeq estimateSize factor, if one sample has 0
Entering edit mode
7.3 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 • 2.8k views
Entering edit mode
7.3 years ago
Martombo ★ 3.0k

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.


Login before adding your answer.

Traffic: 2506 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6