### ========================== ###
### ========================== ###
data(DT_yatesoats)
DT <- DT_yatesoats
head(DT)
# m3 <- mmer(fixed=Y ~ V + N + V:N,
# random = ~ B + B:MP,
# rcov=~units,
# data = DT)
# summary(m3)$varcomp
#
# m3b <- mmec(fixed=Y ~ V + N + V:N,
# random = ~ B + B:MP,
# rcov=~units,
# data = DT)
# summary(m3b)$varcomp
Run the code above in your browser using DataLab