使用单独的数据框将框图添加到基数R中的折线图

我会这么做.

虽然一个类别的三个x变量的平均值可以放在相同的x坐标(1到10)上，但框图应该在不同的x坐标上绘制，以避免彼此覆盖.这可以使用三个x变量上的循环来完成，或者通过将数据转换为长格式，其中只有一个x变量和另一个变量来指示值属于哪个cat/x.我们将采用后一种方法.

# put all x variables in one column
test_data_long <- reshape(test_data, direction="long", varying=list(2:4), 
                          v.names="x", timevar="xn")
# category and variant of x (xn) together for placing boxplots along the x axis
test_data_long$cat2 <- with(test_data_long, cat + (xn - 2)/5)
# reorder so that boxplots are drawn in the correct order
test_data_long <- test_data_long[order(test_data_long$cat2), ]

第一个命令只是将数据转换为具有单个x变量和表示原始索引x(x1,x2,x3)的xn变量的长格式.

然后我们创建cat2变量，它区分xn和cat，并用作箱线图的x坐标.我们希望有x2个值正好在cat(1，2，...，10)，向左x1个值(0.8，1.8，...，9.8)和右侧的x2个值(1.2，2.2，...，10.2)，这是用cat + (xn - 2) / 5.这个间隔是专门针对三个x变量(max(xn) = 3)，通常它可以是cat + (xn - (max(xn) + 1)/2) / (max(xn) + 2).

接下来，我们只需要根据cat2对数据进行重新排序，以便以我们将提供的x个值的相同顺序绘制框图.

现在是密谋本身.在一个情节中有不同的信息层，我喜欢用透明的方式 suppress 不太重要的信息.

xn_col <- palette()[4:2]
box_col <- xn_col[unique(test_data_long$xn)]
# draw boxplots in semi-transparent colors
par(mar=c(4.5, 4.5, .5, .5))
boxplot(x ~ cat2, test_data_long, at=unique(test_data_long$cat2), 
        boxwex=.15, xaxt='n', xlab='cat',
        border=paste0(box_col, '66'), col=paste0(box_col, '11'))
axis(1, at=unique(test_data$cat))

# compute mean values
data_mean <- aggregate(. ~ cat, test_data, mean)
# draw semi-transparent points under the means to make them pop out more
sapply(data_mean[, -1], points, x=data_mean$cat, 
       pch=16, col='#ffffffcc', cex=2)
# draw mean values
for (cat in 1:3) {
  lines(data_mean$cat, data_mean[[paste0('x', cat)]], type='b', 
        col=xn_col[cat], lwd=1.7, cex=1.2)
}