Re: [R] comparing mixed binomial model against the same model without random effect

2011-09-25 Thread lincoln
Thank you very much for answering, I have just tried it and these are the results: > random.model<-glmer(sex~hwp+hcp+(1|colony),family=binomial) Mensajes de aviso perdidos glm.fit: fitted probabilities numerically 0 or 1 occurred > no.random.model<-glm(sex~hwp+hcp,family=binomial) Mensajes de a

Re: [R] comparing mixed binomial model against the same model without random effect

2011-09-24 Thread Billy
Hello Simone, Given that your response variable is binary and, consequently, you should use generalized models, just occurs to me a comparison between a Generalized Linear Model (the model without the random effect) and a Generalized Linear Mixed Model (the model with the random effect). You coul

Re: [R] comparing mixed binomial model against the same model without random effect

2011-09-24 Thread lincoln
Any answer to this? I really need to compare a mixed model with binomial error against the same model without the random effect. I would use anova() but I don't know how to specify both models in order to make them comparable. Thanks for any answer Simone -- View this message in context: http:

[R] comparing mixed binomial model against the same model without random effect

2011-09-22 Thread Simone Santoro
Hi everybody, If I am correct, you can compare a model with random effect with the same model without the random effect by using the nlme function, like this: no.random.model <- gls(Richness ~ NAP * fExp, method = "REML", data = RIKZ) random.model <- lme(Richness ~NAP * fExp, data