####=========================================####
#### For CRAN time limitations most lines in the
#### examples are silenced with one '#' mark,
#### remove them and run the examples
####=========================================####
data(DT_h2)
DT <- DT_h2
head(DT)
####=========================================####
#### fit the mixed model (very heavy model)
####=========================================####
# DT=DT[with(DT, order(Env)), ]
# ans1 <- mmes(y~Env, henderson=TRUE,
# random=~vsm(dsm(Env),ism(Name)) + vsm(dsm(Env),ism(Block)),
# rcov=~vsm(dsm(Env),ism(units)),
# data=DT)
# summary(ans1)$varcomp
#
Run the code above in your browser using DataLab