# Yield
library(agricolae)
data(huasahuasi)
attach(huasahuasi)
YIELD<-huasahuasi$YIELD
attach(YIELD)
market <- y1da + y2da
non_market <- y3da
yield <- market + non_market
model<-strip.plot(block, clon, trt, yield)
comparison<-LSD.test(yield,clon,model$gl.a,model$Ea)
comparison<-LSD.test(yield,trt,model$gl.b,model$Eb)
# simple effects
out<-by(yield,trt, function(x) LSD.test(x, clon, model$gl.c, model$Ec))
out
Run the code above in your browser using DataLab