Examples and tutorials‎ > ‎Volume I‎ > ‎

Code for Surg example

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