Incorporating Subsampling into Bayesian Models for High-Dimensional Spatial Data

by   Sudipto Saha, et al.

Additive spatial statistical models with weakly stationary process assumptions have become standard in spatial statistics. However, one disadvantage of such models is the computation time, which rapidly increases with the number of datapoints. The goal of this article is to apply an existing subsampling strategy to standard spatial additive models and to derive the spatial statistical properties. We call this strategy the “spatial data subset model” approach, which can be applied to big datasets in a computationally feasible way. Our approach has the advantage that one does not require any additional restrictive model assumptions. That is, computational gains increase as model assumptions are removed when using our model framework. This provides one solution to the computational bottlenecks that occur when applying methods such as Kriging to “big data”. We provide several properties of this new spatial data subset model approach in terms of moments, sill, nugget, and range under several sampling designs. The biggest advantage of our approach is that it is scalable to a dataset of any size that can be stored. We present the results of the spatial data subset model approach on simulated datasets, and on a large dataset consists of 150,000 observations of daytime land surface temperatures measured by the MODIS instrument onboard the Terra satellite.


page 13

page 16

page 21

page 40


What is the best predictor that you can compute in five minutes using a given Bayesian hierarchical model?

The goal of this paper is to provide a way for statisticians to answer t...

Kryging: Geostatistical analysis of large-scale datasets using Krylov subspace methods

Analyzing massive spatial datasets using Gaussian process model poses co...

Distributed model building and recursive integration for big spatial data modeling

Motivated by the important need for computationally tractable statistica...

Bayesian Inference for Big Spatial Data Using Non-stationary Spectral Simulation

It is increasingly understood that the assumption of stationarity is unr...

The Matérn Model: A Journey through Statistics, Numerical Analysis and Machine Learning

The Matérn model has been a cornerstone of spatial statistics for more t...

Comparison of new computational methods for geostatistical modelling of malaria

Geostatistical analysis of health data is increasingly used to model spa...

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

With continued advances in Geographic Information Systems and related co...

Please sign up or login with your details

Forgot password? Click here to reset