Hi all,

Thank you for answering my previous questions recently.  I have another one
for you.  I ran the following model:

model<-lmer(A~B*Year+(B|Male),

where B is a within-individual mean-centered continuous environmental
variable, year is a factor with 2 levels, and Male is a factor with 59
levels.

Upon creation of a panel plot  with 'B' on the x axis and 'A' on the y-axis
and each of the 59 males in a different panel, I noticed that the plot
patterns don't correspond to the BLUPs from the prinout coef(model).  For
instance, a male that clearly increases 'A' as 'B' increases according to
the panel plot, decreases 'A' as 'B' increases according to the BLUPs from
coef(model).  This is the case with several of the males.  Am I doing
something wrong?  I feel I shouldn't be seeing this pattern.  Any insight
would be much appreciated.

Thank you,
Jacob

	[[alternative HTML version deleted]]

