# Frequentist coverage and sup-norm convergence rate in Gaussian process regression

arXiv: https://arxiv.org/abs/1708.04753 · HTML: https://mlsanity.com/paper/1708.04753/

## Key facts

- Cited by 17 papers in the map (data through 14 Oct 2026).
- First posted on arXiv 16 Aug 2017.

## Built on by (most cited first, 17 of 17)

- [Fast Statistical Leverage Score Approximation in Kernel Ridge Regression](https://mlsanity.com/paper/2103.05238/) (2021)
- [Improved Regret Bounds for Gaussian Process Upper Confidence Bound in Bayesian Optimization](https://mlsanity.com/paper/2506.01393/) (2025)
- [Likelihood-Free Frequentist Inference: Bridging Classical Statistics and Machine Learning for Reliable Simulator-Based Inference](https://mlsanity.com/paper/2107.03920/) (2021)
- [Empowering parameter-efficient transfer learning by recognizing the kernel structure in self-attention](https://mlsanity.com/paper/2205.03720/) (2022)
- [Sup-Norm Convergence of Deep Neural Network Estimator for Nonparametric Regression by Adversarial Training](https://mlsanity.com/paper/2307.04042/) (2023)
- [Deep Horseshoe Gaussian Processes](https://mlsanity.com/paper/2403.01737/) (2024)
- [Quasi-Bayesian Dual Instrumental Variable Regression](https://mlsanity.com/paper/2106.08750/) (2021)
- [Minimax Analysis for Inverse Risk in Nonparametric Planer Invertible Regression](https://mlsanity.com/paper/2112.00213/) (2021)
- [Memory Efficient And Minimax Distribution Estimation Under Wasserstein Distance Using Bayesian Histograms](https://mlsanity.com/paper/2307.10099/) (2023)
- [Pointwise uncertainty quantification for sparse variational Gaussian process regression with a Brownian motion prior](https://mlsanity.com/paper/2310.00097/) (2023)
- [Trustworthy scientific inference with generative models](https://mlsanity.com/paper/2508.02602/) (2025)
- [Optimal Confidence Band for Kernel Gradient Flow Estimator](https://mlsanity.com/paper/2605.05768/) (2026)
- [Uncertainty quantification for distributed regression](https://mlsanity.com/paper/2105.11425/) (2021)
- [Pure Differential Privacy for Functional Summaries with a Laplace-like Process](https://mlsanity.com/paper/2309.00125/) (2023)
- [Bayesian Optimization for Dynamic Pricing and Learning](https://mlsanity.com/paper/2510.12447/) (2025)
- [Data-driven Learning of Interaction Laws in Multispecies Particle Systems with Gaussian Processes: Convergence Theory and Applications](https://mlsanity.com/paper/2511.02053/) (2025)
- [Nonparametric Goodness-of-fit Testing under Covariate Shift](https://mlsanity.com/paper/2608.04860/) (2026)

## Abstract (arXiv)

Gaussian process (GP) regression is a powerful interpolation technique due to its flexibility in capturing non-linearity. In this paper, we provide a general framework for understanding the frequentist coverage of point-wise and simultaneous Bayesian credible sets in GP regression. As an intermediate result, we develop a Bernstein von-Mises type result under supremum norm in random design GP regression. Identifying both the mean and covariance function of the posterior distribution of the Gaussian process as regularized M-estimators, we show that the sampling distribution of the posterior mean function and the centered posterior distribution can be respectively approximated by two population level GPs. By developing a comparison inequality between two GPs, we provide exact characterization of frequentist coverage probabilities of Bayesian point-wise credible intervals and simultaneous credible bands of the regression function. Our results show that inference based on GP regression tends to be conservative; when the prior is under-smoothed, the resulting credible intervals and bands have minimax-optimal sizes, with their frequentist coverage converging to a non-degenerate value between their nominal level and one. As a byproduct of our theory, we show that the GP regression also yields minimax-optimal posterior contraction rate relative to the supremum norm, which provides a positive evidence to the long standing problem on optimal supremum norm contraction rate in GP regression.
