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
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
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:
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
4 matches
Mail list logo