On importance sampling and independent Metropolis-Hastings with an unbounded weight function

Importance sampling and independent Metropolis-Hastings (IMH) are among the fundamental building blocks of Monte Carlo methods. Both require a proposal distribution that globally approximates the target distribution. The Radon-Nikodym derivative of the target distribution relative to the proposal is called the weight function. Under the weak assumption that the weight is unbounded but has a number of finite moments under the proposal distribution, we obtain new results on the approximation error of importance sampling and of the particle independent Metropolis-Hastings algorithm (PIMH), which includes IMH as a special case. For IMH and PIMH, we show that the common random numbers coupling is maximal. Using that coupling we derive bounds on the total variation distance of a PIMH chain to the target distribution. The bounds are sharp with respect to the number of particles and the number of iterations. Our results allow a formal comparison of the finite-time biases of importance sampling and IMH. We further consider bias removal techniques using couplings of PIMH, and provide conditions under which the resulting unbiased estimators have finite moments. We compare the asymptotic efficiency of regular and unbiased importance sampling estimators as the number of particles goes to infinity.