Tengo un conjunto de datos con muchos ceros que se ve así:
set.seed(1)
x <- c(rlnorm(100),rep(0,50))
hist(x,probability=TRUE,breaks = 25)
Me gustaría dibujar una línea para su densidad, pero la density()
función usa una ventana móvil que calcula valores negativos de x.
lines(density(x), col = 'grey')
Hay density(... from, to)
argumentos, pero estos parecen solo truncar el cálculo, no alterar la ventana para que la densidad en 0 sea consistente con los datos como se puede ver en la siguiente gráfica:
lines(density(x, from = 0), col = 'black')
(si se cambiara la interpolación, esperaría que la línea negra tuviera una densidad mayor a 0 que la línea gris)
¿Existen alternativas a esta función que proporcionarían un mejor cálculo de la densidad en cero?