Hello,

I have a mixed model two-way anova with interaction: y=a+b+a*b+e, where
both a and b are factors, and b and the interaction a*b are random.
I used lme: fm <- lme(y ~ a, random = ~ 1|b/a,method="REML").
I want to extract the variance of the interaction (which should be equal to
(MSAB-MSE)/n   in a balanced model), is it possible?

Thanks,
Tal

	[[alternative HTML version deleted]]

