[R-sig-ME] Multiple independent random effects

Michael Cone coanil at posteo.org
Fri Sep 12 14:59:51 CEST 2014


Mark, I don't think that's possible with lme4/lmer right now.

Michael

Am 12.09.2014 12:12 schrieb Mark Payne:
> Hi,
> 
> I have a mixed-effects model in lme4 like so
> 
> mdl <- lmer(T ~1 + (1|A) + (1|B),...)
> 
> where the factors A and B are being modelled as independent random 
> effects.
> However, there is also heteroscedasticity in the problem, where the
> variance of T depends on a third grouping factor, lets called it C.
> 
> I can fit such a model in the nlme package, using  the
> weights=varIdent(form=~1| C)  argument, but this package doesn't seem 
> to
> easily support independent random effects of the form shown above...
> 
> How can I get the best of both worlds here?
> 
> Mark
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> 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