To give this pot another stir, and to use
mutually accessible data, -logit- and -glm-
with logit link and binomial family give quite
different deviance residuals.
The pattern of -glm-'s makes sense,
but that of -logit-'s is more puzzling.
Ok, after looking in several wrong places, I finally found an explanation
in the Stata Reference Manual G-M, pp. 315-316. It says that "All the
residual and diagnostic statistics calculated by Stata [NOTE: I think it
really means Stata logistic regression and some related routines] are in
terms of covariate patterns, not observations. That is, all observations
with the same covariate patterns are given the same residual and diagnostic
statistics." It says that Hosmer and Lemeshow argue that this is the
better way to do it.
