#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
ic(means, test=1, df=6) # tukey
#intervals
conf=ic(means, test=1, df=6)
#plot intervals
ic.plot(conf)
#more plot parameters
ic.plot(conf, las=2, bg="cornsilk");grid(10)
Run the code above in your browser using DataLab