Mapping the distribution of zero-dose children to assess the performance of vaccine delivery strategies and their relationships with measles incidence in Nigeria.

Geographically precise identification and targeting of populations at risk of vaccine-preventable diseases has gained renewed attention within the global health community over the last few years. District level estimates of vaccination coverage and corresponding zero-dose prevalence constitute a potentially useful evidence base to evaluate the performance of vaccination strategies. These estimates are also valuable for identifying missed communities, hence enabling targeted interventions and better resource allocation. Here, we fit Bayesian geostatistical models to map the routine coverage of the first doses of diphtheria-tetanus-pertussis vaccine (DTP1) and measles-containing vaccine (MCV1) and corresponding zero-dose estimates in Nigeria at 1x1 km resolution and the district level using geospatial data sets. We also map MCV1 coverage before and after the 2019 measles vaccination campaign in the northern states to further explore variations in routine vaccine coverage and to evaluate the effectiveness of both routine immunization (RI) and campaigns in reaching zero-dose children. Additionally, we map the spatial distributions of reported measles cases during 2018 to 2020 and explore their relationships with MCV zero-dose prevalence to highlight the public health implications of varying performance of vaccination strategies across the country. Our analysis revealed strong similarities between the spatial distributions of DTP and MCV zero dose prevalence, with districts with the highest prevalence concentrated mostly in the northwest and the northeast, but also in other areas such as Lagos state and the Federal Capital Territory. Although the 2019 campaign reduced MCV zero-dose prevalence substantially in the north, pockets of vulnerabilities remained in areas that had among the highest prevalence prior to the campaign. Importantly, we found strong correlations between measles case counts and MCV RI zero-dose estimates, which provides a strong indication that measles incidence in the country is mostly affected by RI coverage. Our analyses reveal an urgent and highly significant need to strengthen the country's RI program as a longer-term measure for disease control, whilst ensuring effective campaigns in the short term.

Geographic inequalities in health intervention coverage - mapping the composite coverage index in Peru using geospatial modelling.

BACKGROUND: The composite coverage index (CCI) provides an integrated perspective towards universal health coverage in the context of reproductive, maternal, newborn and child health. Given the sample design of most household surveys does not provide coverage estimates below the first administrative level, approaches for achieving more granular estimates are needed. We used a model-based geostatistical approach to estimate the CCI at multiple resolutions in Peru. METHODS: We generated estimates for the eight indicators on which the CCI is based for the departments, provinces, and areas of 5 × 5 km of Peru using data from two national household surveys carried out in 2018 and 2019 plus geospatial covariates. Bayesian geostatistical models were fit using the INLA-SPDE approach. We assessed model fit using cross-validation at the survey cluster level and by comparing modelled and direct survey estimates at the department-level. RESULTS: CCI coverage in the provinces along the coast was consistently higher than in the remainder of the country. Jungle areas in the north and east presented the lowest coverage levels and the largest gaps between and within provinces. The greatest inequalities were found, unsurprisingly, in the largest provinces where populations are scattered in jungle territory and are difficult to reach. CONCLUSIONS: Our study highlighted provinces with high levels of inequality in CCI coverage indicating areas, mostly low-populated jungle areas, where more attention is needed. We also uncovered other areas, such as the border with Bolivia, where coverage is lower than the coastal provinces and should receive increased efforts. More generally, our results make the case for high-resolution estimates to unveil geographic inequities otherwise hidden by the usual levels of survey representativeness.

Multilevel analysis of predictors of multiple indicators of childhood vaccination in Nigeria.

BACKGROUND: Substantial inequalities exist in childhood vaccination coverage levels. To increase vaccine uptake, factors that predict vaccination coverage in children should be identified and addressed. METHODS: Using data from the 2018 Nigeria Demographic and Health Survey and geospatial data sets, we fitted Bayesian multilevel binomial and multinomial logistic regression models to analyse independent predictors of three vaccination outcomes: receipt of the first dose of Pentavalent vaccine (containing diphtheria-tetanus-pertussis, Hemophilus influenzae type B and Hepatitis B vaccines) (PENTA1) (n = 6059) and receipt of the third dose having received the first (PENTA3/1) (n = 3937) in children aged 12-23 months, and receipt of measles vaccine (MV) (n = 11839) among children aged 12-35 months. RESULTS: Factors associated with vaccination were broadly similar for documented versus recall evidence of vaccination. Based on any evidence of vaccination, we found that health card/document ownership, receipt of vitamin A and maternal educational level were significantly associated with each outcome. Although the coverage of each vaccine dose was higher in urban than rural areas, urban residence was not significant in multivariable analyses that included travel time. Indicators relating to socio-economic status, as well as ethnic group, skilled birth attendance, lower travel time to the nearest health facility and problems seeking health care were significantly associated with both PENTA1 and MV. Maternal religion was related to PENTA1 and PENTA3/1 and maternal age related to MV and PENTA3/1; other significant variables were associated with one outcome each. Substantial residual community level variances in different strata were observed in the fitted models for each outcome. CONCLUSION: Our analysis has highlighted socio-demographic and health care access factors that affect not only beginning but completing the vaccination series in Nigeria. Other factors not measured by the DHS such as health service quality and community attitudes should also be investigated and addressed to tackle inequities in coverage.

