Provides tools to model overdispersed immune cell frequencies using Beta-Binomial generalized linear models.
Originally designed for memory B cells, the framework generalizes to evaluate rare
cellular immune responses across vaccinology, immuno-oncology, and related fields. Includes
method-of-moments dispersion estimation, likelihood ratio testing
across experimental arms, and Monte Carlo simulation to
calculate statistical power and Type I error rates. The methodologies
are directly motivated by the analysis of immunology data in
"Vaccination with mRNA-encoded membrane-anchored HIV envelope trimers
elicited tier 2 neutralizing antibodies in a phase 1 clinical trial"
(Parks et al. 2025)