No problem,

lme1 <- lme(iAUC~Treatment+I(Dose-8),
random=~Treatment+I(Dose-8)|Subject,data=gluc)
summary(lme1)

works fine for me.

Your problem might have been that you had lmer and nlme loaded at the same
time, or it might just have been the spelling mistake (+Does)


Andy.


andydolman@gmail.com


Maybe, since the data are balanced, I could use lme(), but I have no clue
> how to specify this random effect structure there.  I tried
> random=~Treatment+Does|Subject but lme() didn't like that.
>
> Thanks for listening.  I appreciate all you've done already.
>
> Kevin
>
> --
> Kevin E. Thorpe
> Biostatistician/Trialist, Knowledge Translation Program
> Assistant Professor, Dalla Lana School of Public Health
> University of Toronto
> email: kevin.thorpe@utoronto.ca  Tel: 416.864.5776  Fax: 416.864.3016
>
>

	[[alternative HTML version deleted]]


