[R] AIC
juli g. pausas
juli at ceam.es
Tue Sep 11 12:52:50 CEST 2001
Dear R collegues,
I'm trying to understand what's AIC in R (ver. 1.3.1), and I'm getting a
different answer if I look at the AIC(of the fitted model) or the aic in
the summary( of the fitted model). Is this correct? Can somebody explain
me the difference between the two values? or Is the AIC criterion not
appropiated for Poisson models?
R session:
> t1 <- glm(tax ~ areal, family=poisson(link=log))
> t2 <- summary(t1)
> AIC(t1)
[1] 2173.717
> t2[5]
$aic
[1] 35206.69
> t2[12]
$aic
[1] 35206.69
> t1
Call: glm(formula = tax ~ areal, family = poisson(link = log))
Coefficients:
(Intercept) areal
3.19610 0.03029
Degrees of Freedom: 283 Total (i.e. Null); 282 Residual
Null Deviance: 40600
Residual Deviance: 32880 AIC: 35210
> t2
Call:
glm(formula = tax ~ areal, family = poisson(link = log))
Deviance Residuals:
Min 1Q Median 3Q Max
-25.259 -8.768 -1.270 6.653 25.605
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.1960994 0.0362085 88.27 <2e-16 ***
areal 0.0302905 0.0003333 90.88 <2e-16 ***
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 40605 on 283 degrees of freedom
Residual deviance: 32881 on 282 degrees of freedom
AIC: 35207
Number of Fisher Scoring iterations: 3
Thank you very much for your help
Juli
--
Juli G. Pausas
Centro de Estudios Ambientales del Mediterraneo (CEAM)
C/ C.R. Darwin 14, Parc Tecnologic,
46980 Paterna, Valencia, SPAIN
Tel: (+ 34) 96 131 8227; Fax: (+ 34) 96 131 8190
mailto:juli at ceam.es
http://www.gva.es/ceam
GCTE Fire Network - http://www.gva.es/ceam/FireNetwork
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://stat.ethz.ch/pipermail/r-help/attachments/20010911/6365c884/attachment.html
More information about the R-help
mailing list