data(guImmun)
fm1 <-
GLMM(immun ~ kid2p + mom25p + ord + ethn +
momEd + husEd + momWork + rural + pcInd81,
family = binomial, data = guImmun, random = ~1|comm)
summary(fm1)
Run the code above in your browser using DataLab