# NOT RUN {
data(medpar)
med.nb.g <- ml_glm3(los ~ hmo + white,
family = "gNegBinomial",
link = "log",
group = medpar$provnum,
data = medpar)
summary(med.nb.g)
# }
Run the code above in your browser using DataLab