Publications



A Bayesian cohort component projection model to estimate adult populations at the subnational level in data-sparse settings.

M. Alexander, L. Alkema (forthcoming). Demography




Abstract

Accurate estimates of subnational populations are important for policy formulation and monitoring population health indicators. For example, estimates of the number of women of reproductive age are important to understand the population at risk to maternal mortality and unmet need for contraception. However, in many low-income countries, data on population counts and components of population change are limited, and so levels and trends subnationally are unclear. We present a Bayesian constrained cohort component model for the estimation and projection of subnational populations. The model builds on a cohort component projection framework, incorporates census data and estimates from the United Nation's World Population Prospects, and uses characteristic mortality schedules to obtain estimates of population counts and the components of population change, including internal migration. The data required as inputs to the model are minimal and available across a wide range of countries, including most low-income countries. The model is applied to estimate and project populations by county in Kenya for 1979-2019, and validated against the 2019 Kenyan census.

Estimating the stillbirth rate for 195 countries using a Bayesian sparse regression model with temporal smoothing.

Z. Wang, M. Fix, L. Hug, A. Mishra, D. You, H. Blencowe, J. Wakefield, L. Alkema (forthcoming). Annals of Applied Statistics




Abstract

Estimation of stillbirth rates globally is complicated because of the paucity of reliable data from countries where most stillbirths occur. We compiled data and developed a Bayesian hierarchical temporal sparse regression model for estimating stillbirth rates for all countries from 2000 to 2019. The model combines covariates with a temporal smoothing process so that estimates are data-driven in country-periods with high-quality data and deter-mined by covariates for country-periods with limited or no data. Horseshoepriors are used to encourage sparseness. The model adjusts observations with alternative stillbirth definitions and accounts for bias in observations that are subject to non-sampling errors. In-sample goodness of fit and out-of-sample validation results suggest that the model is reasonably well calibrated. The model is used by the UN Inter-agency Group for Child Mortality Estimation to monitor the stillbirth rate for all countries.

Country-specific estimates of unintended pregnancy and abortion incidence: a global comparative analysis of levels in 2015–2019.

J. Bearak, A. Popinchalk, B. Ganatra, A. Moller, Ö. Tunçalp, C. Beavin, L. Kwok, L. Alkema (2022). BMJ Global Health 7:e007151.




Abstract

Introduction Internationally comparable estimates of unintended pregnancy and abortion incidence can illuminate disparities in sexual and reproductive health and autonomy. Country-specific estimates are essential to enable international comparison, and to inform country-level policy and programming.

Methods We developed a Bayesian model which jointly estimated unintended pregnancy and abortion rates using information on contraceptive needs and use, contraceptive method mix, birth rates, the proportions of births from unintended pregnancies and abortion incidence data. Main outcomes were the estimated rates of unintended pregnancy and abortion for 150 countries and territories, reported for the 5-year period 2015–2019, as annual averages per 1000 women aged 15–49 years.

Results Estimated unintended pregnancy rates ranged from 11 (80% uncertainty interval: 9 to 13) in Montenegro to 145 (131 to 159) in Uganda per 1000 women aged 15–49 years. Between-country heterogeneity was substantial in all Sustainable Development Goal (SDG) regions, but was greatest in sub-Saharan Africa. Estimated abortion rates ranged from 5 (5 to 6) in Singapore to 80 (55 to 113) in Georgia. Variation between country estimates was similar in all SDG regions except for Europe and Northern America, where estimated abortion rates were generally lower.

Conclusion The estimates reflect variation in the degree to unintended pregnancy and abortion that are experienced in countries throughout the world. This evidence highlights the importance of investing in access to contraception and comprehensive abortion care, including in regions which may have lower rates of unintended pregnancy or abortion, respectively, as countries may differ substantially from regional averages.

Temporal models for demographic and global health outcomes in multiple populations: Introducing a new framework to review and standardize documentation of model assumptions and facilitate model comparison.

H. Susmann, M. Alexander, L. Alkema (2022). International Statistical Review.




Abstract

There is growing interest in producing estimates of demographic and global health indicators in populations with limited data. Statistical models are needed to combine data from multiple data sources into estimates and projections with uncertainty. Diverse modelling approaches have been applied to this problem, making comparisons between models difficult. We propose a model class, Temporal Models for Multiple Populations (TMMPs), to facilitate both documentation of model assumptions in a standardised way and comparison across models. The class makes a distinction between the process model, which describes latent trends in the indicator interest, and the data model, which describes the data generating process of the observed data. We provide a general notation for the process model that encompasses many popular temporal modelling techniques, and we show how existing models for a variety of indicators can be written using this notation. We end with a discussion of outstanding questions and future directions.

The Global Burden of Disease fertility forecasts: Summary of the approach used and associated statistical concerns.

L. Alkema (2020). Preprint.




Abstract

BACKGROUND The Global Burden of Disease (GBD) project’s forecasts up to 2100 suggest fertility drops will be even greater in sub-Saharan Africa than the UN Population Division (UNPD) has predicted.
OBJECTIVE This reflection summarizes the main assumptions used in the GBD fertility forecasts. I assess the methods used, focusing on high fertility countries and the use of met need for contraceptives as a predictor.
Based on GBD’s forecasting method, I draw two conclusions. Firstly, GBD fertility forecasts are based on unvalidated assumptions about increasing met need for contraception and may overestimate decreases in fertility in countries with low levels of modern contraceptive use. Secondly, the GBD forecast model for fertility is not a causal model for predicting changes. Claims GBD researchers make about the effect of changing access to family planning on fertility are not informative for guiding policy. Based on the GBD validation exercise, I conclude that the GBD study did not check the performance of the method for predicting left-out fertility data. Also the approach used to compare the predictive performance of UNPD and GBD forecasting methods may give the GBD method an inherent advantage.
CONCLUSIONS Communication regarding the GBD method and its findings must avoid causal language and acknowledge the method’s limitations. Future research should examine the performance of the method, especially for countries with low modern contraceptive use.
CONTRIBUTION This paper summarizes the GBD fertility forecasting method and indicates three areas of concern about it and its use.

Estimating misclassification errors in the reporting of maternal mortality in national civil registration vital statistics systems: A Bayesian hierarchical bivariate random walk model to estimate sensitivity and specificity for multiple countries and years with missing data.

E. Peterson, D. Chou, AB. Moller, A. Gemmill, L. Say, L. Alkema (2022). Statistics in Medicine 41:2483–2496.



Abstract

Civil registration vital statistics (CRVS) systems provide data on maternal mortality that can be used for monitoring trends and to inform policies and programs. However, CRVS maternal mortality data may be subject to substantial reporting errors due to misclassification of maternal deaths. Information on misclassification is available for selected countries and periods only. We developed a Bayesian hierarchical bivariate random walk model to estimate sensitivity and specificity for multiple populations and years and used the model to estimate misclassification errors in the reporting of maternal mortality in CRVS systems. The proposed Bayesian misclassification (BMis) model captures differences in sensitivity and specificity across populations and over time, allows for extrapolations to periods with missing data, and includes an exact likelihood function for data provided in aggregated form. Validation exercises using maternal mortality data suggest that BMis is reasonably well calibrated and improves upon the CRVS-adjustment approach used until 2018 by the UN Maternal Mortality Inter-Agency Group (UN-MMEIG) to account for bias in CRVS data resulting from misclassification error. Since 2019, BMis is used by the UN-MMEIG to account for misclassification errors when estimating maternal mortality using CRVS data.

Global, regional and national trends in under-5 mortality 1990-2019 with scenario-based projections to 2030: a systematic analysis by the United Nations Inter-agency Group for Child Mortality Estimation.

D. Sharrow, L. Hug, D. You, L. Alkema, R. Black, S. Cousens, T. Croft, V. Gaigbe-Togbe, P. Gerland, M. Guillot, K. Hill, B. Masquelier, C. Mathers, J. Pedersen, K. Strong, E. Suzuki, J. Wakefield, N. Walker (2022). The Lancet Global Health 10(2): e195–206.




Abstract

Background

The Sustainable Development Goals (SDGs), set in 2015 by the UN General Assembly, call for all countries to reach an under-5 mortality rate (U5MR) of at least as low as 25 deaths per 1000 livebirths and a neonatal mortality rate (NMR) of at least as low as 12 deaths per 1000 livebirths by 2030. We estimated levels and trends in under-5 mortality for 195 countries from 1990 to 2019, and conducted scenario-based projections of the U5MR and NMR from 2020 to 2030 to assess country progress in, and potential for, reaching SDG targets on child survival and the potential under-5 and neonatal deaths over the next decade.

Methods

Levels and trends in under-5 mortality are based on the UN Inter-agency Group for Child Mortality Estimation (UN IGME) database on under-5 mortality, which contains around 18 000 country-year datapoints for 195 countries—nearly 10 000 of those datapoints since 1990. The database includes nationally representative mortality data from vital registration systems, sample registration systems, population censuses, and household surveys. As with previous sets of national UN IGME estimates, a Bayesian B-spline bias-reduction model (B3) that considers the systematic biases associated with the different data source types was fitted to these data to generate estimates of under-5 (age 0–4 years) mortality with uncertainty intervals for 1990–2019 for all countries. Levels and trends in the neonatal mortality rate (0–27 days) are modelled separately as the log ratio of the neonatal mortality rate to the under-5 mortality rate using a Bayesian model. Estimated mortality rates are combined with livebirths data to calculate the number of under-5 and neonatal deaths. To assess the regional and global burden of under-5 deaths in the present decade and progress towards SDG targets, we constructed several scenario-based projections of under-5 mortality from 2020 to 2030 and estimated national, regional, and global under-5 mortality trends up to 2030 for each scenario.

Findings

The global U5MR decreased by 59% (90% uncertainty interval [UI] 56–61) from 93·0 (91·7–94·5) deaths per 1000 livebirths in 1990 to 37·7 (36·1–40·8) in 2019, while the annual number of global under-5 deaths declined from 12·5 (12·3–12·7) million in 1990 to 5·2 (5·0–5·6) million in 2019—a 58% (55–60) reduction. The global NMR decreased by 52% (90% UI 48–55) from 36·6 (35·6–37·8) deaths per 1000 livebirths in 1990, to 17·5 (16·6–19·0) in 2019, and the annual number of global neonatal deaths declined from 5·0 (4·9–5·2) million in 1990, to 2·4 (2·3–2·7) million in 2019, a 51% (47–54) reduction. As of 2019, 122 of 195 countries have achieved the SDG U5MR target, and 20 countries are on track to achieve the target by 2030, while 53 will need to accelerate progress to meet the target by 2030. 116 countries have reached the SDG NMR target with 16 on track, leaving 63 at risk of missing the target. If current trends continue, 48·1 million under-5 deaths are projected to occur between 2020 and 2030, almost half of them projected to occur during the neonatal period. If all countries met the SDG target on under-5 mortality, 11 million under-5 deaths could be averted between 2020 and 2030.

Interpretation

As a result of effective global health initiatives, millions of child deaths have been prevented since 1990. However, the task of ending all preventable child deaths is not done and millions more deaths could be averted by meeting international targets. Geographical and economic variation demonstrate the possibility of even lower mortality rates for children under age 5 years and point to the regions and countries with highest mortality rates and in greatest need of resources and action.

Global estimation and scenario-based projections of sex ratio at birth and missing female births using a Bayesian hierarchical time series mixture model.

F. Chao, P. Gerland, A. Cook, L. Alkema (2021). Annals of Applied Statistics 15(3): 1499-1528. DOI: 10.1214/20-AOAS1436




