# #### =================================== ####
# #### ===== Augmented Block Design 1 ==== ####
# #### =================================== ####
# data(DT_expdesigns)
# DT <- DT_expdesigns
# names(DT)
# data1 <- DT$au1
# head(data1)
# ## response variable: "yield"
# ## check indicator: "entryc" ('nc' for all unreplicated, but personal.name for checks)
# ## blocking factor: "block"
# ## treatments, personal names for replicated and non-replicated: "trt"
# ## check no check indicator: "new"
# mix1 <- mmer(yield~entryc,
# random=~block+trt,
# rcov=~units, tolpar = 1e-6,
# data=data1)
# summary(mix1)$varcomp
#
# mix1b <- mmec(yield~entryc,
# random=~block+trt,
# rcov=~units, tolParConv = 1e-6,
# data=data1)
# summary(mix1b)$varcomp
Run the code above in your browser using DataLab