Estimating under-five mortality in space and time in a developing world context

Asthma hospitalisations and heat exposure in England: a case-crossover study during 2002-2019

BACKGROUND

Previous studies have reported an association between warm temperature and asthma hospitalisation. They have reported different sex-related and age-related vulnerabilities; nevertheless, little is known about how this effect has changed over time and how it varies in space. This study aims to evaluate the association between asthma hospitalisation and warm temperature and investigate vulnerabilities by age, sex, time and space.

METHODS

We retrieved individual-level data on summer asthma hospitalisation at high temporal (daily) and spatial (postcodes) resolutions during 2002-2019 in England from the NHS Digital. Daily mean temperature at 1 km×1 km resolution was retrieved from the UK Met Office. We focused on lag 0-3 days. We employed a case-crossover study design and fitted Bayesian hierarchical Poisson models accounting for possible confounders (rainfall, relative humidity, wind speed and national holidays).

RESULTS

After accounting for confounding, we found an increase of 1.11% (95% credible interval: 0.88% to 1.34%) in the asthma hospitalisation risk for every 1°C increase in the ambient summer temperature. The effect was highest for males aged 16-64 (2.10%, 1.59% to 2.61%) and during the early years of our analysis. We also found evidence of a decreasing linear trend of the effect over time. Populations in Yorkshire and the Humber and East and West Midlands were the most vulnerable.

CONCLUSION

This study provides evidence of an association between warm temperature and hospital admission for asthma. The effect has decreased over time with potential explanations including temporal differences in patterns of heat exposure, adaptive mechanisms, asthma management, lifestyle, comorbidities and occupation.

Impact of ignoring sampling design in the prediction of binary health outcomes through logistic regression: evidence from Malawi demographic and health survey under-five mortality data; 2000-2016

The birth and death rates of a population are among the crucial vital statistics for socio-economic policy planning in any country. Since the under-five mortality rate is one of the indicators for monitoring the health of a population, it requires regular and accurate estimation. The national demographic and health survey data, that are readily available to the puplic, have become a means for answering most health-related questions among African populations, using relevant statistical methods. However, many of such applications tend to ignore survey design effect in the estimations, despite the availability of statistical tools that support the analyses. Little is known about the amount of inaccurate information that is generated when predicting under-five mortality rates. This study estimates and compares the bias encountered when applying unweighted and weighted logistic regression methods to predict under-five mortality rate in Malawi using nationwide survey data. The Malawi demographic and health survey data of 2004, 2010, and 2015-16 were used to determine the bias. The analyses were carried out in R software version 3.6.3 and Stata version 12.0. A logistic regression model that included various bio- and socio-demographic factors concerning the child, mother and households was used to estimate the under-five mortality rate. The results showed that accuracy of predicting the national under-five mortality rate hinges on cluster-weighting of the overall predicted probability of child-deaths, regardless of whether the model was weighted or not. Weighting the model caused small positive and negative changes in various fixed-effect estimates, which diffused the result of weighting in the fitted probabilities of deaths. In turn, there was no difference between the overall predicted mortality rate obtained using the weighted model and that obtained in the unweighted model. We recommend considering survey cluster-weights during the computation of overall predicted probability of events for a binary health outcome. This can be done without worrying about the weights during model fitting, whose aim is prediction of the population parameter.

Progress towards achieving child survival goals in Kenya after devolution: Geospatial analysis with scenario-based projections, 2015-2025

Subnational projections of under-5 mortality (U5M) have increasingly become an essential planning tool to support Sustainable Development Goals (SDGs) agenda and strategies for improving child survival. To support child health policy, planning, and tracking child development goals in Kenya, we projected U5M at units of health decision making. County-specific annual U5M were estimated using a multivariable Bayesian space-time hierarchical model based on intervention coverage from four alternate intervention scale-up scenarios assuming 1) the highest subnational intervention coverage in 2014, 2) projected coverage based on the fastest county-specific rate of change observed in the period between 2003-2014 for each intervention, 3) the projected national coverage based on 2003-2014 trends and 4) the country-specific targets of intervention coverage relative to business as usual (BAU) scenario. We compared the percentage change in U5M based on the four scale-up scenarios relative to BAU and examined the likelihood of reaching SDG 3.2 target of at least 25 deaths/1,000 livebirths by 2022 and 2025. Projections based on 10 factors assuming BAU, showed marginal reductions in U5M across counties with all the counties except Mandera county not achieving the SDG 3.2 target by 2025. Further, substantial reductions in U5M would be achieved based on the various intervention scale-up scenarios, with 63.8% (30), 74.5% (35), 46.8% (22) and 61.7% (29) counties achieving SDG target for scenarios 1,2,3 and 4 respectively by 2025. Scenario 2 yielded the highest reductions of U5M with individual scale-up of access to improved water, recommended treatment of fever and accelerated HIV prevalence reduction showing considerable impact on U5M reduction (≥ 20%) relative to BAU. Our results indicate that sustaining an ambitious intervention scale-up strategy matching the fastest rate observed between 2003-2014 would substantially reduce U5M in Kenya. However, despite this ambitious scale-up scenario, 25% (12 of 47) of the Kenya's counties would still not achieve SDG 3.2 target by 2025.

