Bayesian Modelling of Spatial Data Using Markov Random Fields, With Application to Elemental Composition of Forest Soil |
| |
Authors: | Linda Werner Hartman |
| |
Institution: | (1) Mathematical Statistics, Centre for Mathematical Sciences, Lund University, Box 118, 221 00 Lund, Sweden |
| |
Abstract: | Spatial datasets are common in the environmental sciences. In this study we suggest a hierarchical model for a spatial stochastic field. The main focus of this article is to approximate a stochastic field with a Gaussian Markov Random Field (GMRF) to exploit computational advantages of the Markov field, concerning predictions, etc. The variation of the stochastic field is modelled as a linear trend plus microvariation in the form of a GMRF defined on a lattice. To estimate model parameters we adopt a Bayesian perspective, and use Monte Carlo integration with samples from Markov Chain simulations. Our methods does not demand lattice, or near-lattice data, but are developed for a general spatial data-set, leaving the lattice to be specified by the modeller. The model selection problem that comes with the artificial grid is in this article addressed with cross-validation, but we also suggest other alternatives. From the application of the methods to a data set of elemental composition of forest soil, we obtained predictive distributions at arbitrary locations as well as estimates of model parameters. |
| |
Keywords: | Gaussian Markov random field Markov chain Monte Carlo nonlattice data bilinear interpolation |
本文献已被 SpringerLink 等数据库收录! |
|