## Not run:
# dist <- c(rep(1, 500), rep(2, 250), rep(3, 600))
# Distribution <- factor(dist, 1 : 3, c('Unimodal', 'Bimodal', 'Trimodal'))
# x <- c(rnorm(500, 6, 1),
# rnorm(200, 3, .7), rnorm(50, 7, .4),
# rnorm(200, 2, .7), rnorm(300, 5.5, .4), rnorm(100, 8, .4))
# histboxp(x=x, group=Distribution, sd=TRUE)
# ## End(Not run)
Run the code above in your browser using DataLab