bar.group {agricolae} | R Documentation |
It plots bars of the averages of treatments to compare. It uses the objects generated by a procedure of comparison like LSD, HSD, Kruskall, Waller-Duncan, Friedman or Durbin. It can also display the 'average' value over each bar in a bar chart.
bar.group(x, horiz = FALSE, ...)
x |
Object created by a test of comparison |
horiz |
Horizontal or vertical bars |
... |
Parameters of the function barplot() |
x: data frame formed by 5 columns: name of the bars, height and level of the bar.
x |
Data frame |
horiz |
Logical TRUE or FALSE |
...
Felipe de Meniburu
LSD.test
, HSD.test
, kruskal
, friedman
, durbin.test
, waller.test
# Example 1 library(agricolae) data(sweetpotato) model<-aov(yield~virus,data=sweetpotato) df<-df.residual(model) MSerror<-deviance(model)/df attach(sweetpotato) comparison<- LSD.test(yield, virus,df,MSerror,alpha=0.01,group=TRUE) detach(sweetpotato) #startgraph par(cex=1.5) bar.group(comparison,horiz=TRUE,density=8,col="blue",border="red",xlim=c(0,50)) title(cex.main=0.8,main="Comparison between\ntreatment means",xlab="Yield", ylab="Virus") #endgraph # Example 2 library(agricolae) x <- 1:4 y <- c(0.29, 0.44, 0.09, 0.49) xy <- data.frame(x,y,y) #startgraph par(cex=1.5) bar.group(xy,density=30,angle=90,col="brown",border=FALSE,ylim=c(0,0.6),lwd=2,las=1) #endgraph