Abstract

The sex ratio at birth (SRB) is defined as the ratio of male to female live births. The SRB imbalance in parts of the world over the past several decades is a direct consequence of sex-selective abortion, driven by the coexistence of son preference, readily available technology of prenatal sex determination and fertility decline. Estimation and projection of the degree of SRB imbalance is complicated because of variability in SRB reference levels and because of the uncertainty associated with SRB observations.
We develop Bayesian hierarchical time series mixture models for SRB estimation and scenario-based projections for all countries from 1950 to 2100. We model the SRB regional and national reference levels and the fluctuation around national reference levels. We identify countries at risk of SRB imbalances and model both: (i) the absence or presence of sex ratio transitions in such countries and, if present, (ii) the transition process. The transition model of SRB imbalance captures three stages (increase, stagnation and convergence back to SRB baselines). The model identifies countries with statistical evidence of SRB inflation in a fully Bayesian approach. The scenario-based SRB projections are based on the sex ratio transition model with varying assumptions regarding the occurrence of a sex ratio transition in at-risk countries. Projections are used to quantify the future burden of missing female births due to sex-selective abortions under different scenarios.

Global, regional, and national levels and trends in stillbirths from 2000 to 2019: a systematic assessment.

L. Hug, D. You, H. Blencowe, A. Mishra, Z. Wang, M. Fix, J. Wakefield, A. Moran, V. Gaigbe-Togbe, E. Suzuki, D. Blau, S. Cousens, A. Creanga, T. Croft, K. Hill, K Joseph, S. Maswime, E. McClure, R. Pattinson, J. Pedersen, L. Smith, J. Zeitlin, L. Alkema (2021). The Lancet 398 (10302): 772–85




Abstract

Background: Stillbirths are a major public health issue and a sensitive marker of the quality of care around pregnancy and birth. The UN Global Strategy for Women's, Children's and Adolescents’ Health (2016–30) and the Every Newborn Action Plan (led by UNICEF and WHO) call for an end to preventable stillbirths. A first step to prevent stillbirths is obtaining standardised measurement of stillbirth rates across countries. We estimated stillbirth rates and their trends for 195 countries from 2000 to 2019 and assessed progress over time.
Methods: For a systematic assessment, we created a dataset of 2833 country-year datapoints from 171 countries relevant to stillbirth rates, including data from registration and health information systems, household-based surveys, and population-based studies. After data quality assessment and exclusions, we used 1531 datapoints to estimate country-specific stillbirth rates for 195 countries from 2000 to 2019 using a Bayesian hierarchical temporal sparse regression model, according to a definition of stillbirth of at least 28 weeks’ gestational age. Our model combined covariates with a temporal smoothing process such that estimates were informed by data for country-periods with high quality data, while being based on covariates for country-periods with little or no data on stillbirth rates. Bias and additional uncertainty associated with observations based on alternative stillbirth definitions and source types, and observations that were subject to non-sampling errors, were included in the model. We compared the estimated stillbirth rates and trends to previously reported mortality estimates in children younger than 5 years.
Findings: Globally in 2019, an estimated 2·0 million babies (90% uncertainty interval [UI] 1·9–2·2) were stillborn at 28 weeks or more of gestation, with a global stillbirth rate of 13·9 stillbirths (90% UI 13·5–15·4) per 1000 total births. Stillbirth rates in 2019 varied widely across regions, from 22·8 stillbirths (19·8–27·7) per 1000 total births in west and central Africa to 2·9 (2·7–3·0) in western Europe. After west and central Africa, eastern and southern Africa and south Asia had the second and third highest stillbirth rates in 2019. The global annual rate of reduction in stillbirth rate was estimated at 2·3% (90% UI 1·7–2·7) from 2000 to 2019, which was lower than the 2·9% (2·5–3·2) annual rate of reduction in neonatal mortality rate (for neonates aged less than 28 days) and the 4·3% (3·8–4·7) annual rate of reduction in mortality rate among children aged 1–59 months during the same period. Based on the lower bound of the 90% UIs, 114 countries had an estimated decrease in stillbirth rate since 2000, with four countries having a decrease of at least 50·0%, 28 having a decrease of 25·0–49·9%, 50 having a decrease of 10·0–24·9%, and 32 having a decrease of less than 10·0%. For the remaining 81 countries, we found no decrease in stillbirth rate since 2000. Of these countries, 34 were in sub-Saharan Africa, 16 were in east Asia and the Pacific, and 15 were in Latin America and the Caribbean.
Interpretation: Progress in reducing the rate of stillbirths has been slow compared with decreases in the mortality rate of children younger than 5 years. Accelerated improvements are most needed in the regions and countries with high stillbirth rates, particularly in sub-Saharan Africa. Future prevention of stillbirths needs increased efforts to raise public awareness, improve data collection, assess progress, and understand public health priorities locally, all of which require investment.

Using Family Planning Service Statistics to inform model-based estimates of modern contraceptive prevalence.

N. Cahill, E. Sonneveldt, P. Emmart, J. Williamson, R. Mbu, A. Barrière Fodjo Yetgang, I. Dambula, G. Azambuja, A. Mahumane Govo, B. Joshi, S. Felix, C. Makashaka, V. Ndaruhutse, J. Serucaca, B. Madzima, B. Muzavazi, L. Alkema (2021). Plos One 16(10), e0258304.




Abstract

The annual assessment of Family Planning (FP) indicators, such as the modern contraceptive prevalence rate (mCPR), is a key component of monitoring and evaluating goals of global FP programs and initiatives. To that end, the Family Planning Estimation Model (FPEM) was developed with the aim of producing survey-informed estimates and projections of mCPR and other key FP indictors over time. With large-scale surveys being carried out on average every 3–5 years, data gaps since the most recent survey often exceed one year. As a result, survey-based estimates for the current year from FPEM are often based on projections that carry a larger uncertainty than data informed estimates. In order to bridge recent data gaps we consider the use of a measure, termed Estimated Modern Use (EMU), which has been derived from routinely collected family planning service statistics. However, EMU data come with known limitations, namely measurement errors which result in biases and additional variation with respect to survey-based estimates of mCPR. Here we present a data model for the incorporation of EMU data into FPEM, which accounts for these limitations. Based on known biases, we assume that only changes in EMU can inform FPEM estimates, while also taking inherent variation into account. The addition of this EMU data model to FPEM allows us to provide a secondary data source for informing and reducing uncertainty in current estimates of mCPR. We present model validations using a survey-only model as a baseline comparison and we illustrate the impact of including the EMU data model in FPEM. Results show that the inclusion of EMU data can change point-estimates of mCPR by up to 6.7 percentage points compared to using surveys only. Observed reductions in uncertainty were modest, with the width of uncertainty intervals being reduced by up to 2.7 percentage points.


Projecting sex imbalances at birth at global, regional and national levels from 2021 to 2100: scenario-based Bayesian probabilistic projections of the sex ratio at birth and missing female births based on 3.26 billion birth records.

F. Chao, A. Cook, P. Gerland, C. Guilmoto, L. Alkema (2021). BMJ Global Health (6): e005516.




Abstract

Introduction: Skewed levels of the sex ratio at birth (SRB) due to sex-selective abortions have been observed in several countries since the 1970s. They will lead to long-term sex imbalances in more than one-third of the world’s population with yet unknown social and economic impacts on affected countries. Understanding the potential evolution of sex imbalances at birth is therefore essential for anticipating and planning for changing sex structures across the world.
Methods: We produced probabilistic SRB projections from 2021 to 2100 based on different scenarios of sex ratio transition and assessed their implications in terms of missing female births at global, regional and national levels. Based on a comprehensive SRB database with 3.26 billion birth records, we project the skewed SRB and missing female births with a Bayesian hierarchical time series mixture model. The SRB projections under reference scenario S1 assumed SRB transitions only for countries with strong statistical evidence of SRB inflation, and the more extreme scenario S2 assumed a sex ratio transition for countries at risk of SRB inflation but with no or limited evidence of ongoing inflation.
Results: Under scenario S1, we projected 5.7 (95% uncertainty interval (1.2; 15.3)) million additional missing female births to occur by 2100. Countries affected will be those already affected in the past by imbalanced SRB, such as China and India. If all countries at risk of SRB inflation experience a sex ratio transition as in scenario S2, the projected missing female births increase to 22.1 (12.2; 39.8) million with a sizeable contribution of sub-Saharan Africa.
Conclusion: The scenario-based projections provide important illustrations of the potential burden of future prenatal sex discrimination and the need to monitor SRBs in countries with son preference. Policy planning will be needed in the years to come to minimise future prenatal sex discrimination and its impact on social structures.

Fpemlocal: Estimating family planning indicators in R for a single population of interest. Gates Open Research.

G. Guranich, N. Cahill, L. Alkema (2021). Gates Open Research 5(24).




Abstract

The global Family Planning Estimation model (FPEM) combines a Bayesian hierarchical model with country-specific time trends to yield estimates of contraceptive prevalence and unmet need for family planning for countries worldwide. In this paper, we introduce the R package fpemlocal that carries out the estimation of family planning indicators for a single population, for example, for a single country or smaller area. In this implementation of FPEM, all non-population-specific parameters are fixed at outcomes obtained in a prior global FPEM run. The development of this model was motivated by the demand for computational efficiency, without loss of model accuracy, when estimates and projections from FPEM were needed only for a single country. We present use cases to produce estimates for a single population of women by union status or all women based on package-provided data bases and user-specified data. We also explain how to aggregate estimates across multiple populations. The R package forms the basis of the Track20 Family Planning Estimation Tool to monitor trends in family planning indicators for the FP2020 initiative. Fpemlocal is available from: https://github.com/AlkemaLab/fpemlocal.

Global, regional, and national mortality trends in youth aged 15–24 years between 1990 and 2019: a systematic analysis.

B. Masquelier, L. Hug, D. Sharrow, D. You, C. Mathers, P. Gerland, L. Alkema (2021). The Lancet Global Health. doi: 10.1016/S2214-109X(21)00023-1.




Abstract

Background The global health community is devoting considerable attention to adolescents and young people, but risk of death in this population is poorly measured. We aimed to reconstruct global, regional, and national mortality trends for youths aged 15–24 years between 1990 and 2019.
Methods In this systematic analysis, we used all publicly available data on mortality in the age group 15–24 years for 195 countries, as compiled by the UN Inter-agency Group for Child Mortality Estimation. We used nationally representative vital registration data, estimated the completeness of death registration, and extracted mortality rates from surveys with sibling histories, household deaths reported in censuses, and sample registration systems. We used a Bayesian B-spline bias-reduction model to generate trends in 10q15, the probability that an adolescent aged 15 years would die before reaching age 25 years. This model treats observations of the 10q15 probability as the product of the actual risk of death and an error multiplier that varies depending on the data source. The main outcome that we assessed was the levels of and trends in youth mortality and the global and regional mortality rates from 1990 to 2019.
Findings Globally, the probability of an individual dying between age 15 years and 24 years was 11·2 deaths (90% uncertainty interval [UI] 10·7–12·5) per 1000 youths aged 15 in 2019, which is about 2·5 times less than infant mortality (28·2 deaths [27·2–30·0] by age 1 year per 1000 live births) but is higher than the risk of dying from age 1 to 5 (9·7 deaths [9·1–11·1] per 1000 children aged 1 year). The probability of dying between age 15 years and 24 years declined by 1·4% per year (90% UI 1·1–1·8) between 1990 and 2019, from 17·1 deaths (16·5–18·9) per 1000 in 1990; by contrast with this total decrease of 34% (27–41), under-5 mortality declined by 59% (56–61) in this period. The annual number of deaths declined from 1·7 million (90% UI 1·7–1·9) in 1990 to 1·4 million (1·3–1·5) in 2019. In sub-Saharan Africa, the number of deaths increased by 20·8% from 1990 to 2019. Although 18·3% of the population aged 15–24 years were living in sub-Saharan Africa in 2019, the region accounted for 37·9% (90% UI 34·8–41·9) of all worldwide deaths in youth.
Interpretation It is urgent to accelerate progress in reducing youth mortality. Efforts are particularly needed in sub-Saharan Africa, where the burden of mortality is increasingly concentrated. In the future, a growing number of countries will see youth mortality exceeding under-5 mortality if current trends continue.

