Hierarchical Bayesian Quadrature
Not on the map. Why?
Abstract (arXiv)
Numerical integration is a cornerstone of various scientific computing applications, such as engineering simulations and model evidence computations in probabilistic machine learning. Bayesian Quadrature uses Gaussian process surrogates that explicitly encode structural assumptions about the integrand to obtain integral estimates with quantified uncertainty. These surrogates are predominantly based on stationary covariance functions, which results in model misspecification for integrands exhibiting nonstationary behavior. We tackle this issue through an adaptively growing, tree-based partition of the integration domain into local stationary models. Our method recombines the local integral estimates through a hierarchy of GP conditioning that reintroduces cross-subdomain correlations, while model selection criteria control the tree growth to avoid unnecessary partitioning. The resulting algorithm is simple, requires no MCMC, and adapts its evaluation budget to local integrand complexity. On benchmark integration problems and a model evidence computation for an epidemiological model, Hierarchical Bayesian Quadrature achieves substantial gains over standard Bayesian Quadrature on nonstationary integrands while matching its performance on stationary ones.
Key facts. No citing papers in the map yet (data through 14 Oct 2026). Builds on 18 of its 35 references that are in the map. First posted on arXiv 12 Jul 2026.
Builds on
18 of 35 references.
References · 18
- 2019Scalable Global Optimization via Local Bayesian Optimizationcited by 374
- 2015Distributed Gaussian Processescited by 84
- 2015Probabilistic Numerics and Uncertainty in Computationscited by 60
- 2022Physics-Informed Gaussian Process Regression Generalizes Linear PDE Solverscited by 41
- 2017Non-Stationary Spectral Kernelscited by 38
- 2014Sampling for Inference in Probabilistic Models with Fast Bayesian Quadraturecited by 31
- 2015Probabilistic Integration: A Role in Statistical Computation?cited by 26
- 2019Convergence Guarantees for Adaptive Bayesian Quadrature Methodscited by 22
- 2015Frank-Wolfe Bayesian Quadrature: Probabilistic Integration with Theoretical Guaranteescited by 20
- 2022Fast Bayesian Inference with Batch Bayesian Quadrature via Kernel Recombinationcited by 20
- 2016Nested Kriging predictions for datasets with large number of observationscited by 18
- 2017Convergence Analysis of Deterministic Kernel-Based Quadrature Rules in Misspecified Settingscited by 16
6 more
- 2017Patchwork Kriging for Large-scale Gaussian Process Regressioncited by 9
- 2025A Dictionary of Closed-Form Kernel Mean Embeddingscited by 9
- 2022Multilevel Bayesian Quadraturecited by 8
- 2020Bayesian Probabilistic Numerical Integration with Tree-Based Modelscited by 5
- 2023Easily Computed Marginal Likelihoods from Posterior Simulation Using the THAMES Estimatorcited by 1
- 2026Bayesian Quadraturecited by 1
Built on by
Not cited yet.