Open Access
June 2015 Approximate Bayesian Computation by Modelling Summary Statistics in a Quasi-likelihood Framework
Stefano Cabras, Maria Eugenia Castellanos Nueda, Erlis Ruli
Bayesian Anal. 10(2): 411-439 (June 2015). DOI: 10.1214/14-BA921

Abstract

Approximate Bayesian Computation (ABC) is a useful class of methods for Bayesian inference when the likelihood function is computationally intractable. In practice, the basic ABC algorithm may be inefficient in the presence of discrepancy between prior and posterior. Therefore, more elaborate methods, such as ABC with the Markov chain Monte Carlo algorithm (ABC-MCMC), should be used. However, the elaboration of a proposal density for MCMC is a sensitive issue and very difficult in the ABC setting, where the likelihood is intractable. We discuss an automatic proposal distribution useful for ABC-MCMC algorithms. This proposal is inspired by the theory of quasi-likelihood (QL) functions and is obtained by modelling the distribution of the summary statistics as a function of the parameters. Essentially, given a real-valued vector of summary statistics, we reparametrize the model by means of a regression function of the statistics on parameters, obtained by sampling from the original model in a pilot-run simulation study. The QL theory is well established for a scalar parameter, and it is shown that when the conditional variance of the summary statistic is assumed constant, the QL has a closed-form normal density. This idea of constructing proposal distributions is extended to non constant variance and to real-valued parameter vectors. The method is illustrated by several examples and by an application to a real problem in population genetics.

Citation

Download Citation

Stefano Cabras. Maria Eugenia Castellanos Nueda. Erlis Ruli. "Approximate Bayesian Computation by Modelling Summary Statistics in a Quasi-likelihood Framework." Bayesian Anal. 10 (2) 411 - 439, June 2015. https://doi.org/10.1214/14-BA921

Information

Published: June 2015
First available in Project Euclid: 2 February 2015

zbMATH: 1335.62041
MathSciNet: MR3420888
Digital Object Identifier: 10.1214/14-BA921

Keywords: estimating function , likelihood-free methods , Markov chain Monte Carlo , proposal distribution , pseudo-likelihood

Rights: Copyright © 2015 International Society for Bayesian Analysis

Vol.10 • No. 2 • June 2015
Back to Top