Ben sık sık kernel yoğunluk araziler dağılımları göstermek için kullanın. Bu ve R gibi çok oluşturmak için kolay ve hızlı

set.seed(1)
draws <- rnorm(100)^2
dens <- density(draws)
plot(dens)
#or in one line like this: plot(density(rnorm(100)^2))

Beni bu güzel küçük PDF verir:

75 için 95 yüzdelik gelen PDF altında kalan alan gölge etmek istiyorum. Kolay puan quantile işlevi kullanarak hesaplamak için:

q75 <- quantile(draws, .75)
q95 <- quantile(draws, .95)

Ama nasıl q75 q95 arasındaki alan gölge?

polygon() fonksiyonu ile, yardım sayfasını görmek ve benzer soruları biz de burada olduğuna inanıyorum.

(x,y) gerçek çift almak için quantile değerleri dizini bulmak gerekir.

Düzenleme:Al bakalım:

x1 <- min(which(dens$x >= q75))  
x2 <- max(which(dens$x <  q95))
with(dens, polygon(x=c(x[c(x1,x1:x2,x2)]), y= c(0, y[x1:x2], 0), col="gray"))

Çıkış (JDL tarafından eklenmiştir)