Accelerating delayed-acceptance Markov chain Monte Carlo algorithms
Preprint, 2018

Delayed-acceptance Markov chain Monte Carlo (DA-MCMC) samples from a probability distribution, via a two-stages version of the Metropolis-Hastings algorithm, by combining the target distribution with a "surrogate" (i.e. an approximate and computationally cheaper version) of said distribution. DA-MCMC accelerates MCMC sampling in complex applications, while still targeting the exact distribution. We design a computationally faster DA-MCMC algorithm, which samples from an approximation of the target distribution. As a case study, we also introduce a novel stochastic differential equation model for protein folding data. We consider parameters inference in a Bayesian setting where a surrogate likelihood function is introduced in the delayed-acceptance scheme. In our applications we employ a Gaussian process as a surrogate likelihood, but other options are possible. In our accelerated algorithm the calculations in the "second stage" of the delayed-acceptance scheme are reordered in such as way that we can obtain a significant speed-up in the MCMC sampling, when the evaluation of the likelihood function is computationally intensive. We consider both simulations studies, and the analysis of real protein folding data. Simulation studies for the stochastic Ricker model and the novel stochastic differential equation model for protein-folding data, show that the speed-up is highly problem dependent. The more involved the computations of the likelihood function are, the higher the acceleration becomes when using our algorithm. Inference results for the standard delayed-acceptance algorithm and our approximated version are similar, indicating that our approximated algorithm can return reliable Bayesian inference.

Bayesian inference

pseudo marginal MCMC

protein folding

stochastic differential equation

Gaussian process

Författare

Samuel Wiqvist

Lunds universitet

Umberto Picchini

Göteborgs universitet

Chalmers, Matematiska vetenskaper, Tillämpad matematik och statistik

Julie Lyng Forman

University of Copenhagen

Ämneskategorier

Matematik

Biofysik

Fundament

Grundläggande vetenskaper

Mer information

Senast uppdaterat

2018-09-21