Dear all,
I need to do a linear mixed model and I want to constrain the varcov matrix
of random effects. Let's say I have 2 random effects A and B, I want the
matrix to be [V,1,1,V], i.e. V=VA=VB, and 1=cor(A,B). Both random effects
are factors with the same levels. Is there any easy way to do that in lme,
lmer, mcmcglmm or something else? (I don't care so much about the package to
use as long as I can do that).
Thanks a lot.

-- 
Alexandre Courtiol

http://sites.google.com/site/alexandrecourtiol/home

	[[alternative HTML version deleted]]


