#data3: Kaps and Lamberson (2009): page 347
#Description:
##Latin square design for test four different treatments on hay intake of fattening steers
data(data3)
r<-ea1(data3, design=3)
#plot
contrasts=r[[3]]
contrasts
p.plot(contrasts)
#direct
p.plot(r[3])
# more graphical parameters
p.plot(contrasts, bg="cornsilk", cex=1.5,cex.axis=1.5,
main="P-values of the tukey contrasts",family="sans");grid(10)
data('data7')
r<-ea2(data7,4)
p.plot(r[[5]], bg="cornsilk");grid(10)
par(mfrow=c(1,2))
p.plot(r[[9]][1], xlab="treatment 1", cex=0.5, bg="cornsilk");grid(10)
p.plot(r[[9]][2], xlab="treatment 2", cex=0.5, bg="cornsilk");grid(10)
Run the code above in your browser using DataLab