# NOT RUN {
repart = c(20, 40, 40)
X = data.frame(num = c(rnorm(repart[1], 10, 1), rnorm(repart[2], 11, 1), rnorm(repart[3], 10, 1)),
grp = rep(c("A", "B", "C"), repart))
BoxPlot(X$num, X$grp, data = X, ylab = "num", main = "boxplot with confidence intervals")
# Confidence interval in red with mean in blue.
# }
Run the code above in your browser using DataLab