Erratum: "Estimated Changes in Life Expectancy and Adult Mortality Resulting from Declining PM2.5 Exposures in the Contiguous United States: 1980-2010"

[This corrects the article DOI: https://doi.org/10.1289/EHP507.].

Estimated Changes in Life Expectancy and Adult Mortality Resulting from Declining PM2.5 Exposures in the Contiguous United States: 1980-2010

BACKGROUND

PM2.5 precursor emissions have declined over the course of several decades, following the implementation of local, state, and federal air quality policies. Estimating the corresponding change in population exposure and PM2.5-attributable risk of death prior to the year 2000 is made difficult by the lack of PM2.5 monitoring data.

OBJECTIVES

We used a new technique to estimate historical PM2.5 concentrations, and estimated the effects of changes in PM2.5 population exposures on mortality in adults (age ≥30y), and on life expectancy at birth, in the contiguous United States during 1980-2010.

METHODS

We estimated annual mean county-level PM2.5 concentrations in 1980, 1990, 2000, and 2010 using universal kriging incorporating geographic variables. County-level death rates and national life tables for each year were obtained from the U.S. Census and Centers for Disease Control and Prevention. We used log-linear and nonlinear concentration-response coefficients from previous studies to estimate changes in the numbers of deaths and in life years and life expectancy at birth, attributable to changes in PM2.5.

RESULTS

Between 1980 and 2010, population-weighted PM2.5 exposures fell by about half, and the estimated number of excess deaths declined by about a third. The States of California, Virginia, New Jersey, and Georgia had some of the largest estimated reductions in PM2.5-attributable deaths. Relative to a counterfactual population with exposures held constant at 1980 levels, we estimated that people born in 2050 would experience an ∼1-y increase in life expectancy at birth, and that there would be a cumulative gain of 4.4 million life years among adults ≥30y of age.

CONCLUSIONS

Our estimates suggest that declines in PM2.5 exposures between 1980 and 2010 have benefitted public health. https://doi.org/10.1289/EHP507.

Historical Prediction Modeling Approach for Estimating Long-Term Concentrations of PM2.5 in Cohort Studies before the 1999 Implementation of Widespread Monitoring

INTRODUCTION

Recent cohort studies have used exposure prediction models to estimate the association between long-term residential concentrations of fine particulate matter (PM2.5) and health. Because these prediction models rely on PM2.5 monitoring data, predictions for times before extensive spatial monitoring present a challenge to understanding long-term exposure effects. The U.S. Environmental Protection Agency (EPA) Federal Reference Method (FRM) network for PM2.5 was established in 1999.

OBJECTIVES

We evaluated a novel statistical approach to produce high-quality exposure predictions from 1980 through 2010 in the continental United States for epidemiological applications.

METHODS

We developed spatio-temporal prediction models using geographic predictors and annual average PM2.5 data from 1999 through 2010 from the FRM and the Interagency Monitoring of Protected Visual Environments (IMPROVE) networks. Temporal trends before 1999 were estimated by using a) extrapolation based on PM2.5 data in FRM/IMPROVE, b) PM2.5 sulfate data in the Clean Air Status and Trends Network, and c) visibility data across the Weather Bureau Army Navy network. We validated the models using PM2.5 data collected before 1999 from IMPROVE, California Air Resources Board dichotomous sampler monitoring (CARB dichot), the Children's Health Study (CHS), and the Inhalable Particulate Network (IPN).

RESULTS

In our validation using pre-1999 data, the prediction model performed well across three trend estimation approaches when validated using IMPROVE and CHS data (R2 = 0.84-0.91) with lower R2 values in early years. Model performance using CARB dichot and IPN data was worse (R2 = 0.00-0.85) most likely because of fewer monitoring sites and inconsistent sampling methods.

CONCLUSIONS

Our prediction modeling approach will allow health effects estimation associated with long-term exposures to PM2.5 over extended time periods ≤ 30 years. Citation: Kim SY, Olives C, Sheppard L, Sampson PD, Larson TV, Keller JP, Kaufman JD. 2017. Historical prediction modeling approach for estimating long-term concentrations of PM2.5 in cohort studies before the 1999 implementation of widespread monitoring. Environ Health Perspect 125:38-46; http://dx.doi.org/10.1289/EHP131.

Association between air pollution and coronary artery calcification within six metropolitan areas in the USA (the Multi-Ethnic Study of Atherosclerosis and Air Pollution): a longitudinal cohort study

BACKGROUND

Long-term exposure to fine particulate matter less than 2.5 μm in diameter (PM2.5) and traffic-related air pollutant concentrations are associated with cardiovascular risk. The disease process underlying these associations remains uncertain. We aim to assess association between long-term exposure to ambient air pollution and progression of coronary artery calcium and common carotid artery intima-media thickness.

METHODS

In this prospective 10-year cohort study, we repeatedly measured coronary artery calcium by CT in 6795 participants aged 45-84 years enrolled in the Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air) in six metropolitan areas in the USA. Repeated scans were done for nearly all participants between 2002 and 2005, for a subset of participants between 2005 and 2007, and for half of all participants between 2010 and 2012. Common carotid artery intima-media thickness was measured by ultrasound in all participants at baseline and in 2010-12 for 3459 participants. Residence-specific spatio-temporal pollution concentration models, incorporating community-specific measurements, agency monitoring data, and geographical predictors, estimated concentrations of PM2.5 and nitrogen oxides (NOX) between 1999 and 2012. The primary aim was to examine the association between both progression of coronary artery calcium and mean carotid artery intima-media thickness and long-term exposure to ambient air pollutant concentrations (PM2.5, NOX, and black carbon) between examinations and within the six metropolitan areas, adjusting for baseline age, sex, ethnicity, socioeconomic characteristics, cardiovascular risk factors, site, and CT scanner technology.

