# NOT RUN {
# Hosmer, D. (2013) page 3: Age and coranary Heart Disease (CHD) Status of 20 subjects:
AGE <- c(20, 23, 24, 25, 25, 26, 26, 28, 28, 29, 30, 30, 30, 30, 30, 30, 30, 32, 33, 33)
CHD <- c(0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0)
data <- data.frame (CHD, AGE)
Ela <- lsm(CHD ~ AGE, family = binomial, data)
confint(Ela)
# }
Run the code above in your browser using DataLab