返回
Constructing summary statistics for approximate Bayesian computation: semi-automatic approximate Bayesian computation
DOI:10.1111/j.1467-9868.2011.01010.x.png)
摘要
En 中文
. Many modern statistical applications involve inference for complex stochastic models, where it is easy to simulate from the models, but impossible to calculate likelihoods. Approximate Bayesian computation (ABC) is a method of inference for such models. It replaces calculation of the likelihood by a step which involves simulating artificial data for different parameter values, and comparing summary statistics of the simulated data with summary statistics of the observed data. Here we show how to construct appropriate summary statistics for ABC in a semi-automatic manner. We aim for summary statistics which will enable inference about certain parameters of interest to be as accurate as possible. Theoretical results show that optimal summary statistics are the posterior means of the parameters. Although these cannot be calculated analytically, we use an extra stage of simulation to estimate how the posterior means vary as a function of the data; and we then use these estimates of our summary statistics within ABC. Empirical results show that our approach is a robust method for choosing summary statistics that can result in substantially more accurate ABC analyses than the ad hoc choices of summary statistics that have been proposed in the literature. We also demonstrate advantages over two alternative methods of simulation-based inference.
Keyword:
Indirect inference
Likelihood-free inference
Markov chain Monte Carlo methods
Simulation
Stochastic kinetic networks
期刊
J
IF:
3.6
论文数:
1.5K
被引数:
3.2W
机构
引用论文
Likelihood-Free Inference of Population Structure and Local Adaptation in a Bayesian Hierarchical Model
GENETICS
IF5.1
Efficient Approximate Bayesian Computation Coupled With Markov Chain Monte Carlo Without Likelihood
GENETICS
IF5.1
Model criticism based on likelihood-free inference, with an application to protein network evolution基于无似然推理的模型批评及其在蛋白质网络进化中的应用

