参考GGPLOT的输入数据，并将其用于GEOM中的自定义功能

Question

我正在使用GGPLOT GEOM_VLINE与自定义函数结合使用，以在直方图顶部绘制某些值。

示例函数下面的函数返回三个值的向量（平均值和x SD在平均值以下或高于平均值）。现在，我可以在GEOM_VLINE（Xintercept）中绘制这些值，并在我的图表中查看它们。

#example function
sds_around_the_mean <- function(x, multiplier = 1) {
  mean <- mean(x, na.rm = TRUE)
  sd <- sd(x, na.rm = TRUE)
  
  tibble(low   = mean - multiplier * sd,
         mean  = mean,
         high  = mean + multiplier * sd) %>% 
    pivot_longer(cols = everything()) %>% 
    pull(value)
}

可重复的数据

    #data
set.seed(123)
normal <- tibble(data = rnorm(1000, mean = 100, sd = 5))
outliers <- tibble(data = runif(5, min = 150, max = 200))

df <- bind_rows(lst(normal, outliers), .id = "type")

df %>% 
  ggplot(aes(x = data)) + 
  geom_histogram(bins = 100) + 
  geom_vline(xintercept = sds_around_the_mean(df$data, multiplier = 3),
             linetype = "dashed", color = "red") + 
  geom_vline(xintercept = sds_around_the_mean(df$data, multiplier = 2),
             linetype = "dashed")

问题是，如您所见，我必须在各个地方定义数据$ DF。 当我将任何更改应用于GGPLOT中的原始DF时，这将变得更容易出错，例如在绘制之前过滤离群值。我将不得不在多个地方再次应用相同的更改。

E.g.
df %>% filter(type == "normal")
#also requires 
df$data 
#to be changed to 
df$data[df$type == "normal"] 
#in geom_vline to obtain the correct input values for the xintercept.

因此，我如何首先将DF $数据参数替换为已将其管道的任何内容的列替换为GGPLOT（）？类似于“”。我想操作员。我还尝试了使用GEOM =“ Vline”的Stat_summary来实现这一目标，但没有所需的效果。

Answer 1

您可以将GGPLOT零件包装在卷曲括号中，并在GGPLOT命令中以及计算SDS_AROUND_THE_MEAN时使用。这将使它动态。

df %>% 
  {ggplot(data = ., aes(x = data)) + 
  geom_histogram(bins = 100) + 
  geom_vline(xintercept = sds_around_the_mean(.$data, multiplier = 3),
             linetype = "dashed", color = "red") + 
  geom_vline(xintercept = sds_around_the_mean(.$data, multiplier = 2),
             linetype = "dashed")}