data(Genetics)
options(
contrasts = c(unordered = "contr.SAS", ordered = "contr.poly"))
formula( Genetics )
fm1Gen <- lme( Yield ~ 1, data = Genetics,
random = list(Location = pdCompSymm(~ Family - 1), Block = ~ 1) )
summary( fm1Gen )
VarCorr( fm1Gen ) # compare with Output 4.18, p. 165
anova( fm1Gen )Run the code above in your browser using DataLab