A spatial regression model for the disaggregation of areal unit based data to high-resolution grids with application to vaccination coverage mapping

Reducing spatial heterogeneity in coverage improves the effectiveness of dog vaccination against rabies

Vaccination programs are the mainstay of control for many infectious diseases. Heterogeneous coverage is hypothesised to reduce vaccination effectiveness, but this impact has not been quantified in real systems. We address this gap using fine-scale data from two decades of rabies contact tracing and dog vaccination campaigns in Serengeti district, Tanzania. Using generalised linear mixed models, we find that current local (village-level) dog rabies incidence decreases with increasing recent local vaccination coverage. However, current local incidence is most dependent on recent incidence, both locally and in the wider district, consistent with high population connectivity. Removing the masking effects of prior non-local incidence shows that, for the same average prior vaccination coverage beyond the focal village, more spatial variation increases local incidence. These effects led to outbreaks following years when vaccination campaigns missed many villages, whereas when heterogeneity in coverage was reduced, incidence declined to low levels (<0.4 cases/1,000 dogs annually and no human deaths), such that short vaccination lapses thereafter did not lead to resurgence. We inferred ongoing rabies incursions into the district, suggesting regional connectivity as an important source of residual transmission. Overall, we provide an empirical demonstration of how the same average vaccination coverage can lead to differing outcomes based on its spatial distribution, highlighting the importance of fine-scale monitoring in managing vaccination programs.

Temperature, relative humidity and elderly type 2 diabetes mortality: A spatiotemporal analysis in Shandong, China

BACKGROUND

The mortality of type 2 diabetes mellitus (T2DM) can be affected by environmental factors. However, few studies have explored the effects of environmental factors across diverse regions over time. Given the vulnerability observed in the elderly group in previous research, this research applied Bayesian spatiotemporal models to assess the associations in the elderly group.

METHODS

Data on T2DM death in the elderly group (aged over 60 years old) at the county level were collected from the National Death Surveillance System between 1st January 2013 and 31st December 2019 in Shandong Province, China. A Bayesian spatiotemporal model was employed with the integrated Nested Laplace Approach to explore the associations between socio-environmental factors (i.e., temperatures, relative humidity, the Normalized Difference Vegetation Index (NDVI), particulate matter with a diameter of 2.5 μm or less (PM2.5) and gross domestic product (GDP)) and T2DM mortality.

RESULTS

T2DM mortality in the elderly group was found to be associated with temperature and relative humidity (i.e., temperature: Relative Risk (RR) = 1.41, 95% Credible Interval (CI): 1.27-1.56; relative humidity: RR = 1.05, 95% CI:1.03-1.06), while no significant associations were found with NDVI, PM2.5 and GDP. In winter, significant impacts from temperature (RR = 1.18, 95% CI: 1.06-1.32) and relative humidity (RR = 0.94, 95% CI: 0.89-0.99) were found. Structured and unstructured spatial effects, temporal trends and space-time interactions were considered in the model.

CONCLUSIONS

Higher mean temperatures and relative humidities increased the risk of elderly T2DM mortality in Shandong Province. However, a higher humidity level decreased the T2DM mortality risk in winter in Shandong Province. This research indicated that the spatiotemporal method could be a useful tool to assess the impact of socio-environmental factors on health by combining the spatial and temporal effects.

High-Resolution Geospatial Mapping of Zero-Dose and Underimmunized Children Following Nigeria's 2021 Multiple Indicator Cluster Survey/National Immunization Coverage Survey

BACKGROUND

"Zero-dose" children are those who are without any routine vaccination or are lacking the first dose of the diphtheria, tetanus, and pertussis-containing vaccine. Based on global estimates from the World Health Organization/United Nations Children's Fund in 2022, Nigeria has the highest number of zero-dose children, with >2.3 million unvaccinated.

METHODS

We used data from the 2021 Nigeria Multiple Indicator Cluster Survey/National Immunization Coverage Survey to identify zero-dose and underimmunized children. Geospatial modeling techniques were employed to determine the prevalence of zero-dose children and predict risk areas with underimmunized children at a high resolution (1 × 1 km).

RESULTS

Zero-dose and underimmunized children are more prevalent in socially deprived groups. Univariate and multivariate bayesian analyses showed positive correlations between the prevalence of zero-dose and underimmunized children and factors such as stunting, contraceptive prevalence, and literacy. The prevalence of zero-dose and underimmunized children varies significantly by region and ethnicity, with higher rates observed in the country's northern parts. Significant heterogeneity in the distribution of undervaccinated children was observed.