FINDINGS

In this population, coronary calcium increased on average by 24 Agatston units per year (SD 58), and intima-media thickness by 12 μm per year (10), before adjusting for risk factors or air pollutant exposures. Participant-specific pollutant concentrations averaged over the years 2000-10 ranged from 9.2-22.6 μg PM2.5/m(3) and 7.2-139.2 parts per billion (ppb) NOX. For each 5 μg PM2.5/m(3) increase, coronary calcium progressed by 4.1 Agatston units per year (95% CI 1.4-6.8) and for each 40 ppb NOX coronary calcium progressed by 4.8 Agatston units per year (0.9-8.7). Pollutant exposures were not associated with intima-media thickness change. The estimate for the effect of a 5 μg/m(3) higher long-term exposure to PM2.5 in intima-media thickness was -0.9 μm per year (95% CI -3.0 to 1.3). For 40 ppb higher NOX, the estimate was 0.2 μm per year (-1.9 to 2.4).

INTERPRETATION

Increased concentrations of PM2.5 and traffic-related air pollution within metropolitan areas, in ranges commonly encountered worldwide, are associated with progression in coronary calcification, consistent with acceleration of atherosclerosis. This study supports the case for global efforts of pollution reduction in prevention of cardiovascular diseases.

FUNDING

US Environmental Protection Agency and US National Institutes of Health.

Combining Land-Use Regression and Chemical Transport Modeling in a Spatiotemporal Geostatistical Model for Ozone and PM2.5

Assessments of long-term air pollution exposure in population studies have commonly employed land-use regression (LUR) or chemical transport modeling (CTM) techniques. Attempts to incorporate both approaches in one modeling framework are challenging. We present a novel geostatistical modeling framework, incorporating CTM predictions into a spatiotemporal LUR model with spatial smoothing to estimate spatiotemporal variability of ozone (O3) and particulate matter with diameter less than 2.5 μm (PM2.5) from 2000 to 2008 in the Los Angeles Basin. The observations include over 9 years' data from more than 20 routine monitoring sites and specific monitoring data at over 100 locations to provide more comprehensive spatial coverage of air pollutants. Our composite modeling approach outperforms separate CTM and LUR models in terms of root-mean-square error (RMSE) assessed by 10-fold cross-validation in both temporal and spatial dimensions, with larger improvement in the accuracy of predictions for O3 (RMSE [ppb] for CTM, 6.6; LUR, 4.6; composite, 3.6) than for PM2.5 (RMSE [μg/m(3)] CTM: 13.7, LUR: 3.2, composite: 3.1). Our study highlights the opportunity for future exposure assessment to make use of readily available spatiotemporal modeling methods and auxiliary gridded data that takes chemical reaction processes into account to improve the accuracy of predictions in a single spatiotemporal modeling framework.

Pretreatment with Antioxidants Augments the Acute Arterial Vasoconstriction Caused by Diesel Exhaust Inhalation

RATIONALE

Diesel exhaust inhalation, which is the model traffic-related air pollutant exposure, is associated with vascular dysfunction.

OBJECTIVES

To determine whether healthy subjects exposed to diesel exhaust exhibit acute vasoconstriction and whether this effect could be modified by the use of antioxidants or by common variants in the angiotensin II type 1 receptor (AGTR1) and other candidate genes.

METHODS

In a genotype-stratified, double-blind, four-way crossover study, 21 healthy adult subjects were exposed at rest in a randomized, balanced order to diesel exhaust (200 μg/m(3) particulate matter with an aerodynamic diameter ≤ 2.5 μm [PM2.5]) and filtered air, and to pretreatment with antioxidants (N-acetylcysteine and ascorbate) and placebo. Before and after each exposure, brachial artery diameter (BAd) was assessed using ultrasound. Changes in BAd were compared across pretreatment and exposure sessions. Gene-exposure interactions were evaluated in the AGTR1 A1166C polymorphism, on which recruitment was stratified, and other candidate genes, including TRPV1 and GSTM1.

MEASUREMENTS AND MAIN RESULTS

Compared with filtered air, exposure to diesel exhaust resulted in a significant reduction in BAd (mean, -0.09 mm, 95% confidence interval [CI], -0.01 to -0.17; P = 0.03). Pretreatment with antioxidants augmented diesel exhaust-related vasoconstriction with a mean change in BAd of -0.18 mm (95% CI, -0.28 to -0.07 mm; P = 0.001). Diesel exhaust-related vasoconstriction was primarily observed in the variant alleles of AGTR1 and TRPV1. No association was found between diesel exhaust inhalation and flow-mediated dilation.

CONCLUSIONS

We confirmed that short-term exposure to diesel exhaust in healthy subjects is associated with acute vasoconstriction in a conductance artery and found suggestive evidence of involvement of nociception and renin-angiotensin systems in this effect. Pretreatment with an antioxidant regimen increased vasoconstriction.

A class of non-linear exposure-response models suitable for health impact assessment applicable to large cohort studies of ambient air pollution

