# NOT RUN {
data(example)
mix <- mmer2(Yield~1,
random=~at(Env):Block + us(Env):Name,
rcov=~at(Env):units, data=example)
summary(mix)
## compare with asreml, you'll get the same result
# }
Run the code above in your browser using DataLab