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