[R] GLM different results with the same factors
gaiarrido
gaiarrido at usal.es
Sun Jul 24 17:53:14 CEST 2011
I've read something about this problem, but I don't know how can i avoid this
problem.
Why the order of the factors give different results? I suppose it's because
the order of the factors, i've just changed "lcc" from the first position to
the last in the model, and the significance change completely
> modo<-glm(prevalencia~*lcc*+edadysexo/lcc+edadysexo/mes,binomial)
> anova(modo,test="Chisq")
Df Deviance Resid. Df Resid. Dev P(>|Chi|)
NULL 524 206.97
lcc 1 10.5715 523 196.40 0.001148 **
edadysexo 2 1.0725 521 195.32 0.584929
lcc:edadysexo 2 3.7752 519 191.55 0.151433
edadysexo:mes 9 18.2981 510 173.25 0.031868 *
> mode<-glm(prevalencia~edadysexo/lcc+edadysexo/mes+*lcc*,binomial)
> anova(mode,test="Chisq")
Df Deviance Resid. Df Resid. Dev P(>|Chi|)
NULL 524 206.97
edadysexo 2 9.9165 522 197.05 0.007025 **
lcc 1 1.7275 521 195.32 0.188732
edadysexo:lcc 2 3.7752 519 191.55 0.151433
edadysexo:mes 9 18.2981 510 173.25 0.031868 *
Ijow can i know what's correct? when i test this two factors separately in a
lm (lcc is continuos) and in a chisq.test (edadysexo is categorical) both
are significant, and in the model just one of them is significant.
Thanks very much
-----
Mario Garrido Escudero
PhD student
Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola
Universidad de Salamanca
--
View this message in context: http://r.789695.n4.nabble.com/GLM-different-results-with-the-same-factors-tp3690471p3690471.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list