Click on the links below to access all the ArticlePlus for this article.
Please note that ArticlePlus files may launch a viewer application outside of your web browser.
A review of the literature on the chronic health effects of ambient air pollution suggests that studies using the American Cancer Society (ACS) cohort to assess the relation between particulate air pollution and mortality rank among the most influential and widely cited. The original study1 (a reanalysis that introduced new random-effects methods and spatial analytic techniques2,3) and more recent studies with longer follow up and improved exposure data have all demonstrated air pollution effects on all-cause and cause-specific mortality.4,5 As a result of this robust association and a lack of other studies on the long-term effects, the ACS studies together with the Six-Cities study6 have been important for government regulatory interventions such as the U.S. Environmental Protection Agency's National Air Quality Standard for Fine Particles. The ACS studies have also been used by the World Health Organization as a basis for estimating the burden of mortality attributable to air pollution.7
The assessment of air pollution exposure using only community average concentrations likely underestimates the health burden attributable to elevated concentrations in the vicinity of sources.8,9 Health effects may be larger around sources, and these effects are diminished when using average concentrations for the entire community. Previous ACS studies have relied on between-community exposure contrasts at the scale of a metropolitan area giving all residents of a city the same exposure concentrations. Exposure to air pollution, however, may vary spatially within a city,10–14 and these variations may follow social gradients that influence susceptibility to environmental exposures.3 Residents of poorer neighborhoods may live closer to point sources of industrial pollution or roadways with higher traffic density.15 This exposure misclassification along social gradients may explain the finding of effect modification by educational status in earlier ACS studies.2,16 The spatial correspondence between high exposure and potentially susceptible populations within cities may further bias estimates that rely on central monitors to proxy exposure over wide areas. Theoretically, classic exposure measurement error induced by central monitors may also bias results toward the null.17
Given the potential of the metropolitan scale to bias health effect estimates, we have assessed the association between air pollution and mortality at the within-community or intraurban scale. We sought an urban location with sufficient geographic scope, air pollution data, and enough ACS subjects to test the association. Los Angeles (LA), California, met these selection criteria. The region has high pollution levels, large intraurban gradients in exposure over a wide geographic area, and strong public awareness that air pollution has serious public health consequences.18
We extracted health data from the ACS Cancer Prevention II survey for metropolitan LA at the zip code-area scale (zip codes are used for U.S. mail delivery; average population per zip code in LA is approximately 35,000, with an average area of approximately 22.5 km2). We constructed distribution-weighted centroids using spatial boundary files based on 1980 and 1990 definitions. We were able to assign exposure to 267 zip code areas with a total of 22,905 subjects (5856 deaths based on follow up to 2000). Some subjects reported only postal box addresses and were therefore excluded. These subjects had been enrolled in 1982 along with over one million others as part of the ACS II survey. Similar to earlier ACS analyses, availability of air pollution data and other relevant information led to the subset of study subjects to be used in the health effects research. Although the ACS cohort is not representative of the general population, the cohort allows for internally valid comparisons within large samples of the American population. This study was approved by the Ethics Board of the Ottawa General Hospital, Canada. Subjects had given informed consent at enrollment into the study.
Control for Confounding
We used 44 individual confounders identified in earlier ACS studies of air pollution health effects.4 These variables include lifestyle, dietary, demographic, occupational, and educational factors that may confound the air pollution–mortality association. We had more than 10 variables that measure aspects of smoking. Sensitivity analyses revealed that removal of individual variables had little influence on the estimated pollution coefficients; therefore, to promote comparability with results from earlier studies, we report the results with this standard set of 44 variables.
We also assembled 8 ecologic variables for the zip code areas to control for “contextual” neighborhood confounding. “Contextual” effects occur when individual differences in health outcome are associated with the grouped variables that represent the social, economic, and environmental settings where the individuals live, work, or spend time (eg, poverty or crime rate in a neighborhood).19–22 These contextual effects often operate independently from (or interactively with) the individual-level variables such as smoking. The ecologic variables used represented constructs identified as important in the population health literature and previously tested as potential confounders with the ACS dataset at the metropolitan scale.23,24 These include income, income inequality, education, population size, racial composition (black, white, Hispanic), and unemployment.3 A new variable measuring potential exposure misclassification by the proportion having air conditioning was also tested. Similar variables have been in a metaanalysis of acute effects,25 on the premise that air-conditioned houses are more tightly sealed and have lower penetration of particles indoors. A recent study of personal exposures in LA reported large reductions in penetration of particles for air-conditioned houses.26 This variable adds partial control for the impact of air conditioning, which may relate both to health outcomes (through prevention of heat stroke) and to air pollution (because high air pollution concentrations and lower proportions of air conditioning are related in our study area). We thus expected the proportion of air conditioning in the zip code area to correlate with lower PM exposures and effects. We also computed principal components of all 8 variables to provide maximal control for confounding while avoiding multicollinearity among the ecologic variables.27,28
To derive exposure assessments, we interpolated PM2.5 data from 23 state and local district monitoring stations in the LA basin for the year 2000 using 5 interpolation methods: bicubic splines, 2 ordinary kriging models, universal kriging with a quadratic drift, and a radial basis function multiquadric interpolator. We emphasized kriging interpolation because this stochastic method produces the best linear unbiased estimate of the pollution surface.29 After crossvalidation, we used a combination of universal kriging and multiquadric models. This approach takes advantage of the local detail in the multiquadric surface and the ability to handle trends in the universal surface. We averaged estimated surfaces based on 25-m grid cells. We conducted sensitivity analysis using only the universal estimate and found the results to be similar; therefore, only the findings from the combined model are reported. Sensitivity analyses were also implemented with the kriging variance. Exposure assignments were downweighted with larger errors in exposure estimates in these analyses (ie, weight equal to the inverse of the standard error in the universal kriging estimate).
Although O3 has had few associations in earlier ACS studies using between-city contrasts,1,2,4 exposure to this pollutant is considered a health threat in the LA region, which has some of the highest levels in the United States.18 For O3, we obtained data at 42 sites in and around the LA basin from the California Air Resources Board database. We interpolated 2 surfaces using a universal kriging algorithm: one based on the average of the 4 highest 8-hour concentrations over the year 2000 and another based on the expected peak daily concentration, which is a statistical measure designed to assess the likely exceedance of the 8-hour average at the site based on the previous 3 years (1999–2001). Both measures are used as a basis for either federal or state designation of nonattainment areas. They both capture extreme events, but the expected peak daily concentration provides more stability for estimation of spatial patterns than the 1-year measures based on the 4 highest days. Few studies of chronic effects have found significant ozone effects, although acute effects of a small magnitude have been observed.30 Thus, it seems plausible that an ozone effect would be manifest in those areas most likely to experience exceedances.
Finally, we assessed the impact of traffic by assigning buffers that included zip code-area centroids within either 500 or 1000 meters of a freeway. The U.S. Bureau of the Census feature class codes define freeways as having “limited access,” a numbered assignment, and a speed limit of greater than 50 miles per hour.31 This distance from the zip code-area centroid to the freeway approximated exposure to traffic pollution, which may exert independent effects in addition to pollutants such as PM2.5 and O3 that vary over larger areas.8
Complete residential history information was unavailable for the entire cohort, although we do have information on whether respondents moved between enrollment and 1992 or thereafter (approximately 5633 in LA). Of this group, only 16% moved during follow up, and this diminishes the potential for exposure misclassification resulting from residential mobility.
We used Cox proportional hazards regression for our main analyses of association between air pollution and mortality.32 Because the units of analyses were small zip code areas and previous analyses had indicated spatial autocorrelation in the residual variation of some health effects models, we also developed and used a new spatial random effects Cox model as a crossvalidation of the standard model. We have previously shown that survival experience clusters by community and is spatially autocorrelated between communities.2,3 Lack of statistical control for these factors can bias the estimates of air pollution effects and underestimate associated standard errors.3,33 To characterize the statistical error structure of survival data, novel statistical methodology and computer software have been developed to incorporate spatial clustering at the zip code area. Our model can be expressed mathematically in the form
where hij is the hazard function or instantaneous hazard probability of death for the ith subject in the jth ZCA, whereas s indicates the stratum (defined by sex, race, and age). Here h0, s(t) is the baseline hazard function. The ηj are positive random effects representing the unexplained variation in the response among neighborhoods, in this case zip code areas. Only the moments of the random effects need to be specified within our modeling framework: E(ηj) = 1 and Var(ηj) = τ2. The vector xij represents the known risk factors for the response such as air pollution, smoking habits, and diet. The regression parameter vector is denoted by β. Estimates of the regression vector β, random effects, their variance, and correlation parameter are obtained by methods previously used for random-effects survival models.33 Thiessen polygons, which ensure that all points within the polygon are closer to the centroid of that polygon than to any other centroid, were used to assign first-order nearest neighbor contiguity between the zip code areas. These were derived using ArcView 3.2 (ESRI Corp., Redlands, CA). The standard Moran's I tests of spatial autocorrelation were applied to the random effects.
Figure 1 illustrates the pollution surface used in our main analysis, and the Appendix Figure (available with the online version of this article) illustrates the absolute and relative standard errors of estimation for the interpolated universal kriging surface. Approximately 50% of the modeled surface has errors that are less than 15% of the monitored value, whereas 67% of the surface lies within 20% of the monitored values. For the most part, absolute standard errors for the densely populated areas of the study region are less than 3 μg/m3. Only on the periphery of the study area do errors become large compared with monitored values, but these places have very few of our study subjects. Interestingly, the range of exposure within LA (20 μg/m3) exceeds what we observed in previous studies based on contrasts among 116 cities (16 μg/m3).4
Results for all-cause and cause-specific deaths are reported in Table 1. This table shows the PM2.5 effect with varying levels of control for confounding. Relative risks (RRs) are expressed as 10 μg/m3 exposure contrasts in PM2.5 followed by the 95% confidence interval (95% CI). Using the example of all-cause mortality and with each succeeding stage, including the previous individual-level controls, we find that for PM2.5 alone and controlling just for age, sex, and race, the RR is 1.24 (95% CI = 1.11–1.37), whereas the RR with the 44 individual confounders4 is 1.17 (1.05–1.30). All subsequent results include the 44 individual-level control variables and one or more ecologic variables. For example, with 44 individual variables and the ecologic variable of unemployment, the RR of PM2.5 is 1.15 (1.03–1.28). When we add 4 social factors extracted from the principal component analysis (and accounting for 81% of the total variance in the social variables), the RR is 1.15 (1.03–1.29). Including all ecologic variables associated with mortality in bivariate models reduces the pollution coefficient to RR of 1.11 (0.99–1.25). Finally, for the parsimonious model that includes ecologic confounder variables that both reduce the pollution coefficient and have associations with mortality, the RR is 1.11 (0.99–1.25).
Comparing these results directly with the earlier analyses using between-community contrasts, the health effects are nearly 3 times greater for this analysis (ie, 17% increase compared with 6% in earlier studies in models that control for the 44 individual confounders). With control for neighborhood confounders, effect estimates are still approximately 50% to 90% higher than in previous analyses.
In models with only individual covariates and PM2.5, some residual spatial autocorrelation was present in the random effects from the model clustered on zip code area. We attempted to remove this autocorrelation by fitting a model with a ρ autocorrelation term that used mortality information from nearest neighbors as a predictor of mortality in the ZCA j, but the autocorrelation persisted (results not shown). When contextual socioeconomic status variables were included in the model, however, the Moran's I tests revealed no significant spatial autocorrelation. Table 2 shows the results of the Moran's I test for all-cause and ischemic heart disease mortality. Visual inspection of the random effects, ηj, confirmed the results from the Moran's I testing.
Sensitivity analyses using weighted estimation with weights equal to the inverse of the standard error on the universal kriging exposure model demonstrated that the risk estimates were robust to measurement error in the exposure estimate (results not shown). Point estimates remained elevated.
Our results suggest that the chronic health effects associated with intraurban gradients in exposure to PM2.5 may be even larger than previously reported associations across metropolitan areas. Using the direct comparison to previous ACS studies, we see effects that are nearly 3 times larger than in models relying on between-community exposure contrasts. We also note convincing evidence of specificity in these health effects, with a stronger association between air pollution and ischemic heart disease than for the more general measures of cardiopulmonary deaths or all-cause mortality (Fig. 2 displays the ordering in the risks presented in the tables). These findings concur with recent studies at the metropolitan level, which again demonstrate that ischemic heart disease drives the cardiopulmonary association with air pollution.5
Among the cancer deaths, we also observe ordering in the risks, with decreasing risks as we move from lung cancer to digestive cancer to all cancers. Given that the lung would be most directly affected by air pollution, this finding gives corroborative evidence that the association did not occur by chance.
The larger effects in LA raise the question of whether some underlying aspect of this subcohort differs in characteristics that modify the association between mortality and air pollution. We compared the full cohort with the LA subcohort and found no major differences in attributes likely to modify the air pollution–mortality association, with the exception that the LA cohort was better educated. Based on the findings from the earlier analyses in which subjects with lower education experienced larger health effects,1,2,4 we would expect the effect size in LA to be smaller than in the full cohort. Thus, differences in underlying characteristics appear unlikely to explain the larger effects we observed in LA.
In comparing our results with the earlier national-level ACS studies, we examined the reduction in PM2.5 levels in LA to 50 other metropolitan areas that had data for 1980 and for the year 2000 we used in our study. (See Krewski et al2 for a description of the data.) The mean reduction was 31%, with a range from 0.4% to 59%. LA experienced a reduction of 24.5%, just above the lowest quartile of 23.5%. PM2.5 has therefore declined at a slightly slower rate in LA compared with much of the United States. If we assume that current PM2.5 in LA is at 75.5% (ie, accounting for 24.5% reduction) of the 1980 value and the average metropolitan area is at 69% of the 1980 value, some of the increase in the risks may be attributable to the relatively smaller reductions in LA. We tested this scaling effect by computing the ratio of reductions (0.69/75.5 = 0.914) and multiplied our raw coefficients by this factor before estimating the RR. The RR declines for all-cause mortality with the 44 individual variables to just over 15% and with maximal adjustment for confounders to 10%. Although reduced by up to 1.6%, we conclude that the majority of the increase over previous estimates reported by Pope et al34 is probably not attributable to relative differences in the rate of reduction in ambient air pollution.
The findings for endocrine deaths also reveal another interesting possibility. Chronic air pollution exposure, similar to acute exposures,34 may adversely affect people with diabetes more than the general population. Alternatively, the finding may indicate some uncontrolled confounding because we expect people with type 2 diabetes to live in neighborhoods with poorer social environments. This possibility appears unlikely because of the extensive control we applied for contextual neighborhood variables. This potential problem appears improbable because we see internal validity in the effects of social confounders measured in the zip code areas.
Although the accidental deaths were unexpectedly elevated in men, subsequent analyses revealed that the risks were attributable to deaths in the early years of the cohort before causes of death were coded in detail. As a result, we were unable to assess specific causes for this elevation.
Ozone had few elevated risks in any of our analyses and did not confound the relationship between particles and mortality. This finding agrees with earlier ACS studies indicating that ozone is not associated with elevated mortality risk,2,4 but contradicts studies on nonsmoking Adventists in southern California, where associations between lung cancer in males and ozone exposure were detected.35 Recent national studies have reported elevated acute risks of ozone exposure,30 but risk estimates were small, as would be expected of a study on acute as compared with chronic effects.36
In assessing the association with freeway buffers, point estimates were particularly elevated for lung cancer, endocrine, and digestive mortality. The PM—mortality association remained robust to the freeway buffer, and risk estimates were unchanged when this variable was included in the model. Although imprecision in the freeway exposures resulting from the zip code area assignment of proximity may have biased our results toward the null, we did observe a RR for lung cancer of 1.44, and the other cause-specific mortality metrics indicate that more precise estimation of traffic effects are warranted in future research.
In previous studies based on the ACS cohort, all individuals within the same metropolitan area were assigned the same level of exposure based on the average ambient concentration observed at fixed-site air pollution monitors in that city. We hypothesized that the use of such a broad ecologic indicator of exposure leads to exposure measurement error, which in turn can bias estimates of mortality associated with air pollution exposure. Mallick et al37 analyzed the effect of this source of exposure measurement error based on plausible assumptions about error magnitude in the Six-Cities Study of air pollution and mortality.6 This investigation suggested that the RR of mortality resulting from particulate air pollution may be underestimated by a factor of approximately 2- to 3-fold as a consequence of exposure misclassification, a finding consistent with the present results.
We recognize the possibility of exposure measurement error from using recent exposure models for a cohort enrolled in 1982. There are empiric as well as theoretical reasons that prevent this potential problem from seriously limiting the results. Empirically, other ACS analyses done at the metropolitan scale have found that these more recent exposure estimates predicted mortality with results similar to those based on earlier monitoring data.4 Also, the well-known meteorologic and topographic conditions of LA, along with a dominant on-shore breeze and steep mountains to the north and east, control much of the spatial pattern of pollution in the region. Our results agree with findings of earlier studies on the pattern of spatial variation in PM.38 Although levels may rise and fall in absolute terms, major changes in the spatial patterns within the region over time appear unlikely, and the rank ordering among assigned exposures should be maintained.
From a theoretical perspective, even if spatially heterogeneous changes to pollution levels within a city occurred as a result of new emissions during the follow up, this would lead to larger exposure measurement error, and a bias toward the null would dominate, assuming a classic error structure. With a Berkson error structure, the variance of the dose–response estimate would be inflated.17 In either case, with current exposure models, the health effects likely have a lower probability of false-positive error, and we would expect the measurement error to reduce effect sizes and inflate their variance. High dose–response relationships can be caused by underestimation of concentrations in the high-exposure areas, but for these areas, the monitoring networks tend to be dense and the kriging errors were smaller than in most of the study area. Finally, the findings were robust to weighting for errors in the kriging estimate (ie, eastern parts of the LA region), which decreases the likelihood that elevated risks arise as a result of underestimation in the high-exposure groups.
Although we are unable to reconstruct likely exposures to PM2.5 for our exposure surface, we have assessed the relationship at 51 central monitors between PM2.5 measured in 1980 and those of a period similar to that of our 2000 estimates (ie, 1999–2000). These data were used in previous national studies,4 in which details on their derivation are available. Figure 3 illustrates the regression scatterplot for the 1999–2000 values on the 1980 measurements. The coefficient of determination is approximately 61%, and overall the latter periods are predicted well by the earlier measurements.
In addition, we examined the relationship of historical PM10 data in the LA area with the 2000 PM2.5 estimates used in our analysis. The period of maximal overlap between the sites occurs in 1993, where we had 8 PM10 readings at the same locations as the subsequent PM2.5 measurements. By regressing the 2000 PM2.5 measurements on those PM10 measurements in 1993, we observe an R2 value of 90% (Fig. 4).
In both of these correlation analyses comparing earlier monitoring data with more recent PM2.5 measurements, we found evidence that areas with higher particle concentrations in earlier periods were likely to retain their spatial ranking. Those metropolitan areas likely to be high in 1980 also had a similar tendency in 2000.
Only the Norwegian cohort study has used time windows of exposure.39 In this study, the authors found that timing of the exposure window had little influence on the estimation of health effects; they used exposure windows in the middle of the follow-up period for most of their results. All of the other cohort studies have taken a similar approach to ours and computed the risk based on relatively short-term air pollution monitoring data. A case–control study in Stockholm, Sweden, investigated time windows for lung cancer.40 This study found that windows of exposure 20 years before disease onset were more strongly associated with cancer than later periods. In our study, we found elevated risks of lung and digestive cancers, even with the more recent exposure model. The likely stability in the spatial pattern of exposure in LA probably accounts for this similarity of our findings to the 2 European studies that have used time windows.
Generally, our results agree with recent evidence suggesting that intraurban exposure gradients may be associated with even larger health effects than reported in interurban studies. Hoek et al8 reported a doubling of cardiopulmonary mortality (RR = 1.95; 95% CI 1.09–3.52) for Dutch subjects living near major roads. Canadian cohort studies controlling for medical care utilization and preexisting chronic conditions through record linkage have also uncovered large health effects with proximity to major roads at the intraurban scale.16 Recent results from the cohort in Norway also suggest associations between intraurban gradients in gaseous pollutants and mortality.40 All of these studies have implicated traffic as the source of pollution associated with the larger observed effects. In LA, the proportion of primary particles attributable to traffic is approximately 3.7%, whereas in the rest of the country, it is 0.75%.41 Thus, beyond improved precision in the exposure models, the larger health effects reported here may be partly the result of higher proportions of traffic particulate in LA.
No previous studies have assessed associations based on a continuous exposure model with PM2.5, which limits the use of the estimates for current policy debates that tend to focus on fine particles. In this study, we used PM2.5 with a continuous exposure metric that promotes comparison with previous studies on health effects and contributes to current regulatory debates.
1. Pope CA, Thun MJ, Namboodiri MM, et al. Particulate air pollution as a predictor of mortality in a prospective study of US adults. Am J Respir Crit Care Med
2. Krewski D, Burnet RR, Goldberg MS, et al. Reanalysis of the Harvard Six Cities and the American Cancer Society Study of Particulate Air Pollution and Mortality, Phase II: Sensitivity Analysis
. Cambridge, MA: Health Effects Institute; 2000.
3. Jerrett M, Burnett R, Willis A, et al. Spatial analysis of the air pollution-mortality association in the context of ecologic confounders. J Toxicol Environ Health
4. Pope CA, Burnett RT, Thun MJ, et al. Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution. JAMA
5. Pope CA, Burnett RT, Thurston GD, et al. Cardiovascular mortality and long-term exposure to particulate air pollution. Circulation
6. Dockery D, Pope CA3rd, Xu X, et al. An association between air pollution and mortality in six US cities. N Engl J Med
7. WHO Strategy on Air Quality and Health
. Geneva: World Health Organization; 2001.
8. Hoek G, Brunekreef B, Goldbohm S, et al. Association between mortality and indicators of traffic-related air pollution in The Netherlands: a cohort study. Lancet
9. Peters A, Pope CA. Cardiopulmonary mortality and air pollution. Lancet
10. Briggs D, de Hoogh C, Gulliver J, et al. A regression-based method for mapping traffic-related air pollution: application and testing in four contrasting urban environments. Sci Total Environ
11. Zhu Y, Hinds W, Kim S, et al. Study of ultrafine particulates near a major highway with heavy-duty diesel traffic. Atmos Environ
12. Jerrett M, Burnett RT, Kanaroglou P, et al. A GIS-environmental justice analysis of particulate air pollution in Hamilton, Canada. Environment and Planning
13. Brauer M, Hoek G, van Vliet P, et al. Estimating long-term average particulate air pollution concentrations: application of traffic indicators and geographic information systems. Epidemiology
14. Brunekreef B, Holgate S. Air pollution and health. Lancet
15. O'Neill M, Jerrett M, Kawachi I, et al. Health, wealth and air pollution. Environ Health Perspect
16. Finkelstein M, Jerrett M, DeLuca P, et al. A cohort study of income, air pollution and mortality. Can Med Assoc J
17. Thomas DC, Stram D, Dwyer J. Exposure measurement error: influence on exposure—disease relationships and methods of correction. Annu Rev Public Health
18. Künzli N, McConnell R, Bates D, et al. Breathless in Los Angeles: the exhausting search for clean air. Am J Public Health
19. Diez Roux AV. A glossary for multilevel analysis. J Epidemiol Community Health
20. Pickett KE, Pearl M. Multilevel analyses of neighbourhood socioeconomic context and health outcomes: a critical review. J Epidemiol Community Health
21. Curtis S, Taket A. Health and Societies: Changing Perspectives
. London: Edward Arnold; 1996.
22. Macintyre S, Ellaway A. Ecological approaches: rediscovering the role of the physical and social environments. In: Berkman L, Kawachi I, eds. Social Epidemiology
. Oxford: Oxford University Press; 2000.
23. Evans GW, Kantrowitz E. Socioeconomic status and health: the potential role of environmental risk exposure. Annu Rev Public Health
24. Willis A, Krewski D, Jerrett M, et al. Selection of ecologic covariates in the American Cancer Society study. J Toxicol Environ Health
25. Janssen N, Harssema H, Brunekreef B, et al. Assessment of exposure to traffic related air pollution of children attending schools near motorways. Atmos Environ
26. Meng Q, Turpin B, Korn L, et al. Influence of ambient (outdoor) sources on residential indoor and personal PM2.5
concentrations: analyses of RIOPA data. J Expo Anal Environ Epidemiol
27. Luginaah I, Jerrett M, Elliott S, et al. Health profiles of Hamilton: characterizing neighbourhoods for health investigations. GeoJournal
28. Krieger N, Chen JT, Waterman PD, et al. Geocoding and monitoring of US socioeconomic inequalities in mortality and cancer incidence: does the choice of area-based measure and geographic level matter? The public health disparities geocoding project. Am J Epidemiol
29. Burrough PA, McDonnell RA. Principals of Geographical Information Systems
. Oxford University Press; 1998.
30. Bell ML, McDermott A, Zeger SL, et al. Ozone and short-term mortality in 95 US urban communities, 1987–2000. JAMA
31. First Edition TIGER/Line Technical Documentation: US Bureau of the Census Feature Class Codes (CFCC)
32. Hosmer W, Lemeshow S. Applied Survival Analysis
. New York: Wiley and Sons; 1999.
33. Ma R, Krewski D, Burnett RT. Random effects Cox models: a Poisson regression modelling approach. Biometrika
34. Goldberg MS, Burnett R, Bailar JC III, et al. The association between daily mortality and ambient air particle pollution in Montreal, Quebec. 2. Cause-specific mortality. Environ Res
35. Abbey DE, Nishino N, McDonnell WF, et al. Long-term inhalable particles and other air pollutants related to mortality in nonsmokers. Am J Respir Crit Care Med
36. Künzli N, Medina S, Kaiser R, et al. Assessment of deaths attributable to air pollution: should we use risk estimates based on time series or on cohort studies? Am J Epidemiol
37. Mallick R, Fung F, Krewski D. Adjusting for measurement error in the Cox proportional hazards regression model. J Cancer Epidemiol Prevent
38. Kim BM, Teffera S, Zeldin MD. Characterization of PM2.5
in the South Coast Air Basin of southern California: part 1—spatial variations. J Air Waste Manag Assoc
39. Nafstad P, Lund Håheim L, Wisløff T, et al. Urban air pollution and mortality in a cohort of Norwegian men. Environ Health Perspect
40. Nyberg F, Gustavsson P, Jarup L, et al. Urban air pollution and lung cancer in Stockholm. Epidemiology
41. 1999 National Emission Inventory Documentation and Data—Final Version 3.0. (1999). US Environmental Protection Agency, Technology Transfer Network Clearinghouse for Inventories & Emissions Factors web site. Available at: http://www.epa.gov/ttn/chief/net/1999inventory.html
. Accessed February 1, 2005.