根据组调整ggplot中geom_vline的颜色

Question

因此，我正在为两种虹膜数据的sepal.engength创建密度图。图形如下：

如您所见，我在图表中添加了两个中心趋势的度量（中值和均值）。但是，目前它们会根据度量（即中位数=黑色，平均=灰色）进行颜色。

我想更改颜色，以便setosa/versicolor的含义和中位数具有相同的颜色。此外，我想为setosa/versicolor的平均值/中位数添加一个指示颜色和线型的传奇。因此，最后，我将拥有一个四个部分的传奇：中值setosa，平均setosa，中值versicolor，平均versiocolor。

有人知道该怎么做吗？请参阅下面的复制代码：

library(ggplot2)
library(tidyverse)

iris <- iris %>%
  filter(Species == c("setosa", "versicolor"))

temp <- iris %>% 
  group_by(Species) %>%
  dplyr::summarize(Mean = mean(Sepal.Length, na.rm=TRUE))

temp_2 <- iris %>% 
  group_by(Species) %>%
  dplyr::summarize(Median = median(Sepal.Length, na.rm=TRUE))

ggplot(iris, aes(x = Sepal.Length, fill = factor(Species))) +
  geom_density(alpha = 0.5) +
  theme_minimal() +
  geom_vline(data=temp, aes(xintercept=Mean, color="Mean"),
             linetype="dotted", 
             lwd=1) +
  geom_vline(data=temp_2, aes(xintercept=Median, color="Median"),
             linetype="dashed", 
             lwd=1) +
  scale_color_manual(name = "Statistics", values = c(Median = "black", Mean = "grey50")) + 
  labs(title = "Distribution of Sepal.Length",
       x = "Sepal.Length",
       y = "Density",
       fill = "Species") 

Answer 1

这个比较复杂。您需要将线条的线型和颜色映射到 Species 和表示中位数或平均值的字符串的交互，而不是适当地指定手动比例：

ggplot(iris, aes(x = Sepal.Length, fill = factor(Species))) +
  geom_density(alpha = 0.5) +
  theme_minimal() +
  geom_vline(data=temp, aes(xintercept=Mean, 
                            color = interaction(Species, 'Mean'),
                            linetype = interaction(Species, 'Mean')),
             lwd = 1, key_glyph = draw_key_path) +
  geom_vline(data=temp_2, aes(xintercept=Median, 
                              color = interaction(Species, 'Median'),
                              linetype = interaction(Species, 'Median')),
             lwd=1,  key_glyph = draw_key_path) +
  scale_linetype_manual(values = c(setosa.Mean = 'dotted', 
                                   setosa.Median = 'dashed',
                                   versicolor.Mean = 'dotted',
                                   versicolor.Median = 'dashed'),
                        name = 'Averages') +
  scale_color_manual(values = c(setosa.Mean = "#F8766D", 
                                   setosa.Median = "#F8766D",
                                   versicolor.Mean = "#00BFC4",
                                   versicolor.Median = "#00BFC4"),
                     name = 'Averages') +
  labs(title = "Distribution of Sepal.Length",
       x = "Sepal.Length",
       y = "Density",
       fill = "Species") +
  theme(legend.key.width = unit(15, 'mm'))