# NOT RUN {
data(weiss.incblock)
dat <- weiss.incblock
# True aspect as shown in Weiss and Cox
desplot(yield~col*row, dat, text=gen, shorten='none', cex=.8, out1=block,
main="weiss.incblock", xlab="Soybean yields", aspect=252/96)
# }
# NOT RUN {
# Standard inc block analysis used by Weiss and Cox
require(asreml)
m1 <- asreml(yield ~ gen + block , data=dat)
predict(m1, classify="gen")$predictions$pvals
## gen pred.value std.error est.stat
## G01 24.59 0.8312 Estimable
## G02 26.92 0.8312 Estimable
## G03 32.62 0.8312 Estimable
## G04 26.97 0.8312 Estimable
## G05 26.02 0.8312 Estimable
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab