data(Zambia)
#load map
g = system.file("demodata/zambia.graph", package="INLA")
# add one column for the unstructured spatial effect
Zambia$distr.unstruct = Zambia$district
formula = hazstd ~ f(inla.group(bmi),model="rw2") + f(agc,model="rw2") +
f(district,model="besag",graph=g, param=c(1,0.01)) +
f(distr.unstruct,model="iid",param=c(1,0.01)) +
rcw+edu1+ edu2+ tpr+ sex
mod.zambia = inla(formula,data=Zambia, control.fixed = list(prec.intercept = 0.001, prec = 0.001), verbose=TRUE)