As @smillig mentioned, you can achieve this using ggplot2. The code below reproduces the plot, which you then pretty well warn that is pretty complicated. First download the ggplot2 package and create some data:
library(ggplot2)
dd = data.frame(values=runif(21), type = c("Control", "Treated", "Treated + A"))
:
theme_set(theme_bw())
.
- :
g = ggplot(dd, aes(type, values))
: :
g = g + geom_jitter(aes(pch=type), position=position_jitter(width=0.1))
":" , . . , .
g = g + stat_summary(fun.y = function(i) mean(i),
geom="bar", fill="white", colour="black")
: / :
g = g + stat_summary(
fun.ymax=function(i) mean(i) + qt(0.975, length(i))*sd(i)/length(i),
fun.ymin=function(i) mean(i) - qt(0.975, length(i)) *sd(i)/length(i),
geom="errorbar", width=0.2)
g

- R-
stat_summary , . geom_errorbar geom_bar. - R, question.