pdf("challenger.pdf")

ch0 <- read.fwf("C:\\data\\challenger.dat", width=c(8,8,8,8), col.names=c("flight", "degrees", "numfail", "indfail"))

ch1 <- data.frame(degrees=ch0$degrees, indfail=ch0$indfail)

attach(ch1)

chmod1 <- glm(formula=indfail ~ degrees, family=binomial("logit"))

summary(chmod1)

anova(chmod1, test="Chisq")

dev.off()