data(Scotland)
#load the map
g = system.file("demodata/scotland.graph", package="INLA")
formula = Counts~f(Region,model="besag", graph.file = g,
param=c(0.5,0.0005)) + I(X/10)
mod.scotland = inla(formula,family="poisson",E=E,data=Scotland)