install.packages("TeachingDemos")

source("http://www.r-statistics.com/wp-content/uploads/2011/01/boxplot-with-outlier-label-r.txt")# Load the function

这里有一个例子是这样的：

set.seed(492)

y <- rnorm(2000)

x1 <- sample(letters[1:2], 2000,T)

x2 <- sample(letters[1:2], 2000,T)

lab_y <- sample(letters[1:4], 2000,T)

# plot a boxplot with interactions: