Capturing the time-varying drivers of an epidemic using stochastic dynamical systems.

Epidemics are often modeled using non-linear dynamical systems observed through partial and noisy data. In this paper, we consider stochastic extensions in order to capture unknown influences (changing behaviors, public interventions, seasonal effects, etc.). These models assign diffusion processes to the time-varying parameters, and our inferential procedure is based on a suitably adjusted adaptive particle Markov chain Monte Carlo algorithm. The performance of the proposed computational methods is validated on simulated data and the adopted model is applied to the 2009 H1N1 pandemic in England. In addition to estimating the effective contact rate trajectories, the methodology is applied in real time to provide evidence in related public health decisions. Diffusion-driven susceptible exposed infected retired-type models with age structure are also introduced.

[1]  B Cazelles,et al.  Using the Kalman filter and dynamic models to assess the changing HIV/AIDS epidemic. , 1997, Mathematical biosciences.

[2]  R. May,et al.  Infectious Diseases of Humans: Dynamics and Control , 1991, Annals of Internal Medicine.

[3]  A. Stuart,et al.  Parameter estimation for partially observed hypoelliptic diffusions , 2007, 0710.5442.

[4]  E L Ionides,et al.  Inference for nonlinear dynamical systems , 2006, Proceedings of the National Academy of Sciences.

[5]  Simon Cauchemez,et al.  Likelihood-based estimation of continuous-time epidemic models from time-series data: application to measles transmission in London , 2008, Journal of The Royal Society Interface.

[6]  Paul H. Garthwaite,et al.  Statistical methods for the prospective detection of infectious disease outbreaks: a review , 2012 .

[7]  S. Chib.  Likelihood based inference for diffusion driven state space models , 2006 .

[8]  O. Papaspiliopoulos,et al.  SMC^2: A sequential Monte Carlo algorithm with particle Markov chain Monte Carlo updates , 2011 .

[9]  G. Roberts,et al.  On inference for partially observed nonlinear diffusion models using the Metropolis–Hastings algorithm , 2001 .

[10]  Nicholas G. Polson,et al.  Tracking Flu Epidemics Using Google Flu Trends and Particle Learning , 2009 .

[11]  S. Cauchemez,et al.  Estimating in real time the efficacy of measures to control emerging communicable diseases. , 2006, American journal of epidemiology.

[12]  References , 1971 .

[13]  Charles J. Geyer,et al.  Practical Markov Chain Monte Carlo , 1992 .

[14]  C P Farrington,et al.  Infections with Varying Contact Rates: Application to Varicella , 2004, Biometrics.

[15]  David A. Rasmussen,et al.  Inference for Nonlinear Epidemiological Models Using Genealogies and Time Series , 2011, PLoS Comput. Biol..

[16]  A. Flahault,et al.  Estimating the impact of school closure on influenza transmission from Sentinel data , 2008, Nature.

[17]  Nicolas Chopin,et al.  SMC2: an efficient algorithm for sequential analysis of state space models , 2011, 1101.1528.

[18]  Guy Thomas,et al.  Temporal Variability and Social Heterogeneity in Disease Transmission: The Case of SARS in Hong Kong , 2009, PLoS Comput. Biol..

[19]  K. Kalogeropoulos.  Likelihood-based Inference for a Class of Multivariate Diffusions with Unobserved Paths , 2007 .

[20]  Troy Day,et al.  Mechanistic modelling of the three waves of the 1918 influenza pandemic , 2011, Theoretical Ecology.

[21]  C. Andrieu,et al.  The pseudo-marginal approach for efficient Monte Carlo computations , 2009, 0903.5480.

[22]  B. Øksendal.  Stochastic differential equations : an introduction with applications , 1987 .

[23]  Jeffrey Shaman,et al.  Absolute humidity modulates influenza survival, transmission, and seasonality , 2009, Proceedings of the National Academy of Sciences.

[24]  T. Kurtz.  Approximation of Population Processes , 1987 .

[25]  G. Wahba.  Spline models for observational data , 1990 .

[26]  A. Ghani,et al.  Joint estimation of the basic reproduction number and generation time parameters for infectious disease outbreaks. , 2011, Biostatistics.

[27]  Anthony J McMichael,et al.  Nonstationary Influence of El Niño on the Synchronous Dengue Epidemics in Thailand , 2005, PLoS medicine.

[28]  Jeremy Ginsberg,et al.  Detecting influenza epidemics using search engine query data , 2009, Nature.

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

[30]  William C Hahn,et al.  Oncogenic Transformation by Inhibitor-Sensitive and -Resistant EGFR Mutants , 2005, PLoS medicine.

[31]  A. Doucet,et al.  A Tutorial on Particle Filtering and Smoothing: Fifteen years later , 2008 .

[32]  S. Sarkka,et al.  Application of Girsanov Theorem to Particle Filtering of Discretely Observed Continuous - Time Non-Linear Systems , 2007, 0705.1598.

[33]  M. Jit,et al.  Vaccination against pandemic influenza A/H1N1v in England: a real-time economic evaluation. , 2010, Vaccine.

[34]  S Reza-Paul,et al.  Evaluating large-scale HIV prevention interventions: study design for an integrated mathematical modelling approach , 2007, Sexually Transmitted Infections.

[35]  Aaron A. King,et al.  Time series analysis via mechanistic models , 2008, 0802.0021.

[36]  P. Fearnhead,et al.  Exact and computationally efficient likelihood‐based estimation for discretely observed diffusion processes (with discussion) , 2006 .

[37]  Neil Ferguson,et al.  Capturing human behaviour , 2007, Nature.

[38]  N. Ferguson,et al.  Time lines of infection and disease in human influenza: a review of volunteer challenge studies. , 2008, American journal of epidemiology.

[39]  Darren J. Wilkinson,et al.  Bayesian inference for nonlinear multivariate diffusion models observed with error , 2008, Comput. Stat. Data Anal..

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

[41]  Gareth O. Roberts,et al.  Examples of Adaptive MCMC , 2009 .

[42]  P. Fine,et al.  Measles in England and Wales--I: An analysis of factors underlying seasonal patterns. , 1982, International journal of epidemiology.

[43]  A. Doucet,et al.  Particle Markov chain Monte Carlo methods , 2010 .

[44]  Hao Wang,et al.  Recovering the time-dependent transmission rate from infection data , 2009, 0907.3529.

[45]  S. J. Koopman.  Discussion of `Particle Markov chain Monte Carlo methods – C. Andrieu, A. Doucet and R. Holenstein’ [Review of: Particle Markov chain Monte Carlo methods] , 2010 .

[46]  J. Wallinga,et al.  Different Epidemic Curves for Severe Acute Respiratory Syndrome Reveal Similar Impacts of Control Measures , 2004, American journal of epidemiology.

[47]  N. Andrews,et al.  Incidence of 2009 pandemic influenza A H1N1 infection in England: a cross-sectional serological study , 2010, The Lancet.

[48]  M. Pitt,et al.  Likelihood based inference for diffusion driven models , 2004 .