Bayesian Hierarchical Modeling Approaches for Combining Information From Multiple Data Sources to Produce Annual Estimates of National Immunization Coverage.
4 affiliations
- 1WorldPop, School of Geography and Environmental Science, University of Southampton, Southampton, UK.
- 2World Health Organization, Geneva, Switzerland.
- 3United Nations Children's Fund, New York, USA.
- 4Mathematical Sciences, University of Southampton, Southampton, UK.
Abstract
Estimates of national immunization coverage (ENIC) are crucial for guiding policy and decision-making in national immunization programs and setting the global immunization agenda. World Health Organization (WHO) and UNICEF estimates of national immunization coverage (WUENIC) are produced annually for various vaccine-dose combinations and all WHO Member States using information from multiple data sources and a deterministic computational logic approach. This approach, however, is incapable of characterizing the uncertainties inherent in coverage measurement and estimation. It also provides no statistically principled way of exploiting and accounting for the interdependence in immunization coverage data collected for multiple vaccines, countries, and time points. Here, we develop Bayesian hierarchical modeling approaches for producing accurate ENIC and associated uncertainties. We propose and explore two candidate models: a base single likelihood (BSL) model and an irregular data multiple likelihood (IDML) model, both of which differ in their handling of missing data and characterization of the uncertainties associated with the multiple input data sources. We provide a simulation study that demonstrates a high degree of accuracy of the estimates produced by the proposed models, and which also shows that the IDML model generally outperformed the BSL model in most cases we considered. We apply the methodology to produce coverage estimates for select vaccine-dose combinations for the period 2000-2019. A contributed R package, imcover, implementing the No-U-Turn Sampler in the Stan programming language enhances the utility and reproducibility of the methodology.
Related
- A Sample-Based Approach for Hierarchical Information-Theoretic Compression of Probabilistic Occupancy Grids
- Combining Exponential and Polynomial Convergence Rates for Accelerated Extremum Seeking: A High-Order Lie Bracket Approach
- Asymptotic Theory for Combining Dependent p-Values for Global Hypothesis Testing
- High-order positivity-preserving well-balanced schemes: combining flux globalization with modified Patankar deferred correction methods
- Enforced Gaplessness by Combining Non-Anomalous Symmetries
BibTeX
@article{utazi2026bayesian,
title = {{Bayesian Hierarchical Modeling Approaches for Combining Information From Multiple Data Sources to Produce Annual Estimates of National Immunization Coverage.}},
author = {C Edson Utazi and Warren C Jochem and M Carolina Danovaro-Holliday and Padraic Murphy and Sujit K Sahu and Marta Gacic-Dobo and Andrew J Tatem},
journal = {Biometrical journal. Biometrische Zeitschrift},
year = {2026},
month = oct,
doi = {10.1002/bimj.70189},
url = {https://europepmc.org/article/MED/42839684},
note = {CC BY},
}
RIS
TY - JOUR
TI - Bayesian Hierarchical Modeling Approaches for Combining Information From Multiple Data Sources to Produce Annual Estimates of National Immunization Coverage.
AU - C Edson Utazi
AU - Warren C Jochem
AU - M Carolina Danovaro-Holliday
AU - Padraic Murphy
AU - Sujit K Sahu
AU - Marta Gacic-Dobo
AU - Andrew J Tatem
PY - 2026
DA - 2026/10/01
JO - Biometrical journal. Biometrische Zeitschrift
DO - 10.1002/bimj.70189
UR - https://europepmc.org/article/MED/42839684
AB - Estimates of national immunization coverage (ENIC) are crucial for guiding policy and decision-making in national immunization programs and setting the global immunization agenda. World Health Organization (WHO) and UNICEF estimates of national immunization coverage (WUENIC) are produced annually for various vaccine-dose combinations and all WHO Member States using information from multiple data sources and a deterministic computational logic approach. This approach, however, is incapable of characterizing the uncertainties inherent in coverage measurement and estimation. It also provides no statistically principled way of exploiting and accounting for the interdependence in immunization coverage data collected for multiple vaccines, countries, and time points. Here, we develop Bayesian hierarchical modeling approaches for producing accurate ENIC and associated uncertainties. We propose and explore two candidate models: a base single likelihood (BSL) model and an irregular data multiple likelihood (IDML) model, both of which differ in their handling of missing data and characterization of the uncertainties associated with the multiple input data sources. We provide a simulation study that demonstrates a high degree of accuracy of the estimates produced by the proposed models, and which also shows that the IDML model generally outperformed the BSL model in most cases we considered. We apply the methodology to produce coverage estimates for select vaccine-dose combinations for the period 2000-2019. A contributed R package, imcover, implementing the No-U-Turn Sampler in the Stan programming language enhances the utility and reproducibility of the methodology.
SN - 0323-3847
SN - 1521-4036
ER -