Dear all,

I have a model like: y.m1 <- lmer(y~1+factor(Site)+(1|Species))

so a fixed site effect is crossed with a random effect of  a species.
Does the actual equation yield to:

y(site, species) = mean + SIte1_effect + Site2_effect + ...... SiteN_effect
+ Species_effect + error

ie allowing different Sites to have different variances?

best regars
Nikos

	[[alternative HTML version deleted]]


