# NOT RUN {
####=========================================####
#### AUGMENTED DESIGN EXAMPLE
####=========================================####
data(augment)
head(augment)
####=========================================####
#### fit the mixed model and check summary
####=========================================####
mix1 <- mmer2(TSW ~ Check.Gen, random = ~ Block + Genotype:Check, data=augment)
summary(mix1)
blup <- mix1$u.hat$`Genotype:Check`
####=========================================####
#### get only the blups for unreplicated genotypes
####=========================================####
blup2 <- blup[grep("Check1",rownames(blup)),]
# }
Run the code above in your browser using DataLab