CONCLUSIONS

Nigeria needs to enhance its immunization system and coverage. Geospatial modeling can help deliver vaccines effectively to underserved communities. By adopting this approach, countries can ensure equitable vaccine access and contribute to global vaccination objectives.

High-resolution mapping of essential maternal and child health service coverage in Nigeria: a machine learning approach

BACKGROUND

National-level coverage estimates of maternal and child health (MCH) services mask district-level and community-level geographical inequities. The purpose of this study is to estimate grid-level coverage of essential MCH services in Nigeria using machine learning techniques.

METHODS

Essential MCH services in this study included antenatal care, facility-based delivery, childhood vaccinations and treatments of childhood illnesses. We estimated generalised additive models (GAMs) and gradient boosting regressions (GB) for each essential MCH service using data from five national representative cross-sectional surveys in Nigeria from 2003 to 2018 and geospatial socioeconomic, environmental and physical characteristics. Using the best-performed model for each service, we map predicted coverage at 1 km2 and 5 km2 spatial resolutions in urban and rural areas, respectively.

RESULTS

GAMs consistently outperformed GB models across a range of essential MCH services, demonstrating low systematic prediction errors. High-resolution maps revealed stark geographic disparities in MCH service coverage, especially between rural and urban areas and among different states and service types. Temporal trends indicated an overall increase in MCH service coverage from 2003 to 2018, although with variations by service type and location. Priority areas with lower coverage of both maternal and vaccination services were identified, mostly located in the northern parts of Nigeria.

CONCLUSION

High-resolution spatial estimates can guide geographic prioritisation and help develop better strategies for implementation plans, allowing limited resources to be targeted to areas with lower coverage of essential MCH services.

A scalable approach for short-term disease forecasting in high spatial resolution areal data

Short-term disease forecasting at specific discrete spatial resolutions has become a high-impact decision-support tool in health planning. However, when the number of areas is very large obtaining predictions can be computationally intensive or even unfeasible using standard spatiotemporal models. The purpose of this paper is to provide a method for short-term predictions in high-dimensional areal data based on a newly proposed "divide-and-conquer" approach. We assess the predictive performance of this method and other classical spatiotemporal models in a validation study that uses cancer mortality data for the 7907 municipalities of continental Spain. The new proposal outperforms traditional models in terms of mean absolute error, root mean square error, and interval score when forecasting cancer mortality 1, 2, and 3 years ahead. Models are implemented in a fully Bayesian framework using the well-known integrated nested Laplace estimation technique.

Bayesian transformation model for spatial partly interval-censored data

The transformation model with partly interval-censored data offers a highly flexible modeling framework that can simultaneously support multiple common survival models and a wide variety of censored data types. However, the real data may contain unexplained heterogeneity that cannot be entirely explained by covariates and may be brought on by a variety of unmeasured regional characteristics. Due to this, we introduce the conditionally autoregressive prior into the transformation model with partly interval-censored data and take the spatial frailty into account. An efficient Markov chain Monte Carlo method is proposed to handle the posterior sampling and model inference. The approach is simple to use and does not include any challenging Metropolis steps owing to four-stage data augmentation. Through several simulations, the suggested method's empirical performance is assessed and then the method is used in a leukemia study.

Area-based determinants of outreach vaccination for reaching vulnerable populations: A cross-sectional study in Pakistan

The objective of this study is to gain a comparative understanding of spatial determinants for outreach and clinic vaccination, which is critical for operationalizing efforts and breaking down structural biases; particularly relevant in countries where resources are low, and sub-region variance is high. Leveraging a massive effort to digitize public system reporting by Lady and Community Health Workers (CHWs) with geo-located data on over 4 million public-sector vaccinations from September 2017 through 2019, understanding health service operations in relation to vulnerable spatial determinants were made feasible. Location and type of vaccinations (clinic or outreach) were compared to regional spatial attributes where they were performed. Important spatial attributes were assessed using three modeling approaches (ridge regression, gradient boosting, and a generalized additive model). Consistent predictors for outreach, clinic, and proportion of third dose pentavalent vaccinations by region were identified. Of all Penta-3 vaccination records, 86.3% were performed by outreach efforts. At the tehsil level (fourth-order administrative unit), controlling for child population, population density, proportion of population in urban areas, distance to cities, average maternal education, and other relevant factors, increased poverty was significantly associated with more in-clinic vaccinations (β = 0.077), and lower proportion of outreach vaccinations by region (β = -0.083). Analyses at the union council level (fifth-administrative unit) showed consistent results for the differential importance of poverty for outreach versus clinic vaccination. Relevant predictors for each type of vaccination (outreach vs. in-clinic) show how design of outreach vaccination can effectively augment vaccination efforts beyond healthcare services through clinics. As Pakistan is third among countries with the most unvaccinated and under-vaccinated children, understanding barriers and factors associated with vaccination can be demonstrative for other national and sub-national regions facing challenges and also inform guidelines on supporting CHWs in health systems.

