Abstract
Many vaccines are often administered in multiple doses to boost their effectiveness. In the case of childhood vaccines, the coverage maps of the doses and the differences between these often constitute an evidence base to guide investments in improving access to vaccination services and health system performance in low and middle-income countries. A major problem often encountered when mapping the coverage of multi-dose vaccines is the need to ensure that the coverage maps decrease monotonically with successive doses. That is, for doses i and j , i < j ⇒ pi(s) ≥ pj(s) , where pi(s) is the coverage of dose i at spatial location s . Here, we explore conditional probability (CP) and ratio-based (RB) approaches for mapping pi(s) , embedded within a binomial geostatistical modeling framework, to address this problem. The fully Bayesian model is implemented using the INLA and SPDE approaches. Using a simulation study, we find that both approaches perform comparably for out-of-sample estimation under varying point-level sample size distributions. We apply the methodology to map the coverage of the three doses of diphtheria-tetanus-pertussis vaccine using data from the 2018 Nigeria Demographic and Health Survey. The coverage maps produced using both approaches are almost indistinguishable, although the CP approach yielded more precise estimates on average in this application. We also provide estimates of zero-dose children and the dropout rates between the doses. The methodology is straightforward to implement and can be applied to other vaccines and geographical contexts.