Uso con frecuencia diagramas de densidad del kernel para ilustrar las distribuciones. Estos son fácil y rápida de crear en I de este modo:Sombreando un gráfico de densidad de kernel entre dos puntos.
set.seed(1)
draws <- rnorm(100)^2
dens <- density(draws)
plot(dens)
#or in one line like this: plot(density(rnorm(100)^2))
que me da este pequeño PDF:
me gustaría dar sombra al área bajo la PDF a partir de los 75 hasta percentiles 95. Es fácil de calcular los puntos usando la función quantile
:
q75 <- quantile(draws, .75)
q95 <- quantile(draws, .95)
Pero, ¿cómo hacer que la sombra de la zona comprendida entre q75
y q95
?
pueden hacerle ejemplo de sombreado de la parte exterior de su rango contra el interior de su rango? Gracias. – Milktrader