Geostatistical modelling of the association between malaria and child growth in Africa

Identifying childhood malaria hotspots and risk factors in a Nigerian city using geostatistical modelling approach

Malaria ranks high among prevalent and ravaging infectious diseases in sub-Saharan Africa (SSA). The negative impacts, disease burden, and risk are higher among children and pregnant women as part of the most vulnerable groups to malaria in Nigeria. However, the burden of malaria is not even in space and time. This study explores the spatial variability of malaria prevalence among children under five years (U5) in medium-sized rapidly growing city of Akure, Nigeria using model-based geostatistical modeling (MBG) technique to predict U5 malaria burden at a 100 × 100 m grid, while the parameter estimation was done using Monte Carlo maximum likelihood method. The non-spatial logistic regression model shows that U5 malaria prevalence is significantly influenced by the usage of insecticide-treated nets-ITNs, window protection, and water source. Furthermore, the MBG model shows predicted U5 malaria prevalence in Akure is greater than 35% at certain locations while we were able to ascertain places with U5 prevalence > 10% (i.e. hotspots) using exceedance probability modelling which is a vital tool for policy development. The map provides place-based evidence on the spatial variation of U5 malaria in Akure, and direction on where intensified interventions are crucial for the reduction of U5 malaria burden and improvement of urban health in Akure, Nigeria.

Comparison of new computational methods for spatial modelling of malaria

BACKGROUND

Geostatistical analysis of health data is increasingly used to model spatial variation in malaria prevalence, burden, and other metrics. Traditional inference methods for geostatistical modelling are notoriously computationally intensive, motivating the development of newer, approximate methods for geostatistical analysis or, more broadly, computational modelling of spatial processes. The appeal of faster methods is particularly great as the size of the region and number of spatial locations being modelled increases.

METHODS

This work presents an applied comparison of four proposed 'fast' computational methods for spatial modelling and the software provided to implement them-Integrated Nested Laplace Approximation (INLA), tree boosting with Gaussian processes and mixed effect models (GPBoost), Fixed Rank Kriging (FRK) and Spatial Random Forests (SpRF). The four methods are illustrated by estimating malaria prevalence on two different spatial scales-country and continent. The performance of the four methods is compared on these data in terms of accuracy, computation time, and ease of implementation.

RESULTS

Two of these methods-SpRF and GPBoost-do not scale well as the data size increases, and so are likely to be infeasible for larger-scale analysis problems. The two remaining methods-INLA and FRK-do scale well computationally, however the resulting model fits are very sensitive to the user's modelling assumptions and parameter choices. The binomial observation distribution commonly used for disease prevalence mapping with INLA fails to account for small-scale overdispersion present in the malaria prevalence data, which can lead to poor predictions. Selection of an appropriate alternative such as the Beta-binomial distribution is required to produce a reliable model fit. The small-scale random effect term in FRK overcomes this pitfall, but FRK model estimates are very reliant on providing a sufficient number and appropriate configuration of basis functions. Unfortunately the computation time for FRK increases rapidly with increasing basis resolution.

CONCLUSIONS

INLA and FRK both enable scalable geostatistical modelling of malaria prevalence data. However care must be taken when using both methods to assess the fit of the model to data and plausibility of predictions, in order to select appropriate model assumptions and parameters.

Using model-based geostatistics for assessing the elimination of trachoma

BACKGROUND

Trachoma is the commonest infectious cause of blindness worldwide. Efforts are being made to eliminate trachoma as a public health problem globally. However, as prevalence decreases, it becomes more challenging to precisely predict prevalence. We demonstrate how model-based geostatistics (MBG) can be used as a reliable, efficient, and widely applicable tool to assess the elimination status of trachoma.

METHODS

