This is a test version of Biostars. For the public version, visit https://www.biostars.org.
DESeq estimateSize factor, if one sample has 0

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

1 answer

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.

Log in to answer this question.