Approximate Bayesian inference for large spatial datasets using predictive process models

详细信息	查看全文 | 推荐本文 |

• 作者：Jo Eidsvik^a ; ^{joeid@math.ntnu.no} ; Andrew O. Finley^b ; Sudipto Banerjee^c ; Hå ; vard Rue^a
• 关键词：Approximate Bayesian inference ; Computational statistics ; Gaussian processes ; Geostatistics ; Laplace approximation ; Predictive process model
• 刊名：Computational Statistics and Data Analysis
• 出版年：2012
• 期刊代码：19_01679473
• 类别：cp
• 出版时间：June, 2012
• 卷：56
• 期：6
• 页码：1362-1380
• 文件大小：1186 K

摘要

The challenges of estimating hierarchical spatial models to large datasets are addressed. With the increasing availability of geocoded scientific data, hierarchical models involving spatial processes have become a popular method for carrying out spatial inference. Such models are customarily estimated using Markov chain Monte Carlo algorithms that, while immensely flexible, can become prohibitively expensive. In particular, fitting hierarchical spatial models often involves expensive decompositions of dense matrices whose computational complexity increases in cubic order with the number of spatial locations. Such matrix computations are required in each iteration of the Markov chain Monte Carlo algorithm, rendering them infeasible for large spatial datasets. The computational challenges in analyzing large spatial datasets are considered by merging two recent developments. First, the predictive process model is used as a reduced-rank spatial process, to diminish the dimensionality of the model. Then a computational framework is developed for estimating predictive process models using the integrated nested Laplace approximation. The settings where the first stage likelihood is Gaussian or non-Gaussian are discussed. Issues such as predictions and model comparisons are also discussed. Results are presented for synthetic data and several environmental datasets.