Code for Seeds example

data(Seeds)

formula = r ~ x1*x2+f(plate,model="iid")

mod.seeds = inla(formula,data=Seeds,family="binomial",Ntrials=n)

## improved estimation of the hyperparameters

h.seeds = inla.hyperpar(mod.seeds)