Model Choice in Time Series Studies of Air Pollution and Mortality

Abstract

Summary.Multicity time series studies of particulate matter and mortality and morbidity have provided evidence that daily variation in air pollution levels is associated with daily variation in mortality counts. These findings served as key epidemiological evidence for the recent review of the US national ambient air quality standards for particulate matter. As a result, methodological issues concerning time series analysis of the relationship between air pollution and health have attracted the attention of the scientific community and critics have raised concerns about the adequacy of current model formulations. Time series data on pollution and mortality are generally analysed by using log-linear, Poisson regression models for overdispersed counts with the daily number of deaths as outcome, the (possibly lagged) daily level of pollution as a linear predictor and smooth functions of weather variables and calendar time used to adjust for timevarying confounders. Investigators around the world have used different approaches to adjust for confounding, making it difficult to compare results across studies. To date, the statistical properties of these different approaches have not been comprehensively compared. To address these issues, we quantify and characterize model uncertainty and model choice in adjusting for seasonal and long-term trends in time series models of air pollution and mortality. First, we conduct a simulation study to compare and describe the properties of statistical methods that are commonly used for confounding adjustment. We generate data under several confounding scenarios and systematically compare the performance of the various methods with respect to the mean-squared error of the estimated air pollution coefficient. We find that the bias in the estimates generally decreases with more aggressive smoothing and that model selection methods which optimize prediction may not be suitable for obtaining an estimate with small bias. Second, we apply and compare the modelling approaches with the National Morbidity, Mortality, and Air Pollution Study database which comprises daily time series of several pollutants, weather variables and mortality counts covering the period 1987-2000 for the largest 100 cities in the USA. When applying these approaches to adjusting for seasonal and long-term trends we find that the Study's estimates for the national average effect of PM 10 at lag 1 on mortality vary over approximately a twofold range, with 95% posterior intervals always excluding zero risk.

Effects of temperature on mortality in Hong Kong: a time series analysis

“…The median value of temperature was defined as the baseline temperature (centering value) for calculating the relative risks. Akaike information criterion for quasi-Poisson (AIC-Q) models was used to choose the df (knots) for temperature and lag (Peng et al 2006;Gasparrini et al 2010). It was found that using 5 df for temperature and 4 df for lag produced the best model fit in the present study.…”

Bootstrap model averaging in time series studies of particulate matter air pollution and mortality

“…In many community time series studies on the effect of PM on mortality, an additive Poisson log-linear model is fit to the time series of observed mortality (Kelsall et al, 1997;Daniels et al, 2000;Kim et al, 2003;Peng et al, 2004a). Under this model, the daily mortality counts are modeled as independent Poisson random variables with a time varying mean m t on day t given by…”

Air Pollution Associated with Sumatran Forest Fires and Mortality on the Malay Peninsula

“…Population sensitivity to PM, age distribution, and especially different particle components may strongly affect the exposure-response relationships. In addition, inconsistent with several time series studies [17][18][19], we used the non-parameter Spearman correlation analysis (in SPSS software) for our data analysis. Most prior air pollution time-series studies have been based on the GAM model in S-PLUS.…”