r - 在山脊线图中的所有面板上添加阴影区域

Question

我有一个数据框df：

 set.seed(3467)

 df<- data.frame(method= c(rep("A", 1000), rep("B", 1000), rep("C", 1000)), 
            beta=c(rnorm(1000, mean=0, sd=1),rnorm(1000, mean=2, sd=1.4),rnorm(1000, mean=0, sd=0.5)))

我希望像这样创建一个山脊图：

ggplot(df, aes(x = beta, y = method, color = method, fill = method)) +
geom_density_ridges(size=1)+
scale_fill_manual(values = c("#228b2250", "#0072B250", "#483d8b50")) +
scale_color_manual(values = c("#228b22", "#0072B2", "#483d8b"), guide = "none") +
stat_density_ridges(quantile_lines = TRUE, quantiles = c(0.025, 0.5, 0.975), alpha = 0.6, size=1)

我希望像这样添加一个阴影矩形：

ggplot(df) +
geom_rect(data = data.frame(x = 1),
        xmin = -1, xmax = 1, ymin = -Inf, ymax = Inf,
        alpha = 0.5, fill = "gray") +
geom_density_ridges2(aes(x = beta, fill = method, y = method))

但是，这种尝试是不成功的：

ggplot(df, aes(x = beta, y = method, color = method, fill = method)) +
geom_density_ridges(size=1)+
scale_fill_manual(values = c("#228b2250", "#0072B250", "#483d8b50")) +
scale_color_manual(values = c("#228b22", "#0072B2", "#483d8b"), guide = "none") +
stat_density_ridges(quantile_lines = TRUE, quantiles = c(0.025, 0.5, 0.975), alpha = 0.6, size=1)+
geom_rect(data = data.frame(x = 1),
        xmin = -1, xmax = 1, ymin = -Inf, ymax = Inf,
        alpha = 0.5, fill = "grey") 

该geom_rect()功能似乎不适用于geom_density_ridges()仅geom_density_ridges2().

Answer 1

这是你想要的？这是对您的代码的轻微修改。我认为aes如果某些元素无法识别，我们可以具体指定数据框和每一层。

library(ggplot2)
library(ggridges)

ggplot() +
  geom_rect(data = data.frame(x = 1),
            xmin = -1, xmax = 1, ymin = -Inf, ymax = Inf,
            alpha = 0.5, fill = "gray") +
  geom_density_ridges(data = df, aes(x = beta, y = method, color = method, fill = method),
                      size=1)+
  scale_fill_manual(values = c("#228b2250", "#0072B250", "#483d8b50")) +
  scale_color_manual(values = c("#228b22", "#0072B2", "#483d8b"), guide = "none") +
  stat_density_ridges(data = df, aes(x = beta, y = method, color = method, fill = method),
                      quantile_lines = TRUE, quantiles = c(0.025, 0.5, 0.975), alpha = 0.6, size=1)