We analysed trachoma surveillance data from Brazil, Malawi, and Niger. We developed geostatistical Binomial models to predict trachomatous inflammation-follicular (TF) and trachomatous trichiasis (TT) prevalence. We proposed a general framework to incorporate age and gender in the geostatistical models, whilst accounting for residual spatial and non-spatial variation in prevalence through the use of random effects. We also used predictive probabilities generated by the geostatistical models to quantify the likelihood of having achieved the elimination target in each evaluation unit (EU).

RESULTS

TF and TT prevalence varied considerably by country, with Brazil showing the lowest prevalence and Niger the highest. Brazil and Malawi are highly likely to have met the elimination criteria for TF in each EU, but, for some EUs, there was high uncertainty in relation to the elimination of TT according to the model alone. In Niger, the predicted prevalence varied significantly across EUs, with the probability of having achieved the elimination target ranging from values close to 0% to 100%, for both TF and TT.

CONCLUSIONS

We demonstrated the wide applicability of MBG for trachoma programmes, using data from different epidemiological settings. Unlike the standard trachoma prevalence survey approach, MBG provides a more statistically rigorous way of quantifying uncertainty around the achievement of elimination prevalence targets, through the use of spatial correlation. In addition to the analysis of existing survey data, MBG also provides an approach to identify areas in which more sampling effort is needed to improve EU classification. We advocate MBG as the new standard method for analysing trachoma survey outputs.

Interplay among malnutrition, chemoprevention, and the risk of malaria in young Ugandan children: Longitudinal pharmacodynamic and growth analysis

African children are at risk of malaria and malnutrition. We quantified relationships between malaria and malnutrition among young Ugandan children in a high malaria transmission region. Data were used from a randomized controlled trial where Ugandan HIV-unexposed (n = 393) and HIV-exposed (n = 186) children were randomized to receive no malaria chemoprevention, monthly sulfadoxine-pyrimethamine, daily trimethoprim-sulfamethoxazole, or monthly dihydroartemisinin-piperaquine (DP) from age 6-24 months, and then were followed off chemoprevention until age 36 months. Monthly height and weight, and time of incident malaria episodes were obtained; 89 children who received DP contributed piperaquine (PQ) concentrations. Malaria hazard was modeled using parametric survival analysis adjusted for repeated events, and height and weight were modeled using a Brody growth model. Among 579 children, stunting (height-for-age z-score [ZHA] < -2) was associated with a 17% increased malaria hazard (95% confidence interval [CI] 10-23%) compared with children with a ZHA of zero. DP was associated with a 35% lower malaria hazard (hazard ratio [HR] [95% CI], 0.65 [0.41-0.97]), compared to no chemoprevention. After accounting for PQ levels, stunted children who received DP had 2.1 times the hazard of malaria (HR [95% CI] 2.1 [1.6-3.0]) compared with children with a ZHA of zero who received DP. Each additional malaria episode was associated with a 0.4% reduced growth rate for height. Better dosing regimens are needed to optimize malaria prevention in malnourished populations, but, importantly, malaria chemoprevention may reduce the burden of malnutrition in early childhood.

Modelling chronic malnutrition in Zambia: A Bayesian distributional regression approach

Background

The burden of child under-nutrition still remains a global challenge, with greater severity being faced by low- and middle-income countries, despite the strategies in the Sustainable Development Goals (SDGs). Globally, malnutrition is the one of the most important risk factors associated with illness and death, affecting hundreds of millions of pregnant women and young children. Sub-Saharan Africa is one of the regions in the world struggling with the burden of chronic malnutrition. The 2018 Zambia Demographic and Health Survey (ZDHS) report estimated that 35% of the children under five years of age are stunted. The objective of this study was to analyse the distribution, and associated factors of stunting in Zambia.

Methods

We analysed the relationships between socio-economic, and remote sensed characteristics and anthropometric outcomes in under five children, using Bayesian distributional regression. Georeferenced data was available for 25,852 children from two waves of the ZDHS, 31% observation were from the 2007 and 69% were from the 2013/14. We assessed the linear, non-linear and spatial effects of covariates on the height-for-age z-score.

Results

