Approximate Bayesian computation for spatial SEIR(S) epidemic models.

Approximate Bayesia n Computation (ABC) provides an attractive approach to estimation in complex Bayesian inferential problems for which evaluation of the kernel of the posterior distribution is impossible or computationally expensive. These highly parallelizable techniques have been successfully applied to many fields, particularly in cases where more traditional approaches such as Markov chain Monte Carlo (MCMC) are impractical. In this work, we demonstrate the application of approximate Bayesian inference to spatially heterogeneous Susceptible-Exposed-Infectious-Removed (SEIR) stochastic epidemic models. These models have a tractable posterior distribution, however MCMC techniques nevertheless become computationally infeasible for moderately sized problems. We discuss the practical implementation of these techniques via the open source ABSEIR package for R. The performance of ABC relative to traditional MCMC methods in a small problem is explored under simulation, as well as in the spatially heterogeneous context of the 2014 epidemic of Chikungunya in the Americas.

[1]  G. Gibson,et al.  Estimation of multiple transmission rates for epidemics in heterogeneous populations , 2007, Proceedings of the National Academy of Sciences.

[2]  Mevin B. Hooten,et al.  Assessing North American influenza dynamics with a statistical SIRS model. , 2010, Spatial and spatio-temporal epidemiology.

[3]  S. Cassadou,et al.  Chikungunya in the Americas , 2014, The Lancet.

[4]  Aaron T Porter,et al.  A spatial epidemic model for disease spread over a heterogeneous spatial support. , 2016, Statistics in medicine.

[5]  Jacob J Oleson,et al.  A Path‐Specific SEIR Model for use with General Latent and Infectious Time Distributions , 2013, Biometrics.

[6]  M. Cetron,et al.  Assessing the Origin of and Potential for International Spread of Chikungunya Virus from the Caribbean , 2014, PLoS currents.

[7]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[8]  W. O. Kermack,et al.  A contribution to the mathematical theory of epidemics , 1927 .

[9]  Rob Deardon,et al.  INFERENCE FOR INDIVIDUAL-LEVEL MODELS OF INFECTIOUS DISEASES IN LARGE POPULATIONS. , 2010, Statistica Sinica.

[10]  Dao Nguyen,et al.  Statistical Inference for Partially Observed Markov Processes via the R Package pomp , 2015, 1509.00503.

[11]  Y Dumont,et al.  Mathematical studies on the sterile insect technique for the Chikungunya disease and Aedes albopictus , 2012, Journal of mathematical biology.

[12]  B. Finkenstädt,et al.  Statistical Inference in a Stochastic Epidemic SEIR Model with Control Intervention: Ebola as a Case Study , 2006, Biometrics.

[13]  Aaron T Porter,et al.  An empirically adjusted approach to reproductive number estimation for stochastic compartmental models: A case study of two Ebola outbreaks , 2016, Biometrics.

[14]  V. Colizza,et al.  Local and regional spread of chikungunya fever in the Americas. , 2014, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[15]  Peter Neal,et al.  Forward Simulation Markov Chain Monte Carlo with Applications to Stochastic Epidemic Models , 2015 .

[16]  D. Rubin Bayesianly Justifiable and Relevant Frequency Calculations for the Applied Statistician , 1984 .

[17]  L. Mowatt,et al.  Chikungunya in the Caribbean: An Epidemic in the Making , 2014, Infectious Diseases and Therapy.

[18]  G. Roberts,et al.  Predicting undetected infections during the 2007 foot-and-mouth disease outbreak , 2009, Journal of The Royal Society Interface.

[19]  Arnaud Doucet,et al.  An adaptive sequential Monte Carlo method for approximate Bayesian computation , 2011, Statistics and Computing.

[20]  Y. Dumont,et al.  Vector control for the Chikungunya disease. , 2010, Mathematical biosciences and engineering : MBE.

[21]  Pejman Rohani,et al.  Appropriate Models for the Management of Infectious Diseases , 2005, PLoS medicine.

[22]  Jakub Szuppe,et al.  Boost.Compute: A parallel computing library for C++ based on OpenCL , 2016, IWOCL.

[23]  Libo Sun,et al.  Parameter inference and model selection in deterministic and stochastic dynamical models via approximate Bayesian computation: modeling a wildlife epidemic , 2014, 1409.7715.

[24]  M. Beaumont Approximate Bayesian Computation in Evolution and Ecology , 2010 .

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

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