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.
. list glmdevsum logdevsum in 1 <---- the sum of the dev's for all obs
+-----------------------+
| glmdevsum logdevsum |
|-----------------------|
1. | -9.401649 -9.401649 |
+-----------------------+
. list glmdev logitdev in 1/10, clean
glmdev logitdev
1. -.8137717 -.8137717
2. -.8217416 -.8217416
3. -.8101112 -.8101112
4. -.8225662 -.8225662
5. -.8601973 -.8601973
6. -.8345912 -.8345912
7. -.8181056 -.8181056
8. -.8271667 -.8271667
9. -.8928685 -.8928685
10. -.8135639 -.8135639
Paul
______________________________
Paul F. Visintainer, PhD
Professor and Program Director
Health Quantitative Sciences
School of Public Health
New York Medical College
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
