#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.plot(r, col=gray.colors(4))
# more graphical parameters
means.plot(r, col=c(2,7,3,5), bg="cornsilk", las=1, cex.names=2,
sub="treatments", family="sans");grid(10)
# plot = 2
means.plot(r, plot=2, col="dark green", bg="gray", las=1, cex.names=2,
sub="Treatments", family="Times", ylab="Hay intake")
# plot = 2 decreasing =FALSE
means.plot(r, plot=2, las=1, cex.names=2, col="red",lty=2,pch=20,cex=1.1,
sub="Treatments", family="Times", ylab="Hay intake", decreasing=FALSE, legend=FALSE);grid(10)
# plot=3
means.plot(r, plot=3, las=1, cex.names=2,
sub="Contrasts (Tukey 0.05)", family="Times", ylab="")
# plot=3 alpha=0.10
means.plot(ea1(data3, design=3), plot=3, las=2, cex.names=2,
sub="Contrasts (Tukey 0.10)", family="Times", ylab="", alpha=0.10, bg="cornsilk");grid(10)
Run the code above in your browser using DataLab