means.plot: Plot Means

Description

Plot contrasts of means

Usage

means.plot(data, plot=1, s=1,test=1, legend=TRUE, letters=TRUE,
family="Times", bg="white",cex.names=0.8, cex.text=0.7, cex.legend=1, 
bar.order=2, decreasing=TRUE, alpha=0.05,cex=0.5, pch=19,  ...)

Value

Returns plots and confidence intervals

Arguments

data: output object of ea1 function (see examples)
plot: type of plot 1 = bar plot (default) 2 = means plot 3 = confidence interval of the contrasts
s: s=1 (defalt) plot standard deviation s=2 plot standard error of mean
test: Letters of the post-hoc test 1=Tukey 2=SNK 3=Duncan 4=t 5=Scott-Knott
legend: TRUE = plot p-value of F test FALSE = not plot p-value
letters: TRUE = plot letters FALSE = not plot letters
family: font of plot
bg: background color
cex.names: font size in names of treatments (x-axis)
cex.text: font size in letters and means
cex.legend: font size in legend
bar.order: order of bar or means 1 = order of treatments names 2 = order of the means (default)
decreasing: decreasing bar order (TRUE or FALSE)
alpha: 0.05 (default) is the alpha of confidence intervals
cex: size of points
pch: type of points
...: more plot parameters

Author

Emmanuel Arnhold <emmanuelarnhold@yahoo.com.br>

References

KAPS, M. and LAMBERSON, W. R. Biostatistics for Animal Science: an introductory text. 2nd Edition. CABI Publishing, Wallingford, Oxfordshire, UK, 2009. 504p.

Examples

Run this code


#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