Inference for discretely observed stochastic kinetic networks with applications to epidemic modeling.

We present a new method for Bayesian Markov Chain Monte Carlo-based inference in certain types of stochastic models, suitable for modeling noisy epidemic data. We apply the so-called uniformization representation of a Markov process, in order to efficiently generate appropriate conditional distributions in the Gibbs sampler algorithm. The approach is shown to work well in various data-poor settings, that is, when only partial information about the epidemic process is available, as illustrated on the synthetic data from SIR-type epidemics and the Center for Disease Control and Prevention data from the onset of the H1N1 pandemic in the United States.

[1]  Asger Hobolth,et al.  SIMULATION FROM ENDPOINT-CONDITIONED, CONTINUOUS-TIME MARKOV CHAINS ON A FINITE STATE SPACE, WITH APPLICATIONS TO MOLECULAR EVOLUTION. , 2009, The annals of applied statistics.

[2]  Hervé Philippe,et al.  Uniformization for sampling realizations of Markov processes: applications to Bayesian implementations of codon substitution models , 2008, Bioinform..

[3]  J. Robins,et al.  Transmission Dynamics and Control of Severe Acute Respiratory Syndrome , 2003, Science.

[4]  S. Cornell,et al.  Dynamics of the 2001 UK Foot and Mouth Epidemic: Stochastic Dispersal in a Heterogeneous Landscape , 2001, Science.

[5]  Philip D O'Neill,et al.  A tutorial introduction to Bayesian inference for stochastic epidemic models using Markov chain Monte Carlo methods. , 2002, Mathematical biosciences.

[6]  Matthias Cavassini,et al.  Molecular epidemiology reveals long-term changes in HIV type 1 subtype B transmission in Switzerland. , 2010, The Journal of infectious diseases.

[7]  Marc Lipsitch,et al.  Use of Cumulative Incidence of Novel Influenza A/H1N1 in Foreign Travelers to Estimate Lower Bounds on Cumulative Incidence in Mexico , 2009, PloS one.

[8]  Kiyosi Itô,et al.  Essentials of Stochastic Processes , 2006 .

[9]  R. Durrett.  Essentials of Stochastic Processes , 1999 .

[10]  Alessandro Vespignani,et al.  influenza A(H1N1): a Monte Carlo likelihood analysis based on , 2009 .

[11]  Valeria De Fonzo,et al.  Hidden Markov Models in Bioinformatics , 2007 .

[12]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[13]  Richard J Boys,et al.  Bayesian inference for stochastic epidemic models with time-inhomogeneous removal rates , 2007, Journal of mathematical biology.

[14]  T. Koski.  Hidden Markov Models for Bioinformatics , 2001 .

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

[16]  Simon Cauchemez,et al.  Assessing the severity of the novel influenza A/H1N1 pandemic , 2009, BMJ : British Medical Journal.

[17]  Tom Britton,et al.  Stochastic epidemics in dynamic populations: quasi-stationarity and extinction , 2000, Journal of mathematical biology.

[18]  R. Bellman,et al.  FUNCTIONAL APPROXIMATIONS AND DYNAMIC PROGRAMMING , 1959 .

[19]  P. Fearnhead,et al.  An exact Gibbs sampler for the Markov‐modulated Poisson process , 2006 .

[20]  Joseph Gani,et al.  Stochastic Epidemic Models and Their Statistical Analysis , 2002 .

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

[22]  Darren J. Wilkinson,et al.  Bayesian inference for a discretely observed stochastic kinetic model , 2008, Stat. Comput..

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

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