The effectiveness of regulatory actions designed to improve air quality is often assessed by predicting changes in public health resulting from their implementation. Risk of premature mortality from long-term exposure to ambient air pollution is the single most important contributor to such assessments and is estimated from observational studies generally assuming a log-linear, no-threshold association between ambient concentrations and death. There has been only limited assessment of this assumption in part because of a lack of methods to estimate the shape of the exposure-response function in very large study populations. In this paper, we propose a new class of variable coefficient risk functions capable of capturing a variety of potentially non-linear associations which are suitable for health impact assessment. We construct the class by defining transformations of concentration as the product of either a linear or log-linear function of concentration multiplied by a logistic weighting function. These risk functions can be estimated using hazard regression survival models with currently available computer software and can accommodate large population-based cohorts which are increasingly being used for this purpose. We illustrate our modeling approach with two large cohort studies of long-term concentrations of ambient air pollution and mortality: the American Cancer Society Cancer Prevention Study II (CPS II) cohort and the Canadian Census Health and Environment Cohort (CanCHEC). We then estimate the number of deaths attributable to changes in fine particulate matter concentrations over the 2000 to 2010 time period in both Canada and the USA using both linear and non-linear hazard function models.

Development of Long-term Spatiotemporal Models for Ambient Ozone in Six Metropolitan regions of the United States: The MESA Air Study

BACKGROUND

Current epidemiologic studies rely on simple ozone metrics which may not appropriately capture population ozone exposure. For understanding health effects of long-term ozone exposure in population studies, it is advantageous for exposure estimation to incorporate the complex spatiotemporal pattern of ozone concentrations at fine scales.

OBJECTIVE

To develop a geo-statistical exposure prediction model that predicts fine scale spatiotemporal variations of ambient ozone in six United States metropolitan regions.

METHODS

We developed a modeling framework that estimates temporal trends from regulatory agency and cohort-specific monitoring data from MESA Air measurement campaigns and incorporates land use regression with universal kriging using predictor variables from a large geographic database. The cohort-specific data were measured at home and community locations. The framework was applied in estimating two-week average ozone concentrations from 1999 to 2013 in models of each of the six MESA Air metropolitan regions.

RESULTS

Ozone models perform well in both spatial and temporal dimensions at the agency monitoring sites in terms of prediction accuracy. City-specific leave-one (site)-out cross-validation R² accounting for temporal and spatial variability ranged from 0.65 to 0.88 in the six regions. For predictions at the home sites, the R² is between 0.60 and 0.91 for cross-validation that left out 10% of home sites in turn. The predicted ozone concentrations vary substantially over space and time in all the metropolitan regions.

CONCLUSION

Using the available data, our spatiotemporal models are able to accurately predict long-term ozone concentrations at fine spatial scales in multiple regions. The model predictions will allow for investigation of the long-term health effects of ambient ozone concentrations in future epidemiological studies.

Risk Factors for Long-Term Coronary Artery Calcium Progression in the Multi-Ethnic Study of Atherosclerosis

Background

Coronary artery calcium (CAC) detected by noncontrast cardiac computed tomography scanning is a measure of coronary atherosclerosis burden. Increasing CAC levels have been strongly associated with increased coronary events. Prior studies of cardiovascular disease risk factors and CAC progression have been limited by short follow-up or restricted to patients with advanced disease.

Methods and Results

We examined cardiovascular disease risk factors and CAC progression in a prospective multiethnic cohort study. CAC was measured 1 to 4 times (mean 2.5 scans) over 10 years in 6810 adults without preexisting cardiovascular disease. Mean CAC progression was 23.9 Agatston units/year. An innovative application of mixed-effects models investigated associations between cardiovascular disease risk factors and CAC progression. This approach adjusted for time-varying factors, was flexible with respect to follow-up time and number of observations per participant, and allowed simultaneous control of factors associated with both baseline CAC and CAC progression. Models included age, sex, study site, scanner type, and race/ethnicity. Associations were observed between CAC progression and age (14.2 Agatston units/year per 10 years [95% CI 13.0 to 15.5]), male sex (17.8 Agatston units/year [95% CI 15.3 to 20.3]), hypertension (13.8 Agatston units/year [95% CI 11.2 to 16.5]), diabetes (31.3 Agatston units/year [95% CI 27.4 to 35.3]), and other factors.

Conclusions

CAC progression analyzed over 10 years of follow-up, with a novel analytical approach, demonstrated strong relationships with risk factors for incident cardiovascular events. Longitudinal CAC progression analyzed in this framework can be used to evaluate novel cardiovascular risk factors.

A unified spatiotemporal modeling approach for predicting concentrations of multiple air pollutants in the multi-ethnic study of atherosclerosis and air pollution

BACKGROUND

Cohort studies of the relationship between air pollution exposure and chronic health effects require predictions of exposure over long periods of time.

OBJECTIVES

We developed a unified modeling approach for predicting fine particulate matter, nitrogen dioxide, oxides of nitrogen, and black carbon (as measured by light absorption coefficient) in six U.S. metropolitan regions from 1999 through early 2012 as part of the Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air).

METHODS

We obtained monitoring data from regulatory networks and supplemented those data with study-specific measurements collected from MESA Air community locations and participants' homes. In each region, we applied a spatiotemporal model that included a long-term spatial mean, time trends with spatially varying coefficients, and a spatiotemporal residual. The mean structure was derived from a large set of geographic covariates that was reduced using partial least-squares regression. We estimated time trends from observed time series and used spatial smoothing methods to borrow strength between observations.