Global estimation of unintended pregnancy and abortion using a Bayesian hierarchical random walk model.

J. Bearak, A. Popinchalk, B. Ganatra, A. Moller, Ö. Tunçalp, C. Beavin, L. Kwok, L. Alkema (2020). Preprint.




Abstract

Unintended pregnancy and abortion estimates are needed to inform and motivate investment in global health programmes and policies. Variability in the availability and reliability of data poses challenges for producing estimates. We developed a Bayesian model that simultaneously estimates incidence of unintended pregnancy and abortion for 195 countries and territories. Our modelling strategy was informed by the proximate determinants of fertility with (i) incidence of unintended pregnancy defined by the number of women (grouped by marital and contraceptive use status) and their respective pregnancy rates, and (ii) abortion incidence defined by group-specific pregnancies and propensities to have an abortion. Hierarchical random walk models are used to estimate country-group-period-specific pregnancy rates and propensities to abort.

Unintended Pregnancy and Abortion by Income, Region, and the Legal Status of Abortion: Estimates from a Comprehensive Model for 1990–2019.

J. Bearak, A. Popinchalk, B. Ganatra, A. Moller, Ö. Tunçalp, C. Beavin, L. Kwok, L. Alkema (2020). The Lancet Global Health 8 (9): e1152–61.




Abstract

Background Unintended pregnancy and abortion estimates document trends in sexual and reproductive health and autonomy. These estimates inform and motivate investment in global health programmes and policies. Variability in the availability and reliability of data poses challenges for measuring and monitoring trends in unintended pregnancy and abortion. We developed a new statistical model that jointly estimated unintended pregnancy and abortion that aimed to better inform efforts towards global equity in sexual and reproductive health and rights.
Methods We developed a model that simultaneously estimated incidence of unintended pregnancy and abortion within a Bayesian framework. Data on pregnancy intentions and abortion were compiled from country-based surveys, official statistics, and published studies found through a literature search, and we obtained data on livebirths from the World Population Prospects. We analysed results by World Bank income groups, Sustainable Development Goal regional groupings, and the legal status of abortion.
Findings In 2015–19, there were 121·0 million unintended pregnancies annually (80% uncertainty interval [UI] 112·8–131·5), corresponding to a global rate of 64 unintended pregnancies (UI 60–70) per 1000 women aged 15–49 years. 61% (58–63) of unintended pregnancies ended in abortion (totalling 73·3 million abortions annually [66·7–82·0]), corresponding to a global abortion rate of 39 abortions (36–44) per 1000 women aged 15–49 years. Using World Bank income groups, we found an inverse relationship between unintended pregnancy and income, whereas abortion rates varied non-monotonically across groups. In countries where abortion was restricted, the proportion of unintended pregnancies ending in abortion had increased compared with the proportion for 1990–94, and the unintended pregnancy rates were higher than in countries where abortion was broadly legal.
Interpretation Between 1990–94 and 2015–19, the global unintended pregnancy rate has declined, whereas the proportion of unintended pregnancies ending in abortion has increased. As a result, the global average abortion rate in 2015–19 was roughly equal to the estimates for 1990–94. Our findings suggest that people in high-income countries have better access to sexual and reproductive health care than those in low-income countries. Our findings indicate that individuals seek abortion even in settings where it is restricted. These findings emphasise the importance of ensuring access to the full spectrum of sexual and reproductive health services, including contraception and abortion care, and for additional investment towards equity in health-care services.

What increase in modern contraceptive use is needed in FP2020 countries to reach 75% demand satisfied by 2030? An assessment using the Accelerated Transition Method and Family Planning Estimation Model.

N. Cahill, M. Weinberger, L. Alkema (2020). Gates Open Research 4 (113).




Abstract

Background: Sustainable Development Goal 3.7 aims to ensure universal access to sexual and reproductive health services. One suggested benchmark is to have at least 75% of the demand for contraception satisfied with modern methods (DS) in all countries by 2030. The translation of DS-based targets into targets for the modern contraceptive prevalence rate (mCPR) is needed to make targets actionable.
Methods: We propose the Accelerated Transition (AT) method for determining the mCPR needed to reach demand-satisfied targets by 2030. The starting point for this method is the projection of DS under “business as usual” using the one-country implementation of the Family Planning Estimation Model (FPEMcountry). For countries in which the DS target is projected to be later than 2030, the AT method assumes that meeting the DS target by 2030 requires an acceleration of the contraceptive use transition such that the DS target, and its associated mCPR, will be reached in 2030 as opposed to the later year. The DS-target-associated mCPR becomes the mCPR target for the year 2030.
Results: We apply the AT method to assess progress needed for attaining the 75% DS target for married or in-union women in the world’s poorest countries. For 50 out of 68 countries, we estimate that accelerations are needed, with required mCPR increases ranging from 4.3 to 50.8 percentage points.
Conclusions: The AT method quantifies the acceleration needed – as compared to business as usual projections – for a country to meet a family planning target. The method can be used to determine the mCPR needed to reach demand-satisfied targets.

A systematic assessment of the sex ratio at birth for all countries and estimation of national imbalances and regional reference levels.

F. Chao, A. Cook, P. Gerland, L. Alkema (2019). 9303 Proceedings of the National Academy of Sciences 116(19), 9303-9311.




Abstract

The sex ratio at birth (SRB; ratio of male to female live births) imbalance in parts of the world over the past few decades is a direct consequence of sex-selective abortion, driven by the coexistence of son preference, readily available technology of prenatal sex determination, and fertility decline. Estimation of the degree of SRB imbalance is complicated because of unknown SRB reference levels and because of the uncertainty associated with SRB observations. There are needs for reproducible methods to construct SRB estimates with uncertainty, and to assess SRB inflation due to sex-selective abortion. We compile an extensive database from vital registration systems, censuses and surveys with 10,835 observations, and 16,602 country-years of information from 202 countries. We develop Bayesian methods for SRB estimation for all countries from 1950 to 2017. We model the SRB regional and national reference levels, the fluctuation around national reference levels, and the inflation. The estimated regional reference levels range from 1.031 (95% uncertainty interval [1.027; 1.036]) in sub-Saharan Africa to 1.063 [1.055; 1.072] in southeastern Asia, 1.063 [1.054; 1.072] in eastern Asia, and 1.067 [1.058; 1.077] in Oceania. We identify 12 countries with strong statistical evidence of SRB imbalance during 1970–2017, resulting in 23.1 [19.0; 28.3] million missing female births globally. The majority of those missing female births are in China, with 11.9 [8.5; 15.8] million, and in India, with 10.6 [8.0; 13.6] million.

National, regional, and global levels and trends in neonatal mortality between 1990 and 2015, with scenario-based projections to 2030: a systematic analysis by the United Nations Inter-agency Group for Child Mortality Estimation.

L. Hug, M. Alexander, D. You, L. Alkema (2019). The Lancet Global Health 7(9), e710-e720.




Abstract

Background Reducing neonatal mortality is an essential part of the third Sustainable Development Goal (SDG), to end preventable child deaths. To achieve this aim will require an understanding of the levels of and trends in neonatal mortality. We therefore aimed to estimate the levels of and trends in neonatal mortality by use of a statistical model that can be used to assess progress in the SDG era. With these estimates of neonatal mortality between 1990 and 2017, we then aimed to assess how different targets for neonatal mortality could affect the burden of neonatal mortality from 2018 to 2030.
Methods In this systematic analysis, we used nationally-representative empirical data related to neonatal mortality, including data from vital registration systems, sample registration systems, and household surveys, to estimate country-specific neonatal mortality rates (NMR; the probability of dying during the first 28 days of life) for all countries between 1990 (or the earliest year of available data) and 2017. For our analysis, we used all publicly available data on neonatal mortality from databases compiled annually by the UN Inter-agency Group for Child Mortality Estimation, which were extracted on or before July 31, 2018, for data relating to the period between 1950 and 2017. All nationally representative data were assessed. We used a Bayesian hierarchical penalised B-splines regression model, which allowed for data from different sources to be weighted differently, to account for variable biases and for the uncertainty in NMR to be assessed. The model simultaneously estimated a global association between NMR and under-5 mortality rate and country-specific and time-specific effects, which enabled us to identify countries with an NMR that was higher or lower than expected. Scenario-based projections were made at the county level by use of current levels of and trends in neonatal mortality and historic or annual rates of reduction that would be required to achieve national targets. The main outcome that we assessed was the levels of and trends in neonatal mortality and the global and regional NMRs from 1990 to 2017.
Findings Between 1990 and 2017, the global NMR decreased by 51% (90% uncertainty interval [UI] 46–54), from 36·6 deaths per 1000 livebirths (35·5–37·8) in 1990, to 18·0 deaths per 1000 livebirths (17·0–19·9) in 2017. The estimated number of neonatal deaths during the same period decreased from 5·0 million (4·9 million–5·2 million) to 2·5 million (2·4 million–2·8 million). Annual NMRs vary widely across the world, but west and central Africa and south Asia had the highest NMRs in 2017. All regions have reported reductions in NMRs since 1990, and most regions accelerated progress in reducing neonatal mortality in 2000–17 versus 1990–2000. Between 2018 and 2030, we project that 27·8 million children will die in their first month of life if each country maintains its current rate of reduction in NMR. If each country achieves the SDG neonatal mortality target of 12 deaths per 1000 livebirths or fewer by 2030, we project 22·7 million cumulative neonatal deaths by 2030. More than 60 countries need to accelerate their progress to reach the neonatal mortality SDG target by 2030.
Interpretation Although substantial progress has been made in reducing neonatal mortality since 1990, increased efforts to improve progress are still needed to achieve the SDG target by 2030. Accelerated improvements are most needed in the regions and countries with high NMR, particularly in sub-Saharan Africa and south Asia.

Pregnancies, abortions, and pregnancy intentions: a protocol for modeling and reporting global, regional and country estimates.

J. Bearak, A. Popinchalk, G. Sedgh, B. Ganatra, A. Moller, Ö. Tunçalp, L. Alkema (2019). Reproductive Health 16(1) 36.




Abstract

Background Estimates of pregnancies, abortions and pregnancy intentions can help assess how effectively women and couples are able to fulfil their childbearing aspirations. Abortion incidence estimates are also a necessary foundation for research on the safety of abortions performed and the consequences of unsafe abortion. Furthermore, periodic estimates of these indicators are needed to help inform policy and programmes.
Methods We will develop a Bayesian hierarchical times series model which estimates levels and trends in pregnancy rates, abortion rates, and percentages of pregnancies and births unintended for each five-year period between 1990 and 2019. The model will be informed by data on abortion incidence and the percentage of births or pregnancies that were unintended. We will develop a data classification process to be applied to all available data. Model-based estimates and associated uncertainty will take account of data sparsity and quality. Our proposed approach will advance previous work in two key ways. First, we will estimate pregnancy and abortion rates simultaneously, and model the propensity to abort an unintended pregnancy, as opposed to modeling abortion rates directly as in prior work. Secondly, we will produce estimates that are reproducible at the country level by publishing the data inputs, data classification processes and source code.
Discussion This protocol will form the basis for updated global, regional and national estimates of intended and unintended pregnancy rates, abortion rates, and the percent of unintended pregnancies ending in abortion, from 1990 to 2019.

