r - 使用 ggplot2 的两个不同几何体 'coloured' 和 'filled' 的混合图例

Question

我想将三个类别的核密度图放在一个图和垂直线中，以指示整个分布的均值和中值。使用 ggplot2 绘图相当容易，但我一直在努力使传奇正确，这非常非常困难。

set.seed(1234)
data <- data.frame(value = rgamma(n = 10000, shape = 3, scale = 1),
               type = sample(letters[1:3], size = 10000, replace = TRUE))
data$value[data$type == "b"] <- data$value[data$type == "b"] +
                          rnorm(sum(data$type == "b"), mean = 2)
data$value[data$type == "c"] <- data$value[data$type == "c"] +
                          rnorm(sum(data$type == "c"), mean = 4)

# Let's produce a 'coloured' AND 'filled' density plot
# ('cause I want both the area under the curve and the line to be coloured)
library(ggplot2)
gp <- ggplot(data=data, aes_string(x="value"))
gp <- gp + geom_density(aes_string(fill="type", colour="type"), alpha=0.3)
gp

这给了我没有问题的情节和传说。



现在，我添加了垂直线，情节很好，但是图例真的很丑。

# Now let's add vertical lines to the plot, indicating the mean 
# and median for the whole distribution
vlines <- data.frame(mean_median = c(mean(data$value), median(data$value)),
                 labels = c("Mean", "Median"))

gp <- gp + geom_vline(data=vlines,
                  aes(xintercept=mean_median, colour=labels),
                  size=1.05, linetype="dashed", show_guide=TRUE)
gp

传说与“彩色”混淆了 vlines以及“彩色”和“填充”密度。我确实想要两个图例，其中一个是 vlines 的图例(两个条目)，另一个是密度图颜色的图例。

一种解决方法是获取 geom_density要么着色要么填充，但不能两者兼而有之。它工作得更好，但这不是我想要的(因为我想要密度图的面积和线都具有颜色)。它是这样的。而不是这样做:

gp <- gp + geom_density(aes_string(fill="type", colour="type"), alpha=0.3)

我这样做:

gp <- gp + geom_density(aes_string(fill="type"), alpha=0.3)

结果几乎是我想要的，但我真的希望曲线下的面积和密度图的线都被着色。

我一直试图在整个网络上找到解决方案，但似乎没有任何效果。此类问题有时可以使用 show_guide 解决。几何体，因为您可以覆盖是否为单个几何体绘制图例。我玩过它，但它在我的情况下不起作用，因为我需要来自两个几何体(密度和 vlines)的图例，问题是在其中一个几何体中同时具有填充和颜色，以及另一个几何体中的颜色。

Stack Overflow ( bar and line plot in one chart with a legend under ggplot2 ) 中发布了一个非常相似的问题，但是那里使用的解决方案(使用子集)不适用于我的情况。

我真的很感激任何想法。我已经为此苦苦挣扎了很长一段时间，但找不到解决方案。

Answer 1

很可能是一种更简单的方法来做到这一点，但您可以破解绘图对象，作为上述评论的替代方法。你可以做两个情节，一个是你喜欢的情节，一个是你喜欢的图例，然后交换图例。

# Plot with density area and line coloured but legend not right
p1 <- ggplot(data=data, aes(x=value)) + 
            geom_density(aes(fill=type, colour=type), alpha=0.3 ) + 
            geom_vline(data=vlines, aes(xintercept=mean_median, colour=labels), 
                         linetype="dashed", size=1.5, show_guide=TRUE ) 

g1 <- ggplotGrob(p1)

# Plot with density line not coloured but legend is ok
p2 <- ggplot(data=data, aes(x=value)) + 
            geom_density(aes(fill=type), alpha=0.3 ) + 
            geom_vline(data=vlines, aes(xintercept=mean_median, colour=labels), 
                         linetype="dashed", size=1.5, show_guide=TRUE )  +
            guides(fill = guide_legend(override.aes = list(linetype = 0 ))) 

g2 <- ggplotGrob(p2)


# Add legend of second plot to first plot    
g1$grobs[which(g1$layout$name=="guide-box")] <- 
                              g2$grobs[which(g2$layout$name=="guide-box")]  

grid::grid.newpage()    
grid::grid.draw(g1)