dat <- ivins.herbs
xyplot(herb~nettle|block, dat) # Nettle is primarily established in legumes.
# O'Brien used first 7 species to test gen differences
dat7 <- droplevels(subset(dat, is.element(gen, c('G01','G02','G03','G04','G05','G06','G07'))))
m1 <- lm(herb ~ gen + block, data=dat7)
anova(m1) # gen p-value is .041
friedman.test(herb ~ gen|block, dat7) # gen p-value .056Run the code above in your browser using DataLab