Home > Research > Publications & Outputs > On the efficiency of pseudo-marginal random wal...

Electronic data

  • AOS1278

    Final published version, 603 KB, PDF document

Links

Text available via DOI:

View graph of relations

On the efficiency of pseudo-marginal random walk Metropolis algorithms

Research output: Contribution to Journal/MagazineJournal articlepeer-review

Published
Close
<mark>Journal publication date</mark>01/2015
<mark>Journal</mark>Annals of Statistics
Issue number1
Volume43
Number of pages38
Pages (from-to)238-275
Publication StatusPublished
Early online date9/12/14
<mark>Original language</mark>English

Abstract

We examine the behaviour of the pseudo-marginal random walk Metropolis algorithm, where evaluations of the target density for the accept/reject probability are estimated rather than computed precisely. Under relatively general conditions on the target distribution, we obtain limiting formulae for the acceptance rate and for the expected squared jump distance, as the dimension of the target approaches infinity, under the assumption that the noise in the estimate of the log-target is additive and is independent of the position. For targets with independent and identically distributed components, we also obtain a limiting diffusion for the first component. We then consider the overall efficiency of the algorithm, in terms of both speed of mixing and computational time. Assuming the additive noise is Gaussian and is inversely proportional to the number of unbiased estimates that are used, we prove that the algorithm is optimally efficient when the variance of the noise is approximately 3.3 and the acceptance rate is approximately 7.0%. We also find that the optimal scaling is insensitive to the noise and that the optimal variance of the noise is insensitive to the scaling. The theory is illustrated with a simulation study using the particle random walk Metropolis.