Abstract:Importance sampling (IS) is valuable in reducing the variance of Monte Carlo sam-pling for many areas, including finance, rare event simulation, and Bayesian inference. It is natural and obvious to combine quasi-Monte Carlo (QMC) methods with IS to achieve a faster rate of con-vergence. However, a naive replacement of Monte Carlo with QMC may not work well. This paper investigates the convergence rates of randomized QMC-based IS for estimating integrals with respect to a Gaussian measure, in which the IS measure is a Gaussian or t distribution. We prove that if the target function satisfies the so-called boundary growth condition and the covariance matrix of the IS density has eigenvalues no smaller than 1, then randomized QMC with the Gaussian proposal has a root mean squared error of O(N -1+\epsilon) for arbitrarily small \epsilon > 0. Similar results with a t distribution as the proposal are also established. These sufficient conditions help to assess the effectiveness of IS in QMC. For some particular applications, we find that the Laplace IS, a very general approach to approximate the target function by a quadratic Taylor approximation around its mode, has ei-genvalues smaller than 1, making the resulting integrand less favorable for QMC. From this point of view, when using Gaussian distributions as the IS proposal, a change of measure via Laplace IS may transform a favorable integrand into an unfavorable one for QMC although the variance of Monte Carlo sampling is reduced. We also study the effect of positivization tricks on the error rate when the integrand has mixed sign. If the smooth positivization proposed by Owen and Zhou (J. Amer. Statist. Assoc., 95 (2000), pp. 135--143) is used, the rate O(N -1+\epsilon) is retained. This is not the case if taking the positive and negative parts of the integrand. We also give some examples to verify our propositions and warn against naive replacement of Monte Carlo with QMC under IS proposals. Numerical results suggest that using Laplace IS with t distributions is more robust than that with Gaussian distributions.

A note on the accuracy of adaptive Gauss–Hermite quadrature

The Best Quadrature Formula Based on Hermite Information for the Class KW 2 [ A, B ]

A Review of Error Estimation in Adaptive Quadrature

Theoretical error estimates for computing the matrix logarithm by Padé-type approximants

How sharp are error bounds? --lower bounds on quadrature worst-case errors for analytic functions

How Sharp Are Error Bounds? –Lower Bounds on Quadrature Worst-Case Errors for Analytic Functions–

Monotonicity, bounds and extrapolation of Block-Gauss and Gauss-Radau quadrature for computing $B^T ϕ(A) B$

Error asymptotic expansion for the numerical approximation of logarithmic-kernel integral equations on closed curves

An analysis on stochastic Lanczos quadrature with asymmetric quadrature nodes

On Average-Case Error Bounds for Kernel-Based Bayesian Quadrature

A new class of quadrature rules for estimating the error in Gauss quadrature

Improving the Accuracy of Exponentially Converging Quadratures

ON THE ERROR RATE OF IMPORTANCE SAMPLING WITH RANDOMIZED QUASI-MONTE CARLO*

A posteriori error estimates for fourth-order elliptic problems

An Extension of the Exponential-type Error Bounds for Multiquadric and Gaussian Interpolations

Gauss quadrature for integrals and sums

Error bounds for Gauss–Lobatto quadrature of analytic functions on an ellipse

Accurately Solving the Poisson Equation by Combining Multiscale Radial Basis Functions and Gaussian Quadrature

On the one-point quadrature discretization in peridynamics: a novel perspective from Monte Carlo integration

Adaptive sampling-based quadrature rules for efficient Bayesian prediction

Some Benchmark Problems and Basic Ideas on the Accuracy of Boundary Element Analysis