Asymptotic Bounds for Smoothness Parameter Estimates in Gaussian Process Interpolation

From MaRDI portal
Revision as of 06:49, 10 July 2024 by Import240710060729 (talk | contribs) (Created automatically from import240710060729)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

Publication:6062242

DOI10.1137/22M149288XzbMath1525.60045arXiv2203.05400OpenAlexW4389047038MaRDI QIDQ6062242

Toni Karvonen

Publication date: 30 November 2023

Published in: SIAM/ASA Journal on Uncertainty Quantification (Search for Journal in Brave)

Abstract: It is common to model a deterministic response function, such as the output of a computer experiment, as a Gaussian process with a Mat'ern covariance kernel. The smoothness parameter of a Mat'ern kernel determines many important properties of the model in the large data limit, including the rate of convergence of the conditional mean to the response function. We prove that the maximum likelihood estimate of the smoothness parameter cannot asymptotically undersmooth the truth when the data are obtained on a fixed bounded subset of $mathbb{R}^d$. That is, if the data-generating response function has Sobolev smoothness $

u_0 + d/2$, then the smoothness parameter estimate cannot be asymptotically less than $ u_0 + d/2$. The lower bound is sharp. Additionally, we show that maximum likelihood estimation finds the "correct" smoothness for a class of compactly supported self-similar functions. We also consider cross-validation and prove an asymptotic lower bound $

u_0$, which however is unlikely to be sharp. The results are based on approximation theory in Sobolev spaces and some general theorems that restrict the set of values that the parameter estimators can take.


Full work available at URL: https://arxiv.org/abs/2203.05400





Cites Work






This page was built for publication: Asymptotic Bounds for Smoothness Parameter Estimates in Gaussian Process Interpolation