Skip to main content
Log in

Scaling Analysis of Delayed Rejection MCMC Methods

  • Published:
Methodology and Computing in Applied Probability Aims and scope Submit manuscript

Abstract

In this paper, we study the asymptotic efficiency of the delayed rejection strategy. In particular, the efficiency of the delayed rejection Metropolis–Hastings algorithm is compared to that of the regular Metropolis algorithm. To allow for a fair comparison, the study is carried under optimal mixing conditions for each of these algorithms. After introducing optimal scaling results for the delayed rejection (DR) algorithm, we outline the fact that the second proposal after the first rejection is discarded, with a probability tending to 1 as the dimension of the target density increases. To overcome this drawback, a modification of the delayed rejection algorithm is proposed, in which the direction of the different proposals is fixed once for all, and the Metropolis–Hastings accept-reject mechanism is used to select a proper scaling along the search direction. It is shown that this strategy significantly outperforms the original DR and Metropolis algorithms, especially when the dimension becomes large. We include numerical studies to validate these conclusions.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Similar content being viewed by others

References

  • Bédard M (2007) Weak convergence of Metropolis algorithms for non-i.i.d. target distributions. Ann Appl Probab 17:1222–1244

    Article  MathSciNet  MATH  Google Scholar 

  • Bédard M, Rosenthal JS (2008) Optimal scaling of Metropolis algorithms: heading toward general target distributions. Can J Stat 36:483–503

    Article  MATH  Google Scholar 

  • Bédard M, Douc R, Moulines E (2012) Scaling analysis of multiple-try MCMC methods. Stoch Process Their Appl 122:758–786

    Article  MATH  Google Scholar 

  • Beskos A, Roberts G, Stuart A (2009) Optimal scalings for local Metropolis–Hastings chains on nonproduct targets in high dimensions. Ann Appl Probab 19:863–898

    Article  MathSciNet  MATH  Google Scholar 

  • Billingsley P (1999) Convergence of probability measures, 2nd edn. Wiley Series in Probability and Statistics: Probability and Statistics. Wiley, New York. A Wiley-Interscience Publication

    Book  MATH  Google Scholar 

  • Craiu RV, Lemieux C (2007) Acceleration of the multiple-try Metropolis algorithm using antithetic and stratified sampling. Stat Comput 17:109–120

    Article  MathSciNet  Google Scholar 

  • Gelfand SB, Mitter S (1991) Weak convergence of Markov chain sampling methods and annealing algorithms to diffusions. J Optim Theory Appl 68:483–498

    Article  MathSciNet  MATH  Google Scholar 

  • Green PJ, Mira A (2001) Delayed rejection in reversible jump Metropolis–Hastings. Biometrika 88:1035–1053

    Article  MathSciNet  MATH  Google Scholar 

  • Haario H, Laine M, Mira A, Saksman E (2006) DRAM: efficient adaptive MCMC. Stat Comput 16:339–354

    Article  MathSciNet  Google Scholar 

  • Harkness MA, Green PJ (2000) Parallel chains, delayed rejection and reversible jump MCMC for object recognition. In: British machine vision conference

  • Mattingly J, Pillai N, Stuart A (2012) Diffusion limits of random walk Metropolis in high dimensions. Ann Appl Probab 22:881–930

    Article  MathSciNet  MATH  Google Scholar 

  • Mira A (2001a) On Metropolis–Hastings algorithms with delayed rejection. Metron LIX:231–241

    MathSciNet  Google Scholar 

  • Mira A (2001b) Ordering and improving the performance of Monte Carlo Markov chains. Stat Sci 16:340–350

    Article  MathSciNet  MATH  Google Scholar 

  • Raggi D (2005) Adaptive MCMC for inference on affine stochastic volatility models with jumps. Econom J 8:235–250

    Article  MathSciNet  MATH  Google Scholar 

  • Robert CP, Casella G (2004) Monte Carlo statistical methods, 2nd edn. Springer, New York

    Book  MATH  Google Scholar 

  • Roberts GO, Rosenthal JS (1998) Markov chain Monte Carlo: some practical implications of theoretical results. Can J Stat 26:5–32

    Article  MathSciNet  MATH  Google Scholar 

  • Roberts GO, Rosenthal JS (2001) Optimal scaling for various Metropolis–Hastings algorithms. Stat Sci 16:351–367

    Article  MathSciNet  MATH  Google Scholar 

  • Roberts GO, Gelman A, Gilks WR (1997) Weak convergence and optimal scaling of random walk Metropolis algorithms. Ann Appl Probab 7:110–120

    Article  MathSciNet  MATH  Google Scholar 

  • Tierney LL, Mira AA (1999) Some adaptive Monte Carlo methods for Bayesian inference. Stat Med 18:2507–2515

    Article  Google Scholar 

  • Trias M, Vecchio A, Veitch J (2009) Delayed rejection schemes for efficient Markov-Chain Monte-Carlo sampling of multimodal distributions. ArXiv e-prints

  • Umstätter R, Meyer R, Dupuis R, Veitch J, Woan G, Christensen N (2004) Estimating the parameters of gravitational waves from neutron stars using an adaptive MCMC method. Class Quantum Gravity 21:1655–1675

    Article  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Mylène Bédard.

Additional information

This work of the first author is supported by the National Sciences and Engineering Research Council of Canada.

This work of the second and third authors is supported by the Agence Nationale de la Recherche (ANR, 212, rue de Bercy 75012 Paris) through the 2009-2012 project Big MC.

Rights and permissions

Reprints and permissions

About this article

Cite this article

Bédard, M., Douc, R. & Moulines, E. Scaling Analysis of Delayed Rejection MCMC Methods. Methodol Comput Appl Probab 16, 811–838 (2014). https://doi.org/10.1007/s11009-013-9326-y

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s11009-013-9326-y

Keywords

AMS 2000 Subject Classifications

Navigation