Aircraft manufacturers are understandably interested in the durability of their products.
As one part of a larger investigation, they looked at how metal joint fasteners would hold up under stress. Several fasteners were tested at several different load levels, and the number of fasteners failing out of the total was recorded. We wish to model the probability of failure as a function of load using logistic regression. You should supply the fitted beta values for each model, and the deviance for each model. I would like to see commented code here.