A FLEXIBLE BAYESIAN FRAMEWORK TO ESTIMATE AGE- AND CAUSE-SPECIFIC CHILD MORTALITY OVER TIME FROM SAMPLE REGISTRATION DATA

In order to implement disease-specific interventions in young age groups, policy makers in low- and middle-income countries require timely and accurate estimates of age- and cause-specific child mortality. High-quality data is not available in settings where these interventions are most needed, but there is a push to create sample registration systems that collect detailed mortality information. current methods that estimate mortality from this data employ multistage frameworks without rigorous statistical justification that separately estimate all-cause and cause-specific mortality and are not sufficiently adaptable to capture important features of the data. We propose a flexible Bayesian modeling framework to estimate age- and cause-specific child mortality from sample registration data. We provide a theoretical justification for the framework, explore its properties via simulation, and use it to estimate mortality trends using data from the Maternal and Child Health Surveillance System in China.

Variation in excess all-cause mortality by age, sex, and province during the first wave of the COVID-19 pandemic in Italy

Although previous evidence suggests that the infection fatality rate from COVID-19 varies by age and sex, and that transmission intensity varies geographically within countries, no study has yet explored the age-sex-space distribution of excess mortality associated with the COVID pandemic. By applying the principles of small-area estimation to existing model formulations for excess mortality, this study develops a novel method for assessing excess mortality across small populations and assesses the pattern of COVID excess mortality by province, year, week, age group, and sex in Italy from March through May 2020. We estimate that 53,200 excess deaths occurred across Italy during this time period, compared to just 35,500 deaths where COVID-19 was registered as the underlying cause of death. Out of the total excess mortality burden, 97% of excess deaths occurred among adults over age 60, and 68% of excess deaths were concentrated among adults over age 80. The burden of excess mortality was unevenly distributed across the country, with just three of Italy’s 107 provinces accounting for 32% of all excess mortality. This method for estimating excess mortality can be adapted to other countries where COVID-19 diagnostic capacity is still insufficient, and could be incorporated into public health rapid response systems.

What have we learned about mortality patterns over the past 25 years?

In this paper, I examine progress in the field of mortality over the past 25 years. I argue that we have been most successful in taking advantage of an increasingly data-rich environment to improve aggregate mortality models and test pre-existing theories. Less progress has been made in relating our estimates of mortality risk at the individual level to broader mortality patterns at the population level while appropriately accounting for contextual differences and compositional change. Overall, I find that the field of mortality continues to be highly visible in demographic journals, including Population Studies. However much of what is published today in field journals could just as easily appear in neighbouring disciplinary journals, as disciplinary boundaries are shrinking.

Decomposing County-Level Working-Age Mortality Trends in the United States Between 1999–2001 and 2015–2017

Studies have documented significant geographic divergence in U.S. mortality in recent decades. However, few studies have examined the extent to which county-level trends in mortality can be explained by national, state, and metropolitan-level trends, and which county-specific factors contribute to remaining variation. Combining vital statistics data on deaths and Census data with time-varying county-level contextual characteristics, we use a spatially explicit Bayesian hierarchical model to analyze the associations between working-age mortality, state, metropolitan status and county-level socioeconomic conditions, family characteristics, labor market conditions, health behaviors, and population characteristics between 2000 and 2017. Additionally, we employ a Shapley decomposition to illustrate the additive contributions of each changing county-level characteristic to the observed mortality change in U.S. counties between 1999–2001 and 2015–2017 over and above national, state, and metropolitan–nonmetropolitan mortality trends. Mortality trends varied by state and metropolitan status as did the contribution of county-level characteristics. Metropolitan status predicted more of the county-level variance in mortality than state of residence. Of the county-level characteristics, changes in percent college-graduates, smoking prevalence and the percent of foreign-born population contributed to a decline in all-cause mortality over this period, whereas increasing levels of poverty, unemployment, and single-parent families and declines manufacturing employment slowed down these improvements, and in many nonmetropolitan areas were large enough to overpower the positive contributions of the protective factors.

Methods for Small Area Population Forecasts: State-of-the-Art and Research Needs

