y <- rbinom(100, 1, 0.5)
x <- ifelse(y == 1, rnorm(length(y[y==1]), 100, 20),
rnorm(length(y[y==0]), 50, 20))
mod <- glm(y ~ x, family=binomial)
summary(mod)
mysummary(mod, exp=TRUE, robust=TRUE)
Run the code above in your browser using DataLab