[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