pdf("wintwheat.pdf") ww0 <- read.fwf("C:\\data\\wintwheat.dat", width=c(8,8,8,8), col.names=c("idnum", "variety", "yield", "moist")) ww1 <- data.frame(yield=ww0$yield, moist=ww0$moist, variety=factor(ww0$variety)) attach(ww1) library(nlme) yield.rcr1 <- lme(yield ~ moist , random =~1+moist | variety) summary(yield.rcr1) ranef(yield.rcr1) dev.off()