Small area population forecasts are widely used by government and business for a variety of planning, research and policy purposes, and often influence major investment decisions. Yet, the toolbox of small area population forecasting methods and techniques is modest relative to that for national and large subnational regional forecasting. In this paper, we assess the current state of small area population forecasting, and suggest areas for further research. The paper provides a review of the literature on small area population forecasting methods published over the period 2001–2020. The key themes covered by the review are extrapolative and comparative methods, simplified cohort-component methods, model averaging and combining, incorporating socioeconomic variables and spatial relationships, ‘downscaling’ and disaggregation approaches, linking population with housing, estimating and projecting small area component input data, microsimulation, machine learning, and forecast uncertainty. Several avenues for further research are then suggested, including more work on model averaging and combining, developing new forecasting methods for situations which current models cannot handle, quantifying uncertainty, exploring methodologies such as machine learning and spatial statistics, creating user-friendly tools for practitioners, and understanding more about how forecasts are used.

Space-time modeling of child mortality at the Admin-2 level in a low and middle income countries context

The Sustainable Development Goals call for a total reduction of preventable child mortality before 2030. Further, the goals state the desirability to have subnational mortality estimates. Estimates at this level are required for health interventions at the subnational level. In a low and middle income countries context, the data on mortality typically consist of household surveys, which are carried out with a stratified, cluster design, and census microsamples. Most household surveys collect full birth history (FBH) data on birth and death dates of a mother's children, but censuses collect summary birth history (SBH) data which consist only of the number of children born and the number that died. In previous work, direct (survey-weighted) estimates with associated variances were derived from FBH data and smoothed in space and time. Unfortunately, the FBH data from household surveys are usually not sufficiently abundant to obtain yearly estimates at the Admin-2 level (at which interventions are often made). In this paper we describe four extensions to previous work: (i) combining SBH data with FBH data, (ii) modeling on a yearly scale, to combine data on a yearly scale with data at coarser time scales, (iii) adjusting direct estimates in Admin-2 areas where we do not observe any deaths due to small sample sizes, (iv) acknowledge differences in data sources by modeling potential bias arising from the various data sources. The methods are illustrated using household survey and census data from Kenya and Malawi, to produce mortality estimates from 1980 to the time of the most recent survey, and predictions to 2020.

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.

Pointless spatial modeling

The analysis of area-level aggregated summary data is common in many disciplines including epidemiology and the social sciences. Typically, Markov random field spatial models have been employed to acknowledge spatial dependence and allow data-driven smoothing. In the context of an irregular set of areas, these models always have an ad hoc element with respect to the definition of a neighborhood scheme. In this article, we exploit recent theoretical and computational advances to carry out modeling at the continuous spatial level, which induces a spatial model for the discrete areas. This approach also allows reconstruction of the continuous underlying surface, but the interpretation of such surfaces is delicate since it depends on the quality, extent and configuration of the observed data. We focus on models based on stochastic partial differential equations. We also consider the interesting case in which the aggregate data are supplemented with point data. We carry out Bayesian inference and, in the language of generalized linear mixed models, if the link is linear, an efficient implementation of the model is available via integrated nested Laplace approximations. For nonlinear links, we present two approaches: a fully Bayesian implementation using a Hamiltonian Monte Carlo algorithm and an empirical Bayes implementation, that is much faster and is based on Laplace approximations. We examine the properties of the approach using simulation, and then apply the model to the classic Scottish lip cancer data.

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.

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.

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.

Sub national variation and inequalities in under-five mortality in Kenya since 1965

Background

Despite significant declines in under five mortality (U5M) over the last 3 decades, Kenya did not achieve Millennium Development Goal 4 (MDG 4) by 2015. To better understand trends and inequalities in child mortality, analysis of U5M variation at subnational decision making units is required. Here the comprehensive compilation and analysis of birth history data was used to understand spatio-temporal variation, inequalities and progress towards achieving the reductions targets of U5M between 1965 and 2013 and projected to 2015 at decentralized health planning units (counties) in Kenya.

Methods

Ten household surveys and three censuses with data on birth histories undertaken between 1989 and 2014 were assembled. The birth histories were allocated to the respective counties and demographic methods applied to estimate U5M per county by survey. To generate a single U5M estimate for year and county, a Bayesian spatio-temporal Gaussian process regression was fitted accounting for variation in sample size, surveys and demographic methods. Inequalities and the progress in meeting the goals set to reduce U5M were evaluated subnationally.

Results

