Practical Bayesian Modeling and Inference for Massive Spatial Datasets On Modest Computing Environments

Abstract

With continued advances in Geographic Information Systems and related computational technologies, statisticians are often required to analyze very large spatial datasets. This has generated substantial interest over the last decade, already too vast to be summarized here, in scalable methodologies for analyzing large spatial datasets. Scalable spatial process models have been found especially attractive due to their richness and flexibility and, particularly so in the Bayesian paradigm, due to their presence in hierarchical model settings. However, the vast majority of research articles present in this domain have been geared toward innovative theory or more complex model development. Very limited attention has been accorded to approaches for easily implementable scalable hierarchical models for the practicing scientist or spatial analyst. This article is submitted to the Practice section of the journal with the aim of developing massively scalable Bayesian approaches that can rapidly deliver Bayesian inference on spatial process that are practically indistinguishable from inference obtained using more expensive alternatives. A key emphasis is on implementation within very standard (modest) computing environments (e.g., a standard desktop or laptop) using easily available statistical software packages without requiring message-parsing interfaces or parallel programming paradigms. Key insights are offered regarding assumptions and approximations concerning practical efficiency.

20 pages, 4 figures, 2 tables

References in corpus (4)

• A multi-resolution approximation for massive spatial datasets
• Permutation and Grouping Methods for Sharpening Gaussian Process Approximations
• Bayesian Nonstationary Spatial Modeling for Very Large Datasets
• A Case Study Competition Among Methods for Analyzing Large Spatial Data

Cited by in corpus (15)

• A review of predictive uncertainty estimation with machine learning
• Vecchia-Laplace approximations of generalized Gaussian processes for big non-Gaussian spatial data
• Bayesian spatio-temporal models for stream networks
• A Survey of Bayesian Statistical Approaches for Big Data
• Modeling Massive Spatial Datasets Using a Conjugate Bayesian Linear Regression Framework
• High-dimensional Multivariate Geostatistics: A Bayesian Matrix-Normal Approach
• Prior-preconditioned conjugate gradient method for accelerated Gibbs sampling in "large $n$ & large $p$" Bayesian sparse regression
• Fast increased fidelity approximate Gibbs samplers for Bayesian Gaussian process regression
• Deep Feature Gaussian Processes for Single-Scene Aerosol Optical Depth Reconstruction
• Conjugate Nearest Neighbor Gaussian Process Models for Efficient Statistical Interpolation of Large Spatial Data
• Combining Heterogeneous Spatial Datasets with Process-based Spatial Fusion Models: A Unifying Framework
• Spatial prediction of apartment rent using regression-based and machine learning-based approaches with a large dataset
• Applications of Conjugate Gradient in Bayesian computation
• A Note on the comparison of Nearest Neighbor Gaussian Process (NNGP) based models
• Sparse nearest neighbor Cholesky matrices in spatial statistics