[R] Mixed model help!
bates at stat.wisc.edu
Sat Mar 14 14:04:30 CET 2009
On Fri, Mar 13, 2009 at 2:39 PM, Martín Quiroga <mquiroga at 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
The package is called lme4 and the function is called lmer. To be
able to use the lmer function you must first enter
in a script or console session (I'm not sure how this would be done in
There is a mailing list, R-SIG-Mixed-Models at R-project.org specifically
for discussion of mixed models. I am cc:ing that list on this reply.
More information about the R-help