if (FALSE) {
library(dae)
library(asreml)
library(asremlPlus)
## use ?Wheat.dat for data set details
data(Wheat.dat)
# Fit initial model
current.asr <- asreml(yield ~ Rep + WithinColPairs + Variety,
random = ~ Row + Column + units,
residual = ~ ar1(Row):ar1(Column),
data=Wheat.dat)
# Load current fit into an asrtests object
current.asrt <- as.asrtests(current.asr, NULL, NULL)
# check validity of current.asrt
validAsrtests(current.asrt)
}
Run the code above in your browser using DataLab