Modern contraceptive use, unmet need, and demand satisfied among women of reproductive age who are married or in a union in the focus countries of the Family Planning 2020 initiative: a systematic analysis using the Family Planning Estimation Tool.

N. Cahill, E. Sonneveldt, J. Stover, M. Weinberger, J. Williamson, C. Wei, W. Brown, L. Alkema (2018). The Lancet 391 (10123): 870 – 882.


Background

The London Summit on Family Planning in 2012 inspired the Family Planning 2020 (FP2020) initiative and the 120×20 goal of having an additional 120 million women and adolescent girls become users of modern contraceptives in 69 of the world's poorest countries by the year 2020. Working towards achieving 120 × 20 is crucial for ultimately achieving the Sustainable Development Goals of universal access and satisfying demand for reproductive health. Thus, a performance assessment is required to determine countries' progress.

Methods

An updated version of the Family Planning Estimation Tool (FPET) was used to construct estimates and projections of the modern contraceptive prevalence rate (mCPR), unmet need for, and demand satisfied with modern methods of contraception among women of reproductive age who are married or in a union in the focus countries of the FP2020 initiative. We assessed current levels of family planning indicators and changes between 2012 and 2017. A counterfactual analysis was used to assess if recent levels of mCPR exceeded pre-FP2020 expectations.

Findings

In 2017, the mCPR among women of reproductive age who are married or in a union in the FP2020 focus countries was 45·7% (95% uncertainty interval [UI] 42·4–49·1), unmet need for modern methods was 21·6% (19·7–23·9), and the demand satisfied with modern methods was 67·9% (64·4–71·1). Between 2012 and 2017 the number of women of reproductive age who are married or in a union who use modern methods increased by 28·8 million (95% UI 5·8–52·5). At the regional level, Asia has seen the mCPR among women of reproductive age who are married or in a union grow from 51·0% (95% UI 48·5–53·4) to 51·8% (47·3–56·5) between 2012 and 2017, which is slow growth, particularly when compared with a change from 23·9% (22·9–25·0) to 28·5% (26·8–30·2) across Africa. At the country level, based on a counterfactual analysis, we found that 61% of the countries that have made a commitment to FP2020 exceeded pre-FP2020 expectations for modern contraceptive use. Country success stories include rapid increases in Kenya, Mozambique, Malawi, Lesotho, Sierra Leone, Liberia, and Chad relative to what was expected in 2012.

Interpretation

Whereas the estimate of additional users up to 2017 for women of reproductive age who are married or in a union would suggest that the 120 × 20 goal for all women is overly ambitious, the aggregate outcomes mask the diversity in progress at the country level. We identified countries with accelerated progress, that provide inspiration and guidance on how to increase the use of family planning and inform future efforts, especially in countries where progress has been poor.

Funding

The Bill & Melinda Gates Foundation, through grant support to the University of Massachusetts Amherst and Avenir Health.

Global, regional, and subregional trends in unintended pregnancy and its outcomes from 1990 to 2014: estimates from a Bayesian hierarchical model

J. Bearak, A Popinchalk, L Alkema, G Sedgh (2018). The Lancet Global Health 6(4): e380 - e389.


Background

The progress to achieve the fourth Millennium Development Goal in reducing mortality rate in children younger than 5 years since 1990 has been remarkable. However, work remains to be done in the Sustainable Development Goal era. Estimates of under-5 mortality rates at the national level can hide disparities within countries. We assessed disparities in under-5 mortality rates by household economic status in low-income and middle-income countries (LMICs).


Background

Estimates of pregnancy incidence by intention status and outcome indicate how effectively women and couples are able to fulfil their childbearing aspirations, and can be used to monitor the impact of family-planning programmes. We estimate global, regional, and subregional pregnancy rates by intention status and outcome for 1990–2014.

Methods

We developed a Bayesian hierarchical time series model whereby the unintended pregnancy rate is a function of the distribution of women across subgroups defined by marital status and contraceptive need and use, and of the risk of unintended pregnancy in each subgroup. Data included numbers of births and of women estimated by the UN Population Division, recently published abortion incidence estimates, and findings from surveys of women on the percentage of births or pregnancies that were unintended. Some 298 datapoints on the intention status of births or pregnancies were obtained for 105 countries.

Findings

Worldwide, an estimated 44% (90% uncertainty interval [UI] 42–48) of pregnancies were unintended in 2010–14. The unintended pregnancy rate declined by 30% (90% UI 21–39) in developed regions, from 64 (59–81) per 1000 women aged 15–44 years in 1990–94 to 45 (42–56) in 2010–14. In developing regions, the unintended pregnancy rate fell 16% (90% UI 5–24), from 77 (74–88) per 1000 women aged 15–44 years to 65 (62–76). Whereas the decline in the unintended pregnancy rate in developed regions coincided with a declining abortion rate, the decline in developing regions coincided with a declining unintended birth rate. In 2010–14, 59% (90% UI 54–65) of unintended pregnancies ended in abortion in developed regions, as did 55% (52–60) of unintended pregnancies in developing regions.

Interpretation

The unintended pregnancy rate remains substantially higher in developing regions than in developed regions. Sexual and reproductive health services are needed to help women avoid unintended pregnancies, and to ensure healthy outcomes for those who do experience such pregnancies.

Funding

Dutch Ministry of Foreign Affairs and UK Aid from the UK Government.

Global, regional, and national levels and trends in mortality among older children (5-9) and young adolescents (10-14) from 1990 to 2016.

B. Masquelier, L. Hug, D. Sharrow, D. You, D. Hogan, K. Hill, J. Liu, J. Pedersen, L. Alkema (2018). The Lancet Global Health 6(10): e1087-e1099.




Abstract

Background: From 1990 to 2016, the mortality of children younger than 5 years decreased by more than half, and there are plentiful data regarding mortality in this age group through which we can track global progress in reducing the under-5 mortality rate. By contrast, little is known on how the mortality risk among older children (5-9 years) and young adolescents (10-14 years) has changed in this time. We aimed to estimate levels and trends in mortality of children aged 5-14 years in 195 countries from 1990 to 2016.
Methods: In this analysis of empirical data, we expanded the United Nations Inter-agency Group for Child Mortality Estimation database containing data on children younger than 5 years with 5530 data points regarding children aged 5-14 years. Mortality rates from 1990 to 2016 were obtained from nationally representative birth histories, data on household deaths reported in population censuses, and nationwide systems of civil registration and vital statistics. These data were used in a Bayesian B-spline bias-reduction model to generate smoothed trends with 90% uncertainty intervals, to determine the probability of a child aged 5 years dying before reaching age 15 years.
Findings: Globally, the probability of a child dying between the ages 5 years and 15 years was 7·5 deaths (90% uncertainty interval 7·2-8·3) per 1000 children in 2016, which was less than a fifth of the risk of dying between birth and age 5 years, which was 41 deaths (39-44) per 1000 children. The mortality risk in children aged 5-14 years decreased by 51% (46-54) between 1990 and 2016, despite not being specifically targeted by health interventions. The annual number of deaths in this age group decreased from 1·7 million (1·7 million-1·8 million) to 1 million (0·9 million-1·1 million) in 1990-2016. In 1990-2000, mortality rates in children aged 5-14 years decreased faster than among children aged 0-4 years. However, since 2000, mortality rates in children younger than 5 years have decreased faster than mortality rates in children aged 5-14 years. The annual rate of reduction in mortality among children younger than 5 years has been 4·0% (3·6-4·3) since 2000, versus 2·7% (2·3-3·0) in children aged 5-14 years. Older children and young adolescents in sub-Saharan Africa are disproportionately more likely to die than those in other regions; 55% (51-58) of deaths of children of this age occur in sub-Saharan Africa, despite having only 21% of the global population of children aged 5-14 years. In 2016, 98% (98-99) of all deaths of children aged 5-14 years occurred in low-income and middle-income countries, and seven countries alone accounted for more than half of the total number of deaths of these children.
Interpretation: Increased efforts are required to accelerate reductions in mortality among older children and to ensure that they benefit from health policies and interventions as much as younger children.

Global estimation of neonatal mortality using a Bayesian hierarchical splines regression model

M. Alexander, L. Alkema (2018). Demographic Research 38: 335 – 372.


Background

In recent years, much of the focus in monitoring child mortality has been on assessing changes in the under-5 mortality rate (U5MR). However, as the U5MR decreases, the share of neonatal deaths (within the first month) tends to increase, warranting increased efforts in monitoring the neonatal mortality rate (NMR) in addition to the U5MR.

Objective

Data on neonatal deaths comes from a range of sources across different countries, with the amount of data available and the quality of data varying widely. Our objective in estimating the NMR globally is to combine all data sources available to obtain accurate estimates, be able to project mortality levels, and have some indication of the uncertainty in the estimates and projections.

Methods

We present a new model for estimating the NMR for countries worldwide, using a Bayesian hierarchical model framework.

Contribution

Our modeling approach offers an intuitive way to share information across different countries and time points, and incorporates different sources of error into the estimates. It also improves on previous modeling approaches by allowing for trends observed in NMR to be more driven by the data available, rather than trends in covariates.

National and regional under-5 mortality rate by economic status for low-income and middle-income countries: a systematic assessment.

F. Chao, D. You, J. Pedersen, L. Hug, L. Alkema (2018). The Lancet Global Health 6(5): e535-e547.




Abstract

Background The progress to achieve the fourth Millennium Development Goal in reducing mortality rate in children younger than 5 years since 1990 has been remarkable. However, work remains to be done in the Sustainable Development Goal era. Estimates of under-5 mortality rates at the national level can hide disparities within countries. We assessed disparities in under-5 mortality rates by household economic status in low-income and middle-income countries (LMICs).
Method We estimated country-year-specific under-5 mortality rates by wealth quintile on the basis of household wealth indices for 137 LMICs from 1990 to 2016, using a Bayesian statistical model. We estimated the association between quintile-specific and national-level under-5 mortality rates. We assessed the levels and trends of absolute and relative disparity in under-5 mortality rate between the poorest and richest quintiles, and among all quintiles.
Findings In 2016, for all LMICs (excluding China), the aggregated under-5 mortality rate was 64·6 (90% uncertainty interval [UI] 61·1–70·1) deaths per 1000 livebirths in the poorest households (first quintile), 31·3 (29·5–34·2) deaths per 1000 livebirths in the richest households (fifth quintile), and in between those outcomes for the middle quintiles. Between 1990 and 2016, the largest absolute decline in under-5 mortality rate occurred in the two poorest quintiles: 77·6 (90% UI 71·2–82·6) deaths per 1000 livebirths in the poorest quintile and 77·9 (72·0–82·2) deaths per 1000 livebirths in the second poorest quintile. The difference in under-5 mortality rate between the poorest and richest quintiles decreased significantly by 38·8 (90% UI 32·9–43·8) deaths per 1000 livebirths between 1990 and 2016. The poorest to richest under-5 mortality rate ratio, however, remained similar (2·03 [90% UI 1·94–2·11] in 1990, 1·99 [1·91–2·08] in 2000, and 2·06 [1·92–2·20] in 2016). During 1990–2016, around half of the total under-5 deaths occurred in the poorest two quintiles (48·5% in 1990 and 2000, 49·5% in 2016) and less than a third were in the richest two quintiles (30·4% in 1990, 30·5% in 2000, 29·9% in 2016). For all regions, differences in the under-5 mortality rate between the first and fifth quintiles decreased significantly, ranging from 20·6 (90% UI 15·9–25·1) deaths per 1000 livebirths in eastern Europe and central Asia to 59·5 (48·5–70·4) deaths per 1000 livebirths in south Asia. In 2016, the ratios of under-5 mortality rate in the first quintile to under-5 mortality rate in the fifth quintile were significantly above 2·00 in two regions, with 2·49 (90% UI 2·15–2·87) in east Asia and Pacific (excluding China) and 2·41 (2·05–2·80) in south Asia. Eastern and southern Africa had the smallest ratio in 2016 at 1·62 (90% UI 1·48–1·76). Our model suggested that the expected ratio of under-5 mortality rate in the first quintile to under-5 mortality rate in the fifth quintile increases as national-level under-5 mortality rate decreases.
Interpretation For all LMICs (excluding China) combined, the absolute disparities in under-5 mortality rate between the poorest and richest households have narrowed significantly since 1990, whereas the relative differences have remained stable. To further narrow the rich-and-poor gap in under-5 mortality rate on the relative scale, targeted interventions that focus on the poorest populations are needed.

