Question

如何使geom_boxplot异常值与抖动的geom_points完美叠加？

例如，我希望geom_boxplot中的异常值在抖动后从geom_point显示为“十字线”？

library(ggplot2)
p <- ggplot(mtcars, aes(factor(cyl), mpg)) + 
  geom_boxplot(outlier.shape=10, outlier.size=8)  +
  geom_point(aes(factor(cyl), mpg, color=mpg),  position="jitter", size=4)
p

plot

非常感谢任何协助。

Answer 1

我同意Didzis的观点，那就是完全你想要的目标的解决方案将会相当参与。按字面意思做你的建议需要（我认为）你在ggplot之外进行抖动和异常值计算。如果您对如何突出异常值很灵活，这可能是一个更短的解决方案：

id_outliers <- function(x){
    q <- quantile(x,c(0.25,0.75))
    iqr <- abs(diff(q))
    ifelse((x < q[1] - 1.5*iqr) | (x > q[2] + 1.5*iqr),'Outlier','NotOutlier')
}

mtcars <- ddply(mtcars,
                .(cyl),
                transform,
                out = id_outliers(mpg))

p <- ggplot(mtcars, aes(factor(cyl), mpg)) + 
  geom_boxplot(outlier.colour = NA)  + 
  geom_point(aes(colour = mpg,shape = out),position = "jitter")

Answer 2

这个解决方案将会很长。问题是，使用position="jitter"，您无法获得点的精确坐标，因此需要找到解决方法。

因此，请使用原始图表并使用ggplot_build()保存。第一个数据元素包含有关箱图的信息。我们对列group和outliers感兴趣，因为它显示了ggplot假定为异常值的值。将它们保存为单独的对象。

p <- ggplot(mtcars, aes(factor(cyl), mpg)) + 
                geom_boxplot(outlier.shape=10, outlier.size=8)  +
                geom_point(aes(color=mpg),  position="jitter", size=4)
gg<-ggplot_build(p)
gg$data[[1]]
  ymin lower middle upper ymax         outliers notchupper notchlower x PANEL group weight ymin_final
1 21.4 22.80   26.0 30.40 33.9                    29.62055   22.37945 1     1     1      1       21.4
2 17.8 18.65   19.7 21.00 21.4                    21.10338   18.29662 2     1     2      1       17.8
3 13.3 14.40   15.2 16.25 18.7 10.4, 10.4, 19.2   15.98120   14.41880 3     1     3      1       10.4
  ymax_final  xmin  xmax
1       33.9 0.625 1.375
2       21.4 1.625 2.375
3       19.2 2.625 3.375

xx<-gg$data[[1]][c("group","outliers")]
xx
  group         outliers
1     1                 
2     2                 
3     3 10.4, 10.4, 19.2

现在将group值更改为4,6和8，使其与cyl值相同。

xx$group<-c(4,6,8)

现在将此新数据框与原始mtcars合并，并另存为新数据框。然后应用函数检查mpg级别outliers中是否列出了cyl值的详细信息。这些值（TRUE和FALSE）保存在out列中。

mtcars.new<-merge(mtcars,xx,by.x="cyl",by.y="group")
mtcars.new$out<-apply(mtcars.new,1,function(x) x$mpg %in% x$outliers)

使用新数据框绘制数据。从geom_boxplot()删除异常值。使用列out确定点的形状和大小。 scale_shape_manual()和scale_size_manual()调整外观。

ggplot(mtcars.new, aes(factor(cyl), mpg)) + 
          geom_boxplot(outlier.shape = NA)  +
          geom_point(aes(color=mpg,shape=out,size=out),  position="jitter")+
          scale_shape_manual(values=c(16,10),guide="none")+
          scale_size_manual(values=c(4,8),guide="none")

enter image description here

如何使geom_boxplot异常值与抖动的geom_points“对齐”？

2 个答案: