對兩點之間的核密度圖進(jìn)行陰影處理。我經(jīng)常使用內(nèi)核密度圖來說明發(fā)行版。這些在R中創(chuàng)建起來很容易，而且速度很快，如下所示：set.seed(1)draws <- rnorm(100)^2dens <- density(draws)plot(dens)#or in one line like this: plot(density(rnorm(100)^2))這給了我一個很好的PDF：我想把PDF下的區(qū)域從第75位到第95位數(shù)進(jìn)行遮陽。使用quantile職能：q75 <- quantile(draws, .75)q95 <- quantile(draws, .95)但我該如何遮擋這片區(qū)域q75和q95?
3 回答
MMMHUHU
TA貢獻(xiàn)1834條經(jīng)驗 獲得超8個贊
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"))
- 3 回答
- 0 關(guān)注
- 678 瀏覽
添加回答
舉報
0/150
提交
取消