Hola: On Mon, 7 Dec 2015 15:12:24 +0100 Jes__s Para Fern__ndez <j.para.fernan...@hotmail.com> wrote:
> Buenas, > > Como pudeo calcular el tiempo de vida? Os cuento, tengo una serie de > cuchillas y quiero ver el consumo de las mismas y he pensado en hacer un > estudio por tiempo de vida. No se como hacerlo con R Has una tabla de datos con 4 columnas: 1. Fecha que se empieza a utilizar la cuchilla 2. Fecha de la última revisión de la cuchilla 3. TIEMPO: el tiempo transcurrido entre las dos fechas anteriores 4. ESTADO: estado de la cuchilla cuando se revisó la última vez: 0 = Buen estado 1 = Para tirar Y utiliza el método kaplan-meier: library(survival) KM <- survfit(Surv(DATO$TIEMPO, DATOS$ESTADO) ~ 1) summary(KM) plot(KM) Un ejemplo: > library(survival) > AML.KM <- survfit(Surv(aml$time, aml$status) ~ 1) > summary(AML.KM) Call: survfit(formula = Surv(aml$time, aml$status) ~ 1) time n.risk n.event survival std.err lower 95% CI upper 95% CI 5 23 2 0,9130 0,0588 0,8049 1,000 8 21 2 0,8261 0,0790 0,6848 0,996 9 19 1 0,7826 0,0860 0,6310 0,971 12 18 1 0,7391 0,0916 0,5798 0,942 13 17 1 0,6957 0,0959 0,5309 0,912 18 14 1 0,6460 0,1011 0,4753 0,878 23 13 2 0,5466 0,1073 0,3721 0,803 27 11 1 0,4969 0,1084 0,3240 0,762 30 9 1 0,4417 0,1095 0,2717 0,718 31 8 1 0,3865 0,1089 0,2225 0,671 33 7 1 0,3313 0,1064 0,1765 0,622 34 6 1 0,2761 0,1020 0,1338 0,569 43 5 1 0,2208 0,0954 0,0947 0,515 45 4 1 0,1656 0,0860 0,0598 0,458 48 2 1 0,0828 0,0727 0,0148 0,462 > plot(AML.KM) Dibuja la curva de supervivencia con el intervalo de confianza. Espero que te sea útil. Saludos. > > Gracias > Jes_s > > [[alternative HTML version deleted]] > _______________________________________________ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es