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.

Compressing spectral kernels in Gaussian Process: Enhanced generalization and interpretability

Compressible Spectral Mixture Kernels with Sparse Dependency Structures for Gaussian Processes

A Fast GP Regression Method Using Banded Sparsification of Inverse Covariance

Additive Kernels for Gaussian Process Modeling

Generalized Convolution Spectral Mixture for Multitask Gaussian Processes

A spatial Gaussian process method for hyperspectral remote sensing imagery classification

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

Spectral Mixture Kernels for Multi-Output Gaussian Processes

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

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

Gaussian processes with skewed Laplace spectral mixture kernels for long-term forecasting

On the Identifiability and Interpretability of Gaussian Process Models

When Gaussian Process Meets Big Data: A Review of Scalable GPs

Gaussian Process Regression with Grid Spectral Mixture Kernel: Distributed Learning for Multidimensional Data

Leveraging the Nugget Parameter for Efficient Gaussian Process Modeling

Thin and Deep Gaussian Processes

Composite Gaussian Processes: Scalable Computation and Performance Analysis

Representing Additive Gaussian Processes by Sparse Matrices

A Unifying Perspective on Non-Stationary Kernels for Deeper Gaussian Processes

Global Optimization of Gaussian processes