我正在使用 ggplot2 创建直方图面板,并且我希望能够在每个组的平均值处添加一条垂直线。但是 geom_vline() 对每个面板使用相同的截距(即全局平均值):
require("ggplot2")
# setup some sample data
N <- 1000
cat1 <- sample(c("a","b","c"), N, replace=T)
cat2 <- sample(c("x","y","z"), N, replace=T)
val <- rnorm(N) + as.numeric(factor(cat1)) + as.numeric(factor(cat2))
df <- data.frame(cat1, cat2, val)
# draws a single histogram with vline at mean
qplot(val, data=df, geom="histogram", binwidth=0.2) +
geom_vline(xintercept=mean(val), color="red")
# draws panel of histograms with vlines at global mean
qplot(val, data=df, geom="histogram", binwidth=0.2, facets=cat1~cat2) +
geom_vline(xintercept=mean(val), color="red")
我怎样才能让它使用每个面板的组平均值作为 x 截距?(如果您还可以在具有平均值的行旁边添加文本标签,则可以加分。)