[R] output on glm.nb

[ripley@stats.ox.ac.uk]

glm.nb is part of MASS, although you give no credit.  If you read the book
which this software supports you will find worked examples.

Read ?anova.negbin (more) carefully: it does explain this.

On Thu, 28 Nov 2002, Ronaldo Reis Jr. wrote:

> I make a model and have compared with null model.
>
> anova(mnull,m1)
>
> 	Model theta Resid. df 2xlog-lik   test   df   LR stat.  Pr(chi)
> 1	...   0.39    161     -577.9129          NA   NA        NA
> 2	...   1.30    150     -475.6839  1 vs 2  11   102.229   1.11e10-16
>
> anova(m1)
>
> 	Df  Deviance   Resid. Df   Resid. Dev   P(>|chi|)
> NULL                      161         282.139    9...
> ...     ...    ...        ...           ...       ...
> last    2    8.824        150         138.154     ...
>
> My doubt is:
>
> What is my null deviance? 577.9129 or 282.139?

Only one entry is labelled as a `deviance', so just read more carefully.

> What is 2xlog-lik?

Twice the log-likelihood, of course.

-- 
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._