data(SIMS)
fm1SIMS <- lme(Gain ~ Pretot, data = SIMS,
random = ~ Pretot | Class, control = list(msVerbose = TRUE))
summary(fm1SIMS) # compare to output 7.4, p. 262
VarCorr(fm1SIMS)
anova(fm1SIMS)
intervals( fm1SIMS )Run the code above in your browser using DataLab