Boilerplate code for making dynamite plots in ggplot2.

library("ggplot2")

ggplot(iris) +
  aes(x = Species, y = Sepal.Width, fill = Species) +
  stat_summary(aes(width = .8), geom = "bar", fun.y = "mean") +
  stat_summary(geom = "errorbar", fun.data = "mean_se", width = .2)