ggplot2 içinde her panel için farklı bir kesişme ile dikey bir çizgi ekleyin

Question

Histogramların panellerini oluşturmak için ggplot2 kullanıyorum ve her grubun ortalamasına dikey bir çizgi ekleyebilmek istiyorum.

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") 

Nasıl her panelin grup x-kesişimi olarak ifade kullanmak alabilirsiniz: Ama geom_vline() her bir panel için aynı kesenini (yani küresel ortalama) kullanıyor? (Aynı zamanda, satırın değerine sahip bir metin etiketi ekleyebiliyorsanız, bonus puanları da ekleyebilirsiniz.)

Answer 1

Veri yolu, elden önce ortalama değerlerle oluşturmaktır.

library(reshape) 
dfs <- recast(data.frame(cat1, cat2, val), cat1+cat2~variable, fun.aggregate=mean) 
qplot(val, data=df, geom="histogram", binwidth=0.2, facets=cat1~cat2) + geom_vline(data=dfs, aes(xintercept=val), colour="red") + geom_text(data=dfs, aes(x=val+1, y=1, label=round(val,1)), size=4, colour="red") 

Answer 2

Bu, @ eduardo'nun gerçekten bir yeniden işlenmesi olduğunu tahmin ediyorum, ancak bir satırda.

ggplot(df) + geom_histogram(mapping=aes(x=val)) 
    + geom_vline(data=aggregate(df[3], df[c(1,2)], mean), 
     mapping=aes(xintercept=val), color="red") 
    + facet_grid(cat1~cat2) 

alt text http://www.imagechicken.com/uploads/1264782634003683000.png

ya da (ggplot yazarı, Hadley tarafından require(plyr) bir paket) plyr kullanılarak:

ggplot(df) + geom_histogram(mapping=aes(x=val)) 
    + geom_vline(data=ddply(df, cat1~cat2, numcolwise(mean)), 
     mapping=aes(xintercept=val), color="red") 
    + facet_grid(cat1~cat2) 

O ben, UŞebeke yönleri üzerinde kesik olmadığından emin yetersiz görünüyor emin değilim neden.