[R-sig-ME] [R] Mixed model help!

Douglas Bates 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
> 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-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-sig-mixed-models mailing list