fit <- emil(y~x1+x2+x3, data=incontinence,family=binomial, method="brglmFit")
# display summary of the beta estimates of the model
summary(fit$fit_y)
# for non-ignorability setting of the missing responses
summary(fit$fit_r)
Run the code above in your browser using DataLab