# NOT RUN {
## Run model using SpATS.
myModel <- fitTD(TD = TDHeat05, design = "res.rowcol", traits = "yield",
what = "fixed")
## Create TD object from the fitted model with BLUEs and standard errors.
myTD <- STAtoTD(myModel)
## Add a weight column in the output.
myTDWt <- STAtoTD(myModel, addWt = TRUE)
# }
Run the code above in your browser using DataLab