Examples and tutorials‎ > ‎Volume I‎ > ‎

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)
Comments