Scaling analysis of delayed rejection MCMC methods
نویسندگان
چکیده
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. AMS 2000 subject classifications: Primary 60F05; secondary 65C40.
منابع مشابه
R-package FME : MCMC tests
This vignette tests the Markov chain Monte Carlo (MCMC) implementation of Rpackage FME (Soetaert and Petzoldt 2010). It includes the delayed rejection and adaptive Metropolis algorithm (Haario, Laine, Mira, and Saksman 2006)
متن کاملParallel Chains, Delayed Rejection and Reversible Jump MCMC for Object Recognition
We tackle the problem of object recognition using a Bayesian approach. A marked point process [1] is used as a prior model for the (unknown number of) objects. A sample is generated via Markov chain Monte Carlo (MCMC) techniques using a novel combination of Metropolis-coupled MCMC (MCMCMC) [2] and the Delayed Rejection Algorithm (DRA) [4]. The method is illustrated on some synthetic data contai...
متن کاملStudying Atmospheric Chemistry with MCMC
A model describing chemical reactions in the stratosphere ([1]) is studied with MCMC methods. The model is a large ODE system consisting of 33 components, roughly 150 reactions and 150 reaction rate parameters. Thus, it is a good case study for adaptive MCMC methods designed for high-dimensional problems. In this case, the Delayed Rejection Adaptive Metropolis (DRAM, [2]) is succesfully applied...
متن کاملDRAM: Efficient adaptive MCMC
We propose to combine two quite powerful ideas that have recently appeared in the Markov chain Monte Carlo literature: adaptive Metropolis samplers and delayed rejection. The ergodicity of the resulting non–Markovian sampler is proved, and the efficiency of the combination is demonstrated with various examples. We present situations where the combination outperforms the original methods: adapta...
متن کاملScaling-up Split-Merge MCMC with Locality Sensitive Sampling (LSS)
Split-Merge MCMC (Monte Carlo Markov Chain) is one of the essential and popular variants of MCMC for problems when an MCMC state consists of an unknown number of components. It is well known that state-of-the-art methods for split-merge MCMC do not scale well. Strategies for rapid mixing requires smart and informative proposals to reduce the rejection rate. However, all known smart proposals in...
متن کامل