r - 如何在密度图中添加来自 geom_vline 的 x 截距作为标签?
问题描述
我正在创建一个密度图并对其进行刻面。假设我有 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))
解决方案
一种方法是预先计算方法并使用这些方法来馈送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)