if (FALSE) {
## running this example takes a couple of minutes
#Load the campylobacteriosis data for Germany
data("campyDE")
#Make an sts object from the data.frame
cam.sts <- sts(epoch=campyDE$date,
observed=campyDE$case, state=campyDE$state)
#Define monitoring period
# range <- which(epoch(cam.sts)>=as.Date("2007-01-01"))
# range <- which(epoch(cam.sts)>=as.Date("2011-12-10"))
range <- tail(1:nrow(cam.sts),n=2)
control <- list(range=range, X=NULL, trend=TRUE, season=TRUE,
prior='iid', alpha=0.025, mc.munu=100, mc.y=10,
samplingMethod = "joint")
#Apply the boda algorithm in its simples form, i.e. spline is
#described by iid random effects and no extra covariates
library("INLA") # needs to be attached
cam.boda1 <- boda(cam.sts, control=control)
plot(cam.boda1, xlab='time [weeks]', ylab='No. reported', dx.upperbound=0)
}
Run the code above in your browser using DataLab