Eu tenho um conjunto de dados com muitos zeros que se parece com isso:
set.seed(1)
x <- c(rlnorm(100),rep(0,50))
hist(x,probability=TRUE,breaks = 25)
Gostaria de desenhar uma linha para sua densidade, mas a density()
função usa uma janela em movimento que calcula valores negativos de x.
lines(density(x), col = 'grey')
Há density(... from, to)
argumentos, mas eles parecem apenas truncar o cálculo, não alterar a janela, para que a densidade em 0 seja consistente com os dados, como pode ser visto no seguinte gráfico:
lines(density(x, from = 0), col = 'black')
(se a interpolação fosse alterada, eu esperaria que a linha preta tivesse maior densidade em 0 do que a linha cinza)
Existem alternativas para essa função que proporcionariam um melhor cálculo da densidade em zero?