RESULTS

Prediction accuracy was high for most models, with cross-validation R2 (R2CV) > 0.80 at regulatory and fixed sites for most regions and pollutants. At home sites, overall R2CV ranged from 0.45 to 0.92, and temporally adjusted R2CV ranged from 0.23 to 0.92.

CONCLUSIONS

This novel spatiotemporal modeling approach provides accurate fine-scale predictions in multiple regions for four pollutants. We have generated participant-specific predictions for MESA Air to investigate health effects of long-term air pollution exposures. These successes highlight modeling advances that can be adopted more widely in modern cohort studies.

Reduced-Rank Spatio-Temporal Modeling of Air Pollution Concentrations in the Multi-Ethnic Study of Atherosclerosis and Air Pollution

There is growing evidence in the epidemiologic literature of the relationship between air pollution and adverse health outcomes. Prediction of individual air pollution exposure in the Environmental Protection Agency (EPA) funded Multi-Ethnic Study of Atheroscelerosis and Air Pollution (MESA Air) study relies on a flexible spatio-temporal prediction model that integrates land-use regression with kriging to account for spatial dependence in pollutant concentrations. Temporal variability is captured using temporal trends estimated via modified singular value decomposition and temporally varying spatial residuals. This model utilizes monitoring data from existing regulatory networks and supplementary MESA Air monitoring data to predict concentrations for individual cohort members. In general, spatio-temporal models are limited in their efficacy for large data sets due to computational intractability. We develop reduced-rank versions of the MESA Air spatio-temporal model. To do so, we apply low-rank kriging to account for spatial variation in the mean process and discuss the limitations of this approach. As an alternative, we represent spatial variation using thin plate regression splines. We compare the performance of the outlined models using EPA and MESA Air monitoring data for predicting concentrations of oxides of nitrogen (NO _x )-a pollutant of primary interest in MESA Air-in the Los Angeles metropolitan area via cross-validated R². Our findings suggest that use of reduced-rank models can improve computational efficiency in certain cases. Low-rank kriging and thin plate regression splines were competitive across the formulations considered, although TPRS appeared to be more robust in some settings.

An integrated risk function for estimating the global burden of disease attributable to ambient fine particulate matter exposure

BACKGROUND

Estimating the burden of disease attributable to long-term exposure to fine particulate matter (PM2.5) in ambient air requires knowledge of both the shape and magnitude of the relative risk (RR) function. However, adequate direct evidence to identify the shape of the mortality RR functions at the high ambient concentrations observed in many places in the world is lacking.

OBJECTIVE

We developed RR functions over the entire global exposure range for causes of mortality in adults: ischemic heart disease (IHD), cerebrovascular disease (stroke), chronic obstructive pulmonary disease (COPD), and lung cancer (LC). We also developed RR functions for the incidence of acute lower respiratory infection (ALRI) that can be used to estimate mortality and lost-years of healthy life in children < 5 years of age.

METHODS

We fit an integrated exposure-response (IER) model by integrating available RR information from studies of ambient air pollution (AAP), second hand tobacco smoke, household solid cooking fuel, and active smoking (AS). AS exposures were converted to estimated annual PM2.5 exposure equivalents using inhaled doses of particle mass. We derived population attributable fractions (PAFs) for every country based on estimated worldwide ambient PM2.5 concentrations.

RESULTS

The IER model was a superior predictor of RR compared with seven other forms previously used in burden assessments. The percent PAF attributable to AAP exposure varied among countries from 2 to 41 for IHD, 1 to 43 for stroke, < 1 to 21 for COPD, < 1 to 25 for LC, and < 1 to 38 for ALRI.

CONCLUSIONS

We developed a fine particulate mass-based RR model that covered the global range of exposure by integrating RR information from different combustion types that generate emissions of particulate matter. The model can be updated as new RR information becomes available.

Estimation after classification using lot quality assurance sampling: corrections for curtailed sampling with application to evaluating polio vaccination campaigns

OBJECTIVES

To assess the bias incurred when curtailment of Lot Quality Assurance Sampling (LQAS) is ignored, to present unbiased estimators, to consider the impact of cluster sampling by simulation and to apply our method to published polio immunization data from Nigeria.

METHODS

We present estimators of coverage when using two kinds of curtailed LQAS strategies: semicurtailed and curtailed. We study the proposed estimators with independent and clustered data using three field-tested LQAS designs for assessing polio vaccination coverage, with samples of size 60 and decision rules of 9, 21 and 33, and compare them to biased maximum likelihood estimators. Lastly, we present estimates of polio vaccination coverage from previously published data in 20 local government authorities (LGAs) from five Nigerian states.

RESULTS

Simulations illustrate substantial bias if one ignores the curtailed sampling design. Proposed estimators show no bias. Clustering does not affect the bias of these estimators. Across simulations, standard errors show signs of inflation as clustering increases. Neither sampling strategy nor LQAS design influences estimates of polio vaccination coverage in 20 Nigerian LGAs. When coverage is low, semicurtailed LQAS strategies considerably reduces the sample size required to make a decision. Curtailed LQAS designs further reduce the sample size when coverage is high.

CONCLUSIONS

Results presented dispel the misconception that curtailed LQAS data are unsuitable for estimation. These findings augment the utility of LQAS as a tool for monitoring vaccination efforts by demonstrating that unbiased estimation using curtailed designs is not only possible but these designs also reduce the sample size.

The effect of clustering on lot quality assurance sampling: a probabilistic model to calculate sample sizes for quality assessments

