8 votos

¿Estimación de la densidad con una distribución truncada?

Tengo unos datos que están claramente truncados por la izquierda. Deseo ajustarlo con una estimación de la densidad que lo maneje de alguna manera en lugar de intentar suavizarlo.

¿Qué métodos conocidos (como siempre, en R) pueden abordar esta cuestión?

Código de ejemplo:

set.seed(1341)
x <- c(runif(30, 0, 0.01), rnorm(100,3))
hist(x, br = 10, freq = F)
lines(density(x), col = 3, lwd = 3)

enter image description here

Gracias :)

6voto

Eero Puntos 1612

El paquete logspline para R tiene la función oldlogspline que estimará densidades utilizando una mezcla de datos observados y censurados.

5voto

J Lundqvist Puntos 36

En densidad también tiene una función from para indicar el lado más a la izquierda "de la malla en la que debe estimarse la densidad". Siguiendo con el ejemplo anterior:

lines(density(x, from = 0), col = 4, lwd = 3)

Sin embargo, como puede ver, se trata exactamente de la misma distribución sin el botón from como en el caso anterior. Sólo empieza desde 0, eso es todo.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X