Global, regional and sub-regional classification of abortions by safety: Estimates for 2010-14

B. Ganatra, C. Gerdts, C. Rossier, R. Johnson, Ö. Tunçalp, A. Assifi, G. Sedgh, S. Singh, A. Bankole, A. Popinchalk, J. Bearak, Z. Kang, L. Alkema (2017). The Lancet


Background

Global estimates of unsafe abortions have been produced for 1995, 2003, and 2008. However, reconceptualisation of the framework and methods for estimating abortion safety is needed owing to the increased availability of simple methods for safe abortion (eg, medical abortion), the increasingly widespread use of misoprostol outside formal health systems in contexts where abortion is legally restricted, and the need to account for the multiple factors that affect abortion safety.

Methods

We used all available empirical data on abortion methods, providers, and settings, and factors affecting safety as covariates within a Bayesian hierarchical model to estimate the global, regional, and subregional distributions of abortion by safety categories. We used a three-tiered categorisation based on the WHO definition of unsafe abortion and WHO guidelines on safe abortion to categorise abortions as safe or unsafe and to further divide unsafe abortions into two categories of less safe and least safe.

Findings

Of the 55· 7 million abortions that occurred worldwide each year between 2010-14, we estimated that 30·6 million (54·9%, 90% uncertainty interval 49·9-59·4) were safe, 17·1 million (30·7%, 25·5-35·6) were less safe, and 8·0 million (14·4%, 11·5-18·1) were least safe. Thus, 25·1 million (45·1%, 40·6-50·1) abortions each year between 2010 and 2014 were unsafe, with 24·3 million (97%) of these in developing countries. The proportion of unsafe abortions was significantly higher in developing countries than developed countries (49·5% vs 12·5%). When grouped by the legal status of abortion, the proportion of unsafe abortions was significantly higher in countries with highly restrictive abortion laws than in those with less restrictive laws.

Interpretation

Increased efforts are needed, especially in developing countries, to ensure access to safe abortion. The paucity of empirical data is a limitation of these findings. Improved in-country data for health services and innovative research to address these gaps are needed to improve future estimates.

Funding

UNDP/UNFPA/UNICEF/WHO/World Bank Special Programme of Research, Development and Research Training in Human Reproduction; David and Lucile Packard Foundation; UK Aid from the UK Government; Dutch Ministry of Foreign Affairs; Norwegian Agency for Development Cooperation.

A Bayesian approach to the global estimation of maternal mortality

L. Alkema, S. Zhang, D. Chou, A. Gemmill, A.B. Moller, D.M. Fat, L. Say, C.D. Mathers, D. Hogan (2017). The Annals of Applied Statistics 11(3): 1245 – 1274.


Summary

The maternal mortality ratio (MMR) is defined as the number of maternal deaths in a population per 100,000 live births. Country-specific MMR estimates are published on a regular basis by the United Nations Maternal Mortality Estimation Inter-agency Group (UN MMEIG) to track progress in reducing maternal deaths and to evaluate regional and national performance related to Millennium Development Goal (MDG) 5, which calls for a 75% reduction in the MMR between 1990 and 2015.

Until 2014, the UN MMEIG used a multilevel regression model for producing estimates for countries without sufficient data from vital registration systems. While this model worked well in the past to assess MMR levels for countries with limited data, it was deemed unsatisfactory for final MDG 5 reporting for countries where longer time series of observations had become available because by construction, estimated trends in the MMR were covariate-driven only and did not necessarily track data-driven trends.

We developed a Bayesian maternal mortality estimation model, which extends upon the UN MMEIG multilevel regression model. The new model assesses data-driven trends through the inclusion of an ARIMA time series model that captures accelerations and decelerations in the rate of change in the MMR. Varying reporting and data quality issues are accounted for in source-specific data models. The revised model provides data-driven estimates of MMR levels and trends and will be used for MDG 5 reporting for all countries.

Setting ambitious yet achievable targets using probabilistic projections: meeting demand for family planning

V. Kantorova, J.R. New, A. Biddlecom, L. Alkema (2017). Studies in Family Planning 48(3): 223 – 233.


Abstract

In 2015, governments adopted 17 internationally agreed goals to ensure progress and well-being in the economic, social, and environmental dimensions of sustainable development. These new goals present a challenge for countries to set empirical targets that are ambitious yet achievable and that can account for different starting points and rates of progress. We used probabilistic projections of family planning indicators, based on a global data set and Bayesian hierarchical modeling, to generate illustrative targets at the country level. Targets were defined as the percentage of demand for family planning satisfied with modern contraceptive methods where a country has at least a 10 percent chance of reaching the target by 2030. National targets for 2030 ranged from below 50 percent of demand satisfied with modern contraceptives (for three countries in Africa) to above 90 percent (for 41 countries from all major areas of the world). The probabilistic approach also identified countries for which a global fixed target value of 75 percent demand satisfied was either unambitious or has little chance of achievement. We present the web-based Family Planning Estimation Tool (FPET) enabling national decision makers to compute and assess targets for meeting family planning demand.

Subnational Rates and Trends in Contraceptive Prevalence and Unmet Need for Family Planning from 1990 to 2020: An Analysis for All 29 States in India

J.R. New, N. Cahill, J. Stover, Y.P. Gupta, L. Alkema (2017). The Lancet Global Health 5(3): e350-e358.


Background

Improving access to reproductive health services and commodities is central to development. Efforts to assess progress on this front have been largely focused on national estimates, but such analyses can mask local disparities. We assessed progress in reproductive health services subnationally in India.

Findings

There is a large amount of heterogeneity in India, with a difference of up to 55·1 percentage points (95% uncertainty interval 46·4-62·1) in modern contraceptive use in 2015 between subregions. States such as Andhra Pradesh, with 92·7% (90·9-94·2) demand satisfied with modern methods, are performing well above the national average (71·8%, 56·7-83·6), whereas Manipur, with 26·8% (16·7-38·5) of demand satisfied, and Meghalaya, with 45·0% (40·1-50·0), consistently lag behind the rest of the country. Manipur and Meghalaya require the highest percentage increase in modern contraceptive use to achieve 75% demand satisfied with modern methods by 2030. In terms of absolute numbers, Uttar Pradesh requires the greatest increase, needing 9·2 million (5·5-12·6 million) additional users of modern contraception by 2030 to meet the target of 75%.

Interpretation

The demand for family planning among the states and union territories in India is highly diverse. Greatest attention is needed in Uttar Pradesh, Manipur, and Meghalaya to meet UN targets. The analysis can be generalised to other countries as well as other subpopulations.

Funding

Avenir Health through a grant from the Bill & Melinda Gates Foundation.

Abortion incidence between 1990 and 2014: global, regional, and subregional levels and trends

G. Sedgh, J. Bearak, S. Singh, A. Bankole, A. Popinchalk, B. Ganatra, C. Rossier, C. Gerdts, Ö. Tunçalp, R. Johnson, H.B. Johnston, L. Alkema (2016). The Lancet 388(10041): 258 – 267.


Background

Information about the incidence of induced abortion is needed to motivate and inform efforts to help women avoid unintended pregnancies and to monitor progress toward that end. We estimate subregional, regional, and global levels and trends in abortion incidence for 1990 to 2014, and abortion rates in subgroups of women. We use the results to estimate the proportion of pregnancies that end in abortion and examine whether abortion rates vary in countries grouped by the legal status of abortion.

Methods

We requested abortion data from government agencies and compiled data from international sources and nationally representative studies. With data for 1069 country-years, we estimated incidence using a Bayesian hierarchical time series model whereby the overall abortion rate is a function of the modelled rates in subgroups of women of reproductive age defined by their marital status and contraceptive need and use, and the sizes of these subgroups.

Findings

We estimated that 35 abortions (90% uncertainty interval [UI] 33 to 44) occurred annually per 1000 women aged 15–44 years worldwide in 2010–14, which was 5 points less than 40 (39–48) in 1990–94 (90% UI for decline −11 to 0). Because of population growth, the annual number of abortions worldwide increased by 5·9 million (90% UI −1·3 to 15·4), from 50·4 million in 1990–94 (48·6 to 59·9) to 56·3 million (52·4 to 70·0) in 2010–14. In the developed world, the abortion rate declined 19 points (–26 to −14), from 46 (41 to 59) to 27 (24 to 37). In the developing world, we found a non-significant 2 point decline (90% UI −9 to 4) in the rate from 39 (37 to 47) to 37 (34 to 46). Some 25% (90% UI 23 to 29) of pregnancies ended in abortion in 2010–14. Globally, 73% (90% UI 59 to 82) of abortions were obtained by married women in 2010–14 compared with 27% (18 to 41) obtained by unmarried women. We did not observe an association between the abortion rates for 2010–14 and the grounds under which abortion is legally allowed.

Interpretation

Abortion rates have declined significantly since 1990 in the developed world but not in the developing world. Ensuring access to sexual and reproductive health care could help millions of women avoid unintended pregnancies and ensure access to safe abortion.

Funding

UK Government, Dutch Ministry of Foreign Affairs, Norwegian Agency for Development Cooperation, The David and Lucile Packard Foundation, UNDP/UNFPA/UNICEF/WHO/World Bank Special Programme of Research, Development and Research Training in Human Reproduction.

Global, regional, and national levels and trends in maternal mortality between 1990 and 2015, with scenario-based projections to 2030: a systematic analysis by the UN Inter-agency Group for Maternal Mortality Estimation

L. Alkema, D. Chou, D. Hogan, S. Zhang, A.B. Moller, A. Gemmill, D.M. Fat, T. Boerma, M. Temmerman, C.D. Mathers, L. Say (2016). The Lancet 387(10017): 462 – 474.


Background

Millennium Development Goal 5 calls for a 75% reduction in the maternal mortality ratio (MMR) between 1990 and 2015. We estimated levels and trends in maternal mortality for 183 countries to assess progress made. Based on MMR estimates for 2015, we constructed projections to show the requirements for the Sustainable Development Goal (SDG) of less than 70 maternal deaths per 100 000 livebirths globally by 2030.

Methods

