data(DT_h2)
DT <- DT_h2
head(DT)
# \donttest{
DT=DT[with(DT, order(Env)), ]
Z <- with(DT, smm(Env))
for(i in 1:ncol(Z)){DT[,colnames(Z)[i]] <- Z[,i]}
csdiagFormula <- paste0( "y ~ Env + (", paste(colnames(Z), collapse = "+"), "|| Name)")
blockFormula <- paste0( "(0+", paste(colnames(Z), collapse = "+"), "|| Block)")
form <- paste(csdiagFormula , blockFormula, sep="+")
ans1b <- lmebreed(as.formula(form),
data=DT)
vc <- VarCorr(ans1b); print(vc,comp=c("Variance"))
# }
Run the code above in your browser using DataLab