Assessing the characteristics of un- and under-vaccinated children in low- and middle-income countries: A multi-level cross-sectional study.

Achieving equity in vaccination coverage has been a critical priority within the global health community. Despite increased efforts recently, certain populations still have a high proportion of un- and under-vaccinated children in many low- and middle-income countries (LMICs). These populations are often assumed to reside in remote-rural areas, urban slums and conflict-affected areas. Here, we investigate the effects of these key community-level factors, alongside a wide range of other individual, household and community level factors, on vaccination coverage. Using geospatial datasets, including cross-sectional data from the most recent Demographic and Health Surveys conducted between 2008 and 2018 in nine LMICs, we fitted Bayesian multi-level binary logistic regression models to determine key community-level and other factors significantly associated with non- and under-vaccination. We analyzed the odds of receipt of the first doses of diphtheria-tetanus-pertussis (DTP1) vaccine and measles-containing vaccine (MCV1), and receipt of all three recommended DTP doses (DTP3) independently, in children aged 12-23 months. In bivariate analyses, we found that remoteness increased the odds of non- and under-vaccination in nearly all the study countries. We also found evidence that living in conflict and urban slum areas reduced the odds of vaccination, but not in most cases as expected. However, the odds of vaccination were more likely to be lower in urban slums than formal urban areas. Our multivariate analyses revealed that the key community variables-remoteness, conflict and urban slum-were sometimes associated with non- and under-vaccination, but they were not frequently predictors of these outcomes after controlling for other factors. Individual and household factors such as maternal utilization of health services, maternal education and ethnicity, were more common predictors of vaccination. Reaching the Immunisation Agenda 2030 target of reducing the number of zero-dose children by 50% by 2030 will require country tailored analyses and strategies to identify and reach missed communities with reliable immunisation services.

Estimates of the number and distribution of zero-dose and under-immunised children across remote-rural, urban, and conflict-affected settings in low and middle-income countries.

While there has been great success in increasing the coverage of new childhood vaccines globally, expanding routine immunization to reliably reach all children and communities has proven more challenging in many low- and middle-income countries. Achieving this requires vaccination strategies and interventions that identify and target those unvaccinated, guided by the most current and detailed data regarding their size and spatial distribution. Through the integration and harmonisation of a range of geospatial data sets, including population, vaccination coverage, travel-time, settlement type, and conflict locations. We estimated the numbers of children un- or under-vaccinated for measles and diphtheria-tetanus-pertussis, within remote-rural, urban, and conflict-affected locations. We explored how these numbers vary both nationally and sub-nationally, and assessed what proportions of children these categories captured, for 99 lower- and middle-income countries, for which data was available. We found that substantial heterogeneities exist both between and within countries. Of the total 14,030,486 children unvaccinated for DTP1, over 11% (1,656,757) of un- or under-vaccinated children were in remote-rural areas, more than 28% (2,849,671 and 1,129,915) in urban and peri-urban areas, and up to 60% in other settings, with nearly 40% found to be within 1-hour of the nearest town or city (though outside of urban/peri-urban areas). Of the total number of those unvaccinated, we estimated between 6% and 15% (826,976 to 2,068,785) to be in conflict-affected locations, based on either broad or narrow definitions of conflict. Our estimates provide insights into the inequalities in vaccination coverage, with the distributions of those unvaccinated varying significantly by country, region, and district. We demonstrate the need for further inquiry and characterisation of those unvaccinated, the thresholds used to define these, and for more country-specific and targeted approaches to defining such populations in the strategies and interventions used to reach them.

A review of geospatial methods for population estimation and their use in constructing reproductive, maternal, newborn, child and adolescent health service indicators.