From regional to parcel scale: A high-resolution map of cover crops across Europe combining satellite data with statistical surveys

The reformed Common Agricultural Policy of 2023-2027 aims to promote a more sustainable and fair agricultural system in the European Union. Among the proposed measures, the incentivized adoption of cover crops to cover the soil during winter provides numerous benefits such as improved soil structure and reduced nutrient leaching and erosion. Despite this recognized importance, the availability of spatial data on cover crops is scarce. The increasing availability of field parcel declarations in the European Union has not yet filled this data gap due to its insufficient information content, limited public availability and a lack of standardization at continental scale. At present, the best information available is regionally aggregated survey data, which although indicative, hinders the development of spatially accurate studies. In this work, we propose a statistical model relating Sentinel-1 data to the existence of cover crops at the 100-m spatial resolution over the entirety of the European Union and United Kingdom and estimate its parameters using the spatially aggregated survey data. To validate the method in a spatially-explicit way, predictions were compared against farmers' registered declarations in France, where the adoption of cover crops is widespread. The results indicate a good agreement between predictions and parcel-level data. When interpreted as a binary classifier, the model yielded an Area Under the Curve (AUC) of 0.74 for the whole country. When the country was divided into five regions for the evaluation of regional biases, the AUC values were 0.77, 0.75, 0.74, 0.70, and 0.65 for the North, Center, West, East, and South regions respectively. Despite limitations such as the lack of data for validation outside France, and the non-standardized nomenclature for cover crops among Member States, this work constitutes the first effort to obtain a relevant cover crop map at a European scale for researchers and practitioners.

Generalized Geographically Weighted Regression Model within a Modularized Bayesian Framework

Geographically weighted regression (GWR) models handle geographical dependence through a spatially varying coefficient model and have been widely used in applied science, but its general Bayesian extension is unclear because it involves a weighted log-likelihood which does not imply a probability distribution on data. We present a Bayesian GWR model and show that its essence is dealing with partial misspecification of the model. Current modularized Bayesian inference models accommodate partial misspecification from a single component of the model. We extend these models to handle partial misspecification in more than one component of the model, as required for our Bayesian GWR model. Information from the various spatial locations is manipulated via a geographically weighted kernel and the optimal manipulation is chosen according to a Kullback-Leibler (KL) divergence. We justify the model via an information risk minimization approach and show the consistency of the proposed estimator in terms of a geographically weighted KL divergence.

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.

Bayesian disease mapping: Past, present, and future

On the occasion of the Spatial Statistics' 10th Anniversary, I reflect on the past and present of Bayesian disease mapping and look into its future. I focus on some key developments of models, and on recent evolution of multivariate and adaptive Gaussian Markov random fields and their impact and importance in disease mapping. I reflect on Bayesian disease mapping as a subject of spatial statistics that has advanced to date, and continues to grow, in scope and complexity alongside increasing needs of analytic tools for contemporary health science research, such as spatial epidemiology, population and public health, and medicine. I illustrate (potential) utility and impact of some of the disease mapping models and methods for analysing and monitoring communicable disease such as the COVID-19 infection risks during an ongoing pandemic.

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 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 = 622) was 75.5 min (Range: 62.9 - 94.5 min) and over 87% of the population >18 years reside within 1 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.

Spatiotemporal high-resolution prediction and mapping: methodology and application to dengue disease

