两点ggplot2之间的阴影（刻面）密度图

Question

我正在尝试创建一个阴影两点之间密度图曲线下面积的图。这是我到目前为止所做的：

df.ind <- data.frame(x=rnorm(1000))
df.group <- data.frame(group = factor(1:20))
df.group$gr.mean <- rnorm(20)
df.group$width <- 1 + abs(rnorm(20))
df.group$upper <- df.group$gr.mean + df.group$width
df.group$lower <- df.group$gr.mean - df.group$width

然后，我可以创建一个ggplot，其中我使用来自x的{​​{1}}密度，并将其绘制20个不同的时间（按每个组）。然后我覆盖垂直线。我想要做的是遮蔽线条之间的区域。

df.ind

我在这里发现了一个类似的问题：ggplot2 shade area under density curve by group和此处：Shading a kernel density plot between two points.

但我的数据来自两个不同的ggplot(df.group) + stat_density(data=df.ind, mapping=aes(x), geom="line", position="dodge") + facet_wrap(~group) + geom_vline(aes(xintercept=lower)) + geom_vline(aes(xintercept=upper)) 个对象。因此，我不能使用他们用来聚合数据的聪明技巧......

Answer 1

就像我在评论中提到的那样，您将不得不自己专门为每个面板创建数据。这是一种方式

# calculate density
dx<-density(df.ind$x)

#define areas for each group
mm <- do.call(rbind, Map(function(g, l,u) {
    data.frame(group=g, x=dx$x, y=dx$y, below=dx$x<l, above= dx$x > u)
}, df.group$group, df.group$lower, df.group$upper))

#plot data
p2<-ggplot(mm) + 
   geom_area(aes(x, y, fill=interaction(below,above))) + 
   geom_vline(data=df.group, aes(xintercept=lower)) +
   geom_vline(data=df.group, aes(xintercept=upper)) +
   facet_wrap(~group)