Biological and Environmental Research - Earth and Environmental System Sciences
Earth and Environmental System Modeling

Bayesian Inference for High-Dimensional Nonstationary Gaussian Processes

TitleBayesian Inference for High-Dimensional Nonstationary Gaussian Processes
Publication TypeJournal Article
Year of Publication2020
JournalJournal of Statistical Computation and Simulation
Pages1-27
Abstract / Summary

In spite of the diverse literature on nonstationary spatial modelling and approximate Gaussian process (GP) methods, there are no general approaches for conducting fully Bayesian inference for moderately sized nonstationary spatial data sets on a personal laptop. For statisticians and data scientists who wish to conduct posterior inference and prediction with appropriate uncertainty quantification, the lack of such approaches and software is a limitation. Here, we develop methodology for implementing formal Bayesian inference for a general class of nonstationary GPs. Our novel approach uses pre-existing frameworks for characterizing nonstationarity in a new way while utilizing via modern GP likelihood approximations. Posterior sampling is implemented using flexible MCMC methods, with nonstationary posterior prediction conducted as a post-processing step. We demonstrate our novel methods on two data sets, ranging from several hundred to several thousand locations. All of our methods are implemented in the freely available BayesNSGP software package for R.

URLhttp://dx.doi.org/10.1080/00949655.2020.1792472
DOI10.1080/00949655.2020.1792472
Journal: Journal of Statistical Computation and Simulation
Year of Publication: 2020
Pages: 1-27
Publication Date: 07/2020

In spite of the diverse literature on nonstationary spatial modelling and approximate Gaussian process (GP) methods, there are no general approaches for conducting fully Bayesian inference for moderately sized nonstationary spatial data sets on a personal laptop. For statisticians and data scientists who wish to conduct posterior inference and prediction with appropriate uncertainty quantification, the lack of such approaches and software is a limitation. Here, we develop methodology for implementing formal Bayesian inference for a general class of nonstationary GPs. Our novel approach uses pre-existing frameworks for characterizing nonstationarity in a new way while utilizing via modern GP likelihood approximations. Posterior sampling is implemented using flexible MCMC methods, with nonstationary posterior prediction conducted as a post-processing step. We demonstrate our novel methods on two data sets, ranging from several hundred to several thousand locations. All of our methods are implemented in the freely available BayesNSGP software package for R.

DOI: 10.1080/00949655.2020.1792472
Citation:
Risser, M, and D Turek.  2020.  "Bayesian Inference for High-Dimensional Nonstationary Gaussian Processes."  Journal of Statistical Computation and Simulation 1-27.  https://doi.org/10.1080/00949655.2020.1792472.