Stunting decreased between 2007 and 2013/14 from a mean z-score of 1.59 (credible interval (CI): -1.63; -1.55) to -1.47 (CI: -1.49; -1.44). We found a strong non-linear relationship for the education of the mother and the wealth of the household on the height-for-age z-score. Moreover, increasing levels of maternal education above the eighth grade were associated with a reduced variation of stunting. Our study finds that remote sensed covariates alone explain little of the variation of the height-for-age z-score, which highlights the importance to collect socio-economic characteristics, and to control for socio-economic characteristics of the individual and the household.

Conclusions

While stunting still remains unacceptably high in Zambia with remarkable regional inequalities, the decline is lagging behind goal two of the SDGs. This emphasises the need for policies that help to reduce the share of chronic malnourished children within Zambia.

Global, regional and national epidemiology and prevalence of child stunting, wasting and underweight in low- and middle-income countries, 2006-2018

In 2016, undernutrition, as manifested in childhood stunting, wasting, and underweight were estimated to cause over 1.0 million deaths, 3.9% of years of life lost, and 3.8% of disability-adjusted life years globally. The objective of this study is to estimate the prevalence of undernutrition in low- and middle-income countries (LMICs) using the 2006-2018 cross-sectional nationally representative demographic and health surveys (DHS) data and to explore the sources of regional variations. Anthropometric measurements of children 0-59 months of age from DHS in 62 LMICs worldwide were used. Complete information was available for height-for-age (n = 624,734), weight-for-height (n = 625,230) and weight-for-age (n = 626,130). Random-effects models were fit to estimate the pooled prevalence of stunting, wasting, and underweight. Sources of heterogeneity in the prevalence estimates were explored through subgroup meta-analyses and meta-regression using generalized linear mixed-effects models. Human development index (a country-specific composite index based on life expectancy, literacy, access to education and per capita gross domestic product) and the United Nations region were explored as potential sources of variation in undernutrition. The overall prevalence was 29.1% (95% CI 26.7%, 31.6%) for stunting, 6.3% (95% CI 4.6%, 8.2%) for wasting, and 13.7% (95% CI 10.9%, 16.9%) for underweight. Subgroup analyses suggested that Western Africa, Southern Asia, and Southeastern Asia had a substantially higher estimated prevalence of undernutrition than global average estimates. In multivariable meta-regression, a combination of human development index and United Nations region (a proxy for geographical variation) explained 54%, 56%, and 66% of the variation in stunting, wasting, and underweight prevalence, respectively. Our findings demonstrate that regional, subregional, and country disparities in undernutrition remain, and the residual gaps to close towards achieving the second sustainable development goal-ending undernutrition by 2030.

Modelling children's anthropometric status using Bayesian distributional regression merging socio-economic and remote sensed data from South Asia and sub-Saharan Africa

A history of insufficient nutritional intake is reflected by low anthropometric measures and can lead to growth failures, limited mental development, poor health outcomes and a higher risk of dying. Children below five years are among those most vulnerable and, while improvements in the share of children affected by insufficient nutritional intake has been observed, both sub-Saharan Africa and South Asia have a disproportionately high share of growth failures and large disparities at national and sub-national levels. In this study, we use a Bayesian distributional regression approach to develop models for the standard anthropometric measures, stunting and wasting. This approach allows us to model both the mean and the standard deviation of the underlying response distribution. Accordingly, the whole distribution of the anthropometric measures can be evaluated. This is of particular importance, considering the fact that (severe) growth failures of children are defined having a z-score below -2 (-3), emphasising the need to extend the analysis beyond the conditional mean. In addition, we merge individual data taken from the Demographic and Health Surveys with remote sensed data for a large sample of 38 countries located in sub-Saharan Africa and South Asia for the period 1990-2016, in order to combine individual and household specific characteristics with geophysical and environmental characteristics, and to allow for a comparison over time. Our results show besides gender differences across space, and strong non-linear effects of included socio-economic characteristics, in particular for maternal education and the wealth of the household that, surprisingly, in the presence of socio-economic characteristics, remote sensed data does not contribute to variations in growth failures, and including a pure spatial effect excluding remote sensed data leads to even better results. Further, while all regions showed improvements towards the target of the Sustainable Development Goals (SDGs), our analysis identifies hotspots of growth failures at sub-national levels within India, Nigeria, Niger, and Madagascar, emphasising the need to accelerate progress to reach the target set by the SDGs.

