[R-es] Tiempo de vida
Griera-yandex
griera en yandex.com
Lun Dic 7 17:47:45 CET 2015
Hola:
On Mon, 7 Dec 2015 16:34:14 +0100
Jesús Para Fernández <j.para.fernandez en hotmail.com> wrote:
> Los datos no son de desgaste de cuchilla, sino de consumo de las mismas.
>
> Por ello tengo los datos de la siguiente forma:
>
> Unidades cambiadas Fecha
>
>
> En unidades cambiadas, suele ser una y en fecha el dia que se hizo el cmabio.
No seria:
TIEMPO: Fecha que se coloca una cuchilla nueva hasta que se cambia.
ESTADO: Siempre 1 (siempre se cambia): no existen censuras.
y un registro para cada cuchilla.
Serviría?
Saludos.
>
> Con eso no se muy bien como estructurar los datos para hacer el análisis.
>
> Gracias
> Jesús
>
>
>
> > Date: Mon, 7 Dec 2015 16:27:18 +0100
> > From: griera en yandex.com
> > To: j.para.fernandez en hotmail.com
> > CC: r-help-es en r-project.org
> > Subject: Re: [R-es] Tiempo de vida
> >
> > 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