source("http://192.38.117.59/~linearpredictors/datafiles/readSurgery.R")

## remove surgery type "orthopedic"
surgery2 <- subset(surgery, surgtype %in% 2:3)
surgery2$surgtype <- surgery2$surgtype[drop = TRUE]

logistic <- glm(complication ~ surgtype -1,
                family = binomial,
                data = surgery2)

## estimates and standard deviations
summary(logistic)

## confidence intervals for the parameters
confint(logistic)

## CI for the estimated probabilities
CIlogp <- confint(logistic)
exp(CIlogp)/(1+exp(CIlogp))