#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
#means
means=r[[2]]
means
m.plot(means, col=gray.colors(4))
#direct
m.plot(r[[2]], col=gray.colors(4))
# more graphical parameters
m.plot(means, col=c(2,7,3,5), bg="white", las=1, cex.text=1,main="Tukey (0.05)",
family="sans", bar.order=2, decreasing=FALSE);grid(10)
data('data7')
r<-ea2(data7,4)
m.plot(r[[4]], col=c(2,7,3,5), las=1, bg="cornsilk");grid(10)
par(mfrow=c(1,2))
m.plot(r[[8]][1], test="scott_knott",xlab="treatment 1",col=c(2,7,3,5),
las=2, bg="cornsilk",bar.order=2, decreasing=FALSE);grid(10)
m.plot(r[[8]][2], test="scott_knott",xlab="treatment 2", col=c(2,7,3,5),
las=2, bg="cornsilk", bar.order=2, decreasing=FALSE);grid(10)
Run the code above in your browser using DataLab