<!doctype html public "-//w3c//dtd html 4.0 transitional//en">
<html>
Dear R-collegues,
<br>I'm getting an error message (Error in round) when summarising a glm.nb
model, and when using anova.negbin (in R 1.3.1 for windows):
<p><font face="Courier New,Courier"><font size=-1>> m.nb <- glm.nb(tax
~ areal)</font></font>
<br><font face="Courier New,Courier"><font size=-1>> m.bn</font></font>
<br><font face="Courier New,Courier"><font size=-1>Call: glm.nb(formula
= tax ~ areal, init.theta = 5.08829537115498,
link = log)</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.03146
0.03182</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:
353.3</font></font>
<br><font face="Courier New,Courier"><font size=-1>Residual Deviance: 293.3
AIC: -2334000</font></font>
<p><font face="Courier New,Courier"><font size=-1>> summary(m.nb)</font></font>
<p><font face="Courier New,Courier"><font size=-1>Call:</font></font>
<br><font face="Courier New,Courier"><font size=-1>glm.nb(formula = tax
~ areal, init.theta = 5.08829537115498,</font></font>
<br><font face="Courier New,Courier"><font size=-1> link
= log)</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>-2.7107 -0.8313
-0.1159 0.5866 1.9634</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.031460
0.445165 6.810 9.78e-12 ***</font></font>
<br><font face="Courier New,Courier"><font size=-1>areal
0.031819 0.004151 7.666 1.78e-14 ***</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>
<p><font face="Courier New,Courier"><font size=-1>(Dispersion parameter
for Negative Binomial(5.0883) family taken to be 1)</font></font>
<p><font face="Courier New,Courier"><font size=-1> Null
deviance: 353.34 on 283 degrees of freedom</font></font>
<br><font face="Courier New,Courier"><font size=-1>Residual deviance: 293.34
on 282 degrees of freedom</font></font>
<br><font face="Courier New,Courier"><font size=-1>AIC: -2334394</font></font>
<p><font face="Courier New,Courier"><font size=-1>Number of Fisher Scoring
iterations: 1</font></font>
<p><font face="Courier New,Courier"><font size=-1>Correlation of Coefficients:</font></font>
<br><font face="Courier New,Courier"><font size=-1>
(Intercept)</font></font>
<br><font face="Courier New,Courier"><font size=-1>areal
-0.9982</font></font>
<p><font face="Courier New,Courier"><font size=-1>Error in round(x$theta,
dp, nsmall = dp) :</font></font>
<br><font face="Courier New,Courier"><font size=-1>
unused argument(s) (nsmall ...)</font></font>
<p><font face="Courier New,Courier"><font size=-1>> anova.negbin(m.nb)</font></font>
<br><font face="Courier New,Courier"><font size=-1>Error in get(x, envir,
mode, inherits) : invalid first argument</font></font>
<br><font face="Courier New,Courier"><font size=-1>In addition: Warning
message:</font></font>
<br><font face="Courier New,Courier"><font size=-1>tests made without re-estimating
theta in: anova.negbin(m.nb)</font></font>
<br><font face="Courier New,Courier"><font size=-1>></font></font>
<p>Any suggestion of why I'm getting this error?
<br>Many thanks
<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>