4
上位10%がローカライズされている領域をシェードしたい。私はちょうどこのプロットをプロットするために切り捨て点65を任意に設定します。それは私が探していることです...すべてのデータセットについてです。密度プロットの上位10%の切り捨て点を見つけよう
xf <- rnorm(40000, 50, 10);
plot(density(xf),xlim=c(0,100), main = paste(names(xf), "distribution"))
dens <- density(xf)
x1 <- min(which(dens$x >= 65)) # I want identify this point such that
# the shaded region includes top 10%
x2 <- max(which(dens$x < max(dens$x)))
with(dens, polygon(x=c(x[c(x1,x1:x2,x2)]), y= c(0, y[x1:x2], 0), col="green"))
abline(v= mean(traitF2), col = "black", lty = 1, lwd =2)
が、これは偉大なコメント答えとして、それを改善してください答えていないようだがと...助けてくれてありがとう。私の評判の一部、冗談だけ – jon
なぜこれはdownvotedされましたか?それは質問に答える! –
質問はあなたが探しているものを知っていれば私にとっては完璧な質問なので、私の短い答えです... – Niclas