<!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 &lt;- glm(tax ~
areal, family=poisson(link=log))</font></font>
<br><font face="Courier New,Courier"><font size=-1>> t2 &lt;- 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:&nbsp; 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)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
areal</font></font>
<br><font face="Courier New,Courier"><font size=-1>&nbsp;&nbsp;&nbsp; 3.19610&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
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);&nbsp; 282 Residual</font></font>
<br><font face="Courier New,Courier"><font size=-1>Null Deviance:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
40600</font></font>
<br><font face="Courier New,Courier"><font size=-1>Residual Deviance: 32880&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
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>&nbsp;&nbsp;&nbsp; Min&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
1Q&nbsp;&nbsp; Median&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3Q&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
Max</font></font>
<br><font face="Courier New,Courier"><font size=-1>-25.259&nbsp;&nbsp;
-8.768&nbsp;&nbsp; -1.270&nbsp;&nbsp;&nbsp; 6.653&nbsp;&nbsp; 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>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
Estimate Std. Error z value Pr(>|z|)</font></font>
<br><font face="Courier New,Courier"><font size=-1>(Intercept) 3.1960994&nbsp;
0.0362085&nbsp;&nbsp; 88.27&nbsp;&nbsp; &lt;2e-16 ***</font></font>
<br><font face="Courier New,Courier"><font size=-1>areal&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.0302905&nbsp; 0.0003333&nbsp;&nbsp; 90.88&nbsp;&nbsp; &lt;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:&nbsp;
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>&nbsp;&nbsp;&nbsp; Null
deviance: 40605&nbsp; on 283&nbsp; degrees of freedom</font></font>
<br><font face="Courier New,Courier"><font size=-1>Residual deviance: 32881&nbsp;
on 282&nbsp; 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>&nbsp;
<p>Thank you very much for your help
<p>Juli
<br>&nbsp;
<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>&nbsp;</html>