Dengue disease has become a major public health problem. Accurate and precise identification, prediction and mapping of high-risk areas are crucial elements of an effective and efficient early warning system in countering the spread of dengue disease. In this paper, we present the fusion area-cell spatiotemporal generalized geoadditive-Gaussian Markov random field (FGG-GMRF) framework for joint estimation of an area-cell model, involving temporally varying coefficients, spatially and temporally structured and unstructured random effects, and spatiotemporal interaction of the random effects. The spatiotemporal Gaussian field is applied to determine the unobserved relative risk at cell level. It is transformed to a Gaussian Markov random field using the finite element method and the linear stochastic partial differential equation approach to solve the "big n" problem. Sub-area relative risk estimates are obtained as block averages of the cell outcomes within each sub-area boundary. The FGG-GMRF model is estimated by applying Bayesian Integrated Nested Laplace Approximation. In the application to Bandung city, Indonesia, we combine low-resolution area level (district) spatiotemporal data on population at risk and incidence and high-resolution cell level data on weather variables to obtain predictions of relative risk at subdistrict level. The predicted dengue relative risk at subdistrict level suggests significant fine-scale heterogeneities which are not apparent when examining the area level. The relative risk varies considerably across subdistricts and time, with the latter showing an increase in the period January-July and a decrease in the period August-December.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10109-021-00368-0.

Estimating and predicting snakebite risk in the Terai region of Nepal through a high-resolution geospatial and One Health approach

Most efforts to understand snakebite burden in Nepal have been localized to relatively small areas and focused on humans through epidemiological studies. We present the outcomes of a geospatial analysis of the factors influencing snakebite risk in humans and animals, based on both a national-scale multi-cluster random survey and, environmental, climatic, and socio-economic gridded data for the Terai region of Nepal. The resulting Integrated Nested Laplace Approximation models highlight the importance of poverty as a fundamental risk-increasing factor, augmenting the snakebite odds in humans by 63.9 times. For animals, the minimum temperature of the coldest month was the most influential covariate, increasing the snakebite odds 23.4 times. Several risk hotspots were identified along the Terai, helping to visualize at multiple administrative levels the estimated population numbers exposed to different probability risk thresholds in 1 year. These analyses and findings could be replicable in other countries and for other diseases.

Estimating Health over Space and Time: A Review of Spatial Microsimulation Applied to Public Health

There is an ongoing demand for data on population health, for reasons of resource allocation, future planning and crucially to address inequalities in health between people and between populations. Although there are regular sources of data at coarse spatial scales, such as countries or large sub-national units such as states, there is often a lack of good quality health data at the local level. One method to develop reliable estimates of population health outcomes is spatial microsimulation, an approach that has its roots in economic studies. Here, we share a review of this method for estimating health in populations, explaining the different approaches available and examples where the method is applied successfully for creating both static and dynamic populations. Recent notable advances in the method that allow uncertainty to be represented are highlighted, along with the evolving approaches to validation that are an ongoing challenge in small-area estimation. The summary serves as a primer for academics new to the area of research as well as an overview for non-academic researchers who consider using these models for policy evaluations.

Modeling and presentation of vaccination coverage estimates using data from household surveys

It is becoming increasingly popular to produce high-resolution maps of vaccination coverage by fitting Bayesian geostatistical models to data from household surveys. Usually, the surveys adopt a stratified cluster sampling design. We discuss a number of crucial choices with respect to two key aspects of the map production process: the acknowledgement of the survey design in modeling, and the appropriate presentation of estimates and their uncertainties. Specifically, we consider the importance of accounting for urban/rural stratification and cluster-level non-spatial excess variation in survey outcomes, when fitting geostatistical models. We also discuss the trade-off between the geographical scale and precision of model-based estimates, and demonstrate visualization methods for mapping and ranking that emphasize the probabilistic interpretation of results. A novel approach to coverage map presentation is proposed to allow comparison and control of the overall map uncertainty. We use measles vaccination coverage in Nigeria as a motivating example and illustrate the different issues using data from the 2018 Nigeria Demographic and Health Survey.

Harmonizing child mortality data at disparate geographic levels

