Bayesian parameter inference for stochastic SIR epidemic model with hyperbolic diffusion

ABSTRACT This paper is concerned with the Bayesian estimation parameters of the stochastic SIR (Susceptible-Infective-Removed) epidemic model from the trajectory data. Specifically, the data from the count of both infectives and susceptibles is assumed to be available on some time grid as the epidemic progresses. The diffusion approximation of the appropriate jump process is then used to estimate missing data between every pair of observation times. If the time step of imputations is small enough, we derive the posterior distributions of the infection and recovery rates using the Milstein scheme. The paper also presents Markov-chain Monte Carlo (MCMC) simulation that demonstrates that the method provides accurate estimates, as illustrated by the synthetic data from SIR epidemic model and the real data.

[1]  G. Roberts,et al.  Bayesian inference for partially observed stochastic epidemics , 1999 .

[2]  M. Abramowitz,et al.  Handbook of Mathematical Functions With Formulas, Graphs and Mathematical Tables (National Bureau of Standards Applied Mathematics Series No. 55) , 1965 .

[3]  A. Gut Probability: A Graduate Course , 2005 .

[4]  P. O’Neill,et al.  Bayesian inference for stochastic multitype epidemics in structured populations via random graphs , 2005 .

[5]  Alexander Grey,et al.  The Mathematical Theory of Infectious Diseases and Its Applications , 1977 .

[6]  Ola Elerian,et al.  A note on the existence of a closed form conditional transition density for the Milstein scheme , 1998 .

[7]  G. Roberts,et al.  Statistical inference and model selection for the 1861 Hagelloch measles epidemic. , 2004, Biostatistics.

[8]  P. O’Neill,et al.  Bayesian estimation of the basic reproduction number in stochastic epidemic models , 2008 .

[9]  Michael P H Stumpf,et al.  Sensitivity, robustness, and identifiability in stochastic chemical kinetics models , 2011, Proceedings of the National Academy of Sciences.

[10]  Andrew Papanicolaou Numerical Solution of Stochastic Differential Equations with Jumps in Finance , 2013 .

[11]  Philip D O'Neill,et al.  Bayesian inference for stochastic multitype epidemics in structured populations using sample data. , 2009, Biostatistics.

[12]  Boseung Choi,et al.  Inference for discretely observed stochastic kinetic networks with applications to epidemic modeling. , 2012, Biostatistics.

[13]  P. Fearnhead,et al.  Exact filtering for partially observed continuous time models , 2004 .

[14]  Bjørn Eraker MCMC Analysis of Diffusion Models With Application to Finance , 2001 .

[15]  Tewfik Kernane,et al.  Bayesian inference for nonlinear stochastic SIR epidemic model , 2016 .

[16]  A. Pedersen A new approach to maximum likelihood estimation for stochastic differential equations based on discrete observations , 1995 .

[17]  P. O’Neill,et al.  Bayesian inference for epidemics with two levels of mixing , 2005 .

[18]  Nicola Bruti-Liberati Numerical Solution of Stochastic Differential Equations with Jumps in Finance , 2010 .

[19]  A. G. Greenhill,et al.  Handbook of Mathematical Functions with Formulas, Graphs, , 1971 .

[20]  Gavin J. Gibson,et al.  Bayesian inference for stochastic epidemics in closed populations , 2004 .

[21]  Philip D. O'Neill,et al.  Computation of final outcome probabilities for the generalised stochastic epidemic , 2006, Stat. Comput..

[22]  H. Andersson,et al.  Stochastic Epidemic Models and Their Statistical Analysis , 2000 .

[23]  D. Wilkinson,et al.  Bayesian Inference for Stochastic Kinetic Models Using a Diffusion Approximation , 2005, Biometrics.

[24]  C. Viboud,et al.  A Bayesian MCMC approach to study transmission of influenza: application to household longitudinal data , 2004, Statistics in medicine.

[25]  Michael A. Gibson,et al.  Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels , 2000 .

[26]  N G Becker,et al.  On a general stochastic epidemic model. , 1977, Theoretical population biology.

[27]  Gavin J. Gibson,et al.  Estimating parameters in stochastic compartmental models using Markov chain methods , 1998 .