eggs<-c(11,17,16,14,15,12,10,15,19,11,23,20,18,17,27,33,22,26,28)
trt<-c(1,1,1,1,1,2,2,2,2,2,3,3,3,3,4,4,4,4,4)
bplot(y=eggs, x=trt,int="SE",xlab="Treatment",ylab="Mean number of eggs",names.arg=c(1,2,3,4),simlett=TRUE,
lett=c("b","b","b","a"))
Run the code above in your browser using DataLab