A two-component model for counts of infectious diseases

Biostatistics. 2006 Jul;7(3):422-37. doi: 10.1093/biostatistics/kxj016. Epub 2006 Jan 11.

Abstract

We propose a stochastic model for the analysis of time series of disease counts as collected in typical surveillance systems on notifiable infectious diseases. The model is based on a Poisson or negative binomial observation model with two components: a parameter-driven component relates the disease incidence to latent parameters describing endemic seasonal patterns, which are typical for infectious disease surveillance data. An observation-driven or epidemic component is modeled with an autoregression on the number of cases at the previous time points. The autoregressive parameter is allowed to change over time according to a Bayesian changepoint model with unknown number of changepoints. Parameter estimates are obtained through the Bayesian model averaging using Markov chain Monte Carlo techniques. We illustrate our approach through analysis of simulated data and real notification data obtained from the German infectious disease surveillance system, administered by the Robert Koch Institute in Berlin. Software to fit the proposed model can be obtained from http://www.statistik.lmu.de/ approximately mhofmann/twins.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Algorithms
  • Bayes Theorem
  • Communicable Diseases / epidemiology*
  • Disease Outbreaks / statistics & numerical data
  • Germany / epidemiology
  • Hepatitis A / epidemiology*
  • Hepatitis B / epidemiology*
  • Hepatitis B / prevention & control
  • Humans
  • Markov Chains
  • Monte Carlo Method
  • Population Surveillance
  • Seasons
  • Stochastic Processes*
  • Vaccination / statistics & numerical data