data(Surg)
formula = r ~ f(hospital,model="iid",param=c(0.001,0.001))
mod.surg = inla(formula,data=Surg,family="binomial",Ntrials=n)
# improved estimation of the hyperparameters
h.surg = inla.hyperpar(mod.surg)