BACKGROUND

Traditional Lot Quality Assurance Sampling (LQAS) designs assume observations are collected using simple random sampling. Alternatively, randomly sampling clusters of observations and then individuals within clusters reduces costs but decreases the precision of the classifications. In this paper, we develop a general framework for designing the cluster(C)-LQAS system and illustrate the method with the design of data quality assessments for the community health worker program in Rwanda.

RESULTS

To determine sample size and decision rules for C-LQAS, we use the beta-binomial distribution to account for inflated risk of errors introduced by sampling clusters at the first stage. We present general theory and code for sample size calculations.The C-LQAS sample sizes provided in this paper constrain misclassification risks below user-specified limits. Multiple C-LQAS systems meet the specified risk requirements, but numerous considerations, including per-cluster versus per-individual sampling costs, help identify optimal systems for distinct applications.

CONCLUSIONS

We show the utility of C-LQAS for data quality assessments, but the method generalizes to numerous applications. This paper provides the necessary technical detail and supplemental code to support the design of C-LQAS for specific programs.

Comparing the performance of cluster random sampling and integrated threshold mapping for targeting trachoma control, using computer simulation

BACKGROUND

Implementation of trachoma control strategies requires reliable district-level estimates of trachomatous inflammation-follicular (TF), generally collected using the recommended gold-standard cluster randomized surveys (CRS). Integrated Threshold Mapping (ITM) has been proposed as an integrated and cost-effective means of rapidly surveying trachoma in order to classify districts according to treatment thresholds. ITM differs from CRS in a number of important ways, including the use of a school-based sampling platform for children aged 1-9 and a different age distribution of participants. This study uses computerised sampling simulations to compare the performance of these survey designs and evaluate the impact of varying key parameters.

METHODOLOGY/PRINCIPAL FINDINGS

Realistic pseudo gold standard data for 100 districts were generated that maintained the relative risk of disease between important sub-groups and incorporated empirical estimates of disease clustering at the household, village and district level. To simulate the different sampling approaches, 20 clusters were selected from each district, with individuals sampled according to the protocol for ITM and CRS. Results showed that ITM generally under-estimated the true prevalence of TF over a range of epidemiological settings and introduced more district misclassification according to treatment thresholds than did CRS. However, the extent of underestimation and resulting misclassification was found to be dependent on three main factors: (i) the district prevalence of TF; (ii) the relative risk of TF between enrolled and non-enrolled children within clusters; and (iii) the enrollment rate in schools.

CONCLUSIONS/SIGNIFICANCE

Although in some contexts the two methodologies may be equivalent, ITM can introduce a bias-dependent shift as prevalence of TF increases, resulting in a greater risk of misclassification around treatment thresholds. In addition to strengthening the evidence base around choice of trachoma survey methodologies, this study illustrates the use of a simulated approach in addressing operational research questions for trachoma but also other NTDs.

The state of US health, 1990-2010: burden of diseases, injuries, and risk factors

IMPORTANCE

Understanding the major health problems in the United States and how they are changing over time is critical for informing national health policy.

OBJECTIVES

To measure the burden of diseases, injuries, and leading risk factors in the United States from 1990 to 2010 and to compare these measurements with those of the 34 countries in the Organisation for Economic Co-operation and Development (OECD) countries.

DESIGN

We used the systematic analysis of descriptive epidemiology of 291 diseases and injuries, 1160 sequelae of these diseases and injuries, and 67 risk factors or clusters of risk factors from 1990 to 2010 for 187 countries developed for the Global Burden of Disease 2010 Study to describe the health status of the United States and to compare US health outcomes with those of 34 OECD countries. Years of life lost due to premature mortality (YLLs) were computed by multiplying the number of deaths at each age by a reference life expectancy at that age. Years lived with disability (YLDs) were calculated by multiplying prevalence (based on systematic reviews) by the disability weight (based on population-based surveys) for each sequela; disability in this study refers to any short- or long-term loss of health. Disability-adjusted life-years (DALYs) were estimated as the sum of YLDs and YLLs. Deaths and DALYs related to risk factors were based on systematic reviews and meta-analyses of exposure data and relative risks for risk-outcome pairs. Healthy life expectancy (HALE) was used to summarize overall population health, accounting for both length of life and levels of ill health experienced at different ages.

RESULTS

US life expectancy for both sexes combined increased from 75.2 years in 1990 to 78.2 years in 2010; during the same period, HALE increased from 65.8 years to 68.1 years. The diseases and injuries with the largest number of YLLs in 2010 were ischemic heart disease, lung cancer, stroke, chronic obstructive pulmonary disease, and road injury. Age-standardized YLL rates increased for Alzheimer disease, drug use disorders, chronic kidney disease, kidney cancer, and falls. The diseases with the largest number of YLDs in 2010 were low back pain, major depressive disorder, other musculoskeletal disorders, neck pain, and anxiety disorders. As the US population has aged, YLDs have comprised a larger share of DALYs than have YLLs. The leading risk factors related to DALYs were dietary risks, tobacco smoking, high body mass index, high blood pressure, high fasting plasma glucose, physical inactivity, and alcohol use. Among 34 OECD countries between 1990 and 2010, the US rank for the age-standardized death rate changed from 18th to 27th, for the age-standardized YLL rate from 23rd to 28th, for the age-standardized YLD rate from 5th to 6th, for life expectancy at birth from 20th to 27th, and for HALE from 14th to 26th.

