Just as one more example, if you replace MPG with PRICE (a continuous
variable with one observation per value) in Nick's program, the -logit-
command produces the same dev's as -glm-, because each observation is a
unique covariate pattern. So if your model contains at least one
continuous variable, you are likely to get very close agreement among
commands, as well as programs.
Also, if you don't mind results being slightly off in the 4th or 5th
decimal place, I found that something like this works pretty good:
