Emulating a gravity model to infer the spatiotemporal dynamics of an infectious disease

• Activity Number: 401
• Type: Topic Contributed
• Date/Time: Tuesday, August 2, 2011 : 2:00 PM to 3:50 PM
• Sponsor: Section on Bayesian Statistical Science
• Author(s): Roman Jandarov*+ and Murali Haran and Ottar Bjornstad and Bryan Grenfell
• Companies: Penn State University and Penn State University and Penn State University and Princeton University
• Keywords: Gravity Model ; Disease Dynamics ; Measles ; Bayesian Inference ; Gaussian Processes ; SIR

Abstract

We study a metapopulation model for regional measles dynamics that uses a gravity coupling model and a time series susceptible-infected-recovered (T-SIR) model for local dynamics. Standard maximum likelihood or Bayesian inference for this model is infeasible as there are potentially tens of thousands of latent variables in the model and each evaluation of the likelihood is expensive. We develop an efficient discretized MCMC algorithm for Bayesian inference with these expensive likelihood evaluations. However, we find through a simulation study that parameter estimates are biased and simulations at the obtained parameter settings do not explain some important biological characteristics of the data. We propose fitting a Gaussian process (GP) model to forward simulations of the gravity model at a number of parameter settings. Treating this GP model as an approximation ('emulator') for the gravity model, we perform a full Bayesian analysis of a given data set. This approach allows us to conveniently study posterior distributions of the key parameters of the gravity model and has number of advantages over the classic likelihood based inference.