Spatio-temporal modeling of particulate matter concentration through the SPDE approach

AStA Advances in Statistical Analysis - Tập 97 - Trang 109-131 - 2012
Michela Cameletti1, Finn Lindgren2, Daniel Simpson2, Håvard Rue2
1Dip. di Matematica, Statistica, Informatica e Applicazioni, Università di Bergamo, Bergamo, Italy
2Dept. of Mathematical Sciences, Norwegian University of Science and Technology, Trondheim, Norway

Tóm tắt

In this work, we consider a hierarchical spatio-temporal model for particulate matter (PM) concentration in the North-Italian region Piemonte. The model involves a Gaussian Field (GF), affected by a measurement error, and a state process characterized by a first order autoregressive dynamic model and spatially correlated innovations. This kind of model is well discussed and widely used in the air quality literature thanks to its flexibility in modelling the effect of relevant covariates (i.e. meteorological and geographical variables) as well as time and space dependence. However, Bayesian inference—through Markov chain Monte Carlo (MCMC) techniques—can be a challenge due to convergence problems and heavy computational loads. In particular, the computational issue refers to the infeasibility of linear algebra operations involving the big dense covariance matrices which occur when large spatio-temporal datasets are present. The main goal of this work is to present an effective estimating and spatial prediction strategy for the considered spatio-temporal model. This proposal consists in representing a GF with Matérn covariance function as a Gaussian Markov Random Field (GMRF) through the Stochastic Partial Differential Equations (SPDE) approach. The main advantage of moving from a GF to a GMRF stems from the good computational properties that the latter enjoys. In fact, GMRFs are defined by sparse matrices that allow for computationally effective numerical methods. Moreover, when dealing with Bayesian inference for GMRFs, it is possible to adopt the Integrated Nested Laplace Approximation (INLA) algorithm as an alternative to MCMC methods giving rise to additional computational advantages. The implementation of the SPDE approach through the R-library INLA ( www.r-inla.org ) is illustrated with reference to the Piemonte PM data. In particular, providing the step-by-step R-code, we show how it is easy to get prediction and probability of exceedance maps in a reasonable computing time.

Tài liệu tham khảo

Banerjee, S., Carlin, B., Gelfand, A.: Hierarchical Modeling and Analysis for Spatial Data. Monographs on Statistics and Applied Probability. Chapman and Hall, New York (2004)

Banerjee, S., Gelfand, A., Finley, A., Sang, H.: Gaussian predictive process models for large spatial datasets. J. R. Stat. Soc. B 70(4), 825–848 (2008)

Besag, J.: Spatial interaction and the statistical analysis of lattice systems (with discussion). J. R. Stat. Soc. B 36(2), 192–225 (1974)

Cameletti, M., Ignaccolo, R., Bande, S.: Comparing spatio-temporal models for particulate matter in Piemonte. Environmetrics (2011). doi:10.1002/env.1139

Cocchi, D., Greco, F., Trivisano, C.: Hierarchical space-time modelling of PM10 pollution. Atmos. Environ. 41, 532–542 (2007)

Cressie, N.: Statistics for Spatial Data. Wiley, New York (1993)

Cressie, N., Johannesson, G.: Fixed rank kriging for large spatial datasets. J. R. Stat. Soc. B 70, 209–226 (2008)

Cressie, N., Wikle, C.: Statistics for Spatio-Temporal Data. Wiley, New York (2011)

Fassò, A., Finazzi, F.: Maximum likelihood estimation of the dynamic coregionalization model with heterotopic data. Environmetrics 22, 735–748 (2011)

Finardi, S., De Maria, R., D’Allura, A., Cascone, C., Calori, G., Lollobrigida, F.: A deterministic air quality forecasting system for Torino urban area, Italy. Environ. Model. Softw. 23(3), 344–355 (2008)

Furrer, R., Genton, M., Nychka, D.: Covariance tapering for interpolation of large spatial datasets. J. Comput. Graph. Stat. 15(3), 502–523 (2006)

Gelfand, A., Diggle, P., Fuentes, M., Guttorp, P. (eds.): Handbook of Spatial Statistics. Chapman & Hall, New York (2010)

Lindgren, F., Rue, H., Lindström, J.: An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach (with discussion). J. R. Stat. Soc. B 73(4), 423–498 (2011)

R Development Core Team: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna (2011). URL http://www.R-project.org, ISBN 3-900051-07-0

Rue, H., Held, L.: Gaussian Markov Random Fields. Theory and Applications. Chapman & Hall, New York (2005)

Rue, H., Martino, S., Chopin, N.: Approximate Bayesian inference for latent Gaussian model by using integrated nested Laplace approximations (with discussion). J. R. Stat. Soc. B 71, 319–392 (2009)

Sahu, S.: Hierarchical Bayesian Models for Space-time Air Pollution Data. Rao, C. (ed.): Handbook of Statistics—Time Series Analysis, Methods and Applications. Handbook of Statistics, vol. 30. Elsevier, Amsterdam (2011)

Samet, J., Dominici, F., Curriero, F., Coursac, I., Zeger, S.: Fine particulate air pollution and mortality in 20 US cities: 1987–1994. N. Engl. J. Med. 343, 1742–1749 (2000)

Samoli, E., Peng, R., Ramsay, T., Pipikou, M., Touloumi, G., Dominici, F., Burnett, R., Cohen, A., Krewski, D., Samet, J., Katsouyanni, K.: Acute effects of ambient particulate matter on mortality in Europe and North America: results from the APHENA study. Environ. Health Perspect. 116, 1480–1486 (2008)

Thông tin
Thông tin xuất bản

AStA Advances in Statistical Analysis

Tập 97

109-131

Thông tin tác giả