#load the data
data(Munich)
#load the map
g = system.file("demodata/munich.graph", package="INLA")
## Note that here we what to have an estimator of the effect of `year`
## also the for years where we have no observation, therefore we give a
## vector with all possible values assumed by the covariate `year`, that
## is seq(1918,2001)
formula = rent ~ f(location,model="besag",graph.file=g,initial=1,param=c(1,0.01)) +
f(year,model="rw2",values=seq(1918,2001),param=c(1,0.01)) +
f(floor.size,model="rw2",param=c(1,0.01)) +
Gute.Wohnlage + Beste.Wohnlage + Keine.Wwv + Keine.Zh +
Kein.Badkach + Besond.Bad + Gehobene.Kueche +
zim1 + zim2 + zim3 + zim4 + zim5 + zim6 -1
mod = inla(formula,data=Munich, verbose = TRUE,
control.family=list(initial=-1),
control.predictor = list(initial = 12),
control.inla = list(h=1e-4))