Multilevel rejection sampling for approximate Bayesian computation

Abstract

Likelihood-free methods, such as approximate Bayesian computation, are powerful tools for practical inference problems with intractable likelihood functions. Markov chain Monte Carlo and sequential Monte Carlo variants of approximate Bayesian computation can be effective techniques for sampling posterior distributions in an approximate Bayesian computation setting. However, without careful consideration of convergence criteria and selection of proposal kernels, such methods can lead to very biased inference or computationally inefficient sampling. In contrast, rejection sampling for approximate Bayesian computation, despite being computationally intensive, results in independent, identically distributed samples from the approximated posterior. An alternative method is proposed for the acceleration of likelihood-free Bayesian inference that applies multilevel Monte Carlo variance reduction techniques directly to rejection sampling. The resulting method retains the accuracy advantages of rejection sampling while significantly improving the computational efficiency.

References in corpus (3)

• Approximate Bayesian computation scheme for parameter inference and model selection in dynamical systems
• Optimal quantification of contact inhibition in cell populations
• Approximate Bayesian Computation by Modelling Summary Statistics in a Quasi-likelihood Framework

Cited by in corpus (8)

• Simulation and inference algorithms for stochastic biochemical reaction networks: from basic concepts to state-of-the-art
• Rapid Bayesian inference for expensive stochastic models
• Multifidelity Approximate Bayesian Computation
• Multifidelity multilevel Monte Carlo to accelerate approximate Bayesian parameter inference for partially observed stochastic processes
• Multifidelity Approximate Bayesian Computation with Sequential Monte Carlo Parameter Sampling
• Vector operations for accelerating expensive Bayesian computations -- a tutorial guide
• Importance sampling for a robust and efficient multilevel Monte Carlo estimator for stochastic reaction networks
• Multilevel Monte Carlo Variational Inference