data(DT_example)
DT <- DT_example
A <- A_example
head(DT)
## Compound simmetry (CS) model
ans1 <- lmebreed(Yield~Env + (1|Name) + (1|Env:Name),
data=DT)
vc <- VarCorr(ans1); print(vc,comp=c("Variance"))
PEV <- getCi(ans1)
# image(PEV)
Run the code above in your browser using DataLab