I made the following plot in R :
library(MASS) a = rnorm(100, 10, 10) b = rnorm(100, 10, 5) c = rnorm(100, 5, 10) group <- sample( LETTERS[1:4], 100, repla