dat <- byers.apple
library(lattice)
xyplot(diameter ~ time | factor(appleid), data=dat, type=c('p','l'))
# Overall fixed linear trend, plus random intercept/slope deviations
# for each apple. Observations within each apple are correlated.
library(nlme)
m1 <- lme(diameter ~ 1 + time, data=dat,
random = ~ time|appleid, method='ML',
cor = corAR1(0, form=~ time|appleid),
na.action=na.omit)
VarCorr(m1)Run the code above in your browser using DataLab