surveillance (version 1.5-4)

hagelloch: 1861 measles epidemic in the city of Hagelloch, Germany

Description

Data on the 188 cases in the measles outbreak among children in the German city of Hagelloch (near T�{ue}bingen) 1861. The data were originally collected by Dr. Albert Pfeilsticker and re-analysed by Dr. Heike Oesterle.

Usage

data("hagelloch")

encoding

latin1

source

Thanks to Peter J. Neal, University of Manchester, for providing us with these data, which he again became from Niels Becker, Australian National University. To cite the data, the main references are Pfeilsticker (1863) and Oesterle (1992).

References

Pfeilsticker, A. (1863). Beitr�{ae}ge zur Pathologie der Masern mit besonderer Ber�{ue}cksichtgung der statistischen Verh�{ae}ltnisse, M.D. Thesis, Eberhard-Karls Universit�{ae}t, T�{ue}bingen. Available as http://www.archive.org/details/beitrgezurpatho00pfeigoog. Oesterle, H. (1992). Statistische Reanalyse einer Masernepidemie 1861 in Hagelloch, M.D. Thesis, Eberhard-Karls Universit�{ae}at, T�{ue}bingen.

Neal, P. J. and Roberts, G. O (2004). Statistical inference and model selection for the 1861 Hagelloch measles epidemic, Biostatistics 5(2):249-261

See Also

twinSIR, epidata

Examples

data("hagelloch")
  #Data frame with the original data documented in Oesterle (1992)
  head(hagelloch.df)
  #Represented as epidata object
  (s <- summary(hagelloch))

  ### Use data.frame to illustrate the basics

  #Show case locations as in Neal & Roberts (different scaling) using
  #the data.frame (promoted to a SpatialPointsDataFrame)
  coordinates(hagelloch.df) <- c("HNX","HNY")
  plot(hagelloch.df,xlab="x [m]",ylab="x [m]",pch=15,axes=TRUE,
       cex=sqrt(multiplicity(hagelloch.df)))

  #Epicurve
  hist(as.numeric(hagelloch.df$tS),xlab="Time (days)",ylab="Cases",main="")

  ### SIR model information for population & individuals

  plot(s, col=c("green","red","darkgray"))
  stateplot(s, id=c("187"))

  #Show a dynamic illustration of the spread of the infection
    animate(hagelloch,time.spacing=0.1,legend.opts=list(x="topleft"),sleep=1/100)