data(Salm)
mod.salm = inla(y~log(dose+10)+dose+f(rand,model="iid",param=c(0.001,0.001)),
family="poisson",data=Salm)
## improved estimation of the hyperparameters
h.salm = inla.hyperpar(mod.salm)