glm.nb, anova.negbin
On Tue, 25 Sep 2001, juli g. pausas wrote:
Dear R-collegues, 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):
m.nb <- glm.nb(tax ~ areal) m.bn
Call: glm.nb(formula = tax ~ areal, init.theta = 5.08829537115498,
link = log)
Coefficients:
(Intercept) areal
3.03146 0.03182
Degrees of Freedom: 283 Total (i.e. Null); 282 Residual
Null Deviance: 353.3
Residual Deviance: 293.3 AIC: -2334000
summary(m.nb)
Call:
glm.nb(formula = tax ~ areal, init.theta = 5.08829537115498,
link = log)
Deviance Residuals:
Min 1Q Median 3Q Max
-2.7107 -0.8313 -0.1159 0.5866 1.9634
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.031460 0.445165 6.810 9.78e-12 ***
areal 0.031819 0.004151 7.666 1.78e-14 ***
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1
(Dispersion parameter for Negative Binomial(5.0883) family taken to be
1)
Null deviance: 353.34 on 283 degrees of freedom
Residual deviance: 293.34 on 282 degrees of freedom
AIC: -2334394
Number of Fisher Scoring iterations: 1
Correlation of Coefficients:
(Intercept)
areal -0.9982
Error in round(x$theta, dp, nsmall = dp) :
unused argument(s) (nsmall ...)
anova.negbin(m.nb)
Error in get(x, envir, mode, inherits) : invalid first argument In addition: Warning message: tests made without re-estimating theta in: anova.negbin(m.nb)
Any suggestion of why I'm getting this error?
Because there is a copying error in lines
cat("\n Theta: ", format(round(x$theta, dp, nsmall=dp)),
"\n Std. Err.: ", format(round(x$SE.theta, dp,
nsmall=dp)),
"\n")
if(!is.null(x$th.warn))
cat("Warning while fitting theta:", x$th.warn,"\n")
cat("\n 2 x log-likelihood: ", format(round(x$twologlik, 3,
nsmall=dp)), "\n
Move the second ) before the nsmall.
Or just use the debugging tools to find it ....
Many thanks 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
Brian D. Ripley, ripley at stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272860 (secr) Oxford OX1 3TG, UK Fax: +44 1865 272595 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._