There is an increasing focus on reducing inequalities in health outcomes in developing countries. Subnational variation is of particular interest, with geographically-indexed data being used to understand the spatial risk of detrimental outcomes and to identify who is at greatest risk. While some health surveys provide observations with associated geographic coordinates (point data), many others provide data that have their locations masked and instead only report the strata (polygon information) within which the data resides (masked data). How to harmonize these data sources for spatial analysis has been previously considered although only ad hoc methods and comparison of methods is lacking. In this paper, we present a new method for analyzing masked survey data, using a method that is consistent with the data-generating process. In addition, we critique two previously proposed approaches to analyzing masked data and illustrate that they are fundamentally flawed methodologically. To validate our method, we compare our approach with previously formulated solutions in several realistic simulation environments in which the underlying structure of the risk field is known. We simulate samples from spatiotemporal fields in a way that mimics the sampling frame implemented in the most common health surveys in low- and middle-income countries, the Demographic and Health Surveys and Multiple Indicator Cluster Surveys. In simulations, the newly proposed approach outperforms previously proposed approaches in terms of minimizing error while increasing the precision of estimates. The approaches are subsequently compared using child mortality data from the Dominican Republic where our findings are reinforced. The ability to accurately increase precision of child mortality estimates, and health outcomes in general, by leveraging various types of data, improves our ability to implement precision public health initiatives and better understand the landscape of geographic health inequalities.

Model-based spatial-temporal mapping of opisthorchiasis in endemic countries of Southeast Asia

Opisthorchiasis is an overlooked danger to Southeast Asia. High-resolution disease risk maps are critical but have not been available for Southeast Asia. Georeferenced disease data and potential influencing factor data were collected through a systematic review of literatures and open-access databases, respectively. Bayesian spatial-temporal joint models were developed to analyze both point- and area-level disease data, within a logit regression in combination of potential influencing factors and spatial-temporal random effects. The model-based risk mapping identified areas of low, moderate, and high prevalence across the study region. Even though the overall population-adjusted estimated prevalence presented a trend down, a total of 12.39 million (95% Bayesian credible intervals [BCI]: 10.10–15.06) people were estimated to be infected with O. viverrini in 2018 in four major endemic countries (i.e., Thailand, Laos, Cambodia, and Vietnam), highlighting the public health importance of the disease in the study region. The high-resolution risk maps provide valuable information for spatial targeting of opisthorchiasis control interventions.

Spatio-temporal analysis of misaligned burden of disease data using a geo-statistical approach

Data used to estimate the burden of diseases (BOD) are usually sparse, noisy, and heterogeneous. These data are collected from surveys, registries, and systematic reviews that have different areal units, are conducted at different times, and are reported for different age groups. In this study, we developed a Bayesian geo-statistical model to combine aggregated sparse, noisy BOD data from different sources with misaligned areal units. Our model incorporates the correlation of space, time, and age to estimate health indicators for areas with no data or a small number of observations. The model also considers the heterogeneity of data sources and the measurement errors of input data in the final estimates and uncertainty intervals. We applied the model to combine data from nine different sources of body mass index in a national and sub-national BOD study. The cross-validation results confirmed a high out-of-sample predictive ability in sparse and noisy data. The proposed model can be used by other BOD studies especially at the sub-national level when the areal units are subject to misalignment.

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.

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.

Estimating uncertainty in geospatial modelling at multiple spatial resolutions: the pattern of delivery via caesarean section in Tanzania

Visualising maternal and newborn health (MNH) outcomes at fine spatial resolutions is crucial to ensuring the most vulnerable women and children are not left behind in improving health. Disaggregated data on life-saving MNH interventions remain difficult to obtain, however, necessitating the use of Bayesian geostatistical models to map outcomes at small geographical areas. While these methods have improved model parameter estimates and precision among spatially correlated health outcomes and allowed for the quantification of uncertainty, few studies have examined the trade-off between higher spatial resolution modelling and how associated uncertainty propagates. Here, we explored the trade-off between model outcomes and associated uncertainty at increasing spatial resolutions by quantifying the posterior distribution of delivery via caesarean section (c-section) in Tanzania. Overall, in modelling delivery via c-section at multiple spatial resolutions, we demonstrated poverty to be negatively correlated across spatial resolutions, suggesting important disparities in obtaining life-saving obstetric surgery persist across sociodemographic factors. Lastly, we found that while uncertainty increased with higher spatial resolution input, model precision was best approximated at the highest spatial resolution, suggesting an important policy trade-off between identifying concealed spatial heterogeneities in health indicators.

Using models to shape measles control and elimination strategies in low- and middle-income countries: A review of recent applications

