[R] glm.nb, anova.negbin

juli g. pausas juli at ceam.es
Tue Sep 25 18:30:44 CEST 2001


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?
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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://stat.ethz.ch/pipermail/r-help/attachments/20010925/5dbae6ec/attachment.html


More information about the R-help mailing list