On Fri, Mar 13, 2009 at 2:39 PM, Martín Quiroga <mquir...@ssdfe.com.ar> wrote: > Hi everyone! I am a biologist from Argentina and have to solve this problem. > I have an insect population obtained from 10 different nests and need to > know its sex ratio. But as I cannot ensure insects independence I need to > run a model where I can include the variable “nest” as with a random effect. > The response variable has a binomial distribution (males or females). > I’ve been reading for a while and found the MASS and lmer packages that will > allow me to do such a thing with my data. I found the script I should write > “lmer(y~fixed+(time | random), family=binomial)” but, despite of using the R > commander cannot make it work. Can you help me with this? I’ll appreciate > that!
The package is called lme4 and the function is called lmer. To be able to use the lmer function you must first enter library(lme4) in a script or console session (I'm not sure how this would be done in R Commander). There is a mailing list, r-sig-mixed-mod...@r-project.org specifically for discussion of mixed models. I am cc:ing that list on this reply. ______________________________________________ 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.
