r - 对两点之间的核密度图进行着色。

Question

我经常使用核密度图来说明分布。这些可以在 R 中轻松快速地创建，如下所示:

set.seed(1)
draws <- rnorm(100)^2
dens <- density(draws)
plot(dens)
#or in one line like this: plot(density(rnorm(100)^2))

这给了我这个漂亮的小 PDF:

我想对 PDF 下第 75 到 95 个百分位数的区域进行阴影处理。使用 quantile 函数可以轻松计算点:

q75 <- quantile(draws, .75)
q95 <- quantile(draws, .95)

但是如何对 q75 和 q95 之间的区域进行阴影处理？

Answer 1

使用polygon()函数，请参阅它的帮助页面，我相信我们在这里也有类似的问题。

您需要找到分位数值的索引才能获得实际的 (x,y) 对。

编辑:给你:

x1 <- min(which(dens$x >= q75))  
x2 <- max(which(dens$x <  q95))
with(dens, polygon(x=c(x[c(x1,x1:x2,x2)]), y= c(0, y[x1:x2], 0), col="gray"))

输出(由 JDL 添加)