Block Nearest Neighboor Gaussian processes for large datasets

08/18/2019
by   Zaida C. Quiroz, et al.
0

This work develops a valid spatial block-Nearest Neighbor Gaussian process (block-NNGP) for estimation and prediction of location-referenced large spatial datasets. The key idea behind our approach is to subdivide the spatial domain into several blocks which are dependent under some constraints. The cross-blocks capture the large-scale spatial variation, while each block capture the small-scale dependence. The block-NNGP is embeded as a sparsity-inducing prior within a hierarchical modeling framework. Markov chain Monte Carlo (MCMC) algorithms are executed without storing or decomposing large matrices, while the sparse block precision matrix is efficiently computed through parallel computing. We also consider alternate MCMC algorithms through composite sampling for faster computing time, and more reproducible Bayesian inference. The performance of the block-NNGP is illustrated using simulation studies and applications with massive real data, for locations in the order of 10^4.

READ FULL TEXT
POST COMMENT

Comments

There are no comments yet.

Authors

page 17

01/24/2020

R package for Nearest Neighbor Gaussian Process models

This paper describes and illustrates functionality of the spNNGP R packa...
02/26/2021

Sparse Cholesky matrices in spatial statistics

Gaussian Processes (GP) is a staple in the toolkit of a spatial statisti...
07/23/2019

Conjugate Nearest Neighbor Gaussian Process Models for Efficient Statistical Interpolation of Large Spatial Data

A key challenge in spatial statistics is the analysis for massive spatia...
10/02/2020

Improving performances of MCMC for Nearest Neighbor Gaussian Process models with full data augmentation

Even though Nearest Neighbor Gaussian Processes (NNGP) alleviate conside...
03/23/2020

High-dimensional multivariate Geostatistics: A Bayesian Matrix-Normal Approach

Joint modeling of spatially-oriented dependent variables are commonplace...
03/25/2020

Highly Scalable Bayesian Geostatistical Modeling via Meshed Gaussian Processes on Partitioned Domains

We introduce a class of scalable Bayesian hierarchical models for the an...
07/07/2020

Fast Bayesian Estimation of Spatial Count Data Models

Spatial count data models are used to explain and predict the frequency ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.