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