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.

QMC integration for lognormal-parametric, elliptic PDEs: local supports and product weights

Nonasymptotic Convergence Rate of Quasi-Monte Carlo: Applications to Linear Elliptic PDEs with Lognormal Coefficients and Importance Samplings

Circulant embedding with QMC -- analysis for elliptic PDE with lognormal coefficients

Higher order QMC Galerkin discretization for parametric operator equations

Density estimation for elliptic PDE with random input by preintegration and quasi-Monte Carlo methods

Quasi-Monte Carlo and discontinuous Galerkin

QMC integration based on arbitrary (t,m,s)-nets yields optimal convergence rates on several scales of function spaces

Sparse-grid polynomial interpolation approximation and integration for parametric and stochastic elliptic PDEs with lognormal inputs

Quasi-Monte Carlo for partial differential equations with generalized Gaussian input uncertainty

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

Liberating The Weights

A Universal Median Quasi-Monte Carlo Integration

Quasi-Monte Carlo for unbounded integrands with importance sampling

Quasi-Monte Carlo for Bayesian design of experiment problems governed by parametric PDEs

On the convergence of adaptive Galerkin FEM for parametric PDEs with lognormal coefficients

Finite Element Error Analysis of Elliptic PDEs with Random Coefficients and Its Application to Multilevel Monte Carlo Methods

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

A-posteriori QMC-FEM error estimation for Bayesian inversion and optimal control with entropic risk measure

Goal-oriented adaptive finite element multilevel Monte Carlo with convergence rates

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