CONCLUSIONS AND RELEVANCE

From 1990 to 2010, the United States made substantial progress in improving health. Life expectancy at birth and HALE increased, all-cause death rates at all ages decreased, and age-specific rates of years lived with disability remained stable. However, morbidity and chronic disability now account for nearly half of the US health burden, and improvements in population health in the United States have not kept pace with advances in population health in other wealthy nations.

Choosing a design to fit the situation: how to improve specificity and positive predictive values using Bayesian lot quality assurance sampling

BACKGROUND

Lot Quality Assurance Sampling (LQAS) is a provably useful tool for monitoring health programmes. Although LQAS ensures acceptable Producer and Consumer risks, the literature alleges that the method suffers from poor specificity and positive predictive values (PPVs). We suggest that poor LQAS performance is due, in part, to variation in the true underlying distribution. However, until now the role of the underlying distribution in expected performance has not been adequately examined.

METHODS

We present Bayesian-LQAS (B-LQAS), an approach to incorporating prior information into the choice of the LQAS sample size and decision rule, and explore its properties through a numerical study. Additionally, we analyse vaccination coverage data from UNICEF's State of the World's Children in 1968-1989 and 2008 to exemplify the performance of LQAS and B-LQAS.

RESULTS

Results of our numerical study show that the choice of LQAS sample size and decision rule is sensitive to the distribution of prior information, as well as to individual beliefs about the importance of correct classification. Application of the B-LQAS approach to the UNICEF data improves specificity and PPV in both time periods (1968-1989 and 2008) with minimal reductions in sensitivity and negative predictive value.

CONCLUSIONS

LQAS is shown to be a robust tool that is not necessarily prone to poor specificity and PPV as previously alleged. In situations where prior or historical data are available, B-LQAS can lead to improvements in expected performance.

Years lived with disability (YLDs) for 1160 sequelae of 289 diseases and injuries 1990-2010: a systematic analysis for the Global Burden of Disease Study 2010

BACKGROUND

Non-fatal health outcomes from diseases and injuries are a crucial consideration in the promotion and monitoring of individual and population health. The Global Burden of Disease (GBD) studies done in 1990 and 2000 have been the only studies to quantify non-fatal health outcomes across an exhaustive set of disorders at the global and regional level. Neither effort quantified uncertainty in prevalence or years lived with disability (YLDs).

METHODS

Of the 291 diseases and injuries in the GBD cause list, 289 cause disability. For 1160 sequelae of the 289 diseases and injuries, we undertook a systematic analysis of prevalence, incidence, remission, duration, and excess mortality. Sources included published studies, case notification, population-based cancer registries, other disease registries, antenatal clinic serosurveillance, hospital discharge data, ambulatory care data, household surveys, other surveys, and cohort studies. For most sequelae, we used a Bayesian meta-regression method, DisMod-MR, designed to address key limitations in descriptive epidemiological data, including missing data, inconsistency, and large methodological variation between data sources. For some disorders, we used natural history models, geospatial models, back-calculation models (models calculating incidence from population mortality rates and case fatality), or registration completeness models (models adjusting for incomplete registration with health-system access and other covariates). Disability weights for 220 unique health states were used to capture the severity of health loss. YLDs by cause at age, sex, country, and year levels were adjusted for comorbidity with simulation methods. We included uncertainty estimates at all stages of the analysis.

FINDINGS

Global prevalence for all ages combined in 2010 across the 1160 sequelae ranged from fewer than one case per 1 million people to 350,000 cases per 1 million people. Prevalence and severity of health loss were weakly correlated (correlation coefficient -0·37). In 2010, there were 777 million YLDs from all causes, up from 583 million in 1990. The main contributors to global YLDs were mental and behavioural disorders, musculoskeletal disorders, and diabetes or endocrine diseases. The leading specific causes of YLDs were much the same in 2010 as they were in 1990: low back pain, major depressive disorder, iron-deficiency anaemia, neck pain, chronic obstructive pulmonary disease, anxiety disorders, migraine, diabetes, and falls. Age-specific prevalence of YLDs increased with age in all regions and has decreased slightly from 1990 to 2010. Regional patterns of the leading causes of YLDs were more similar compared with years of life lost due to premature mortality. Neglected tropical diseases, HIV/AIDS, tuberculosis, malaria, and anaemia were important causes of YLDs in sub-Saharan Africa.

INTERPRETATION

Rates of YLDs per 100,000 people have remained largely constant over time but rise steadily with age. Population growth and ageing have increased YLD numbers and crude rates over the past two decades. Prevalences of the most common causes of YLDs, such as mental and behavioural disorders and musculoskeletal disorders, have not decreased. Health systems will need to address the needs of the rising numbers of individuals with a range of disorders that largely cause disability but not mortality. Quantification of the burden of non-fatal health outcomes will be crucial to understand how well health systems are responding to these challenges. Effective and affordable strategies to deal with this rising burden are an urgent priority for health systems in most parts of the world.

FUNDING

Bill & Melinda Gates Foundation.

Disability-adjusted life years (DALYs) for 291 diseases and injuries in 21 regions, 1990-2010: a systematic analysis for the Global Burden of Disease Study 2010

BACKGROUND

