# NOT RUN {
data(medpar)
glmp <- glm(los ~ hmo + white + factor(type),
family = poisson, data = medpar)
summary(glmp)
exp(coef(glmp))
ml.p <- ml_glm(los ~ hmo + white + factor(type),
family = "poisson",
link = "log",
data = medpar)
summary(ml.p)
library(MASS)
glmnb <- glm.nb(los ~ hmo + white + factor(type),
data = medpar)
summary(glmnb)
exp(coef(glmnb))
# }
Run the code above in your browser using DataLab