Geostatistical analysis and mapping of malaria risk in children of Mozambique

Malaria remains one of the most prevalent infectious diseases in the tropics and subtropics, and Mozambique is not an exception. To design geographically targeted and effective intervention mechanisms of malaria, an up-to-date map that shows the spatial distribution of malaria is needed. This study analyzed 2018 Mozambique Malaria Indicator Survey using geostatistical methods to: i) explore individual, household, and community-level determinants of malaria in under-five children, ii) prepare a malaria prevalence map in Mozambique, and iii) produce prediction prevalence maps and exceedence probability across the country. The results show the overall weighted prevalence of malaria was 38.9% (N = 4347, with 95% CI: 36.9%–40.8%). Across different provinces of Mozambique, the prevalence of malaria ranges from 1% in Maputo city to 57.3% in Cabo Delgado province. Malaria prevalence was found to be higher in rural areas, increased with child’s age, and decreased with household wealth index and mother’s level of education. Given the high prevalence of childhood malaria observed in Mozambique there is an urgent need for effective public health interventions in malaria hot spot areas. The household determinants of malaria infection that are identified in this study as well as the maps of parasitaemia risk could be used by malaria control program implementers to define priority intervention areas.

Contextual Data in IPUMS DHS: Physical and Social Environment Variables linked to the Demographic and Health Surveys

The Demographic and Health Surveys (DHS) are the most important source of comparative information on the health of women and young children in low- and middle-income countries and are well-suited for studies of the relationship between environmental factors and health. However, barriers have limited the use of the DHS for these purposes. IPUMS DHS, an online data dissemination tool, overcomes these barriers, simplifying comparative analyses with DHS. IPUMS DHS recently incorporated environmental variables that can easily be attached to individual or household records, facilitating the use of DHS data for the study of population and environment issues. We provide a brief introduction to IPUMS DHS, describe the current and anticipated environmental variables and how to use them, and provide an example of the novel research possibilities facilitated by this latest IPUMS DHS development. IPUMS-DHS is available free online at dhs.ipums.org.

A Scoping Review of Spatial Analysis Approaches Using Health Survey Data in Sub-Saharan Africa

Spatial analysis has become an increasingly used analytic approach to describe and analyze spatial characteristics of disease burden, but the depth and coverage of its usage for health surveys data in Sub-Saharan Africa are not well known. The objective of this scoping review was to conduct an evaluation of studies using spatial statistics approaches for national health survey data in the SSA region. An organized literature search for studies related to spatial statistics and national health surveys was conducted through PMC, PubMed/Medline, Scopus, NLM Catalog, and Science Direct electronic databases. Of the 4,193 unique articles identified, 153 were included in the final review. Spatial smoothing and prediction methods were predominant (n = 108), followed by spatial description aggregation (n = 25), and spatial autocorrelation and clustering (n = 19). Bayesian statistics methods and lattice data modelling were predominant (n = 108). Most studies focused on malaria and fever (n = 47) followed by health services coverage (n = 38). Only fifteen studies employed nonstandard spatial analyses (e.g., spatial model assessment, joint spatial modelling, accounting for survey design). We recommend that for future spatial analysis using health survey data in the SSA region, there must be an improve recognition and awareness of the potential dangers of a naïve application of spatial statistical methods. We also recommend a wide range of applications using big health data and the future of data science for health systems to monitor and evaluate impacts that are not well understood at local levels.

The Economic Burden of Malaria: Revisiting the Evidence

