# NOT RUN {
## Fit a model using lme4.
myModel <- fitTD(TD = TDHeat05, traits = "yield", design = "res.rowcol",
engine = "lme4")
## Detect outliers in the standardized residuals of the fitted model.
outliers <- outlierSTA(STA = myModel, traits = "yield")
# }
Run the code above in your browser using DataLab