Sorry, I edited the lmer function...
============================================
Dear list,
I fitted a linear mixed effects models to a set of 41 field trials
with plot-level assessments of x,y, for estimating the linear
regression coefficients ?_0 and ?_1
res1 <- lmer(y ~ x+ (x|trial), data=mydata, REML=F)
I wish to write the model equation for its publication, so this is my first try:
W_ij= (?_0 + u_0i)+ (?_1+ u_1i) x_ij + e_ij
where j subscript represents the j-plot within i-trial, both for y or
x. ?0 and ?1 are the population average intercept and slope; u0i and
u1i are the effect of the i-trial on the intercept and the slope,
respectively, considered as random variables (with mean 0 and
variances ?_u0 and ?_u1 a )
I?m not sure if I?m in the right path... I would really appreciate any guidance.
Juan Edwards
National Institute of Agriculture Technology - Argentina