BACKGROUND: Household survey data are frequently used to measure reproductive, maternal, newborn, child and adolescent health (RMNCAH) service utilisation in low and middle income countries. However, these surveys are typically only undertaken every 5 years and tend to be representative of larger geographical administrative units. Investments in district health management information systems (DHMIS) have increased the capability of countries to collect continuous information on the provision of RMNCAH services at health facilities. However, reliable and recent data on population distributions and demographics at subnational levels necessary to construct RMNCAH coverage indicators are often missing. One solution is to use spatially disaggregated gridded datasets containing modelled estimates of population counts. Here, we provide an overview of various approaches to the production of gridded demographic datasets and outline their potential and their limitations. Further, we show how gridded population estimates can be used as alternative denominators to produce RMNCAH coverage metrics in combination with data from DHMIS, using childhood vaccination as examples. METHODS: We constructed indicators on the percentage of children one year old for diphtheria, pertussis and tetanus vaccine dose 3 (DTP3) and measles vaccine dose (MCV1) in Zambia and Nigeria at district levels. For the numerators, information on vaccines doses was obtained from each country's respective DHMIS. For the denominators, the number of children was obtained from 3 different sources including national population projections and aggregated gridded estimates derived using top-down and bottom-up geospatial methods. RESULTS: In Zambia, vaccination estimates utilising the bottom-up approach to population estimation substantially reduced the number of districts with > 100% coverage of DTP3 and MCV1 compared to estimates using population projection and the top-down method. In Nigeria, results were mixed with bottom-up estimates having a higher number of districts > 100% and estimates using population projections performing better particularly in the South. CONCLUSIONS: Gridded demographic data utilising traditional and novel data sources obtained from remote sensing offer new potential in the absence of up to date census information in the estimation of RMNCAH indicators. However, the usefulness of gridded demographic data is dependent on several factors including the availability and detail of input data.

District-level estimation of vaccination coverage: Discrete vs continuous spatial models.

Health and development indicators (HDIs) such as vaccination coverage are regularly measured in many low- and middle-income countries using household surveys, often due to the unreliability or incompleteness of routine data collection systems. Recently, the development of model-based approaches for producing subnational estimates of HDIs using survey data, particularly cluster-level data, has been an active area of research. This is mostly driven by the increasing demand for estimates at certain administrative levels, for example, districts, at which many development goals are set and evaluated. In this study, we explore spatial modeling approaches for producing district-level estimates of vaccination coverage. Specifically, we compare discrete spatial smoothing models which directly model district-level data with continuous Gaussian process (GP) models that utilize geolocated cluster-level data. We adopt a fully Bayesian framework, implemented using the INLA and SPDE approaches. We compare the predictive performance of the models by analyzing vaccination coverage using data from two Demographic and Health Surveys (DHS), namely the 2014 Kenya DHS and the 2015-16 Malawi DHS. We find that the continuous GP models performed well, offering a credible alternative to traditional discrete spatial smoothing models. Our analysis also revealed that accounting for between-cluster variation in the continuous GP models did not have any real effect on the district-level estimates. Our results provide guidance to practitioners on the reliability of these model-based approaches for producing estimates of vaccination coverage and other HDIs.

Geospatial estimation of reproductive, maternal, newborn and child health indicators: a systematic review of methodological aspects of studies based on household surveys.

BACKGROUND: Geospatial approaches are increasingly used to produce fine spatial scale estimates of reproductive, maternal, newborn and child health (RMNCH) indicators in low- and middle-income countries (LMICs). This study aims to describe important methodological aspects and specificities of geospatial approaches applied to RMNCH coverage and impact outcomes and enable non-specialist readers to critically evaluate and interpret these studies. METHODS: Two independent searches were carried out using Medline, Web of Science, Scopus, SCIELO and LILACS electronic databases. Studies based on survey data using geospatial approaches on RMNCH in LMICs were considered eligible. Studies whose outcomes were not measures of occurrence were excluded. RESULTS: We identified 82 studies focused on over 30 different RMNCH outcomes. Bayesian hierarchical models were the predominant modeling approach found in 62 studies. 5 × 5 km estimates were the most common resolution and the main source of information was Demographic and Health Surveys. Model validation was under reported, with the out-of-sample method being reported in only 56% of the studies and 13% of the studies did not present a single validation metric. Uncertainty assessment and reporting lacked standardization, and more than a quarter of the studies failed to report any uncertainty measure. CONCLUSIONS: The field of geospatial estimation focused on RMNCH outcomes is clearly expanding. However, despite the adoption of a standardized conceptual modeling framework for generating finer spatial scale estimates, methodological aspects such as model validation and uncertainty demand further attention as they are both essential in assisting the reader to evaluate the estimates that are being presented.

Spatial inequalities in skilled attendance at birth in Ghana: a multilevel analysis integrating health facility databases with household survey data.

