options(
contrasts = c(unordered = "contr.SAS", ordered = "contr.poly"))
data(PBIB)
formula( PBIB )
fm1PBIB <-
lme( response ~ Treatment, data = PBIB, random = ~ 1 | Block)
summary( fm1PBIB ) # compare with output 1.7 pp. 24-25
VarCorr( fm1PBIB )
anova( fm1PBIB )Run the code above in your browser using DataLab