¿Sabe alguien si hay algún código disponible en R (u otro lenguaje) de código abierto que implemente la estimación no paramétrica de la función de densidad del tiempo de fallo sujeta a una restricción de monotonicidad en la función de densidad del tiempo de fallo (ejemplo: densidad decreciente) utilizando datos de tiempo de fallo censurados a la derecha?
Respuesta
¿Demasiados anuncios?La parte censurada a la derecha sólo le permitirá estimar una probabilidad de cola derecha, pero obviamente no le permite estimar la forma de la densidad en la cola. Así que, en realidad, todo lo que hay que hacer es ignorar los datos censurados y obtener una estimación de la densidad dentro de la parte no censurada de la distribución. No estoy seguro de por qué querrías imponer la restricción de que la densidad debe ser decreciente, pero si quieres hacer esto podrías hacerlo usando el scdensity package
(véase la documentación aquí ).