data(est)
f_fit <- emforbeta(survive ~ Fetoprtn + Antigen + Jaundice + Age,
data = est,
family = binomial, method = "glm.fit")
summary(f_fit$mfit)
f_fit <- emforbeta(survive ~ Fetoprtn + Antigen + Jaundice + Age,
data = est,
family = binomial, method = "brglmFit")
summary(f_fit$mfit)
Run the code above in your browser using DataLab