Fast Bayesian inference of Block Nearest Neighbor Gaussian process for large data

2019/08/18 by Quiroz, Zaida C., Prates, Marcos O., Dey, Dipak K. +1
#FOS: Computer and information sciences #Methodology (stat.ME)

paper · doi:10.48550/arxiv.1908.06437

Abstract

This paper presents the development of a spatial block-Nearest Neighbor Gaussian process (block-NNGP) for location-referenced large spatial data. The key idea behind this approach is to divide the spatial domain into several blocks which are dependent under some constraints. The cross-blocks capture the large-scale spatial dependence, while each block captures the small-scale spatial dependence. The resulting block-NNGP enjoys Markov properties reflected on its sparse precision matrix. It is embedded as a prior within the class of latent Gaussian models, thus Bayesian inference is obtained using the integrated nested Laplace approximation (INLA). The performance of the block-NNGP is illustrated on simulated examples and massive real data for locations in the order of 10⁴.

Related

• Conjugate Nearest Neighbor Gaussian Process Models for Efficient Statistical Interpolation of Large Spatial Data
• Hierarchical Nearest-Neighbor Gaussian Process Models for Large Geostatistical Datasets
• Efficient algorithms for Bayesian Nearest Neighbor Gaussian Processes
• High-Dimensional Bayesian Geostatistics
• INLA or MCMC? A Tutorial and Comparative Evaluation for Spatial Prediction in log-Gaussian Cox Processes