violin and CIs

gistfile1.txt

df1 <- data.frame(y = rnorm(100), x = letters[1:2])
df2 <- data.frame(aggregate(. ~ x, df1, mean),
                  upperCI = aggregate(. ~ x, df1, function(x) quantile(x, 0.025))[, 2],
                  lowerCI = aggregate(. ~ x, df1, function(x) quantile(x, 0.975))[, 2])

ggplot() +
  geom_violin(data = df1, aes(x = x, y = y)) +
  geom_errorbar(data = df2, aes(x = x, ymax = upperCI, ymin = lowerCI), width = 0.5) +
  geom_point(data = df2, aes(x, y))