有条件地填充同一组内的密度图

Question

我想创建一个带有凸脊的密度图，其中颜色以超过特定阈值为条件。

目的是说明收入分配情况，同时突出显示超过一定阈值（每周1368澳元）的收入。我在下面提供了我的尝试的简图。

ggplot(subset(vic, totalweekslost > 0 & nwe %in% 100:3000), 
       aes(nwe, fill = nwe > 1368)) + 
  geom_density()

它不仅改变了密度图中单个组中颜色的颜色，而且还按照自己的比例创建了一个新组。

我想要的是使颜色变化在同一组中，并说明清楚的切点。我可以看到我的代码有什么问题，但是我不知道如何创建自己想要的代码。

Answer 1

也许您可以获取geom_density()的值并使用geom_area()。
参考here
这应该可以解决问题：

library(ggplot2)

# Build some data
vic = data.frame(totalweekslost = 100:3000, nwe = 100:3000)

#Yours code
df = subset(vic, totalweekslost > 0 & nwe %in% 100:3000 )    

# Creating density plot
p = ggplot(df, 
           aes(x = nwe)
          ) + 
  geom_density()

# Getting the values of plot
d = ggplot_build(p)$data[[1]]

# Building shaded area
p = p +  geom_area(data = subset(d, x>1368), aes(x=x,y=y), fill = "#00AFBB", alpha = 0.5)

plot(p)