Portal de Periódicos da CAPES

National Kriging Exposure Estimation: Liao et al. Respond

2007; National Institute of Environmental Health Sciences; Volume: 115; Issue: 7 Linguagem: Francês

10.1289/ehp.10205r

1552-9924

Duanping Liao, Donna J. Peuquet, Hung-Mo Lin, Yinkang Duan, Eric A. Whitsel, Richard L. Smith, Gerardo Heiss,

Climate Change and Health Impacts

Vol. 115, No. 7 PerspectivesOpen AccessNational Kriging Exposure Estimation: Liao et al. Respond Duanping Liao, Donna J. Peuquet, Hung-Mo Lin, and Yinkang Duan Eric A. Whitsel, Richard L. Smith, and Gerardo Heiss Duanping Liao , Donna J. Peuquet , Hung-Mo Lin , Yinkang Duan , Eric A. Whitsel , Richard L. Smith , and Gerardo Heiss Published:1 July 2007https://doi.org/10.1289/ehp.10205RCited by:8AboutSectionsPDF ToolsDownload CitationsTrack Citations ShareShare onFacebookTwitterLinked InReddit Szpiro et al. suggest that our findings Liao et al. (2006) do not adequately support using national-scale, log-normal ordinary kriging to estimate daily mean concentrations of PM10 (particulate matter with aerodynamic diameter ≤ 10 μm) at unmonitored locations in the contiguous United States. They posit that the absence of the cross-validation SE prevents evaluating the validity of kriging estimation, as we implemented in this context, and the comparability of both regional-versus national-scale kriging and manually modified versus semiautomated, default-calculated semivariograms.Little literature is available on the use of kriging methods to estimate daily air pollution data for large population–based multi-center epidemiologic studies. The four studies cited by Szpiro et al. (Dockery et al. 1993; Jerrett et al. 2005; Miller et al. 2007; Pope et al. 1995) all used cohort analyses for which only long-term average exposures are required, and only one of those (Jerrett et al. 2005) actually involved interpolation methods at all, although the study was restricted to a single city. In contrast, our objective was to create an interpolated daily pollutant concentration database for a multisite population-based epidemiologic study.The cross-validation mean square error (MSE) mentioned by Szpiro et al. is also termed the "root-mean-square (RMS) prediction error," which is the empirical SE based on the mean square of the predictions, as opposed to SE, the mathematical formula for the RMS prediction error. RMS and SE, both are available from ArcView (ESRI Inc., Redlands, CA), are often considered jointly as an alternative measure (to RMSS) of the validity of spatial analysis.The average RMS and SE from 366 daily PM10 spherical model cross-validations based on year 2000 PM10 data were 19.48 and 16.19 μg/m3, respectively, from the log-normal regular kriging, and 26.43 and 25.60 μg/m3, respectively from a ordinary kriging. The validity of the model is supported by RMSS alone (≈1), by the similarity of RMS and SE, and by SPE (≈ 0). Additionally, the average daily SD of PM10 measured at the monitor locations was 27.20 μg/m3. Comparing SD with the kriging-RMS provides a measure of the reduction in error due to interpolation. If RMS is less than the SD, then the kriging approach has some benefit, compared with using long-run averages. From both ordinary and log-normal kriging, especially for the latter, we see a notable reduction in RMS compared with SD. Meanwhile, substantial variability remains, suggesting that kriging error should be taken into account when using the kriged values.Szpiro et al. also implicitly criticize our use of daily kriging when the objective was to interpolate daily data. Spatial–temporal models have potentially greater power than a 1-day-at-a-time spatial analysis but are not easy to apply in practice, with large datasets and many missing values.Regional kriging could be superior to national kriging if the spatial dependence parameters (range, sill, and nugget) vary substantially from region to region, in which case a national kriging model could result in misspecified covariances. However, regional kriging also uses fewer data points to estimate those parameters and could result in greater errors. We would welcome theoretical or empirical studies that could cast further light on this trade-off. However, as far as our article (Liao et al. 2006) is concerned, our main purpose was to note that the national kriging method appears to be competitive when assessed by overall RMS error. We compared the results of regional- and national-scale kriging on a small set (17%) of days when the largest number of monitors (≥ 400) were reporting data—a scenario heavily favoring regional spatial interpolation strategies. On the remaining days when only 120–400 monitors were reporting data, regional kriging was inherently problematic given the restricted availability of monitors within regions. Szpiro et al. suggest that the problems of interpolation near the boundary could be solved by "overlapping," but this is only one of the issues encountered using regional-kriging: One would still need to decide how to consistently define the regions, considering the number of available data points that change substantially from day to day, to achieve a meaningful reduction in RMS error.Based on the 12 "optimal" days in 2000, the average RMS and SE were 12.68 and 12.82 μg/m3, respectively, from the national scale kriging, compared with 12.22 and 12.49 μg/m3, respectively, from regional-scale kriging (Liao et al. 2006). These results, together with RMSS and SPE we reported, support our conclusion that national kriging performs comparably to regional kriging even when restricted to optimal days.Szpiro et al. correctly note that it is possible to improve the RMSS values by manual adjustment. However, typically we found that when one of the validity measures (RMSS, PE, or SPE) was improved by manual adjustment, other measures became worse. It is difficult to manually adjust models to improve all cross-validation parameters simultaneously. Manually adjusting daily semivariogrms is not feasible when kriging over 10 years. Moreover, the predicted SE at unmeasured locations was uniformly lower in automatically fit models.Szpiro et al. are correct that cross-validation may not be representative of the performance at participant address locations, although it is unclear what alternative methods they would like us to use. The ability to do semi-automatic cross-validations was a major attraction of ArcView and, despite limitations, is the best tool we know for validating spatial predictions.The semiautomated kriging approach presents considerable advantages in estimating daily residential-level pollutant concentrations in large cohorts over long periods. Our proposed method (Liao et al. 2006) used log-normal kriging based on a spherical model to interpolate daily data on a national scale, and the weighted least squares method of parameter estimation without manual adjustment. We believe that the cross-validation statistics, presented in our article and amplified here, provide adequate support for these recommendations against reasonable alternatives that we considered.ReferencesDockery D, Pope C, Xu X, Spengler J, Ware J, Fay Met al.. 1993. An association between air pollution and mortality in six US cities. N Engl J Med 329:1753-1759doi:10.1056/NEJM1993120932924018179653. Crossref, Medline, Google ScholarJerrett M, Burnett R, Ma R, Pope C, Krewski D, Newbold Ket al.. 2005. Spatial analysis of air pollution and mortality in Los Angeles. Epidemiology 16:727-736doi:10.1097/01.ede.0000181630.15826.7d16222161. Crossref, Medline, Google ScholarLiao D, Peuquet DJ, Duan Y, Whitsel EA, Dou J, Smith RLet al.. 2006. GIS approaches for the estimation of residential-level ambient PM concentrations. Environ Health Perspect 114:1374-1380doi:10.1289/ehp.9169[Online 8 June 2006]16966091. Link, Google ScholarMiller K, Siscovick D, Sheppard L, Shepherd K, Sullivan J, Anderson Get al.. 2007. Long-term exposure to air pollution and incidence of cardiovascular events in women. N Engl J Med 356:447-458doi:10.1056/NEJMoa05440917267905. Crossref, Medline, Google ScholarPope C, Thun M, Namboodiri M, Dockery D, Evans J, Speizer Fet al.. 1995. Particulate air pollution as a predictor of mortality in a prospective study of US adults. Am J Respir Crit Care Med 151:669-6747881654. Crossref, Medline, Google ScholarFiguresReferencesRelatedDetailsCited byGondalia R, Holliday K, Baldassari A, Justice A, Stewart J, Liao D, Yanosky J, Engel S, Jordahl K, Bhatti P, Horvath S, Assimes T, Pankow J, Demerath E, Guan W, Fornage M, Bressler J, North K, Conneely K, Li Y, Hou L, Baccarelli A and Whitsel E (2020) Leukocyte Traits and Exposure to Ambient Particulate Matter Air Pollution in the Women's Health Initiative and Atherosclerosis Risk in Communities Study, Environmental Health Perspectives, 128:1, Online publication date: 1-Jan-2020.Gondalia R, Avery C, Napier M, Méndez-Giráldez R, Stewart J, Sitlani C, Li Y, Wilhelmsen K, Duan Q, Roach J, North K, Reiner A, Zhang Z, Tinker L, Yanosky J, Liao D and Whitsel E (2017) Genome-wide Association Study of Susceptibility to Particulate Matter–Associated QT Prolongation, Environmental Health Perspectives, 125:6, Online publication date: 23-Jun-2017. Penfold B, Hafner H and Brown S (2013) Predicting areas of high diesel particulate matter emissions in Phoenix, Arizona, using spatial analysis techniques, Journal of Environmental Engineering and Science, 10.1680/jees.2013.0050, 8:4, (493-504), Online publication date: 1-Nov-2013. Sivaraman V, Carrapetta J, Hu K and Luxan B (2013) HazeWatch: A participatory sensor system for monitoring air pollution in Sydney 2013 IEEE 38th Conference on Local Computer Networks Workshops (LCN Workshops), 10.1109/LCNW.2013.6758498, 978-1-4799-0540-9, (56-64) Penfold B, Hafner H and Brown S (2010) Predicting areas of high diesel particulate matter emissions in Phoenix, Arizona, using spatial analysis techniquesA paper submitted to the Journal of Environmental Engineering and Science., Canadian Journal of Civil Engineering, 10.1139/L09-162, 37:2, (354-365), Online publication date: 1-Feb-2010. Whitsel E, Quibrera P, Christ S, Liao D, Prineas R, Anderson G and Heiss G (2009) Heart Rate Variability, Ambient Particulate Matter Air Pollution, and Glucose Homeostasis: The Environmental Epidemiology of Arrhythmogenesis in the Women's Health Initiative, American Journal of Epidemiology, 10.1093/aje/kwn400, 169:6, (693-703), Online publication date: 15-Mar-2009., Online publication date: 15-Mar-2009. Zhang Z, Whitsel E, Quibrera P, Smith R, Liao D, Anderson G and Prineas R (2009) Ambient Fine Particulate Matter Exposure and Myocardial Ischemia in the Environmental Epidemiology of Arrhythmogenesis in the Women's Health Initiative (EEAWHI) Study, Environmental Health Perspectives, 117:5, (751-756), Online publication date: 1-May-2009. Liao D, Whitsel E, Duan Y, Lin H, Quibrera P, Smith R, Peuquet D, Prineas R, Zhang Z and Anderson G (2008) Ambient Particulate Air Pollution and Ectopy—The Environmental Epidemiology of Arrhythmogenesis in Women's Health Initiative Study, 1999–2004, Journal of Toxicology and Environmental Health, Part A, 10.1080/15287390802445483, 72:1, (30-38), Online publication date: 7-Nov-2008. Vol. 115, No. 7 July 2007Metrics About Article Metrics Publication History Originally published1 July 2007Published in print1 July 2007 Financial disclosuresPDF download License information EHP is an open-access journal published with support from the National Institute of Environmental Health Sciences, National Institutes of Health. All content is public domain unless otherwise noted. Note to readers with disabilities EHP strives to ensure that all journal content is accessible to all readers. However, some figures and Supplemental Material published in EHP articles may not conform to 508 standards due to the complexity of the information being presented. If you need assistance accessing journal content, please contact [email protected]. Our staff will work with you to assess and meet your accessibility needs within 3 working days.