gistfile1.txt

# パッケージ読込
library(dplyr)
library(ggplot2)

# 水準ごとの平均値を算出
xx <- iris %>% 
  group_by(Species) %>% 
  summarise(a=mean(Sepal.Width))

# ボックスプロット(ベース)作成
p <- ggplot(iris, aes(x=Species, y=Sepal.Width))+
  geom_boxplot()

# 横線を水準ごとに引く
p + mapply(function(i,j) geom_segment(x=i-0.4, y=j, xend=i+0.4, yend=j, colour="red"),
         1:nrow(xx), xx$a)


#mapplyなしでプロット(1本だけ。やってることは上と一緒)
p + geom_segment(x=1-0.4,y=xx$a[1],xend=1+0.4,yend=xx$a[1],colour="red")