在箱线图上使用 stat_summary 到 label 中线

Question

我有一个 function ，其中我正在尝试使用stat_summary()到 plot 中位数的值，就在geom_boxplot()的中线上方。 我已经减少了我的问题并创建了一个玩具示例来简化但保留上下文。

library(ggplot2)

set.seed(20191120)

dat <- data.frame(var   = sample(c("a", "b"), 
                           50, 
                           replace = TRUE),
                   value = rpois(50, 5))

lims <- c(0, 10)

myplot <- function(DATA, YLIMS) {
  ggplot(data  = DATA,
         aes(x = var)) +
    geom_boxplot(aes(y = value),
                 outlier.shape = NA,
                 coef = 0) +
    stat_summary(aes(y = ifelse(value > (YLIMS[2]*0.9),   # if median in top10% of plot window
                                (value - (YLIMS[2]/10)),  # put it below bar
                                (value + (YLIMS[2]/10))), # else put it above
                     label = round(..y.., 2)), #round(median(value), 2))
                 fun.y = median, 
                 geom = "text") +
    coord_cartesian(ylim = YLIMS)
}

myplot(dat, lims)

我的实际地块有几个方面，各种范围，一些中位数位于范围的顶部或底部。 如您所见，我已经排除了胡须和异常值。 这就是 YLIMS 参数用于放大和关注框并排除未使用的 plot 空间的地方。 我已经将这些 YLIMS 值也用于 position 和 label 的 +/- 10% 范围内，效果非常好。

我尝试使用..y.. 值来获取stat_summary(aes())的 label 参数的中值，但它取而代之的是新值。 正如您从 plot 中看到的那样，我们希望两个标签都是“5”，但它们是“6”，因为添加了 10 的 10%。

我还尝试重新计算中位数（如您所见已注释掉），但这需要所有数据的简单中位数，并且不能控制分组/方面/等。

我知道如何重构我的代码以计算为数据中的 y 标签和位置创建值，或者通过聚合和使用箱线图的标识，但我想知道是否有一种方法可以像我的尝试即将完成。

Answer 1

解决这个问题的关键不是试图调整值，而是使用postion=position_nudge()选项来移动位置。

library(ggplot2)
set.seed(20191120)

dat <- data.frame(var   = sample(c("a", "b"), 50, replace = TRUE),
                  value = rpois(50, 5))
lims <- c(0, 10)

myplot <- function(DATA, YLIMS) {
  ggplot(data  = DATA, aes(x = var)) +
    geom_boxplot(aes(y = value), outlier.shape = NA, coef = 0) +
    stat_summary(aes(y = value , label = round(..y.., 2)), 
                 fun.y = median, geom = "text", 
                 position=position_nudge(y = ifelse(value > (YLIMS[2]*0.9), #if median in top 10% of plot window
                                                            (-YLIMS[2]/10), #put it below bar
                                                            (YLIMS[2]/10)), x = 0)) +
    coord_cartesian(ylim = YLIMS)
}

myplot(dat, lims)