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.

Fabry-Perot Imaging Spectroscopy of Starburst and AGN Winds

Monte-Carlo Simulation of a Stochastic Differential Equation

A new algorithm for Monte Carlo simulation of Ising spin systems

Relieving the fermionic and the dynamical sign problem: Multilevel Blocking Monte Carlo simulations

Monte Carlo simulations for Ising spins with spin greater than 1/2 applied to the square and triangular lattices with antiferromagnetic interactions and comparing results using Kawasaki and Glauber dynamics

Complex-valued Bayesian parameter estimation via Markov chain Monte Carlo

Criticality of the random-site Ising model: Metropolis, Swendsen-Wang and Wolff Monte Carlo algorithms

Monte Carlo Simulation of Critical Dynamic Spin Systems

Monte Carlo Simulation of Anisotropic Ising Model Using Metropolis and Wolff Algorithm

Monte Carlo simulation for Barndorff-Nielsen and Shephard model under change of measure

Non-equilibrium quantum Monte Carlo algorithm for stabilizer Renyi entropy in spin systems

Post-measurement Quantum Monte Carlo

Applications of the Monte Carlo Method in Statistical Physics

Geometric allocation approaches in Markov chain Monte Carlo

Interacting Markov chain Monte Carlo methods for solving nonlinear measure-valued equations

Fluctuations of Interacting Markov Chain Monte Carlo Methods

Extreme Long-time Dynamic Monte Carlo Simulations

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

Combinatorial Dehn surgery on cubed and Haken 3-manifolds

Milstein schemes and antithetic multilevel Monte Carlo sampling for delay McKean-Vlasov equations and interacting particle systems

The three-dimensional randomly dilute Ising model: Monte Carlo results