On 01/06/12 02:52, Mark Leeds wrote:
just one other thing about the AIC issue:
there is a line in glm.fit which is the following:
aic = aic(y, n, mu, weights, dev) + 2 * rank
but I couldn't find the function aic so I couldn't investigate further. It
looks suspicious though because it seems to me like
it should be
aic = -2*likelihood + 2 * rank
if anyone could help me find the aic function it's appreciated.
Dear Mark,
The aic() function depends on the family that you are using.
E.g.:
binomial()$aic
gaussian()$aic
cheers,
Rolf
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.