Replacing Normalizations with Interval AssumptionsImproves the Rigor and Robustness of DifferentialExpression and Differential Abundance Analyses
McGovern, K. C.; Silverman, J. D.
Show abstract
Standard methods for differential expression and differential abundance analysis rely on normalization to address sample-to-sample variation in sequencing depth. However, normalizations imply strict, unrealistic assumptions about the unmeasured scale of biological systems (e.g., microbial load or total cellular transcription). This introduces bias that can lead to false positives and false negatives. To overcome these limitations, we suggest replacing normalizations with interval assumptions. This approach allows researchers to explicitly define plausible lower and upper bounds on the unmeasured biological systems scale, making these assumptions more realistic, transparent, and flexible than those imposed by traditional normalizations. Compared to recent alternatives like scale models and sensitivity analyses, interval assumptions are easier to use, resulting in potentially reduced false positives and false negatives, and have stronger guarantees of Type-I error control. We make interval assumptions accessible by introducing a modified version of ALDEx2 as a publicly available software package. Through simulations and real data studies, we show these methods can reduce false positives and false negatives compared to normalization-based tools.
Matching journals
The top 2 journals account for 50% of the predicted probability mass.
Similar papers in this journal
Similar papers in this journal
- A negative binomial latent factor model for paired microbiome sequencing data 96%
- Feature selection and causal analysis for microbiome studies in the presence of confounding using standardization 95%
- Statistical Analysis of Variability in TnSeq Data Across Conditions Using Zero-Inflated Negative Binomial Regression 95%
"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.