if (FALSE) {
library(agridat)
data(versteeg.maize)
dat <- versteeg.maize
dat <- transform(dat, farm=as.factor(farm),
gen=as.factor(gen), system=as.factor(system))
# anova of Versteeg page 115
# plots are different sizes, so we have to manually
# specify the error terms for testing
m1 <- aov(yield ~ 1 + farm + gen + system + gen:system +
Error(gen:farm + system:farm), data=dat)
summary(m1)
libs(dplyr)
d1 <- summarize(dat, mn=mean(yield), .by=farm)
dat <- left_join(dat, d1)
# stability analysis, figure 1 of Versteeg, page 119
libs(lattice)
xyplot(yield~mn, data=dat, groups=gen, type=c('p','r'))
}
Run the code above in your browser using DataLab