OBJECTIVE: This study aimed at using survey data to predict skilled attendance at birth (SBA) across Ghana from healthcare quality and health facility accessibility. METHODS: Through a cross-sectional, observational study, we used a random intercept mixed effects multilevel logistic modelling approach to estimate the odds of having SBA and then applied model estimates to spatial layers to assess the probability of SBA at high-spatial resolution across Ghana. We combined data from the Demographic and Health Survey (DHS), routine birth registers, a service provision assessment of emergency obstetric care services, gridded population estimates and modelled travel time to health facilities. RESULTS: Within an hour's travel, 97.1% of women sampled in the DHS could access any health facility, 96.6% could reach a facility providing birthing services, and 86.2% could reach a secondary hospital. After controlling for characteristics of individual women, living in an urban area and close proximity to a health facility with high-quality services were significant positive determinants of SBA uptake. The estimated variance suggests significant effects of cluster and region on SBA as 7.1% of the residual variation in the propensity to use SBA is attributed to unobserved regional characteristics and 16.5% between clusters within regions. CONCLUSION: Given the expansion of primary care facilities in Ghana, this study suggests that higher quality healthcare services, as opposed to closer proximity of facilities to women, is needed to widen SBA uptake and improve maternal health.

OBJECTIF: Cette étude visait à utiliser les données d'enquête pour prédire l'assistance qualifiée à l'accouchement (AQA) à travers le Ghana à partir de la qualité des soins de santé et de l'accessibilité des établissements de santé. MÉTHODES: Grâce à une étude observationnelle transversale, nous avons utilisé une approche de modélisation logistique à multiniveau à effets mixtes d'interception aléatoire pour estimer les chances d'avoir une AQA, puis avons appliqué des estimations de modèle aux couches spatiales pour évaluer la probabilité d'AQA avec une résolution spatiale élevée à travers le Ghana. Nous avons combiné les données de l'Enquête démographique et de santé (EDS), les registres de naissance de routine, une évaluation de la prestation des services de soins obstétricaux d'urgence, des estimations démographiques quadrillées et un temps de trajet modélisé vers les établissements de santé. RÉSULTATS: En moins d'une heure de trajet, 97,1% des femmes échantillonnées dans l'EDS pouvaient accéder à un établissement de santé, 96,6% pouvaient atteindre un établissement fournissant des services d'accouchement et 86,2% pouvaient atteindre un hôpital secondaire. Après avoir ajusté pour les caractéristiques de chaque femme, le fait de vivre dans une zone urbaine et à proximité d'un établissement de santé offrant des services de haute qualité étaient des déterminants positifs significatifs de l'adoption de l'AQA. La variance estimée suggère des effets significatifs de regroupement et de la région sur l'AQA, car 7,1% de la variation résiduelle de la propension à utiliser l'AQA est attribuée à des caractéristiques régionales non observées et 16,5% entre les regroupements au sein des régions. CONCLUSION: Compte tenu de l'expansion des établissements de soins primaires au Ghana, cette étude suggère que des services de santé de meilleure qualité, par opposition à une plus grande proximité des établissements aux femmes, sont nécessaires pour élargir le recours à l'AQA et améliorer la santé maternelle.

Geospatial variation in measles vaccine coverage through routine and campaign strategies in Nigeria: Analysis of recent household surveys.

Measles vaccination campaigns are conducted regularly in many low- and middle-income countries to boost measles control efforts and accelerate progress towards elimination. National and sometimes first-level administrative division campaign coverage may be estimated through post-campaign coverage surveys (PCCS). However, these large-area estimates mask significant geographic inequities in coverage at more granular levels. Here, we undertake a geospatial analysis of the Nigeria 2017-18 PCCS data to produce coverage estimates at 1 × 1 km resolution and the district level using binomial spatial regression models built on a suite of geospatial covariates and implemented in a Bayesian framework via the INLA-SPDE approach. We investigate the individual and combined performance of the campaign and routine immunization (RI) by mapping various indicators of coverage for children aged 9-59 months. Additionally, we compare estimated coverage before the campaign at 1 × 1 km and the district level with predicted coverage maps produced using other surveys conducted in 2013 and 2016-17. Coverage during the campaign was generally higher and more homogeneous than RI coverage but geospatial differences in the campaign's reach of previously unvaccinated children are shown. Persistent areas of low coverage highlight the need for improved RI performance. The results can help to guide the conduct of future campaigns, improve vaccination monitoring and measles elimination efforts. Moreover, the approaches used here can be readily extended to other countries.

Mapping vaccination coverage to explore the effects of delivery mechanisms and inform vaccination strategies.