Nationally, U5M reduced by 61·6%, from 141·7 (121·6–164·0) in 1965 to 54·5 (44·6–65·5) in 2013. The declining U5M was uneven ranging between 19 and 80% across the counties with some years when rates increased. By 2000, 25 counties had achieved the World Summit for Children goals. However, as of 2015, no county had achieved MDG 4. There was a striking decline in the levels of inequality between counties over time, however, disparities persist. By 2013 there persists a 3·8 times difference between predicted U5M rates when comparing counties with the highest U5M rates against those with the lowest U5M rates.

Conclusion

Kenya has made huge progress in child survival since independence. However, U5M remains high and heterogeneous with substantial differences between counties. Better use of the current resources through focused allocation is required to achieve further reductions, reduce inequalities and increase the likelihood of achieving Sustainable Development Goal 3·2 on U5M by 2030.

Electronic supplementary material

The online version of this article (10.1186/s12889-019-6474-1) contains supplementary material, which is available to authorized users.

Changes in the spatial distribution of the under-five mortality rate: Small-area analysis of 122 DHS surveys in 262 subregions of 35 countries in Africa

The under-five mortality rate (U5MR) is a critical and widely available population health indicator. Both the MDGs and SDGs define targets for improvement in the U5MR, and the SDGs require spatial disaggregation of indicators. We estimate trends in the U5MR for Admin-1 subnational areas using 122 DHS surveys in 35 countries in Africa and assess progress toward the MDG target reductions for each subnational region and each country as a whole. In each country, direct weighted estimates of the U5MR from each survey are calculated and combined into a single estimate for each Admin-1 region across five-year periods. Our method fully accounts for the sample design of each survey. The region-time-specific estimates are smoothed using a Bayesian, space-time model that produces more precise estimates (when compared to the direct estimates) at a one-year scale that are consistent with each other in both space and time. The resulting estimated distributions of the U5MR are summarized and used to assess subnational progress toward the MDG 4 target of two-thirds reduction in the U5MR during 1990–2015. Our space-time modeling approach is tractable and can be readily applied to a large collection of sample survey data. Subnational, regional spatial heterogeneity in the levels and trends in the U5MR vary considerably across Africa. There is no generalizable pattern between spatial heterogeneity and levels or trends in the U5MR. Subnational, small-area estimates of the U5MR: (i) identify subnational regions where interventions are still necessary and those where improvement is well under way; and (ii) countries where there is very little spatial variation and others where there are important differences between subregions in both levels and trends. More work is necessary to improve both the data sources and methods necessary to adequately measure subnational progress toward the SDG child survival targets.

Development and validation of a new method for indirect estimation of neonatal, infant, and child mortality trends using summary birth histories

BACKGROUND

The addition of neonatal (NN) mortality targets in the Sustainable Development Goals highlights the increased need for age-specific quantification of mortality trends, detail that is not provided by summary birth histories (SBHs). Several methods exist to indirectly estimate trends in under-5 mortality from SBHs; however, efforts to monitor mortality trends in important age groups such as the first month and first year of life have yet to utilize the vast amount of SBH data available from household surveys and censuses.

METHODS AND FINDINGS

We analyzed 243 Demographic and Health Surveys (DHS) from 76 countries, which collected both complete and SBHs from 8.5 million children from 2.3 million mothers to develop a new empirically based method to indirectly estimate time trends in age-specific mortality. We used complete birth history (CBH) data to train a discrete hazards generalized additive model in order to predict individual hazard functions for children based on individual-, mother-, and country-year-level covariates. Individual-level predictions were aggregated over time by assigning probability weights to potential birth years from mothers from SBH data. Age-specific estimates were evaluated in three ways: using cross-validation, using an external database of an additional 243 non-DHS census and survey data sources, and comparing overall under-5 mortality to existing indirect methods. Our model was able to closely approximate trends in age-specific child mortality. Depending on age, the model was able to explain between 80% and 95% of the variation in the validation data. Bias was close to zero in every age, with median relative errors spanning from 0.96 to 1.09. For trends in all under-5s, performance was comparable to the methods used for the Global Burden of Disease (GBD) study and significantly better than the standard indirect (Brass) method, especially in the 5 years preceding a survey. For the 15 years preceding surveys, the new method and GBD methods could explain more than 95% of the variation in the validation data for under-5s, whereas the standard indirect variants tested could only explain up to 88%. External validation using census and survey data found close agreement with concurrent direct estimates of mortality in the NN and infant age groups. As a predictive method based on empirical data, one limitation is that potential issues in these training data could be reflected in the resulting application of the method out of sample.

CONCLUSIONS

This new method for estimating child mortality produces results that are comparable to current best methods for indirect estimation of under-5 mortality while additionally producing age-specific estimates. Use of such methods allows researchers to utilize a massive amount of SBH data for estimation of trends in NN and infant mortality. Systematic application of these methods could further improve the evidence base for monitoring of trends and inequalities in age-specific child mortality.