This function is a wrap around for barplot() that will add error bars. This code was written by Dr. Mike Hammond at the U of A.
eb(dv, fac, ...)
eb(mtcars$mpg, as.factor(mtcars$cyl), ylab = "Miles per gallon", xlab = "cyl")
library(untidydata)
data(spirantization)
eb(spirantization$cIntensity, as.factor(spirantization$position))