Bayesian penalized spline models for the analysis of spatio‐temporal count data

In recent years, the availability of infectious disease counts in time and space has increased, and consequently, there has been renewed interest in model formulation for such data. In this paper, we describe a model that was motivated by the need to analyze hand, foot, and mouth disease surveillance data in China. The data are aggregated by geographical areas and by week, with the aims of the analysis being to gain insight into the space-time dynamics and to make short-term predictions, which will aid in the implementation of public health campaigns in those areas with a large predicted disease burden. The model we develop decomposes disease-risk into marginal spatial and temporal components and a space-time interaction piece. The latter is the crucial element, and we use a tensor product spline model with a Markov random field prior on the coefficients of the basis functions. The model can be formulated as a Gaussian Markov random field and so fast computation can be carried out using the integrated nested Laplace approximation approach. A simulation study shows that the model can pick up complex space-time structure and our analysis of hand, foot, and mouth disease data in the central north region of China provides new insights into the dynamics of the disease.

[1]  L. Fahrmeir,et al.  High dimensional structured additive regression models: Bayesian regularization, smoothing and predictive performance , 2011 .

[2]  Haavard Rue,et al.  Think continuous: Markovian Gaussian models in spatial statistics , 2011, 1110.6796.

[3]  Bradley P. Carlin,et al.  Hierarchical Spatio-Temporal Mapping of Disease Rates , 1997 .

[4]  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.

[5]  Ying C. MacNab,et al.  Spline smoothing in Bayesian disease mapping , 2007 .

[6]  Y. MacNab,et al.  Spatio‐temporal modelling of rates for the construction of disease maps , 2002, Statistics in medicine.

[7]  Y. Xia,et al.  Measles Metapopulation Dynamics: A Gravity Model for Epidemiological Coupling and Dynamics , 2004, The American Naturalist.

[8]  J. Wakefield,et al.  Bayesian inference for generalized linear mixed models. , 2010, Biostatistics.

[9]  A. Riebler,et al.  Bayesian bivariate meta‐analysis of diagnostic test studies using integrated nested Laplace approximations , 2010, Statistics in medicine.

[10]  J. Wakefield,et al.  Hand, Foot, and Mouth Disease in China: Patterns of Spread and Transmissibility , 2011, Epidemiology.

[11]  Leonhard Held,et al.  Modeling seasonality in space‐time infectious disease surveillance data , 2012, Biometrical journal. Biometrische Zeitschrift.

[12]  H. Rue,et al.  An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach , 2011 .

[13]  M. Höhle,et al.  Assessment of varicella vaccine effectiveness in Germany: a time-series approach , 2010, Epidemiology and Infection.

[14]  Dae-Jin Lee,et al.  P-spline ANOVA-type interaction models for spatio-temporal smoothing , 2011 .

[15]  María Dolores Ugarte,et al.  Spatio‐temporal modeling of mortality risks using penalized splines , 2009 .

[16]  M. D. Ugarte,et al.  On fitting spatio-temporal disease mapping models using approximate Bayesian inference , 2014, Statistical methods in medical research.

[17]  Sylvia Richardson,et al.  Bayesian spatio-temporal analysis of joint patterns of male and female lung cancer risks in Yorkshire (UK) , 2006, Statistical methods in medical research.

[18]  M. Cameletti,et al.  Spatial and Spatio-temporal Bayesian Models with R - INLA , 2015 .

[19]  D. Clayton,et al.  Bayesian analysis of space-time variation in disease risk. , 1995, Statistics in medicine.

[20]  Leonhard Held,et al.  Using integrated nested Laplace approximations for the evaluation of veterinary surveillance data from Switzerland: a case‐study , 2011 .

[21]  Leonhard Held,et al.  Posterior and Cross-validatory Predictive Checks: A Comparison of MCMC and INLA , 2010 .

[22]  L Held,et al.  Predictive assessment of a non‐linear random effects model for multivariate time series of infectious disease counts , 2011, Statistics in medicine.

[23]  H. Rue,et al.  Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations , 2009 .

[24]  Noel Cressie,et al.  Hierarchical statistical modelling of influenza epidemic dynamics in space and time , 2002, Statistics in medicine.

[25]  Haavard Rue,et al.  Estimation and extrapolation of time trends in registry data—Borrowing strength from related populations , 2011, 1108.0606.

[26]  Nicola G. Best,et al.  A shared component model for detecting joint and selective clustering of two diseases , 2001 .

[27]  L. Held,et al.  Multivariate modelling of infectious disease surveillance data , 2008, Statistics in medicine.

[28]  Leonhard Held,et al.  Spatio‐temporal disease mapping using INLA , 2011 .

[29]  C. McCulloch,et al.  Generalized Linear Mixed Models , 2005 .

[30]  L. Fahrmeir,et al.  Structured Additive Regression for Categorical Space–Time Data: A Mixed Model Approach , 2006 .

[31]  H. Rue,et al.  In order to make spatial statistics computationally feasible, we need to forget about the covariance function , 2012 .

[32]  Jon Wakefield,et al.  Disease mapping and spatial regression with count data. , 2007, Biostatistics.

[33]  A. Cook,et al.  Inference in Epidemic Models without Likelihoods , 2009 .

[34]  D. Ruppert.  Selecting the Number of Knots for Penalized Splines , 2002 .

[35]  Murali Haran,et al.  Emulating a gravity model to infer the spatiotemporal dynamics of an infectious disease , 2011, 1110.6451.

[36]  R. Assunção,et al.  Diffusion and prediction of Leishmaniasis in a large metropolitan area in Brazil with a Bayesian space–time model , 2001, Statistics in medicine.

[37]  L. Pettit.  The Conditional Predictive Ordinate for the Normal Distribution , 1990 .

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

[39]  B P Carlin,et al.  Spatio-temporal models with errors in covariates: mapping Ohio lung cancer mortality. , 1998, Statistics in medicine.

[40]  L. Held,et al.  Heterogeneity in vaccination coverage explains the size and occurrence of measles epidemics in German surveillance data , 2010, Epidemiology and Infection.

[41]  Sylvia Richardson,et al.  A hierarchical model for space–time surveillance data on meningococcal disease incidence , 2003 .

[42]  Leonhard Held,et al.  A statistical framework for the analysis of multivariate infectious disease surveillance counts , 2005 .

[43]  A. F. Militino,et al.  Evaluating space‐time models for short‐term cancer mortality risk predictions in small areas , 2014, Biometrical journal. Biometrische Zeitschrift.

[44]  C. Tong,et al.  Global epidemiology of Enterovirus 71 , 2009 .

[45]  L Knorr-Held,et al.  Bayesian modelling of inseparable space-time variation in disease risk. , 2000, Statistics in medicine.

[46]  L. Fahrmeir,et al.  A Mixed Model Approach for Geoadditive Hazard Regression , 2007 .

[47]  J. Besag,et al.  Bayesian image restoration, with two applications in spatial statistics , 1991 .

[48]  Gianluca Baio,et al.  Spatial and spatio-temporal models with R-INLA. , 2013, Spatial and spatio-temporal epidemiology.

[49]  Edward L. Ionides,et al.  Plug-and-play inference for disease dynamics: measles in large and small populations as a case study , 2009, Journal of The Royal Society Interface.

[50]  David Welch,et al.  Approximate Bayesian computation scheme for parameter inference and model selection in dynamical systems , 2009, Journal of The Royal Society Interface.

[51]  Ying C MacNab,et al.  Regression B‐spline smoothing in Bayesian disease mapping: with an application to patient safety surveillance , 2007, Statistics in medicine.