Differential Expression from Counts: Dispersion, Fit and Shrinkage
Sequencing gives counts, not measurements, and counts carry a mean-variance relationship that a normal model gets wrong. This stop fits a count model, borrows information across genes to estimate dispersion, and reads the result honestly.
What a learner can do afterwards
- Explain why read counts need a negative binomial rather than a normal or Poisson model
- Say what library-size normalisation corrects and what composition effects it cannot fix
- Interpret shrunken fold changes and explain why low-count genes move the most
1 · Read
You treat sequencing output as counts, not measurements. Two samples differ in total reads, and replicates vary far more than coin flip sampling predicts. A normal model gets this wrong, and a Poisson model underfits the extra variance.
You scale counts before you compare anything. A sample sequenced twice as deeply shows roughly twice the counts for every gene, so size factors rescale using the bulk of moderately expressed genes. But scaling has limits: a few huge transcripts can distort every share, and then the constant majority assumption breaks. Diagnostic plots show when count distributions warp after scaling; suspect composition effects and distrust comparisons until the failure is fixed.
You see a tenfold rise built from three reads versus zero and feel excited. Then you check the counts and the uncertainty, and the excitement fades. The pipeline models counts with a negative binomial, borrows strength across genes to estimate dispersion, and shrinks your flimsy fold change toward zero.
You read shrunken fold changes as conservative statements of what the data support. Strongly expressed genes barely move, while low count genes move the most because their raw estimates carry the least information. A shrunken value near zero with wide uncertainty means no trustworthy evidence of change.
You scale for depth, you model counts with a negative binomial plus borrowed dispersion, and you trust shrunken changes over raw ratios.
2 · Watch
Take it off screen
Where it sits
8 questions wait behind this lesson, each with its answer explained. Every answer feeds the sky: stars light as they are learned, and dim when it is time to come back.