[R-sig-ME] lmer with binomial distribution of random effects

Viechtbauer Wolfgang (STAT) wolfgang.viechtbauer at maastrichtuniversity.nl
Tue Mar 10 09:45:07 CET 2015


The lcmm package looks promising:

http://finzi.psych.upenn.edu/R/library/lcmm/html/lcmm-package.html

And it seems there may be relevant stuff in DPpackage:

http://finzi.psych.upenn.edu/R/library/DPpackage/html/00Index.html

Best,
Wolfgang

> -----Original Message-----
> From: R-sig-mixed-models [mailto:r-sig-mixed-models-bounces at r-
> project.org] On Behalf Of Ben Bolker
> Sent: Monday, March 09, 2015 02:16
> To: r-sig-mixed-models at r-project.org
> Subject: Re: [R-sig-ME] lmer with binomial distribution of random effects
> 
> El Kihal, Siham <Siham.ElKihal at ...> writes:
> 
> >
> > Dear lmer() friends,
> >
> > I am trying to estimate a model with a random
> > intercept, and 2 random slopes.
> > I believe that my betas (slopes) do not follow
> >a normal distribution, but rather a bimodal distribution.
> > The reason for this that there are two possible
> > mechanisms that influence the evolution of this variable,
> > one with a negative influence and another one with a
> >  positive influence. This is why I need to use a bimodal
> > distribution for my slopes to avoid the fact that
> > both effects right now cancel out.
> >
> > Does anyone of you has already done this or has
> > an idea how to concretely implement this using lmer()?
> 
>   This sounds like a latent mixture model problem.  lme4 doesn't
> do this; you *might* be able to implement an expectation-maximization
> wrapper around lme4 that would do it, but it wouldn't be entirely
> trivial.  If I had to do this I would probably turn to JAGS/BUGS.
> Looking forward to other answers from the list ...
> 
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models



More information about the R-sig-mixed-models mailing list