Dear all,
I am new to R and would like your help with lme formula for partially
crossed random effect in a random-intercept, random-slope model.
In the longitudinal data I have, each subject was tested at 5 different
occasions. The standardized tests were administered by 3 different
examiners, with 2 of them present at all occasions, and the 3rd one
administering tests only on the last two occasions. The subjects  were
randomly assigned to the examiners. There was some attrition over time. I
would like to account for random effect of subject & examiner if needed.

I tested the following models:
model1<-lme(score~time*covariate,random=~time|subject,method="REML",na.action=na.omit,data=dat)
model2<-lme(score~time*covariate,random=list(examiner=~1,subject=~time),method="REML",na.action=na.omit,data=dat)

anova(model1,model2) gives p<0.05 with better model fit for model2.

I would like to know if model2 is the correct way to specify the
partially-crossed random effect in the data I described.

thanks!

	[[alternative HTML version deleted]]