Measuring disease and injury burden in populations requires a composite metric that captures both premature mortality and the prevalence and severity of ill-health. The 1990 Global Burden of Disease study proposed disability-adjusted life years (DALYs) to measure disease burden. No comprehensive update of disease burden worldwide incorporating a systematic reassessment of disease and injury-specific epidemiology has been done since the 1990 study. We aimed to calculate disease burden worldwide and for 21 regions for 1990, 2005, and 2010 with methods to enable meaningful comparisons over time.

METHODS

We calculated DALYs as the sum of years of life lost (YLLs) and years lived with disability (YLDs). DALYs were calculated for 291 causes, 20 age groups, both sexes, and for 187 countries, and aggregated to regional and global estimates of disease burden for three points in time with strictly comparable definitions and methods. YLLs were calculated from age-sex-country-time-specific estimates of mortality by cause, with death by standardised lost life expectancy at each age. YLDs were calculated as prevalence of 1160 disabling sequelae, by age, sex, and cause, and weighted by new disability weights for each health state. Neither YLLs nor YLDs were age-weighted or discounted. Uncertainty around cause-specific DALYs was calculated incorporating uncertainty in levels of all-cause mortality, cause-specific mortality, prevalence, and disability weights.

FINDINGS

Global DALYs remained stable from 1990 (2·503 billion) to 2010 (2·490 billion). Crude DALYs per 1000 decreased by 23% (472 per 1000 to 361 per 1000). An important shift has occurred in DALY composition with the contribution of deaths and disability among children (younger than 5 years of age) declining from 41% of global DALYs in 1990 to 25% in 2010. YLLs typically account for about half of disease burden in more developed regions (high-income Asia Pacific, western Europe, high-income North America, and Australasia), rising to over 80% of DALYs in sub-Saharan Africa. In 1990, 47% of DALYs worldwide were from communicable, maternal, neonatal, and nutritional disorders, 43% from non-communicable diseases, and 10% from injuries. By 2010, this had shifted to 35%, 54%, and 11%, respectively. Ischaemic heart disease was the leading cause of DALYs worldwide in 2010 (up from fourth rank in 1990, increasing by 29%), followed by lower respiratory infections (top rank in 1990; 44% decline in DALYs), stroke (fifth in 1990; 19% increase), diarrhoeal diseases (second in 1990; 51% decrease), and HIV/AIDS (33rd in 1990; 351% increase). Major depressive disorder increased from 15th to 11th rank (37% increase) and road injury from 12th to 10th rank (34% increase). Substantial heterogeneity exists in rankings of leading causes of disease burden among regions.

INTERPRETATION

Global disease burden has continued to shift away from communicable to non-communicable diseases and from premature death to years lived with disability. In sub-Saharan Africa, however, many communicable, maternal, neonatal, and nutritional disorders remain the dominant causes of disease burden. The rising burden from mental and behavioural disorders, musculoskeletal disorders, and diabetes will impose new challenges on health systems. Regional heterogeneity highlights the importance of understanding local burden of disease and setting goals and targets for the post-2015 agenda taking such patterns into account. Because of improved definitions, methods, and data, these results for 1990 and 2010 supersede all previously published Global Burden of Disease results.

FUNDING

Bill & Melinda Gates Foundation.

A comparative risk assessment of burden of disease and injury attributable to 67 risk factors and risk factor clusters in 21 regions, 1990-2010: a systematic analysis for the Global Burden of Disease Study 2010

BACKGROUND

Quantification of the disease burden caused by different risks informs prevention by providing an account of health loss different to that provided by a disease-by-disease analysis. No complete revision of global disease burden caused by risk factors has been done since a comparative risk assessment in 2000, and no previous analysis has assessed changes in burden attributable to risk factors over time.

METHODS

We estimated deaths and disability-adjusted life years (DALYs; sum of years lived with disability [YLD] and years of life lost [YLL]) attributable to the independent effects of 67 risk factors and clusters of risk factors for 21 regions in 1990 and 2010. We estimated exposure distributions for each year, region, sex, and age group, and relative risks per unit of exposure by systematically reviewing and synthesising published and unpublished data. We used these estimates, together with estimates of cause-specific deaths and DALYs from the Global Burden of Disease Study 2010, to calculate the burden attributable to each risk factor exposure compared with the theoretical-minimum-risk exposure. We incorporated uncertainty in disease burden, relative risks, and exposures into our estimates of attributable burden.

FINDINGS

In 2010, the three leading risk factors for global disease burden were high blood pressure (7·0% [95% uncertainty interval 6·2-7·7] of global DALYs), tobacco smoking including second-hand smoke (6·3% [5·5-7·0]), and alcohol use (5·5% [5·0-5·9]). In 1990, the leading risks were childhood underweight (7·9% [6·8-9·4]), household air pollution from solid fuels (HAP; 7·0% [5·6-8·3]), and tobacco smoking including second-hand smoke (6·1% [5·4-6·8]). Dietary risk factors and physical inactivity collectively accounted for 10·0% (95% UI 9·2-10·8) of global DALYs in 2010, with the most prominent dietary risks being diets low in fruits and those high in sodium. Several risks that primarily affect childhood communicable diseases, including unimproved water and sanitation and childhood micronutrient deficiencies, fell in rank between 1990 and 2010, with unimproved water and sanitation accounting for 0·9% (0·4-1·6) of global DALYs in 2010. However, in most of sub-Saharan Africa childhood underweight, HAP, and non-exclusive and discontinued breastfeeding were the leading risks in 2010, while HAP was the leading risk in south Asia. The leading risk factor in Eastern Europe, most of Latin America, and southern sub-Saharan Africa in 2010 was alcohol use; in most of Asia, North Africa and Middle East, and central Europe it was high blood pressure. Despite declines, tobacco smoking including second-hand smoke remained the leading risk in high-income north America and western Europe. High body-mass index has increased globally and it is the leading risk in Australasia and southern Latin America, and also ranks high in other high-income regions, North Africa and Middle East, and Oceania.

