Я использую 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-перехват? (Бонусные баллы, если вы также можете добавить текстовую метку в строке со значением среднего.)