[R] significance of glm complete model
Ronaldo Reis Jr.
chrysopa at terra.com.br
Thu Oct 25 13:49:20 CEST 2001
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list