# Beta-binomial distribution

The **beta-binomial distribution** is a discrete probability distribution that describes the number of successes in a fixed number of binomial trials when the probability of success is itself random and drawn from a beta distribution. It is the standard model for binomial count data that show more variability than a single binomial distribution allows, a property called overdispersion.<sup>[1](https://docs.scipy.org/doc/scipy/tutorial/stats/discrete_betabinom.html)</sup>

| Key fact | Detail |
|---|---|
| Definition | A binomial distribution whose success probability p follows a beta distribution<sup>[1](https://docs.scipy.org/doc/scipy/tutorial/stats/discrete_betabinom.html)</sup> |
| Support | x = 0, 1, 2, ..., n, with shape parameters alpha, beta > 0<sup>[2](https://www.itl.nist.gov/div898/software/dataplot/refman2/auxillar/bbnpdf.htm)</sup> |
| Probability mass function | p(x, alpha, beta, n) = B(alpha + n − x, x + beta) / [(n + 1) B(n − x + 1, x + 1) B(alpha, beta)], where B is the complete beta function<sup>[2](https://www.itl.nist.gov/div898/software/dataplot/refman2/auxillar/bbnpdf.htm)</sup> |
| Origin | Derived by J. G. Skellam in 1948 by regarding the binomial probability of success as variable between sets of trials<sup>[3](https://ecommons.cornell.edu/bitstream/handle/1813/8678/TR000795.pdf?sequence=1)</sup> |
| Key parameters in practice | Two estimable parameters: the overall mean proportion and the intraclass correlation, which measures clustering among trials in the same group<sup>[4](https://academic.oup.com/jaoac/article/106/6/1629/7224238)</sup> |
| Typical uses | Microbiome relative-abundance modeling, validation studies of binary analytical methods, and rare-events meta-analysis<sup>[5](https://pmc.ncbi.nlm.nih.gov/articles/PMC7514055/)</sup><sup> • </sup><sup>[4](https://academic.oup.com/jaoac/article/106/6/1629/7224238)</sup><sup> • </sup><sup>[6](https://onlinelibrary.wiley.com/doi/full/10.1002/jrsm.1662)</sup> |
| Software | Available as `betabinom` in SciPy, defined for 0 ≤ k ≤ n<sup>[1](https://docs.scipy.org/doc/scipy/tutorial/stats/discrete_betabinom.html)</sup> |

## How the distribution arises

The construction is a two-stage randomization. First, a success probability p is drawn from a beta distribution, a continuous distribution on the interval from 0 to 1 with two shape parameters, often written alpha and beta, whose density involves the beta function B(alpha, beta) = ∫₀¹ t^(alpha−1)(1−t)^(beta−1) dt.<sup>[7](https://www.itl.nist.gov/div898/handbook/eda/section3/eda366h.htm)</sup> Second, given that p, the number of successes in n trials follows an ordinary binomial distribution. Averaging the binomial over the beta distribution of p produces the beta-binomial.<sup>[1](https://docs.scipy.org/doc/scipy/tutorial/stats/discrete_betabinom.html)</sup>

The resulting mass function is expressed with beta functions rather than binomial coefficients:<sup>[2](https://www.itl.nist.gov/div898/software/dataplot/refman2/auxillar/bbnpdf.htm)</sup>

\[ p(x, \alpha, \beta, n) = \frac{B(\alpha + n - x,\; x + \beta)}{(n+1)\, B(n - x + 1,\; x + 1)\, B(\alpha, \beta)}, \qquad x = 0, 1, \ldots, n, \quad \alpha, \beta > 0. \]

When the beta shape parameters are large, p is concentrated near its mean and the beta-binomial approaches an ordinary binomial distribution. When they are small, p varies widely between groups of trials, and the counts become more spread out than a binomial allows. This extra spread is exactly what the model is used to capture.

## History

The distribution is associated with the statistician J. G. Skellam, whose 1948 paper, *A Probability Distribution Derived from the Binomial Distribution by Regarding the Probability of Success as Variable Between the Sets of Trials*, derived expressions for the distribution of counts when p varies between sets of trials.<sup>[3](https://ecommons.cornell.edu/bitstream/handle/1813/8678/TR000795.pdf?sequence=1)</sup> Later surveys note that Skellam was not the first to consider the model, and a substantial literature on the beta-binomial model and its extensions accumulated through the following decades, including work on inference for the intraclass correlation in beta-binomial analysis of proportions.<sup>[8](https://www.researchgate.net/publication/238874723_A_Markov_Chain_Model_of_Extrabinomial_Variation)</sup><sup> • </sup><sup>[9](https://journals.sagepub.com/doi/10.3102/10769986006001003)</sup>

## Overdispersion and the intraclass correlation

The practical meaning of the beta-binomial's extra parameters is clearest through the intraclass correlation (ICC), the correlation between two individual binary observations from the same group. In the beta-binomial parameterization used for rare-events meta-analysis, the counts in each study follow Betabin(n, mu, psi), and the corresponding correlation parameter rho has this intraclass interpretation.<sup>[6](https://onlinelibrary.wiley.com/doi/full/10.1002/jrsm.1662)</sup>

In validation studies of analytical methods with binary responses, the beta-binomial is described as the simplest and most commonly used model of this family. Each collaborator's true success fraction (the true probability of detection) is presumed drawn at random from a beta distribution, and that collaborator's test-portion results then follow a binomial distribution. The model has only two estimable parameters: the overall mean proportion of "1"s across collaborators, called LPOD, and the common intraclass correlation. If the ICC estimate exceeds a specified performance criterion, the study results lack reproducibility and cannot be generalized to other sets of collaborators.<sup>[4](https://academic.oup.com/jaoac/article/106/6/1629/7224238)</sup>

## Applications

**Microbiome data.** [High-throughput sequencing](https://www.edgepedia.org/high-throughput-sequencing) produces count data for microbial taxa whose variability exceeds binomial expectations, and beta-binomial regression has been proposed specifically for modeling relative abundances and dysbiosis. The model supports tests not only for differential relative abundance between conditions but also for differential variability.<sup>[5](https://pmc.ncbi.nlm.nih.gov/articles/PMC7514055/)</sup><sup> • </sup><sup>[10](https://arxiv.org/abs/1902.02776)</sup> This line of work underlies methods such as CORNCOB, count regression for correlated observations with the beta-binomial, used for differential-abundance analysis with covariate adjustments and repeated measures.<sup>[11](https://doi.org/10.1038/s41592-023-02092-7)</sup> Extensions handle features of sequencing data directly, including zero-inflated beta-binomial models for microbiome count data.<sup>[12](https://pmc.ncbi.nlm.nih.gov/articles/PMC6124506/)</sup>

**Rare-events meta-analysis.** When combining studies in which the outcome of interest is uncommon, estimating a pooled effect from a small number of studies is difficult. In the Bayesian beta-binomial approach, the event counts in each study follow a beta-binomial distribution, and the model estimated by maximum likelihood performs well compared with alternative random-effects meta-analysis methods for rare events.<sup>[6](https://onlinelibrary.wiley.com/doi/full/10.1002/jrsm.1662)</sup>

**Method validation.** The beta-binomial serves as the reference distributional model in multi-laboratory validation studies of qualitative analytical methods, where the two-parameter LPOD-plus-ICC structure directly matches the quantities the study is designed to estimate.<sup>[4](https://academic.oup.com/jaoac/article/106/6/1629/7224238)</sup>

## Related models

The beta-binomial sits within a family of models for extrabinomial variation. The beta-negative binomial distribution, also called the generalized Waring distribution, is a related discrete distribution computed by the same statistical software tradition.<sup>[2](https://www.itl.nist.gov/div898/software/dataplot/refman2/auxillar/bbnpdf.htm)</sup> For microbiome compositional data, alternatives to beta-binomial regression include two-part beta regression models that separately capture zero inflation and skewness, and ratio-based pairwise methods that trade some of the beta-binomial's structure for computational speed in pairwise comparisons.<sup>[13](https://pubmed.ncbi.nlm.nih.gov/30036363/)</sup><sup> • </sup><sup>[14](https://pmc.ncbi.nlm.nih.gov/articles/PMC9853100/)</sup>

## References

1. [Beta-Binomial Distribution, SciPy v1.18.0 Manual](https://docs.scipy.org/doc/scipy/tutorial/stats/discrete_betabinom.html)
2. [Beta-Binomial Probability Mass Function, NIST/SEMATECH e-Handcase Dataplot reference](https://www.itl.nist.gov/div898/software/dataplot/refman2/auxillar/bbnpdf.htm)
3. [Skellam, J. G. (1948). A Probability Distribution Derived from the Binomial Distribution by Regarding the Probability of Success as Variable Between the Sets of Trials](https://ecommons.cornell.edu/bitstream/handle/1813/8678/TR000795.pdf?sequence=1)
4. [Beta-Binomial Statistical Model for Validation Studies of Analytes with a Binary Response, Journal of AOAC INTERNATIONAL](https://academic.oup.com/jaoac/article/106/6/1629/7224238)
5. [Modeling Microbial Abundances and Dysbiosis with Beta-Binomial Regression, PMC](https://pmc.ncbi.nlm.nih.gov/articles/PMC7514055/)
6. [Jansen et al. (2023). Rare events meta-analysis using the Bayesian beta-binomial model, Research Synthesis Methods](https://onlinelibrary.wiley.com/doi/full/10.1002/jrsm.1662)
7. [1.3.6.6.17. Beta Distribution, NIST/SEMATECH e-Handbook of Statistical Methods](https://www.itl.nist.gov/div898/handbook/eda/section3/eda366h.htm)
8. [A Markov Chain Model of Extrabinomial Variation](https://www.researchgate.net/publication/238874723_A_Markov_Chain_Model_of_Extrabinomial_Variation)
9. [Wilcox, R. R. (1981). A Review of the Beta-Binomial Model and its Extensions, Journal of Educational Statistics](https://journals.sagepub.com/doi/10.3102/10769986006001003)
10. [Modeling microbial abundances and dysbiosis with beta-binomial regression, arXiv](https://arxiv.org/abs/1902.02776)
11. [Multigroup analysis of compositions of microbiomes with covariate adjustments and repeated measures, Nature Methods](https://doi.org/10.1038/s41592-023-02092-7)
12. [A Zero-inflated Beta-binomial Model for Microbiome Data Analysis, PMC](https://pmc.ncbi.nlm.nih.gov/articles/PMC6124506/)
13. [A marginalized two-part Beta regression model for microbiome compositional data, PubMed](https://pubmed.ncbi.nlm.nih.gov/30036363/)
14. [Pairwise ratio-based differential abundance analysis of infant microbiome 16S sequencing data, PMC](https://pmc.ncbi.nlm.nih.gov/articles/PMC9853100/)

---
*Topic: Encyclopedia › Physical world and mathematics › Mathematics and statistics › Statistics and probability › Probability theory › Probability distributions › Compound, infinitely divisible, and convolved distributions › Compound distributions*

*Initially written Sep 17, 2026 · Reviewed: Sep 30, 2026; Oct 15, 2026 · Edited: Sep 30, 2026; Oct 15, 2026 · Last review: Oct 15, 2026*

*Copyright 2026 EdgeChat AI, a subsidiary of Biostate AI.*

License: Edgepedia Community License 1.0, https://www.edgepedia.org/edgepedia/license