DT <- DT_example
A <- A_example
DT2 <- fillData(DT , toBalanceSplit="Env", toBalanceFill="Name" )
head(DT2)
# \donttest{
## Compound simmetry (CS) model
ans2 <- lmebreed(Yield~Env+ (0+Env|Name),
relmat = list(Name=A),
rotation = TRUE,
data=DT2)
vc <- VarCorr(ans2); print(vc,comp=c("Variance"))
BLUP <- ranef(ans2, condVar=TRUE)$Name
SEs <- attr(BLUP, which="postVar")[,,]
# }
Run the code above in your browser using DataLab