使用ggplot绘制每个组的平均值facet_wrap

Question

有问题的数据集是钻石（在dplyr中）。我正试图通过使用facet_wrap剪切来突破价格的直方图。我也希望每个地块都有中线。 这就是我所做的 -

by_cut <- group_by(diamonds, cut) 

med <- by_cut %>% 
    summarise(medn = median(price)) 

diam <- diamonds %>% 
    mutate(med_cut = ifelse(cut == med$cut[1], med$medn[1], 
         ifelse(cut == med$cut[2], med$medn[2], 
          ifelse(cut == med$cut[3], med$medn[3], 
            ifelse(cut == med$cut[4], med$medn[4], med$medn[5]))))) 

diam %>% 
    ggplot(aes(price)) + 
    geom_histogram(binwidth = 100) + 
    facet_wrap(~ cut, scales = "free_y") + 
    geom_vline(aes(xintercept= med_cut), colour='red') 

，并获得所需的图（如下图，这是我想要的..） -

不过，我敢肯定这是不是理想的方式实现结果，因此我想知道什么是实现最终情节的最佳方式？

基于@yeedle评论

这里比以前更好的尝试：

diam <- diamonds %>% 
    group_by(cut) %>% 
    mutate(medn = median(as.numeric(price))) 

diam %>% 
    ggplot(aes(price)) + 
    geom_histogram(binwidth = 100) + 
    facet_wrap(~ cut, scales = "free_y") + 
    geom_vline(aes(xintercept= medn, group= cut), colour='red') 

我们能否进一步改善？

Answer 1

diamonds %>% 
    group_by(cut) %>% 
    mutate(med_price = median(as.numeric(price))) %>% 
    ggplot(aes(price)) + 
    geom_histogram(binwidth = 100) + 
    geom_vline(aes(xintercept= med_price), colour='red') + 
    facet_wrap(~ cut, scales = "free_y")