Hi, Is possible to anova function on R calculate the significance of complete glm model? I explain better. I have a glm model:
model1 <- glm(nsp~caules+mo,family=poisson)
I make an anova of this:
anova(model1,test="F")
Analysis of Deviance Table
Model: poisson, link: log
Response: nsp
Terms added sequentially (first to last)
Df Deviance Resid. Df Resid. Dev F Pr(>F)
NULL 6 40.315
caules 1 23.436 5 16.879 23.4360 1.291e-06 ***
mo 1 0.545 4 16.333 0.5455 0.4601808
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1
It OK, but I need to know the significance of model, I know how to calculate,
is not the problem, the question is: Exist a function that calculate this?
The deviance table that I need is:
Df Deviance Resid. Df Resid. Dev F Pr(>F)
NULL 6 40.315
caules 1 23.436 5 16.879 23.4360 1.291e-06 ***
mo 1 0.545 4 16.333 0.5455 0.4601808
Total 2 23.981 4 16.333 11.9905 0.000??? ***
This is my doubt. If the function anova dont make this, is possible to
implement this in function, or is a wrong idea?
Thanks for all :)
Ronaldo
Steele's Plagiarism of Somebody's Philosophy: Everybody should believe in something -- I believe I'll have another drink. -- | //|\\ [*****************************][*******************] || ( ? ? ) [Ronaldo Reis J?nior ][PentiumIII-600 ] | V [ESALQ/USP-Entomologia, CP-09 ][HD: 30 + 10 Gb ] || / l \ [13418-900 Piracicaba - SP ][RAM: 128 Mb ] | /(lin)\ [Fone: 19-429-4199 r.229 ][Video: SiS620-8Mb ] ||/(linux)\ [chrysopa at insecta.ufv.br ][Modem: Pctel-onboar] |/ (linux) \[ICQ#: 5692561 ][SO: CL 7.0 (2.2.19)] || ( x ) [*****************************][*******************] ||| _/ \_Powered by Conectiva Linux 7.0 D+:) | Lxuser#: 205366 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._