<!doctype html public "-//w3c//dtd html 4.0 transitional//en">
<html>
Dear R collegues,
<br>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?
<p>R session:
<p><font face="Courier New,Courier"><font size=-1>> t1 <- glm(tax ~
areal, family=poisson(link=log))</font></font>
<br><font face="Courier New,Courier"><font size=-1>> t2 <- summary(t1)</font></font>
<br><font face="Courier New,Courier"><font size=-1>> AIC(t1)</font></font>
<br><font face="Courier New,Courier"><font size=-1>[1] 2173.717</font></font><font face="Courier New,Courier"><font size=-1></font></font>
<p><font face="Courier New,Courier"><font size=-1>> t2[5]</font></font>
<br><font face="Courier New,Courier"><font size=-1>$aic</font></font>
<br><font face="Courier New,Courier"><font size=-1>[1] 35206.69</font></font><font face="Courier New,Courier"><font size=-1></font></font>
<p><font face="Courier New,Courier"><font size=-1>> t2[12]</font></font>
<br><font face="Courier New,Courier"><font size=-1>$aic</font></font>
<br><font face="Courier New,Courier"><font size=-1>[1] 35206.69</font></font><font face="Courier New,Courier"><font size=-1></font></font>
<p><font face="Courier New,Courier"><font size=-1>> t1</font></font>
<br><font face="Courier New,Courier"><font size=-1>Call: glm(formula
= tax ~ areal, family = poisson(link = log))</font></font><font face="Courier New,Courier"><font size=-1></font></font>
<p><font face="Courier New,Courier"><font size=-1>Coefficients:</font></font>
<br><font face="Courier New,Courier"><font size=-1>(Intercept)
areal</font></font>
<br><font face="Courier New,Courier"><font size=-1> 3.19610
0.03029</font></font><font face="Courier New,Courier"><font size=-1></font></font>
<p><font face="Courier New,Courier"><font size=-1>Degrees of Freedom: 283
Total (i.e. Null); 282 Residual</font></font>
<br><font face="Courier New,Courier"><font size=-1>Null Deviance:
40600</font></font>
<br><font face="Courier New,Courier"><font size=-1>Residual Deviance: 32880
AIC: 35210</font></font><font face="Courier New,Courier"><font size=-1></font></font>
<p><font face="Courier New,Courier"><font size=-1>> t2</font></font>
<br><font face="Courier New,Courier"><font size=-1>Call:</font></font>
<br><font face="Courier New,Courier"><font size=-1>glm(formula = tax ~
areal, family = poisson(link = log))</font></font><font face="Courier New,Courier"><font size=-1></font></font>
<p><font face="Courier New,Courier"><font size=-1>Deviance Residuals:</font></font>
<br><font face="Courier New,Courier"><font size=-1> Min
1Q Median 3Q
Max</font></font>
<br><font face="Courier New,Courier"><font size=-1>-25.259
-8.768 -1.270 6.653 25.605</font></font><font face="Courier New,Courier"><font size=-1></font></font>
<p><font face="Courier New,Courier"><font size=-1>Coefficients:</font></font>
<br><font face="Courier New,Courier"><font size=-1>
Estimate Std. Error z value Pr(>|z|)</font></font>
<br><font face="Courier New,Courier"><font size=-1>(Intercept) 3.1960994
0.0362085 88.27 <2e-16 ***</font></font>
<br><font face="Courier New,Courier"><font size=-1>areal
0.0302905 0.0003333 90.88 <2e-16 ***</font></font>
<br><font face="Courier New,Courier"><font size=-1>---</font></font>
<br><font face="Courier New,Courier"><font size=-1>Signif. codes:
0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1</font></font><font face="Courier New,Courier"><font size=-1></font></font>
<p><font face="Courier New,Courier"><font size=-1>(Dispersion parameter
for poisson family taken to be 1)</font></font><font face="Courier New,Courier"><font size=-1></font></font>
<p><font face="Courier New,Courier"><font size=-1> Null
deviance: 40605 on 283 degrees of freedom</font></font>
<br><font face="Courier New,Courier"><font size=-1>Residual deviance: 32881
on 282 degrees of freedom</font></font>
<br><font face="Courier New,Courier"><font size=-1>AIC: 35207</font></font><font face="Courier New,Courier"><font size=-1></font></font>
<p><font face="Courier New,Courier"><font size=-1>Number of Fisher Scoring
iterations: 3</font></font>
<br>
<p>Thank you very much for your help
<p>Juli
<br>
<p>--
<br>Juli G. Pausas
<br>Centro de Estudios Ambientales del Mediterraneo (CEAM)
<br>C/ C.R. Darwin 14, Parc Tecnologic,
<br>46980 Paterna, Valencia, SPAIN
<br>Tel: (+ 34) 96 131 8227; Fax: (+ 34) 96 131 8190
<br><A HREF="mailto:juli@ceam.es">mailto:juli@ceam.es</A>
<br><A HREF="http://www.gva.es/ceam">http://www.gva.es/ceam</A>
<p>GCTE Fire Network - <A HREF="http://www.gva.es/ceam/FireNetwork">http://www.gva.es/ceam/FireNetwork</A>
<br> </html>