Approximate Bayesian Computation Algorithms for Estimating Network Model Parameters

Studies on Approximate Bayesian Computation (ABC) replacing the intractable likelihood function in evaluation of the posterior distribution have been developed for several years. However, their field of application has to date essentially been limited to inference in population genetics. Here, we propose to extend this approach to estimating the structure of transmission networks of viruses in human populations. In particular, we are interested in estimating the transmission parameters under four very general network structures: random, Watts-Strogatz, Barabasi-Albert and an extension that incorporates aging. Estimation was evaluated under three approaches, based on ABC, ABC-Markov chain Monte Carlo (ABC-MCMC) and ABC-Sequential Monte Carlo (ABC-SMC) samplers. We show that ABC-SMC samplers outperform both ABC and ABC-MCMC, achieving high accuracy and low variance in simulations. This approach paves the way to estimating parameters of real transmission networks of transmissible diseases.

[1]  David Welch,et al.  Approximate Bayesian computation scheme for parameter inference and model selection in dynamical systems , 2009, Journal of The Royal Society Interface.

[2]  O. François,et al.  Approximate Bayesian Computation (ABC) in practice. , 2010, Trends in ecology & evolution.

[3]  Alan M. Frieze,et al.  Random graphs , 2006, SODA '06.

[4]  Duncan J. Watts,et al.  Collective dynamics of ‘small-world’ networks , 1998, Nature.

[5]  Albert,et al.  Emergence of scaling in random networks , 1999, Science.

[6]  M. Kenward,et al.  An Introduction to the Bootstrap , 2007 .

[7]  Gábor Csárdi,et al.  The igraph software package for complex network research , 2006 .

[8]  Paul Marjoram,et al.  Markov chain Monte Carlo without likelihoods , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[9]  Tim Hesterberg,et al.  Monte Carlo Strategies in Scientific Computing , 2002, Technometrics.

[10]  Daniel Wegmann,et al.  Bayesian Computation and Model Selection Without Likelihoods , 2010, Genetics.

[11]  D. Balding,et al.  Approximate Bayesian computation in population genetics. , 2002, Genetics.

[12]  M. Gutmann,et al.  Approximate Bayesian Computation , 2019, Annual Review of Statistics and Its Application.

[13]  Mark M. Tanaka,et al.  Sequential Monte Carlo without likelihoods , 2007, Proceedings of the National Academy of Sciences.

[14]  Olivier François,et al.  Non-linear regression models for Approximate Bayesian Computation , 2008, Stat. Comput..

[15]  M. Feldman,et al.  Population growth of human Y chromosomes: a study of Y chromosome microsatellites. , 1999, Molecular biology and evolution.

[16]  ANTHONY WREN,et al.  Population Growth , 1972, Nature.

[17]  P. Bickel,et al.  A nonparametric view of network models and Newman–Girvan and other modularities , 2009, Proceedings of the National Academy of Sciences.

[18]  H. Haario,et al.  An adaptive Metropolis algorithm , 2001 .