INTERPRETATION

Worldwide, the contribution of different risk factors to disease burden has changed substantially, with a shift away from risks for communicable diseases in children towards those for non-communicable diseases in adults. These changes are related to the ageing population, decreased mortality among children younger than 5 years, changes in cause-of-death composition, and changes in risk factor exposures. New evidence has led to changes in the magnitude of key risks including unimproved water and sanitation, vitamin A and zinc deficiencies, and ambient particulate matter pollution. The extent to which the epidemiological shift has occurred and what the leading risks currently are varies greatly across regions. In much of sub-Saharan Africa, the leading risks are still those associated with poverty and those that affect children.

FUNDING

Bill & Melinda Gates Foundation.

Multiple category-lot quality assurance sampling: a new classification system with application to schistosomiasis control

Background

Originally a binary classifier, Lot Quality Assurance Sampling (LQAS) has proven to be a useful tool for classification of the prevalence of Schistosoma mansoni into multiple categories (≤10%, >10 and <50%, ≥50%), and semi-curtailed sampling has been shown to effectively reduce the number of observations needed to reach a decision. To date the statistical underpinnings for Multiple Category-LQAS (MC-LQAS) have not received full treatment. We explore the analytical properties of MC-LQAS, and validate its use for the classification of S. mansoni prevalence in multiple settings in East Africa.

Methodology

We outline MC-LQAS design principles and formulae for operating characteristic curves. In addition, we derive the average sample number for MC-LQAS when utilizing semi-curtailed sampling and introduce curtailed sampling in this setting. We also assess the performance of MC-LQAS designs with maximum sample sizes of n = 15 and n = 25 via a weighted kappa-statistic using S. mansoni data collected in 388 schools from four studies in East Africa.

Principle Findings

Overall performance of MC-LQAS classification was high (kappa-statistic of 0.87). In three of the studies, the kappa-statistic for a design with n = 15 was greater than 0.75. In the fourth study, where these designs performed poorly (kappa-statistic less than 0.50), the majority of observations fell in regions where potential error is known to be high. Employment of semi-curtailed and curtailed sampling further reduced the sample size by as many as 0.5 and 3.5 observations per school, respectively, without increasing classification error.

Conclusion/Significance

This work provides the needed analytics to understand the properties of MC-LQAS for assessing the prevalance of S. mansoni and shows that in most settings a sample size of 15 children provides a reliable classification of schools.

The control of schistosomiasis calls for rapid and reliable classification tools. This study evaluates the performance of one such tool, Lot Quality Assurance Sampling (LQAS) for assessing the prevalence of S. mansoni in African schoolchildren. We outline the design considerations and introduce novel sequential sampling plans for Multiple Category-LQAS. We use data from 388 schools in Uganda, Kenya, and Tanzania to assess the performance of LQAS as a tool for classification of S. mansoni infection into one of three classes: ≤10% >10 and <50%, ≥50%. Our findings suggest that an LQAS-based multiple classification system performs as well as the World Health Organization recommended methods at a fraction of the sampling effort. Our work validates LQAS as a rapid assessment tool and extends it to allow investigators to apply the method to control other infectious diseases.

Bayes-LQAS: classifying the prevalence of global acute malnutrition

Lot Quality Assurance Sampling (LQAS) applications in health have generally relied on frequentist interpretations for statistical validity. Yet health professionals often seek statements about the probability distribution of unknown parameters to answer questions of interest. The frequentist paradigm does not pretend to yield such information, although a Bayesian formulation might. This is the source of an error made in a recent paper published in this journal. Many applications lend themselves to a Bayesian treatment, and would benefit from such considerations in their design. We discuss Bayes-LQAS (B-LQAS), which allows for incorporation of prior information into the LQAS classification procedure, and thus shows how to correct the aforementioned error. Further, we pay special attention to the formulation of Bayes Operating Characteristic Curves and the use of prior information to improve survey designs. As a motivating example, we discuss the classification of Global Acute Malnutrition prevalence and draw parallels between the Bayes and classical classifications schemes. We also illustrate the impact of informative and non-informative priors on the survey design. Results indicate that using a Bayesian approach allows the incorporation of expert information and/or historical data and is thus potentially a valuable tool for making accurate and precise classifications.

Cluster designs to assess the prevalence of acute malnutrition by lot quality assurance sampling: a validation study by computer simulation

Traditional lot quality assurance sampling (LQAS) methods require simple random sampling to guarantee valid results. However, cluster sampling has been proposed to reduce the number of random starting points. This study uses simulations to examine the classification error of two such designs, a 67x3 (67 clusters of three observations) and a 33x6 (33 clusters of six observations) sampling scheme to assess the prevalence of global acute malnutrition (GAM). Further, we explore the use of a 67x3 sequential sampling scheme for LQAS classification of GAM prevalence. Results indicate that, for independent clusters with moderate intracluster correlation for the GAM outcome, the three sampling designs maintain approximate validity for LQAS analysis. Sequential sampling can substantially reduce the average sample size that is required for data collection. The presence of intercluster correlation can impact dramatically the classification error that is associated with LQAS analysis.