R: valores de geom_densidad en el eje Y

Question

¿Por qué geom_density me muestra valores superiores a 1 en el gráfico de densidad? ¿Cómo se puede cambiar a una fracción?

Y mi código utilizado para generar la trama

ggplot(data = input2, aes(x = r.close)) +
  geom_density(aes(y = ..density.., fill = `Próba`), alpha = 0.3, stat = "density", position = "identity") +

  xlab("y") + ylab("density") +
  theme_bw() +
  theme(plot.title=element_text(size = rel(1.6), face = "bold"),
        legend.position = "bottom",
        legend.background = element_rect(colour = "gray"),
        legend.key = element_rect(fill = "gray90"),
        axis.title = element_text(face = "bold", size = 13))

Answer 1

O bien, puede utilizar simplemente el cálculo ..scaled.. valor stat_density proporciona:

library(ggplot2)

set.seed(1)
vals1 <- rbeta(1000, 0.5, 0.1)
vals2 <- rbeta(1000, 0.25, 0.3)

gg <- ggplot(data.frame(x=c(vals1, vals2),
                        grp=c(rep("a", 1000), rep("b", 1000))))
gg <- gg + geom_density(aes(x=x, y=..scaled.., fill=grp), alpha=1/2)
gg <- gg + theme_bw()
gg

Answer 2

Parece que geom_density() muestra los valores adecuados. El área bajo toda esa curva debe ser 1.

Para obtener una estimación de la probabilidad de ciertos valores, tendría que integrar sobre un intervalo en su eje "y", y que nunca debe ser mayor que 1.