A portion of the economics literature has long debated about the relative importance of historical, institutional, geographical, and health determinants of economic growth. In 2001, Gallup and Sachs quantified the association between malaria and the level and growth of per capita income over the period 1965–1995 in a cross-country regression framework. We took a contemporary look at Gallup and Sachs’ seminal work in the context of significant progress in malaria control achieved globally since 2000. Focusing on the period 2000–2017, we used the latest data available on malaria case incidence and other determinants of economic growth, as well as macro-econometric methods that are now the professional norm. In our preferred specification using a fixed-effects model, a 10% decrease in malaria incidence was associated with an increase in income per capita of nearly 0.3% on average and a 0.11 percentage point faster per capita growth per annum. Greater average income gains were expected among higher burden countries and those with lower income. Growth of industries with the same level of labor intensity was found to be significantly slower in countries with higher malaria incidence. To analyze the causal impact of malaria on economic outcomes, we used malaria treatment failure and pyrethroid-only insecticide resistance as exogeneous instruments in two-stage least squares estimations. Despite several methodological challenges, as expected in these types of analyses, our findings confirm the intrinsic link between malaria and economic growth and underscore the importance of malaria control in the agenda for sustainable development.

New Perspectives for Mapping Global Population Distribution Using World Settlement Footprint Products

In the production of gridded population maps, remotely sensed, human settlement datasets rank among the most important geographical factors to estimate population densities and distributions at regional and global scales. Within this context, the German Aerospace Centre (DLR) has developed a new suite of global layers, which accurately describe the built-up environment and its characteristics at high spatial resolution: (i) the World Settlement Footprint 2015 layer (WSF-2015), a binary settlement mask; and (ii) the experimental World Settlement Footprint Density 2015 layer (WSF-2015-Density), representing the percentage of impervious surface. This research systematically compares the effectiveness of both layers for producing population distribution maps through a dasymetric mapping approach in nine low-, middle-, and highly urbanised countries. Results indicate that the WSF-2015-Density layer can produce population distribution maps with higher qualitative and quantitative accuracies in comparison to the already established binary approach, especially in those countries where a good percentage of building structures have been identified within the rural areas. Moreover, our results suggest that population distribution accuracies could substantially improve through the dynamic preselection of the input layers and the correct parameterisation of the Settlement Size Complexity (SSC) index.

Mapping 123 million neonatal, infant and child deaths between 2000 and 2017

Since 2000, many countries have achieved considerable success in improving child survival, but localized progress remains unclear. To inform efforts towards United Nations Sustainable Development Goal 3.2-to end preventable child deaths by 2030-we need consistently estimated data at the subnational level regarding child mortality rates and trends. Here we quantified, for the period 2000-2017, the subnational variation in mortality rates and number of deaths of neonates, infants and children under 5 years of age within 99 low- and middle-income countries using a geostatistical survival model. We estimated that 32% of children under 5 in these countries lived in districts that had attained rates of 25 or fewer child deaths per 1,000 live births by 2017, and that 58% of child deaths between 2000 and 2017 in these countries could have been averted in the absence of geographical inequality. This study enables the identification of high-mortality clusters, patterns of progress and geographical inequalities to inform appropriate investments and implementations that will help to improve the health of all populations.

Geostatistical analysis and mapping of malaria risk in children under 5 using point-referenced prevalence data in Ghana

BACKGROUND

Malaria remains a major challenge in sub-Saharan Africa and Ghana is not an exception. Effective malaria transmission control requires evidence-based targeting and utilization of resources. Disease risk mapping provides an effective and efficient tool for monitoring transmission and control efforts. The aim of this study is to analyse and map malaria risk in children under 5 years old, with the ultimate goal of identifying areas where control efforts can be targeted.

METHODS

Data collected from the 2016 Ghana demographic and health survey was analyzed. Binomial logistic regression was applied to examine the determinants of malaria risk among children. Model-based geostatistical methods were applied to analyze, predict and map malaria prevalence.

