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

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&nbsp;<-&nbsp;rnorm(100)^2dens&nbsp;<-&nbsp;density(draws)plot(dens)q2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<-&nbsp;2q65&nbsp;&nbsp;&nbsp;&nbsp;<-&nbsp;6.5qn08&nbsp;&nbsp;&nbsp;<-&nbsp;-0.8qn02&nbsp;&nbsp;&nbsp;<-&nbsp;-0.2x1&nbsp;<-&nbsp;min(which(dens$x&nbsp;>=&nbsp;q2))&nbsp;&nbsp;x2&nbsp;<-&nbsp;max(which(dens$x&nbsp;<&nbsp;&nbsp;q65))x3&nbsp;<-&nbsp;min(which(dens$x&nbsp;>=&nbsp;qn08))&nbsp;&nbsp;x4&nbsp;<-&nbsp;max(which(dens$x&nbsp;<&nbsp;&nbsp;qn02))with(dens,&nbsp;polygon(x=c(x[c(x1,x1:x2,x2)]),&nbsp;y=&nbsp;c(0,&nbsp;y[x1:x2],&nbsp;0),&nbsp;col="gray"))with(dens,&nbsp;polygon(x=c(x[c(x3,x3:x4,x4)]),&nbsp;y=&nbsp;c(0,&nbsp;y[x3:x4],&nbsp;0),&nbsp;col="gray"))