## Fit model using SpATS.
modSp <- fitTD(TD = TDHeat05,
design = "res.rowcol",
traits = "yield",
what = "fixed")
## Create TD object from the fitted model with BLUEs and standard errors.
TDRes <- STAtoTD(modSp,
what = c("BLUEs", "seBLUEs"))
## Add a weight column in the output.
TDResWt <- STAtoTD(modSp,
what = c("BLUEs", "seBLUEs"),
addWt = TRUE)
Run the code above in your browser using DataLab