dat <- cochran.bib
m1 = aov(yield ~ gen + Error(loc), data=dat)
summary(m1)
require(nlme)
m2 = lme(yield ~ -1 + gen, data=dat, random=~1|loc)
require(agricolae)
BIB.test(dat$loc, dat$gen, y=dat$yield)Run the code above in your browser using DataLab