π€ AI Summary
Existing Gaussian process regression (GPR) uncertainty quantification methods rely on fixed input locations, posterior variance scaling, or hyperparameter tuning, limiting their ability to characterize global extremal behavior and yielding poorly robust upper/lower bounds on unseen data. This work proposes the first chain-structured GPR uncertainty quantification framework that requires neither prespecified input points nor variance scaling. By leveraging kernel-specific analysis (e.g., RBF, MatΓ©rn) and partition-diameter-driven local geometric modeling, it delivers globally valid extremal bounds and adaptive local uncertainty measures. Theoretically, its bound tightness surpasses that of analytic relaxation approaches. Empirical evaluation on synthetic and real-world benchmarks demonstrates significant improvements over state-of-the-art baselines in bound tightness, robustness to distributional shift, and generalization across diverse tasks.
π Abstract
Gaussian process regression (GPR) is a popular nonparametric Bayesian method that provides predictive uncertainty estimates and is widely used in safety-critical applications. While prior research has introduced various uncertainty bounds, most existing approaches require access to specific input features and rely on posterior mean and variance estimates or tuning hyperparameters. These limitations hinder robustness and fail to capture the model's global behavior in expectation. To address these limitations, we propose a chaining-based framework for estimating upper and lower bounds on the expected extreme values over unseen data, without requiring access to specific input locations. We provide kernel-specific refinements for commonly used kernels such as RBF and Mat'ern, in which our bounds are tighter than generic constructions. We further improve numerical tightness by avoiding analytical relaxations. In addition to global estimation, we also develop a novel method for local uncertainty quantification at specified inputs. This approach leverages chaining geometry through partition diameters, adapting to local structure without relying on posterior variance scaling. Our experimental results validate the theoretical findings and demonstrate that our method outperforms existing approaches on both synthetic and real-world datasets.