Abstract:The Gaussian process (GP) is a widely used probabilistic machine learning method for stochastic function approximation, stochastic modeling, and analyzing real-world measurements of nonlinear processes. Unlike many other machine learning methods, GPs include an implicit characterization of uncertainty, making them extremely useful across many areas of science, technology, and engineering. Traditional implementations of GPs involve stationary kernels (also termed covariance functions) that limit their flexibility and exact methods for inference that prevent application to data sets with more than about ten thousand points. Modern approaches to address stationarity assumptions generally fail to accommodate large data sets, while all attempts to address scalability focus on approximating the Gaussian likelihood, which can involve subjectivity and lead to inaccuracies. In this work, we explicitly derive an alternative kernel that can discover and encode both sparsity and nonstationarity. We embed the kernel within a fully Bayesian GP model and leverage high-performance computing resources to enable the analysis of massive data sets. We demonstrate the favorable performance of our novel kernel relative to existing exact and approximate GP methods across a variety of synthetic data examples. Furthermore, we conduct space-time prediction based on more than one million measurements of daily maximum temperature and verify that our results outperform state-of-the-art methods in the Earth sciences. More broadly, having access to exact GPs that use ultra-scalable, sparsity-discovering, nonstationary kernels allows GP methods to truly compete with a wide variety of machine learning methods.

Sparse gaussian processes using backward elimination

A Fast GP Regression Method Using Banded Sparsification of Inverse Covariance

Recursive estimation for sparse Gaussian process regression

Representing Additive Gaussian Processes by Sparse Matrices

EigenGP: Gaussian Processes with Sparse Data-Dependent Eigenfunctions

Adaptive Sparse Gaussian Process

EigenGP: Gaussian Process Models with Adaptive Eigenfunctions.

Sparse Gaussian processes for solving nonlinear PDEs

Sparse Gaussian Process Based on Hat Basis Functions

Greedy Forward Selection Algorithms to Sparse Gaussian Process Regression

Kernel Multigrid: Accelerate Back-fitting via Sparse Gaussian Process Regression

Sparsity-Aware Distributed Learning for Gaussian Processes with Linear Multiple Kernel

Sparse Gaussian Process regression model based on ℓ 1/2 regularization

Efficient Gaussian process based on BFGS updating and logdet approximation

Nonparametric Bayesian Mixed-effect Model: a Sparse Gaussian Process Approach

Learning robot inverse dynamics using sparse online Gaussian process with forgetting mechanism

Variable Sigma Gaussian Processes: an Expectation Propagation Perspective

Compactly-supported nonstationary kernels for computing exact Gaussian processes on big data

Learning gradients with gaussian processes

Scaling Gaussian Process Regression with Derivatives

Sparse Gaussian Processes with Manifold-Preserving Graph Reduction