# NOT RUN {
data(DT_yatesoats)
DT <- DT_yatesoats
m3 <- mmer(fixed=Y ~ V + N + V:N,
random = ~ B + B:MP,
rcov=~units,
data = DT)
summary(m3)$varcomp
p0 <- predict.mmer(object=m3, classify = "N")
p0$pvals
p0 <- predict.mmer(object=m3, classify = "V")
p0$pvals
p0 <- predict.mmer(object=m3, classify = "B")
p0$pvals
# }
Run the code above in your browser using DataLab