Description Usage Format Source Examples
Weekly number of reported campylobacteriosis cases in Germany 2002-2011 together with the corresponding absolute humidity (in g/m^3) that week. The absolute humidity was computed according to the procedure by Dengler (1997) using the means of representative weather station data from the German Climate service.
1 |
A data.frame
containing the following columns
date
Date
instance containing the monday of the
reporting week.
case
Number of reported cases that week.
state
Boolean indicating whether there is external knowledge about an outbreak that week
hum
Mean absolute humidity (in g/m^3) of that week as measured by a single representative weather station.
l1.hum
-l5.hum
Lagged version (lagged by 1-5) of
the hum
covariate.
Boolean indicating whether the reporting week corresponds to the first two weeks of the year (TRUE) or not (FALSE). Note: The first week of a year is here defined as the first reporting week, which has its corresponding monday within new year.
Boolean indicating whether the reporting week
corresponds to the last two weeks of the year (TRUE) or not
(FALSE). Note: This are the first two weeks before the
newyears
weeks.
Boolean indicating whether the reporting week corresponds to the W21-W30 period of increased gastroenteritis awareness during the O104:H4 STEC outbreak.
The data on campylobacteriosis cases are queried from the Survstat@RKI database of the German Robert Koch Institute (http://www3.rki.de/SurvStat/).
Data for the computation of absolute humidity were obtained from the German Climiate Service (Deutscher Wetterdienst), Climate data of Germany, available at http://www.dwd.de.
A complete data description and an analysis of the data can be found in:
Bayesian model algorithm for monitoring reported cases of campylobacteriosis in Germany (2013), Manitz J and Höhle M, Biometrical Journal, 55(4), pp. 509 526.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | #Load the data
data("campyDE")
#O104 period is W21-W30 in 2011
stopifnot(all(campyDE$O104period == (
(campyDE$date >= as.Date("2011-05-23")) &
(campyDE$date < as.Date("2011-07-31"))
)))
#Make an sts object from the data.frame
cam.sts <- new("sts",epoch=as.numeric(campyDE$date),observed=campyDE$case,
state=campyDE$state, epochAsDate=TRUE)
#Plot the result
plot(cam.sts)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.