# \donttest{
## Fit model using lme4.
lmeMod <- fitTD(TD = TDHeat05,
design = "ibd",
traits = "yield",
what = "fixed",
engine = "lme4")
## Summarize results.
summary(lmeMod)
## Create base plots of the results.
plot(lmeMod)
## Create a pdf report summarizing results.
report(lmeMod,
outfile = tempfile(fileext = ".pdf"),
what = "fixed")
## Fit model using SpATS.
spaMod <- fitTD(TD = TDHeat05,
design = "res.rowcol",
traits = "yield",
what = "fixed")
summary(spaMod)
## Create spatial plots of the results.
plot(spaMod, plotType = "spatial")
## Fit model using SpATS.
## Specify the number of segments to use in PSANOVA
spaMod2 <- fitTD(TD = TDHeat05,
design = "res.rowcol",
traits = "yield",
what = "fixed",
control = list(nSeg = c(13, 26)))
## Create a pdf report summarizing results.
report(spaMod, outfile = tempfile(fileext = ".pdf"), what = "fixed")
## Fit model using asreml.
if (requireNamespace("asreml", quietly = TRUE)) {
asrMod <- fitTD(TD = TDHeat05,
design = "res.rowcol",
traits = "yield",
what = "fixed",
engine = "asreml")
summary(asrMod)
report(asrMod,
outfile = tempfile(fileext = ".pdf"),
what = "fixed")
}
# }
Run the code above in your browser using DataLab