RESULTS

There is a significant association of malaria prevalence with area of residence (rural/urban), age, indoor residual spray use, social economic status and mother's education level. Overall, parasitaemia prevalence among children under 5 years old for the year 2016 is low albeit characterized by "hotspots" in specific areas.

CONCLUSION

The risk maps indicate the spatial heterogeneity of malaria prevalence. The high resolution maps can serve as an effective tool in the identification of locations that require targeted interventions by programme implementers; this is key and relevant for reducing malaria burden in Ghana.

Spatiotemporal characterisation and risk factor analysis of malaria outbreak in Cabo Verde in 2017

Background

Cabo Verde is a country that has been in the pre-elimination stage of malaria since the year 2000. The country is still reporting cases, particularly in the capital of Praia, where more than 50% of the national population live. This study aims to examine the spatial and temporal epidemiological profile of malaria across the country during the 2017 outbreak and to analyse the risk factors, which may have influenced the trend in malaria cases.

Methods

Longitudinal data collected from all malaria cases in Cabo Verde for the year 2017 were used in this study. The epidemiological characteristics of the cases were analysed. Local and spatial clusters of malaria from Praia were detected by applying the Cluster and Outlier Analysis (Anselin Local Moran's I) to determine the spatial clustering pattern. We then used the Pearson correlation coefficient to analyse the relationship between malaria cases and meteorological variables to identify underlying drivers.

Results

In 2017, 446 cases of malaria were reported in Cabo Verde with the peak of cases in October. These cases were primarily Plasmodium falciparum infections. Of these cases, 423 were indigenous infections recorded in Praia, while 23 were imported malaria cases from different African countries. One case of P. vivax infection was imported from Brazil. Spatial autocorrelation analysis revealed a cluster of high-high malaria cases in the centre of the city. Malaria case occurrence has a very weak correlation (r = 0.16) with breeding site location. Most of the cases (69.9%, R ² = 0.699) were explained by the local environmental condition, with temperature being the primary risk factor followed by relative humidity. A moderately positive relationship was noted with the total pluviometry, while wind speed had a strong negative influence on malaria infections.

Conclusions

In Cabo Verde, malaria remains a serious public health issue, especially in Praia. The high number of cases recorded in 2017 demonstrates the fragility of the situation and the challenges to eliminating indigenous malaria cases and preventing imported cases. Mosquito breeding sites have been the main risk factor, while temperature and precipitation were positively associated with malaria infection. In light of this study, there is an urgent need to reinforce control strategies to achieve the elimination goal in the country.

malariaAtlas: an R interface to global malariometric data hosted by the Malaria Atlas Project

Background

The Malaria Atlas Project (MAP) has worked to assemble and maintain a global open-access database of spatial malariometric data for over a decade. This data spans various formats and topics, including: geo-located surveys of malaria parasite rate; global administrative boundary shapefiles; and global and regional rasters representing the distribution of malaria and associated illnesses, blood disorders, and intervention coverage. MAP has recently released malariaAtlas, an R package providing a direct interface to MAP’s routinely-updated malariometric databases and research outputs.

Methods and results

The current paper reviews the functionality available in malariaAtlas and highlights its utility for spatial epidemiological analysis of malaria. malariaAtlas enables users to freely download, visualise and analyse global malariometric data within R. Currently available data types include: malaria parasite rate and vector occurrence point data; subnational administrative boundary shapefiles; and a large suite of rasters covering a diverse range of metrics related to malaria research. malariaAtlas is here used in two mock analyses to illustrate how this data may be incorporated into a standard R workflow for spatial analysis.

Conclusions

malariaAtlas is the first open-access R-interface to malariometric data, providing a new and reproducible means of accessing such data within a freely available and commonly used statistical software environment. In this way, the malariaAtlas package aims to contribute to the environment of data-sharing within the malaria research community.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2500-5) contains supplementary material, which is available to authorized users.