# NOT RUN {
data(kang.peanut)
dat <- kang.peanut
# Table 5.1 of Kang et al. (Chapter 5 of Saxton)
require(reshape2)
round(acast(dat, env~gen, value.var='yield', fun=mean),2)
# GGE biplot of Kang, p. 82.
require(gge)
m1 <- gge(yield~gen*env, data=dat, scale=FALSE)
biplot(m1, flip=c(1,1), main="kang.peanut - GGE biplot")
# }
Run the code above in your browser using DataLab