Bayesian Analysis of Geostatistical Models With an Auxiliary Lattice

Jincheol Park, Faming Liang

Research output: Contribution to journalArticlepeer-review

4 Scopus citations

Abstract

The Gaussian geostatistical model has been widely used for modeling spatial data. However, this model suffers from a severe difficulty in computation: it requires users to invert a large covariance matrix. This is infeasible when the number of observations is large. In this article, we propose an auxiliary lattice-based approach for tackling this difficulty. By introducing an auxiliary lattice to the space of observations and defining a Gaussian Markov random field on the auxiliary lattice, our model completely avoids the requirement of matrix inversion. It is remarkable that the computational complexity of our method is only O(n), where n is the number of observations. Hence, our method can be applied to very large datasets with reasonable computational (CPU) times. The numerical results indicate that our model can approximate Gaussian random fields very well in terms of predictions, even for those with long correlation lengths. For real data examples, our model can generally outperform conventional Gaussian random field models in both prediction errors and CPU times. Supplemental materials for the article are available online. © 2012 American Statistical Association, Institute of Mathematical Statistics, and Interface Foundation of North America.
Original languageEnglish (US)
Pages (from-to)453-475
Number of pages23
JournalJournal of Computational and Graphical Statistics
Volume21
Issue number2
DOIs
StatePublished - Jun 14 2012
Externally publishedYes

Fingerprint Dive into the research topics of 'Bayesian Analysis of Geostatistical Models With an Auxiliary Lattice'. Together they form a unique fingerprint.

Cite this