Categories
Nevin Manimala Statistics

Pseudobulk with proper offsets has the same statistical properties as generalized linear mixed models in single-cell case-control studies

Bioinformatics. 2024 Aug 8:btae498. doi: 10.1093/bioinformatics/btae498. Online ahead of print.

ABSTRACT

MOTIVATION: Generalized linear mixed models (GLMMs), such as the negative-binomial or Poisson linear mixed model, are widely applied to single-cell RNA sequencing data to compare transcript expression between different conditions determined at the subject level. However, the model is computationally intensive, and its relative statistical performance to pseudobulk approaches is poorly understood.

RESULTS: We propose offset-pseudobulk as a lightweight alternative to GLMMs. We prove that a count-based pseudobulk equipped with a proper offset variable has the same statistical properties as GLMMs in terms of both point estimates and standard errors. We confirm our findings using simulations based on real data. Offset-pseudobulk is substantially faster (>x10) and numerically more stable than GLMMs.

AVAILABILITY: Offset pseudobulk can be easily implemented in any generalized linear model software by tweaking a few options. The codes can be found at https://github.com/hanbin973/pseudobulk_is_mm.

SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

PMID:39115884 | DOI:10.1093/bioinformatics/btae498

By Nevin Manimala

Portfolio Website for Nevin Manimala