Displaying 1 result from an estimated 1 matches for "fbeach".
Did you mean:
beach
2011 Sep 22
2
comparing mixed binomial model against the same model without random effect
...e 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 = RIKZ,
random = ~1 | fBeach, method = "REML")
anova(no.random.model,random.model)
But, nlme is valid only for the gaussian family, isn't it? In my case I have a mixed model with binomial family, like this:
random.model <- lme(sex ~hwp+hcp, data = mydata,
random = ~1 | colony, method = "REML&...