对两点之间的核密度图进行阴影处理。

Question

对两点之间的核密度图进行阴影处理。我经常使用内核密度图来说明发行版。这些在R中创建起来很容易，而且速度很快，如下所示：set.seed(1)draws <- rnorm(100)^2dens <- 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

扩大的解决办法：

如果您想要对两个尾部进行阴影(复制和粘贴Dirk的代码)并使用已知的x值：

set.seed(1)draws <- rnorm(100)^2dens <- density(draws)plot(dens)q2     <- 2q65    <- 6.5qn08   <- -0.8qn02   <- -0.2x1 <- min(which(dens$x >= q2))  x2 <- max(which(dens$x <  q65))x3 <- min(which(dens$x >= qn08))  x4 <- max(which(dens$x <  qn02))with(dens, polygon(x=c(x[c(x1,x1:x2,x2)]), y= c(0, y[x1:x2], 0), col="gray"))with(dens, polygon(x=c(x[c(x3,x3:x4,x4)]), y= c(0, y[x3:x4], 0), col="gray"))