We updated the UN Maternal Mortality Estimation Inter-Agency Group (MMEIG) database with more than 200 additional records (vital statistics from civil registration systems, surveys, studies, or reports). We generated estimates of maternal mortality and related indicators with 80% uncertainty intervals (UIs) using a Bayesian model. The model combines the rate of change implied by a multilevel regression model with a time-series model to capture data-driven changes in country-specific MMRs, and includes a data model to adjust for systematic and random errors associated with different data sources.

Results

We had data for 171 of 183 countries. The global MMR fell from 385 deaths per 100 000 livebirths (80% UI 359-427) in 1990, to 216 (207-249) in 2015, corresponding to a relative decline of 43 9% (34 0-48 7), with 303 000 (291 000-349 000) maternal deaths worldwide in 2015. Regional progress in reducing the MMR since 1990 ranged from an annual rate of reduction of 1.8% (0.0-3.1) in the Caribbean to 5.0% (4.0-6.0) in eastern Asia. Regional MMRs for 2015 ranged from 12 deaths per 100 000 livebirths (11-14) for high-income regions to 546 (511-65 2) for sub-Saharan Africa. Accelerated progress will be needed to achieve the SDG goal; countries will need to reduce their MMRs at an annual rate of reduction of at least 7.5%.

Interpretation

Despite global progress in reducing maternal mortality, immediate action is needed to meet the ambitious SDG 2030 target, and ultimately eliminate preventable maternal mortality. Although the rates of reduction that are needed to achieve country-specific SDG targets are ambitious for most high mortality countries, countries that made a concerted effort to reduce maternal mortality between 2000 and 2010 provide inspiration and guidance on how to accomplish the acceleration necessary to substantially reduce preventable maternal deaths.

Global, regional, and national levels and trends in under-5 mortality between 1990 and 2015, with scenario-based projections to 2030: a systematic analysis by the UN Inter-agency Group for Child Mortality Estimation

D. You, L. Hug, S. Ejdemyr, P. Idele, D. Hogan, C. Mathers, P. Gerland, J.R. New, L. Alkema (2015). The Lancet 386(10010): 2275–2286.


Background

In 2000, world leaders agreed on the Millennium Development Goals (MDGs). MDG 4 called for a two-thirds reduction in the under-5 mortality rate between 1990 and 2015. We aimed to estimate levels and trends in under-5 mortality for 195 countries from 1990 to 2015 to assess MDG 4 achievement and then intended to project how various post-2015 targets and observed rates of change will affect the burden of under-5 deaths from 2016 to 2030.

Methods

We updated the UN Inter-agency Group for Child Mortality Estimation (UN IGME) database with 5700 country-year datapoints. As of July, 2015, the database contains about 17 000 country-year datapoints for mortality of children younger than 5 years for 195 countries, and includes all available nationally-representative data from vital registration systems, population censuses, household surveys, and sample registration systems. We used these data to generate estimates, with uncertainty intervals, of under-5 (age 0-4 years) mortality using a Bayesian B-spline bias-reduction model (B3 model). This model includes a data model to adjust for systematic biases associated with different types of data sources. To provide insights into the global and regional burden of under-5 deaths associated with post-2015 targets, we constructed five scenario-based projections for under-5 mortality from 2016 to 2030 and estimated national, regional, and global under-5 mortality rates up to 2030 for each scenario.

Results

The global under-5 mortality rate has fallen from 90·6 deaths per 1000 livebirths (90% uncertainty interval 89·3-92·2) in 1990 to 42·5 (40·9-45·6) in 2015. During the same period, the annual number of under-5 deaths worldwide dropped from 12·7 million (12·6 million-13·0 million) to 5·9 million (5·7 million-6·4 million). The global under-5 mortality rate reduced by 53% (50-55%) in the past 25 years and therefore missed the MDG 4 target. Based on point estimates, two regions-east Asia and the Pacific, and Latin America and the Caribbean-achieved the MDG 4 target. 62 countries achieved the MDG 4 target, of which 24 were low-income and lower-middle income countries. Between 2016 and 2030, 94·4 million children are projected to die before the age of 5 years if the 2015 mortality rate remains constant in each country, and 68·8 million would die if each country continues to reduce its mortality rate at the pace estimated from 2000 to 2015. If all countries achieve the Sustainable Development Goal of an under-5 mortality rate of 25 or fewer deaths per 1000 livebirths by 2030, we project 56·0 million deaths by 2030. About two-thirds of all sub-Saharan African countries need to accelerate progress to achieve this target.

Interpretation

Despite substantial progress in reducing child mortality, concerted efforts remain necessary to avoid preventable under-5 deaths in the coming years and to accelerate progress in improving child survival further. Urgent actions are needed most in the regions and countries with high under-5 mortality rates, particularly those in sub-Saharan Africa and south Asia.

Funding No funding

The United Nations Probabilistic Population Projections: An Introduction to Demographic Forecasting with Uncertainty.

L. Alkema, P. Gerland, A. Raftery, J. Wilmoth (2015). Foresight: The International Journal of Applied Forecasting 37: 19–24.


The United Nations publishes projections of populations around the world and breaks these down by age and sex. Traditionally, they are produced with standard demographic methods based on assumptions about future fertility rates, survival probabilities, and migration counts. Such projections, however, were not accompanied by formal statements of uncertainty expressed in probabilistic terms. In July 2014 the UN for the first time issued official probabilistic population projections for all countries to 2100. These projections quantify uncertainty associated with future fertility and mortality trends worldwide.

This review article summarizes the probabilistic population projection methods and presents forecasts for population growth over the rest of this century.

Global estimation of child mortality using a Bayesian B-spline bias-reduction method

L. Alkema, J.R. New (2014). The Annals of Applied Statistics 8(4): 2122–2149.


World fertpopulation stabilization unlikely this century

P. Gerland, A. Raftery, H. Ševčíková, N. Li, D. Gu, T. Spoorenberg, L. Alkema, B. Fosdick, J. Chunn, N. Lalic, G. Bay, T. Buettner, G.K. Heilig, J. Wilmoth (2014). Science 346(6206): 234–237.


Abstract

The United Nations recently released population projections based on data until 2012 and a Bayesian probabilistic methodology. Analysis of these data reveals that, contrary to previous literature, world population is unlikely to stop growing this century. There is an 80% probability that world population, now 7.2 billion, will increase to between 9.6 and 12.3 billion in 2100. This uncertainty is much smaller than the range from the traditional UN high and low variants. Much of the increase is expected to happen in Africa, in part due to higher fertility and a recent slowdown in the pace of fertility decline. Also, the ratio of working age people to older people is likely to decline substantially in all countries, even those that currently have young populations.

National, regional, and global sex ratios of infant, child, and under-5 mortality and identification of countries with outlying ratios: a systematic assessment

L. Alkema, F. Chao, D. You, J. Pedersen, C.C. Sawyer (2014). The Lancet Global Health 2(9): e521–e530.


Background

Under natural circumstances, the sex ratio of male to female mortality up to the age of 5 years is greater than one but sex discrimination can change sex ratios. The estimation of mortality by sex and identification of countries with outlying levels is challenging because of issues with data availability and quality, and because sex ratios might vary naturally based on differences in mortality levels and associated cause of death distributions.

Methods

For this systematic analysis, we estimated country-specific mortality sex ratios for infants, children aged 1–4 years, and children under the age of 5 years (under 5s) for all countries from 1990 (or the earliest year of data collection) to 2012 using a Bayesian hierarchical time series model, accounting for various data quality issues and assessing the uncertainty in sex ratios. We simultaneously estimated the global relation between sex ratios and mortality levels and constructed estimates of expected and excess female mortality rates to identify countries with outlying sex ratios.

Findings

Global sex ratios in 2012 were 1·13 (90% uncertainty interval 1·12–1·15) for infants, 0·95 (0·93–0·97) for children aged 1–5 years, and 1·08 (1·07–1·09) for under 5s, an increase since 1990 of 0·01 (−0·01 to 0·02) for infants, 0·04 (0·02 to 0·06) for children aged 1–4 years, and 0·02 (0·01 to 0·04) for under 5s. Levels and trends varied across regions and countries. Sex ratios were lowest in southern Asia for 1990 and 2012 for all age groups. Highest sex ratios were seen in developed regions and the Caucasus and central Asia region. Decreasing mortality was associated with increasing sex ratios, except at very low infant mortality, where sex ratios decreased with total mortality. For 2012, we identified 15 countries with outlying under-5 sex ratios, of which ten countries had female mortality higher than expected (Afghanistan, Bahrain, Bangladesh, China, Egypt, India, Iran, Jordan, Nepal, and Pakistan). Although excess female mortality has decreased since 1990 for the vast majority of countries with outlying sex ratios, the ratios of estimated to expected female mortality did not change substantially for most countries, and worsened for India.

Interpretation

Important differences exist between boys and girls with respect to survival up to the age of 5 years. Survival chances tend to improve more rapidly for girls compared with boys as total mortality decreases, with a reversal of this trend at very low infant mortality. For many countries, sex ratios follow this pattern but important exceptions exist. An explanation needs to be sought for selected countries with outlying sex ratios and action should be undertaken if sex discrimination is present.

Funding

The National University of Singapore and the United Nations Children's Fund (UNICEF).

Child Mortality Estimation 2013: An overview of updates in estimation methods by the United Nations Inter-agency Group for Child Mortality Estimation

L. Alkema, J.R. New, J. Pedersen, D. You, on behalf of the members of the UN Inter-agency Group for Child Mortality Estimation and its Technical Advisory Group (2014). PLOS ONE 9(7): e101112.


Global causes of maternal deaths: A WHO systematic analysis

L. Say, D. Chou, A. Gemmill, Ö. Tunçalp, A. B. Moller, J. Daniels, A. M. Gülmezoglu, M. Temmerman, L. Alkema (2014). The Lancet Global Health 2(6): e323–e3332.


Background

Data for the causes of maternal deaths are needed to inform policies to improve maternal health. We developed and analysed global, regional, and subregional estimates of the causes of maternal death during 2003-09, with a novel method, updating the previous WHO systematic review.

Methods

We searched specialised and general bibliographic databases for articles published between between Jan 1, 2003, and Dec 31, 2012, for research data, with no language restrictions, and the WHO mortality database for vital registration data. On the basis of prespecified inclusion criteria, we analysed causes of maternal death from datasets. We aggregated country level estimates to report estimates of causes of death by Millennium Development Goal regions and worldwide, for main and subcauses of death categories with a Bayesian hierarchical model.

Findings

We identified 23 eligible studies (published 2003-12). We included 417 datasets from 115 countries comprising 60 799 deaths in the analysis. About 73% (1 771 000 of 2 443 000) of all maternal deaths between 2003 and 2009 were due to direct obstetric causes and deaths due to indirect causes accounted for 27·5% (672 000, 95% UI 19·7-37·5) of all deaths. Haemorrhage accounted for 27·1% (661 000, 19·9-36·2), hypertensive disorders 14·0% (343 000, 11·1-17·4), and sepsis 10·7% (261 000, 5·9-18·6) of maternal deaths. The rest of deaths were due to abortion (7·9% [193 000], 4·7-13·2), embolism (3·2% [78 000], 1·8-5·5), and all other direct causes of death (9·6% [235 000], 6·5-14·3). Regional estimates varied substantially.

Interpretation

Between 2003 and 2009, haemorrhage, hypertensive disorders, and sepsis were responsible for more than half of maternal deaths worldwide. More than a quarter of deaths were attributable to indirect causes. These analyses should inform the prioritisation of health policies, programmes, and funding to reduce maternal deaths at regional and global levels. Further efforts are needed to improve the availability and quality of data related to maternal mortality.

Forecasting the burden of type 2 diabetes in Singapore using a demographic epidemiological model of Singapore

