Methods
Poisson regression models incorporating natural spline smoothing functions were used to adjust for seasonality and other time-varying covariates that might confound the association between air pollution and mortality. Effect estimates were determined for each city and then for the cities combined using a random effects method.