Re: [R] anova of glm output

2010-09-07 Thread Peng, C

The ANOVA adds the factors only in the order given in the model formula from
left to right. You may try

drop1(out, test=Chisq)


-- 
View this message in context: 
http://r.789695.n4.nabble.com/anova-of-glm-output-tp2528336p2530620.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] anova of glm output

2010-09-06 Thread Dieter Menne


francogrex wrote:
 
 out - glm(response~Var1+Var2+Var3..,family=binomial,data=mydata)
 summary(out)
 stepAIC(out)
 anova(out, test='Chisq')
 I understand that stepAIC is used to select the model with the lowest AIC
 (the best model) but can someone explain what is the purpose of doing the
 anova: anova(out, test='Chisq')? What extra information does it bring?
 Thanks
 

out is of class glm (and lm, if that matters). So we have anova working on
a glm class, which is documented in anova.glm.

You should be aware that this anova depends on the order of you Varx, so be
cautious.

Dieter


-- 
View this message in context: 
http://r.789695.n4.nabble.com/anova-of-glm-output-tp2528336p2528354.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] anova of glm output

2010-09-06 Thread francogrex

Hi, this is more related to understanding some statistics while using R; I've
see such output in a paper:
out - glm(response~Var1+Var2+Var3..,family=binomial,data=mydata)
summary(out)
stepAIC(out)
anova(out, test='Chisq')
I understand that stepAIC is used to select the model with the lowest AIC
(the best model) but can someone explain what is the purpose of doing the
anova: anova(out, test='Chisq')? What extra information does it bring?
Thanks
-- 
View this message in context: 
http://r.789695.n4.nabble.com/anova-of-glm-output-tp2528336p2528336.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.