Abstract:We study a random sampling technique to approximate integrals $\int_{[0,1]^s}f(\mathbf{x})\,\mathrm{d}\mathbf{x}$ by averaging the function at some sampling points. We focus on cases where the integrand is smooth, which is a problem which occurs in statistics. The convergence rate of the approximation error depends on the smoothness of the function $f$ and the sampling technique. For instance, Monte Carlo (MC) sampling yields a convergence of the root mean square error (RMSE) of order $N^{-1/2}$ (where $N$ is the number of samples) for functions $f$ with finite variance. Randomized QMC (RQMC), a combination of MC and quasi-Monte Carlo (QMC), achieves a RMSE of order $N^{-3/2+\varepsilon}$ under the stronger assumption that the integrand has bounded variation. A combination of RQMC with local antithetic sampling achieves a convergence of the RMSE of order $N^{-3/2-1/s+\varepsilon}$ (where $s\ge1$ is the dimension) for functions with mixed partial derivatives up to order two. Additional smoothness of the integrand does not improve the rate of convergence of these algorithms in general. On the other hand, it is known that without additional smoothness of the integrand it is not possible to improve the convergence rate. This paper introduces a new RQMC algorithm, for which we prove that it achieves a convergence of the root mean square error (RMSE) of order $N^{-\alpha-1/2+\varepsilon}$ provided the integrand satisfies the strong assumption that it has square integrable partial mixed derivatives up to order $\alpha>1$ in each variable. Known lower bounds on the RMSE show that this rate of convergence cannot be improved in general for integrands with this smoothness. We provide numerical examples for which the RMSE converges approximately with order $N^{-5/2}$ and $N^{-7/2}$, in accordance with the theoretical upper bound.

Optimal convergence rates of MCMC integration for functions with unbounded second moment

Achieving High Convergence Rates by Quasi-Monte Carlo and Importance Sampling for Unbounded Integrands

Convergence Speed and Approximation Accuracy of Numerical MCMC

Monte Carlo convergence rates for kth moments in Banach spaces

Almost sure convergence rates of adaptive increasingly rare Markov chain Monte Carlo

Multilevel Monte Carlo methods for positivity-preserving approximations of the Heston 3/2-model

Monte Carlo convergence rates for $k$th moments in Banach spaces

Quasi-Monte Carlo for unbounded integrands with importance sampling

Speeding up Monte Carlo Integration: Control Neighbors for Optimal Convergence

Convergence of Uniformity Criteria and the Application in Numerical Integration

Multilevel Monte Carlo with Numerical Smoothing for Robust and Efficient Computation of Probabilities and Densities

Linear Monte Carlo quadrature with optimal confidence intervals

Antithetic multilevel Monte Carlo method for approximations of SDEs with non-globally Lipschitz continuous coefficients

On the Error Rate of Conditional Quasi-Monte Carlo for Discontinuous Functions.

A remark on the rates of convergence for integrated volatility estimation in the presence of jumps

On the Optimal Rate for the Convergence Problem in Mean Field Control

Higher order scrambled digital nets achieve the optimal rate of the root mean square error for smooth integrands

Asymptotic bias of inexact Markov Chain Monte Carlo methods in high dimension

An invitation to adaptive Markov chain Monte Carlo convergence theory

Contraction Rate Estimates of Stochastic Gradient Kinetic Langevin Integrators