data(challenger)
## Not run:
# ch <- glm(cbind(Damage, 6-Damage) ~ Temp, family = binomial, data = challenger)
# with(challenger, plot(Temp, Damage/6))
# lines(challenger$Temp, predict(ch, type = "response") )
# abline(v=32, col = "red", lwd = 2)## temp when challenger launched## End(Not run)
Run the code above in your browser using DataLab