DATA = rnorm(25,3)
# PLOT BOXPLOT first
par(mar=c(3,3,.5,0.5), oma=c(0,0,0,0), mgp=c(2,1,0))
boxplot(DATA, outpch = NA)
# PLOT POINTS
stripchart(DATA,vertical = TRUE, method = "jitter", pch = 21, col = "maroon", bg = "bisque", add = TRUE)