## Run a single trait analysis using SpATS.
modSp <- fitTD(TD = TDHeat05,
design = "res.rowcol",
traits = "yield")
## Create base plots.
plot(modSp,
what = "fixed",
plotType = "base")
## Create spatial plots.
plot(modSp,
what = "fixed",
plotType = "spatial")
## Create spatial plots showing the spatial trend as percentage.
plot(modSp,
what = "fixed",
plotType = "spatial",
spaTrend = "percentage")
Run the code above in your browser using DataLab