A Bayesian latent process spatiotemporal regression model for areal count data

Abstract 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

Geospatial modeling of child mortality across 27 countries in Sub-Saharan Africa

Preventable mortality of children has been targeted as one of the UN’s Sustainable Development Goals for the 2015-30 period. Global decreases in child mortality (4q1) have been seen, although sub-Saharan Africa remains an area of concern, with child mortality rates remaining high relative to global averages or even increasing in some cases. Furthermore, the spatial distribution of child mortality in sub-Saharan Africa is highly heterogeneous. Thus, research that identifies primary risk factors and protective measures in the geographic context of sub-Saharan Africa is needed. In this study, household survey data collected by The Demographic and Health Surveys (DHS) Program aggregated at DHS sub-national area scale are used to evaluate the spatial distribution of child mortality (age 1 to 4) across 27 sub-Saharan Africa countries in relation to a number of demographic and health indicators collected in the DHS surveys. In addition, this report controls for spatial variation in potential environmental drivers of child mortality by modeling it against a suite of geospatial datasets. These datasets vary across the study area in an autoregressive spatial model that accounts for the spatial autocorrelation present in the data. This study shows that socio-demographic factors such as birth interval, stunting, access to health facilities and literacy, along with geospatial factors such as prevalence of Plasmodium falciparum malaria, variety of ethnic groups, mean temperature, and intensity of lights at night can explain up to 60% of the variance in child mortality across 255 DHS sub-national areas in the 27 countries. Additionally, three regions - Western, Central, and Eastern Africa - have markedly different mortality rates. By identifying the relative importance of policy-relevant socio-demographic and environmental factors, this study highlights priorities for research and programs targeting child mortality over the next decade. <br/

Bayesian hierarchical modelling approaches for combining information from multiple data sources to produce annual estimates of national immunization coverage

Estimates of national immunization coverage are crucial for guiding policy and decision-making in national immunization programs and setting the global immunization agenda. 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 estimates of national immunization coverage and their associated uncertainties. We propose and explore two candidate models: a balanced data single likelihood (BDSL) 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 is the better model. We apply the methodology to produce coverage estimates for select vaccine-dose combinations for the period 2000-2019. A contributed R package {\tt imcover} implementing the No-U-Turn Sampler (NUTS) in the Stan programming language enhances the utility and reproducibility of the methodology.Comment: 31 pages (main), 4 figure

Mapping the timeliness of routine childhood vaccination in The Gambia: A spatial modelling study.

INTRODUCTION: Timeliness of routine vaccination shapes childhood infection risk and thus is an important public health metric. Estimates of indicators of the timeliness of vaccination are usually produced at the national or regional level, which may conceal epidemiologically relevant local heterogeneities and makeitdifficultto identify pockets of vulnerabilities that could benefit from targeted interventions. Here, we demonstrate the utility of geospatial modelling techniques in generating high-resolution maps of the prevalence of delayed childhood vaccination in The Gambia. To guide local immunisation policy and prioritize key interventions, we also identified the districts with a combination of high estimated prevalence and a significant population of affected infants. METHODS: We used the birth dose of the hepatitis-B vaccine (HepB0), third-dose of the pentavalent vaccine (PENTA3), and the first dose of measles-containing vaccine (MCV1) as examples to map delayed vaccination nationally at a resolution of 1 × 1-km2 pixel. We utilized cluster-level childhood vaccination data from The Gambia 2019-20 Demographic and Health Survey. We adopted a fully Bayesian geostatistical model incorporating publicly available geospatial covariates to aid predictive accuracy. The model was implemented using the integrated nested Laplace approximation-stochastic partial differential equation (INLA-SPDE) approach. RESULTS: We found significant subnational heterogeneity in delayed HepB0, PENTA3 and MCV1 vaccinations. Specificdistricts in the central and eastern regions of The Gambia consistentlyexhibited the highest prevalence of delayed vaccination, while the coastal districts showed alower prevalence forallthree vaccines. We also found that districts in the eastern, central, as well as in coastal parts of The Gambia had a combination of high estimated prevalence of delayed HepB0, PENTA3 and MCV1 and a significant population of affected infants. CONCLUSIONS: Our approach provides decision-makers with a valuable tool to better understand local patterns of untimely childhood vaccination and identify districts where strengthening vaccine delivery systems could have the greatest impact

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

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

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 undervaccinated 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

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

BackgroundGeospatial 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.MethodsTwo 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.ResultsWe 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.ConclusionsThe 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

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

BackgroundThe 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.MethodsUsing 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.ResultsMeasles 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.ConclusionThe 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

Unmet need for COVID-19 vaccination coverage in Kenya

COVID-19 has impacted the health and livelihoods of billions of people since it emerged in 2019. Vaccination for COVID-19 is a critical intervention that is being rolled out globally to end the pandemic. Understanding the spatial inequalities in vaccination coverage and access to vaccination centres is important for planning this intervention nationally. Here, COVID-19 vaccination data, representing the number of people given at least one dose of vaccine, a list of the approved vaccination sites, population data and ancillary GIS data were used to assess vaccination coverage, using Kenya as an example. Firstly, physical access was modelled using travel time to estimate the proportion of population within 1&#xA0;hour of a vaccination site. Secondly, a Bayesian conditional autoregressive (CAR) model was used to estimate the COVID-19 vaccination coverage and the same framework used to forecast coverage rates for the first quarter of 2022. Nationally, the average travel time to a designated COVID-19 vaccination site (n&#xA0;=&#xA0;622) was 75.5&#xA0;min (Range: 62.9 - 94.5&#xA0;min) and over 87% of the population&#xA0;&gt;18&#xA0;years reside within 1&#xA0;hour to a vaccination site. The COVID-19 vaccination coverage in December 2021 was 16.70% (95% CI: 16.66 - 16.74) - 4.4 million people and was forecasted to be 30.75% (95% CI: 25.04 - 36.96) - 8.1 million people by the end of March 2022. Approximately 21 million adults were still unvaccinated in December 2021 and, in the absence of accelerated vaccine uptake, over 17.2 million adults may not be vaccinated by end March 2022 nationally. Our results highlight geographic inequalities at sub-national level and are important in targeting and improving vaccination coverage in hard-to-reach populations. Similar mapping efforts could help other countries identify and increase vaccination coverage for such populations