Hi,

As far as I know such heterogenity of residual variances in lmer is not
possible. However you can easily fit such model in MCMCglmm using
rcov=~us(Factor):units where factor is a categorical variable in which
levels residual var shoul vary. Should you want to restrict residual
covariance to 0 use idh() instead of us().

Cheers,
Sz.

	[[alternative HTML version deleted]]


