[R-es] Tiempo de vida
Griera-yandex
griera en yandex.com
Lun Dic 7 16:27:18 CET 2015
Hola:
On Mon, 7 Dec 2015 15:12:24 +0100
Jes__s Para Fern__ndez <j.para.fernandez en 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]]
>
Más información sobre la lista de distribución R-help-es