After many decades of vaccination, measles epidemiology varies greatly between and within countries. National immunization programs are therefore encouraged to conduct regular situation analyses and to leverage models to adapt interventions to local needs. Here, we review applications of models to develop locally tailored interventions to support control and elimination efforts. In general, statistical and semi-mechanistic transmission models can be used to synthesize information from vaccination coverage, measles incidence, demographic, and/or serological data, offering a means to estimate the spatial and age-specific distribution of measles susceptibility. These estimates complete the picture provided by vaccination coverage alone, by accounting for natural immunity. Dynamic transmission models can then be used to evaluate the relative impact of candidate interventions for measles control and elimination and the expected future epidemiology. In most countries, models predict substantial numbers of susceptible individuals outside the age range of routine vaccination, which affects outbreak risk and necessitates additional intervention to achieve elimination. More effective use of models to inform both vaccination program planning and evaluation requires the development of training to enhance broader understanding of models and where feasible, building capacity for modelling in-country, pipelines for rapid evaluation of model predictions using surveillance data, and clear protocols for incorporating model results into decision-making.

Downscaling livestock census data using multivariate predictive models: Sensitivity to modifiable areal unit problem

The analysis of census data aggregated by administrative units introduces a statistical bias known as the modifiable areal unit problem (MAUP). Previous researches have mostly assessed the effect of MAUP on upscaling models. The present study contributes to clarify the effects of MAUP on the downscaling methodologies, highlighting how a priori choices of scales and shapes could influence the results. We aggregated chicken and duck fine-resolution census in Thailand, using three administrative census levels in regular and irregular shapes. We then disaggregated the data within the Gridded Livestock of the World analytical framework, sampling predictors in two different ways. A sensitivity analysis on Pearson's r correlation statistics and RMSE was carried out to understand how size and shapes of the response variables affect the goodness-of-fit and downscaling performances. We showed that scale, rather than shapes and sampling methods, affected downscaling precision, suggesting that training the model using the finest administrative level available is preferable. Moreover, datasets showing non-homogeneous distribution but instead spatial clustering seemed less affected by MAUP, yielding higher Pearson's r values and lower RMSE compared to a more spatially homogenous dataset. Implementing aggregation sensitivity analysis in spatial studies could help to interpret complex results and disseminate robust products.

Modelling the Wealth Index of Demographic and Health Surveys within Cities Using Very High-Resolution Remotely Sensed Information

A systematic and precise understanding of urban socio-economic spatial inequalities in developing regions is needed to address global sustainability goals. At the intra-urban scale, access to detailed databases (i.e., a census) is often a difficult exercise. Geolocated surveys such as the Demographic and Health Surveys (DHS) are a rich alternative source of such information but can be challenging to interpolate at such a fine scale due to their spatial displacement, survey design and the lack of very high-resolution (VHR) predictor variables in these regions. In this paper, we employ satellite-derived VHR land-use/land-cover (LULC) datasets and couple them with the DHS Wealth Index (WI), a robust household wealth indicator, in order to provide city-scale wealth maps. We undertake several modelling approaches using a random forest regressor as the underlying algorithm and predict in several geographic administrative scales. We validate against an exhaustive census database available for the city of Dakar, Senegal. Our results show that the WI was modelled to a satisfactory degree when compared against census data even at very fine resolutions. These findings might assist local authorities and stakeholders in rigorous evidence-based decision making and facilitate the allocation of resources towards the most disadvantaged populations. Good practices for further developments are discussed with the aim of upscaling these findings at the global scale.

Geospatial Disaggregation of Population Data in Supporting SDG Assessments: A Case Study from Deqing County, China

Quantitative assessments and dynamic monitoring of indicators based on fine-scale population data are necessary to support the implementation of the United Nations (UN) 2030 Agenda and to comprehensively achieve its 17 Sustainable Development Goals (SDGs). However, most population data are collected by administrative units, and it is difficult to reflect true distribution and uniformity in space. To solve this problem, based on fine building information, a geospatial disaggregation method of population data for supporting SDG assessments is presented in this paper. First, Deqing County in China, which was divided into residential areas and nonresidential areas according to the idea of dasymetric mapping, was selected as the study area. Then, the town administrative areas were taken as control units, building area and number of floors were used as weighting factors to establish the disaggregation model, and population data with a resolution of 30 m in Deqing County in 2016 were obtained. After analyzing the statistical population of 160 villages and the disaggregation results, we found that the global average accuracy was 87.08%. Finally, by using the disaggregation population data, indicators 3.8.1, 4.a.1, and 9.1.1 were selected to conduct an accessibility analysis and a buffer analysis in a quantitative assessment of the SDGs. The results showed that the SDG measurement and assessment results based on the disaggregated population data were more accurate and effective than the results obtained using the traditional method.