Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
st: AIC and BIC in Poisson and GLM
From
"Lachenbruch, Peter" <[email protected]>
To
"'[email protected]'" <[email protected]>
Subject
st: AIC and BIC in Poisson and GLM
Date
Mon, 25 Apr 2011 09:21:08 -0700
I was developing a simple class example for poisson regression using the poisson command and the glm command. The log likelihood was the same in both regressions; the coefficients were the same. When I looked at the AIC and BIC reported by glm I got
AIC = 7.920031
Log likelihood = -33.60015344 BIC = 2.922026
For the estat command with poisson I got
. estat ic
-----------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+---------------------------------------------------------------
full | 10 -495.0676 -33.60015 6 79.20031 81.01582
-----------------------------------------------------------------------------
Note: N=Obs used in calculating BIC; see [R] BIC note
I would expect these to be the same, but they ain't. I suspect there may be a normalizing constant lurking around here somewhere.
I don't want to fix this unless it's truly a bug; but I would like to be able to explain this to my students.
Tony
Peter A. Lachenbruch
Department of Public Health
Oregon State University
Corvallis, OR 97330
Phone: 541-737-3832
FAX: 541-737-4001
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/