Median randomized quasi-Monte Carlo for nested integration
Abstract
Nested integrals ∫ f(∫ g(y,x)dx)dy where y∈[0,1]d1 and x∈[0,1]d2 frequently arise in Bayesian experimental design, medical decision making, and computational finance. Statistical estimators for nested integrals are more computationally demanding than for standard integrals. Recent studies demonstrate that the median randomized quasi-Monte Carlo (mRQMC) method achieves super-polynomial convergence rates for analytic integrands. We extend mRQMC methods to nested integrations and propose a novel double-loop mRQMC estimator with super-polynomial convergence rate if both the inner and outer integrands are analytic. Established estimators for nested integrals, such as double-loop or multilevel Monte Carlo, are by contrast limited to algebraic convergence rates. However, expected information gain (EIG) applications from experimental design typically do not fulfill this condition due to a lack of regularity in the outer integrand. We address this through a theoretically-grounded importance sampling regularization scheme. Although super-polynomial convergence rates remain achievable for the inner integral, only the algebraic rate α+3/2−ε for any ε>0 can be theoretically guaranteed for the outer integral in the EIG, where α∈ℕ∪{0} denotes the order of the mixed derivatives of the regularized outer integrand. Our approach shows substantial improvements over double-loop randomized quasi-Monte Carlo methods based on the sample mean in moderately high-dimensional settings for a simple polynomial example. Moreover, we demonstrate improved convergence behavior of the median-based estimator for the more challenging setting of EIG estimation.
Related
- Optimal Monte Carlo integration on matrix manifolds
- Slice Monte Carlo Integration
- Quasi-Monte Carlo based meshfree methods for parabolic equations on parameterized Riemannian manifolds
- Faster Hamiltonian Monte Carlo by learning leapfrog scale: an offline randomized solution
- Beyond Exchangeability: Distribution-Shift-Aware Integration of External Control Data in Randomized Trials
BibTeX
@misc{bartuska2026median,
title = {{Median randomized quasi-Monte Carlo for nested integration}},
author = {Arved Bartuska and Yang Liu and Raúl Tempone},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.12643},
eprint = {2610.12643},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.12643},
}
RIS
TY - UNPB
TI - Median randomized quasi-Monte Carlo for nested integration
AU - Arved Bartuska
AU - Yang Liu
AU - Raúl Tempone
PY - 2026
DA - 2026/10/08
PB - arXiv
DO - 10.48550/arxiv.2610.12643
UR - https://arxiv.org/abs/2610.12643
AB - Nested integrals ∫ f(∫ g(y,x)dx)dy where y∈[0,1]^d_1 and x∈[0,1]^d_2 frequently arise in Bayesian experimental design, medical decision making, and computational finance. Statistical estimators for nested integrals are more computationally demanding than for standard integrals. Recent studies demonstrate that the median randomized quasi-Monte Carlo (mRQMC) method achieves super-polynomial convergence rates for analytic integrands. We extend mRQMC methods to nested integrations and propose a novel double-loop mRQMC estimator with super-polynomial convergence rate if both the inner and outer integrands are analytic. Established estimators for nested integrals, such as double-loop or multilevel Monte Carlo, are by contrast limited to algebraic convergence rates. However, expected information gain (EIG) applications from experimental design typically do not fulfill this condition due to a lack of regularity in the outer integrand. We address this through a theoretically-grounded importance sampling regularization scheme. Although super-polynomial convergence rates remain achievable for the inner integral, only the algebraic rate α+3/2−ε for any ε>0 can be theoretically guaranteed for the outer integral in the EIG, where α∈ℕ∪{0} denotes the order of the mixed derivatives of the regularized outer integrand. Our approach shows substantial improvements over double-loop randomized quasi-Monte Carlo methods based on the sample mean in moderately high-dimensional settings for a simple polynomial example. Moreover, we demonstrate improved convergence behavior of the median-based estimator for the more challenging setting of EIG estimation.
ER -