T.P. Phan, L. Alkema, E. Tai, K. Tan, Q. Yang, W. Lim, Y. Teo, C. Cheng, X. Wang, T. Wong, K. Chia, A. Cook (2014). BMJ Open Diabetes Research and Care 2(1): e000012.


Objective

Singapore is a microcosm of Asia as a whole, and its rapidly ageing, increasingly sedentary population heralds the chronic health problems other Asian countries are starting to face and will likely face in the decades ahead. Forecasting the changing burden of chronic diseases such as type 2 diabetes in Singapore is vital to plan the resources needed and motivate preventive efforts.

Methods

This paper describes an individual-level simulation model that uses evidence synthesis from multiple data streams-national statistics, national health surveys, and four cohort studies, and known risk factors-aging, obesity, ethnicity, and genetics-to forecast the prevalence of type 2 diabetes in Singapore. This comprises submodels for mortality, fertility, migration, body mass index trajectories, genetics, and workforce participation, parameterized using Markov chain Monte Carlo methods, and permits forecasts by ethnicity and employment status.

Results

We forecast that the obesity prevalence will quadruple from 4.3% in 1990 to 15.9% in 2050, while the prevalence of type 2 diabetes (diagnosed and undiagnosed) among Singapore adults aged 18-69 will double from 7.3% in 1990 to 15% in 2050, that ethnic Indians and Malays will bear a disproportionate burden compared with the Chinese majority, and that the number of patients with diabetes in the workforce will grow markedly.

Conclusion

If the recent rise in obesity prevalence continues, the lifetime risk of type 2 diabetes in Singapore will be one in two by 2050 with concomitant implications for greater healthcare expenditure, productivity losses, and the targeting of health promotion programmes.

Bayesian population projections for the United Nations

A. Raftery, L. Alkema, P. Gerland (2014). Statistical Science 29(1): 56–68.


Abstract

The United Nations regularly publishes projections of the populations of all the world’s countries broken down by age and sex. These projections are the de facto standard and are widely used by international organizations, governments and researchers. Like almost all other population projections, they are produced using the standard deterministic cohort-component projection method and do not yield statements of uncertainty. We describe a Bayesian method for producing probabilistic population projections for most countries that the United Nations could use. It has at its core Bayesian hierarchical models for the total fertility rate and life expectancy at birth. We illustrate the method and show how it can be extended to address concerns about the UN’s current assumptions about the long-term distribution of fertility. The method is implemented in the R packages bayesTFR, bayesLife, bayesPop and bayesDem.

How informative are vital registration data for estimating maternal mortality? A Bayesian analysis of WHO adjustment data and parameters

F. Chao, L. Alkema (2014). Statistics and Public Policy 1(1): 6–14.


Summary

Monitoring maternal mortality is challenging due to fragmented data of varying quality. The maternal mortality estimates published by the WHO in 2012 included data adjustment parameters to account for these data quality issues, but there was a discrepancy between the WHO assumption about, and the observed variability in, misclassification errors in vital registration (VR) observations. We developed a Bayesian hierarchical time series model to estimate the extent of VR misclassification errors and to provide a plausible assessment of the uncertainty associated with VR observations for countries with and without external information on VR adjustment parameters. The resulting Bayesian distribution for VR adjustments was more comparable to the observed biases than the WHO expert distribution and the model allows for estimation of VR adjustment values for any period of interest for countries with partial information on such adjustments. We also illustrated that a fully Bayesian modeling approach for estimating maternal mortality can provide more data-driven insights into maternal mortality estimates and data adjustment parameters. However, given the paucity of, and the issues with, maternal mortality data, validation of modeling assumptions and fin

Millennium Development Goals national targets are moving targets and the results will not be known until well after the 5 deadline of 2015.

M. Oestergaard, L. Alkema, J. Lawn (Editorial, 2013). International Journal of Epidemiology 42(3): 645–647.


Summary

The Millennium Development Goal 4 (MDG-4) was agreed upon at the United Nations Millennium Summit in September 2000 to accelerate national and international efforts to reduce child mortality, and improve development by setting explicit targets. However, there is very little awareness among global health data experts, let alone policymakers, that in reality, the Millennium Development Goals related to mortality outcomes are moving targets that may change every year.

The MDG-4 target is the under-five mortality rate (U5MR, deaths in children less than 5 years of age per 1000 live births) in 2015 given a two-thirds reduction in the rate compared with 1990. The target is set for each country of the world and at the global level and is determined by the estimated baseline mortality levels for year 1990.

We explain why MDG-4 targets are moving and summarize indicators of progress, quantify how much national and global targets have moved and outline key policy implications and considerations for tracking country progress. Our aim is to strengthen the process of evaluating progress and to inform the debate on a next generation of goals that are likely to succeed when the MDG period ends in 2015. We illustrate changes in targets by using MDG4 predictions from the United Nations Inter-agency Group on Child Mortality Estimation (UN-IGME),1,2 but the same principles and implications apply to any burden of disease prediction exercise, for example for MDG-5.3-5

The UN-IGME, comprising the United Nations Children’s Fund, World Health Organization, United Nations Population Division and the World Bank, was formed in 2004 with a brief to produce annual updates of child mortality levels and trends for all 195 UN member states in order to facilitate consistent tracking for MDG-4.6 UN-IGME has identified the annual rate of change as the indicator of choice for tracking progress towards the MDG-4 target. A twothirds reduction in U5MR between 1990 and 2015 amounts to an annual rate of change of 4.4%, so if a country’s annual rate of change is 4.4% or greater, the country is said to be on target. With lower rates than this, the country is off target. Both the progress between 1990 and today and the required change going forward to reach the MDG target are changing every year as these are determined by mortality levels for both 1990 and today. Paradoxically, this means that the required change today to reach the MDG target may differ from previous years even if the 1990 level and thus the MDG-4 target stay the same; or it may mean that the required change fluctuates between years even if the mortality level in recent years stays constant but the estimated 1990 level changes.

National, regional and global rates and trends in contraceptive prevalence and unmet need for family planning between 1990 and 2015: a systematic and comprehensive analysis

L. Alkema, V. Kantorova, C. Menozzi, A. Biddlecom (2013). The Lancet 381(9878): 1642–1652.


Background

Expansion of access to contraception and reduction of unmet n eed for family planning are key components to improve reproductive health, but scarce data and variability in data sources create diffi culties in monitoring of progress for these outcomes. We estimated and projected indicators of contraceptive prevalence and unmet need for family planning from 1990 to 2015.

Methods

We obtained data from nationally representative surveys, for women aged 15–49 years who were married or in a union. Estimates were based on 930 observations of contraceptive prevalence between 1950 and 2011 from 194 countries or areas, and 306 observations of unmet need for family planning from 111 countries or areas. We used a Bayesian hierarchical model combined with country-specifi c time trends to yield estimates of these indicators and uncertainty assessments. The model accounted for diff erences by data source, sample population, and contraceptive methods included in the measure.

Findings

Worldwide, contraceptive prevalence increased from 54·8% (95% uncertainty interval 52·3–57·1) in 1990, to 63·3% (60·4–66·0) in 2010, and unmet need for family planning decreased from 15·4% (14·1–16·9) in 1990, to 12·3% (10·9–13·9) in 2010. Almost all subregions, except for those where contraceptive prevalence was already high in 1990, had an increase in contraceptive prevalence and a decrease in unmet need for family planning between 1990 and 2010, although the pace of change over time varied between countries and subregions. In 2010, 146 million (130–166 million) women worldwide aged 15–49 years who were married or in a union had an unmet need for family planning. The absolute number of married women who either use contraception or who have an unmet need for family planning is projected to grow from 900 million (876–922 million) in 2010 to 962 million (927–992 million) in 2015, and will increase in most developing countries.

Interpretation

Trends in contraceptive prevalence and unmet need for family planning, and the projected growth in the number of potential contraceptive users indicate that increased investment is necessary to meet demand for contraceptive methods and improve reproductive health worldwide.

Funding

United Nations Population Division and National University of Singapore.

Progress toward global reduction in under-5 mortality: A bootstrap analysis of uncertainty in Millennium Development Goal 4 estimates

L. Alkema, J.R. New (2012). PLOS Medicine 9(12): e1001355.


Background

Millennium Development Goal 4 calls for an annual rate of reduction (ARR) of the under-five mortality rate (U5MR) of 4.4% between 1990 and 2015. Progress is measured through the point estimates of the United Nations Inter-agency Group for Child Mortality Estimation (UN IGME). To facilitate evidence-based conclusions about progress toward the goal, we assessed the uncertainty in the estimates arising from sampling errors and biases in data series and the inferior quality of specific data series.

Methods and Findings

We implemented a bootstrap procedure to construct 90% uncertainty intervals (UIs) for the U5MR and ARR to complement the UN IGME estimates. We constructed the bounds for all countries without a generalized HIV epidemic, where a standard estimation approach is carried out (174 countries). In the bootstrap procedure, potential biases in levels and trends of data series of different source types were accounted for. There is considerable uncertainty about the U5MR, particularly for high mortality countries and in recent years. Among 86 countries with a U5MR of at least 40 deaths per 1,000 live births in 1990, the median width of the UI, relative to the U5MR level, was 19% for 1990 and 48% for 2011, with the increase in uncertainty due to more limited data availability. The median absolute width of the 90% UI for the ARR from 1990 to 2011 was 2.2%. Although the ARR point estimate for all high mortality countries was greater than zero, for eight of them uncertainty included the possibility of no improvement between 1990 and 2011. For 13 countries, it is deemed likely that the ARR from 1990 to 2011 exceeded 4.4%.

Conclusion

In light of the upcoming evaluation of Millennium Development Goal 4 in 2015, uncertainty assessments need to be taken into account to avoid unwarranted conclusions about countries' progress based on limited data.

Child Mortality Estimation: a comparison of UN-IGME and IHME estimates of levels and trends in under-5 mortality rates and deaths

L. Alkema, D. You (2012). PLOS Medicine 9(8): e1001288.


Background

Millennium Development Goal 4 calls for a reduction in the under-five mortality rate (U5MR) by two-thirds between 1990 and 2015. In 2011, estimates were published by the United Nations Inter-agency Group for Child Mortality Estimation (UN IGME) and the Institute for Health Metrics and Evaluation (IHME). The difference in the U5MR estimates produced by the two research groups was more than 10% and corresponded to more than ten deaths per 1,000 live births for 10% of all countries in 1990 and 20% of all countries in 2010, which can lead to conflicting conclusions with respect to countries' progress. To understand what caused the differences in estimates, we summarised differences in underlying data and modelling approaches used by the two groups, and analysed their effects.

Methods and Findings

UN IGME and IHME estimation approaches differ with respect to the construction of databases and the pre-processing of data, trend fitting procedures, inclusion and exclusion of data series, and additional adjustment procedures. Large differences in U5MR estimates between the UN IGME and the IHME exist in countries with conflicts or civil unrest, countries with high HIV prevalence, and countries where the underlying data used to derive the estimates were different, especially if the exclusion of data series differed between the two research groups. A decomposition of the differences showed that differences in estimates due to using different data (inclusion of data series and pre-processing of data) are on average larger than the differences due to using different trend fitting methods.

Conclusions

Substantial country-specific differences between UN IGME and IHME estimates for U5MR and the number of under-five deaths exist because of various differences in data and modelling assumptions used. Often differences are illustrative of the lack of reliable data and likely to decrease as more data become available. Improved transparency on methods and data used will help to improve understanding about the drivers of the differences.

Monitoring progress towards Millennium Development Goal 4: A call for improved validation of under-5 mortality rate estimates