The success of vaccination programs depends largely on the mechanisms used in vaccine delivery. National immunization programs offer childhood vaccines through fixed and outreach services within the health system and often, additional supplementary immunization activities (SIAs) are undertaken to fill gaps and boost coverage. Here, we map predicted coverage at 1 × 1 km spatial resolution in five low- and middle-income countries to identify areas that are under-vaccinated via each delivery method using Demographic and Health Surveys data. We compare estimates of the coverage of the third dose of diphtheria-tetanus-pertussis-containing vaccine (DTP3), which is typically delivered through routine immunization (RI), with those of measles-containing vaccine (MCV) for which SIAs are also undertaken. We find that SIAs have boosted MCV coverage in some places, but not in others, particularly where RI had been deficient, as depicted by DTP coverage. The modelling approaches outlined here can help to guide geographical prioritization and strategy design.

A Bayesian latent process spatiotemporal regression model for areal count data.

Model-based approaches for the analysis of areal count data are commonplace in spatiotemporal analysis. In Bayesian hierarchical models, a latent process is incorporated in the mean function to account for dependence in space and time. Typically, the latent process is modelled using a conditional autoregressive (CAR) prior. The aim of this paper is to offer an alternative approach to CAR-based priors for modelling the latent process. The proposed approach is based on a spatiotemporal generalization of a latent process Poisson regression model developed in a time series setting. Spatiotemporal dependence in the autoregressive model for the latent process is modelled through its transition matrix, with a structured covariance matrix specified for its error term. The proposed model and its parameterizations are fitted in a Bayesian framework implemented via MCMC techniques. Our findings based on real-life examples show that the proposed approach is at least as effective as CAR-based models.

Geographic coverage of demographic surveillance systems for characterising the drivers of childhood mortality in sub-Saharan Africa.

A major focus of international health and development goals is the reduction of mortality rates in children under 5 years of age. Achieving this requires understanding the drivers of mortality and how they vary geographically to facilitate the targeting and prioritisation of appropriate interventions. Much of our knowledge on the causes of, and trends in, childhood mortality come from longitudinal demographic surveillance sites, with a renewed focus recently on the establishment and growth of networks of sites from which standardised outputs can facilitate broader understanding of processes. To ensure that the collective outputs from surveillance sites can be used to derive a comprehensive understanding and monitoring system for driving policy on tackling childhood mortality, confidence is needed that existing and planned networks of sites are providing a reliable and representative picture of the geographical variation in factors associated with mortality. Here, we assembled subnational data on childhood mortality as well as key factors known to be associated with it from household surveys in 27 sub-Saharan African countries. We then mapped the locations of existing longitudinal demographic surveillance sites to assess the extent of current coverage of the range of factors, identifying where gaps exist. The results highlight regions with unique combinations of factors associated with childhood mortality that are poorly represented by the current distribution of sites, such as southern Mali, central Nigeria and southern Zambia. Finally, we determined where the establishment of new surveillance systems could improve coverage.

High resolution age-structured mapping of childhood vaccination coverage in low and middle income countries.

BACKGROUND: The expansion of childhood vaccination programs in low and middle income countries has been a substantial public health success story. Indicators of the performance of intervention programmes such as coverage levels and numbers covered are typically measured through national statistics or at the scale of large regions due to survey design, administrative convenience or operational limitations. These mask heterogeneities and 'coldspots' of low coverage that may allow diseases to persist, even if overall coverage is high. Hence, to decrease inequities and accelerate progress towards disease elimination goals, fine-scale variation in coverage should be better characterized. METHODS: Using measles as an example, cluster-level Demographic and Health Surveys (DHS) data were used to map vaccination coverage at 1 km spatial resolution in Cambodia, Mozambique and Nigeria for varying age-group categories of children under five years, using Bayesian geostatistical techniques built on a suite of publicly available geospatial covariates and implemented via Markov Chain Monte Carlo (MCMC) methods. RESULTS: Measles vaccination coverage was found to be strongly predicted by just 4-5 covariates in geostatistical models, with remoteness consistently selected as a key variable. The output 1 × 1 km maps revealed significant heterogeneities within the three countries that were not captured using province-level summaries. Integration with population data showed that at the time of the surveys, few districts attained the 80% coverage, that is one component of the WHO Global Vaccine Action Plan 2020 targets. CONCLUSION: The elimination of vaccine-preventable diseases requires a strong evidence base to guide strategies and inform efficient use of limited resources. The approaches outlined here provide a route to moving beyond large area summaries of vaccination coverage that mask epidemiologically-important heterogeneities to detailed maps that capture subnational vulnerabilities. The output datasets are built on open data and methods, and in flexible format that can be aggregated to more operationally-relevant administrative unit levels.