如何在密度图中添加 geom_vline 的 x 截距作为标签？

Question

我正在创建一个密度图并将其分面。假设我有 4 个变量组，使用该变量对其进行分面将生成 4 个密度图。

对于每个密度图，我放置了一条代表平均值的垂直线。

然而，人们需要仔细观察 x 轴和垂直线之间的交点，才能大致了解平均值是多少。

我想要创建的是每个密度图，我还想将其平均值显示为绘图区域中的标签。

下面的示例代码

x <- rnorm(n = 100, mean = 10, sd = 1)
y <- rnorm(n = 100, mean = 20, sd = 1)
z <- rnorm(n = 100, mean = 40, sd = 1)



df <- as_tibble(cbind(
  c(x,y,z), 
  c(rep('x',length(x)), rep('y',length(y)), rep('z',length(z))),
  c(rep('a',length(x)/2), rep('b',length(x)/2))))

df$V1 <- as.numeric(df$V1)

df <- df %>% group_by(V2, V3) %>%
  summarise(mumean = mean(V1)) %>%
  right_join(df)


df %>%
  ggplot(aes(x = V1, color = V2)) + 
  geom_density(aes(fill = V2)) + facet_grid(V3 ~ V2) + theme_bw() +
  geom_vline(data = df, aes(xintercept = mumean))

Answer 1

一种方法是预先计算平均值并使用它们来提供geom_vline和geom_text：

library(dplyr)
iris_means <- iris %>%
  group_by(Species) %>%
  summarize(mean = mean(Sepal.Length))

ggplot(iris, aes(Sepal.Length)) +
  geom_density() +
  geom_vline(data = iris_means, aes(xintercept = mean)) +
  geom_text(data = iris_means, aes(x = mean, label = mean), 
             y = 0.1, angle = 90, vjust = -0.2) +
  facet_wrap(~Species)