L. Alkema, M. Wong, P. Seah (2012). Statistics, Politics, and Policy 3(2): Art. 2.


Background

With 4 years until 2015, it is essential to monitor progress towards Millennium Development Goals (MDGs) 4 and 5. Although estimates of maternal and child mortality were published in 2010, an update of estimates is timely in view of additional data sources that have become available and new methods developed. Our aim was to update previous estimates of maternal and child mortality using better data and more robust methods to provide the best available evidence for tracking progress on MDGs 4 and 5.

Methods

We update the analyses of the progress towards MDGs 4 and 5 from 2010 with additional surveys, censuses, vital registration, and verbal autopsy data. For children, we estimate early neonatal (0-6 days), late neonatal (7-28 days), postneonatal (29-364 days), childhood (ages 1-4 years), and under-5 mortality. We use an improved model for estimating mortality by age under 5 years. For maternal mortality, our updated analysis includes greater than 1000 additional site-years of data. We tested a large set of alternative models for maternal mortality; we used an ensemble model based on the models with the best out-of-sample predictive validity to generate new estimates from 1990 to 2011.

Findings

Under-5 deaths have continued to decline, reaching 7·2 million in 2011 of which 2·2 million were early neonatal, 0·7 million late neonatal, 2·1 million postneonatal, and 2·2 million during childhood (ages 1-4 years). Comparing rates of decline from 1990 to 2000 with 2000 to 2011 shows that 106 countries have accelerated declines in the child mortality rate in the past decade. Maternal mortality has also continued to decline from 409,100 (uncertainty interval 382,900-437,900) in 1990 to 273,500 (256,300-291,700) deaths in 2011. We estimate that 56,100 maternal deaths in 2011 were HIV-related deaths during pregnancy. Based on recent trends in developing countries, 31 countries will achieve MDG 4, 13 countries MDG 5, and nine countries will achieve both.

Interpretations

Even though progress on reducing maternal and child mortality in most countries is accelerating, most developing countries will take many years past 2015 to achieve the targets of the MDGs 4 and 5. Similarly, although there continues to be progress on maternal mortality the pace is slow, without any overall evidence of acceleration. Immediate concerted action is needed for a large number of countries to achieve MDG 4 and MDG 5.

Funding

Bill & Melinda Gates Foundation.

Estimating trends in the total fertility rate with uncertainty using imperfect data: Examples from West Africa

L. Alkema, A. Raftery, P. Gerland, S. Clark, F. Pelletier (2012). Demographic Research 26(15): 331–362.


Background

Estimating the total fertility rate is challenging for many developing countries because of limited data and varying data quality. A standardized, reproducible approach to produce estimates that include an uncertainty assessment is desired.

Methods

We develop a method to estimate and assess uncertainty in the total fertility rate over time, based on multiple imperfect observations from different data sources, including surveys and censuses. We take account of measurement error in observations by decomposing it into bias and variance, and assess both by linear regression on a variety of data quality covariates. We estimate the total fertility rate using a local smoother, and assess uncertainty using the weighted likelihood bootstrap.

Results

We apply our method to data from seven countries in West Africa and construct estimates and uncertainty intervals for the total fertility rate. Based on cross-validation exercises, we find that accounting for differences in data quality between observations gives better calibrated confidence intervals and reduces bias.

Conclusion

When working with multiple imperfect observations from different data sources to estimate the total fertility rate, or demographic indicators in general, potential biases and differences in error variance should be taken into account to improve the estimates and their uncertainty assessment.

Estimating the under-5 mortality rate using a Bayesian hierarchical time series model

L. Alkema, W.L. Ann (2011). PLOS ONE 6(9): e23954.


Background

Millennium Development Goal 4 calls for a reduction in the under-five mortality rate by two-thirds between 1990 and 2015, which corresponds to an annual rate of decline of 4.4%. The United Nations Inter-Agency Group for Child Mortality Estimation estimates under-five mortality in every country to measure progress. For the majority of countries, the estimates within a country are based on the assumption of a piece-wise constant rate of decline.

Methods and Findings

This paper proposes an alternative method to estimate under-five mortality, such that the underlying rate of change is allowed to vary smoothly over time using a time series model. Information about the average rate of decline and changes therein is exchanged between countries using a Bayesian hierarchical model. Cross-validation exercises suggest that the proposed model provides credible bounds for the under-five mortality rate that are reasonably well calibrated during the observation period. The alternative estimates suggest smoother trends in under-five mortality and give new insights into changes in the rate of decline within countries.

Conclusions

The proposed model offers an alternative modeling approach for obtaining estimates of under-five mortality which removes the restriction of a piece-wise linear rate of decline and introduces hierarchy to exchange information between countries. The newly proposed estimates of the rate of decline in under-5 mortality and the uncertainty assessments would help to monitor progress towards Millennium Development Goal 4.

Probabilistic projections of the total fertility rate for all countries.

L. Alkema, A. Raftery, P. Gerland, S. Clark, F. Pelletier, T. Buettner, G. Heilig (2011). Demography 48(3): 815–839


Abstract

We describe a Bayesian projection model to produce country-specific projections of the total fertility rate (TFR) for all countries. The model decomposes the evolution of TFR into three phases: pre-transition high fertility, the fertility transition, and post-transition low fertility. The model for the fertility decline builds on the United Nations Population Division's current deterministic projection methodology, which assumes that fertility will eventually fall below replacement level. It models the decline in TFR as the sum of two logistic functions that depend on the current TFR level, and a random term. A Bayesian hierarchical model is used to project future TFR based on both the country's TFR history and the pattern of all countries. It is estimated from United Nations estimates of past TFR in all countries using a Markov chain Monte Carlo algorithm. The post-transition low fertility phase is modeled using an autoregressive model, in which long-term TFR projections converge toward and oscillate around replacement level. The method is evaluated using out-of-sample projections for the period since 1980 and the period since 1995, and is found to be well calibrated.

BayesTFR: An R package for probabilistic projections of the total fertility rate

H. Ševčíková, L. Alkema, A. Raftery (2011). Journal of Statistical Software 43: 1–29.


Abstract

The bayesTFR package for R provides a set of functions to produce probabilistic projections of the total fertility rate (TFR) for all countries. In the model, a random walk with drift is used to project the TFR during the fertility transition, using a Bayesian hierarchical model to estimate the parameters of the drift term. The TFR is modeled with a first order autoregressive process during the post-transition phase. The computationally intensive part of the projection model is a Markov chain Monte Carlo algorithm for estimating the parameters of the drift term. This article summarizes the projection model and describes the basic steps to generate probabilistic projections, as well as other functionalities such as projecting aggregate outcomes and dealing with missing values.

A Bayesian approach to uncertainty analysis of sexually transmitted infection models

L. F. Johnson, L. Alkema, R. E. Dorrington (2010). Sexually Transmitted Infections 86: 169–174.


Objectives

To propose a Bayesian approach to uncertainty analysis of sexually transmitted infection (STI) models, that can be used to quantify uncertainty in model assessments of policy options, estimate regional STI prevalence from sentinel surveillance data and make inferences about STI transmission and natural history parameters.

Methods

Prior distributions are specified to represent uncertainty regarding STI parameters. A likelihood function is defined using a hierarchical approach that takes account of variation between study populations, variation in diagnostic accuracy as well as random binomial variation. The method is illustrated using a model of syphilis, gonorrhoea, chlamydial infection and trichomoniasis in South Africa.

Results

Model estimates of STI prevalence are in good agreement with observations. Out-of-sample projections and cross-validations also show that the model is reasonably well calibrated. Model predictions of the impact of interventions are subject to significant uncertainty: the predicted reductions in the prevalence of syphilis by 2020, as a result of doubling the rate of health seeking, increasing the proportion of private practitioners using syndromic management protocols, and screening all pregnant women for syphilis, are 43% (95% CI: 3–77%), 9% (95% CI: 1–19%) and 6% (95% CI: 4–7%) respectively.

Conclusions

This study extends uncertainty analysis techniques for fitted HIV/AIDS models to models that are fitted to other STI prevalence data. There is significant uncertainty regarding the relative effectiveness of different STI control strategies. The proposed technique is reasonable for estimating uncertainty in past STI prevalence levels and for projections of future STI prevalence.

Bayesian melding for estimating uncertainty in national HIV prevalence estimates

L. Alkema, A. Raftery, T. Brown (2008). Sexually Transmitted Infections 84 (Suppl I): i11–i16.


Objective

To construct confidence intervals for HIV prevalence in countries with generalised epidemics.

Methods

In the Bayesian melding approach, a sample of country-specific epidemic curves describing HIV prevalence over time is derived based on time series of antenatal clinic prevalence data and general information on the parameters that describe the HIV epidemic. The prevalence trends at antenatal clinics are calibrated to population-based HIV prevalence estimates from national surveys. For countries without population based estimates, a general calibration method is developed. Based on the sample of calibrated epidemic curves, we derive annual 95% confidence intervals for HIV prevalence. The curve that best represents the data at antenatal clinics and population-based surveys, as well as general information about the epidemic, is chosen to represent the best estimates and predictions.

Results

We present results for urban areas in Haiti and Namibia to illustrate the estimates and confidence intervals that are derived with the methodology.

Progress and challenges in modelling country-level HIV/AIDS epidemics: the UNAIDS Estimation and Projection Package 2007

T. Brown, J. Salomon, L. Alkema, A. Raftery, E. Gouws (2008) Sexually Transmitted Infections 84 (Suppl I): i5–i11.


Abstract

The UNAIDS Estimation and Projection Package (EPP) was developed to aid in country-level estimation and short-term projection of HIV/AIDS epidemics. This paper describes advances reflected in the most recent update of this tool (EPP 2007), and identifies key issues that remain to be addressed in future versions. The major change to EPP 2007 is the addition of uncertainty estimation for generalised epidemics using the technique of Bayesian melding, but many additional changes have been made to improve the user interface and efficiency of the package. This paper describes the interface for uncertainty analysis, changes to the user interface for calibration procedures and other user interface changes to improve EPP’s utility in different settings. While formal uncertainty assessment remains an unresolved challenge in low-level and concentrated epidemics, the Bayesian melding approach has been applied to provide analysts in these settings with a visual depiction of the range of models that may be consistent with their data. In fitting the model to countries with longer-running epidemics in sub-Saharan Africa, a number of limitations have been identified in the current model with respect to accommodating behaviour change and accurately replicating certain observed epidemic patterns. This paper discusses these issues along with their implications for future changes to EPP and to the underlying UNAIDS Reference Group model.

Probabilistic projections of HIV prevalence using Bayesian melding

L. Alkema, A. Raftery, S. Clark (2007) The Annals of Applied Statistics 1(1): 229–248


Abstract

The Joint United Nations Programme on HIV/AIDS (UNAIDS) has developed the Estimation and Projection Package (EPP) for making national estimates and short-term projections of HIV prevalence based on observed prevalence trends at antenatal clinics. Assessing the uncertainty about its estimates and projections is important for informed policy decision making, and we propose the use of Bayesian melding for this purpose. Prevalence data and other information about the EPP model's input parameters are used to derive a probabilistic HIV prevalence projection, namely a probability distribution over a set of future prevalence trajectories. We relate antenatal clinic prevalence to population prevalence and account for variability between clinics using a random effects model. Predictive intervals for clinic prevalence are derived for checking the model. We discuss predictions given by the EPP model and the results of the Bayesian melding procedure for Uganda, where prevalence peaked at around 28% in 1990; the 95% prediction interval for 2010 ranges from 2% to 7%.