r - 带状图中按组的点颜色 (col) 和填充颜色 (bg)

Question

我正在尝试重现以下图像 image http://www.davidzeleny.net/wiki/lib/exe/fetch.php/vizualizace:figures:boxplots-jitter-rdbu-colors.png?cache=

我使用的代码大致是这样的:

library(RColorBrewer) 
library(reshape2)

a=rnorm(100, mean=1)
b=rnorm(100, mean=0, sd=1)
ab=data.frame(a,b)
melt=melt(ab)
bpColor=brewer.pal(4, 'RdBu')

boxplot(melt$value ~ melt$variable, notch=T, col=c(bpColor[1], bpColor[4]), outline=F, varwidth=T)
stripchart(melt$value ~ melt$variable, add=T, vertical=T, pch=21,
         bg=bpColor[2:3][melt$variable], method='jitter', jitter=0.02)

除了带状图点的颜色外，我从中得到的几乎相同

my_image http://is.muni.cz/de/256262/Rplot.png

我应该如何编辑我的代码以重现正确的着色？
我认为

bg=bpColor[2:3][melt$variable]

可以完成这项工作，但是我得到了这个输出，如果我要删除 [] 括号，我会得到两种颜色，但在组内混合。提前感谢您的帮助。

Answer 1

不是最优雅的方式，但是嘿，它起作用了

boxplot(melt$value ~ melt$variable, notch=T, col=c(bpColor[1], bpColor[4]), outline=F, varwidth=T)
stripchart(melt[melt$variable == "a", "value"] ~ melt[melt$variable == "a", "variable"], add=T, vertical=T, pch=21, bg=c(bpColor[2]), method='jitter', jitter=0.02)
stripchart(melt[melt$variable == "b", "value"] ~ melt[melt$variable == "b", "variable"], add=T, vertical=T, pch=21, bg=c(bpColor[3]), method='jitter', jitter=0.02)