[R-es] pregunta

[Javier Marcuzzi]

Estimado Eric

http://cran.r-project.org/doc/contrib/Epicalc_Book.pdf

El resto tendría que leerlo, pero ahora no puedo.

Javier Marcuzzi

El 04/10/2013 04:59 p.m., neo escribió:
> Hola Jose, y podrias enviar algun link al libro o enviar esa parte
> escaneada ?? para minimizar la posibilidad de que haya algo en la
> transcripcion del script ...
>
> Slds, eric.
>
>
> On Fri 04 Oct 2013 01:59:33 PM CLT, daniel wrote:
>> José,
>>
>> La función survfit espera como argumento una fórmula. Al pie del help de la
>> función survfit dice:
>>
>> "Versiones anteriores del código también permiten la especificación de una
>> sola curva al omitir el lado derecho de la fórmula, es decir, ~ 1. El
>> manejo de este caso requiere algunas manipulaciones poco comunes y
>> frágiles, este caso ya no es compatible." (sepan disculpar la traducción)
>>
>> Desconozco el libro que citas y no se si tiene sentido con los datos, pero
>> agregando el lado derecho a la fórmula se encuentra un resultado.
>>
>> survfit(surv.marr ~ 1)
>>
>> Repito no se si tiene sentido, quizás tú lo encuentres o alguien más pueda
>> opinar sobre este uso en este caso.
>>
>> Daniel Merino
>>
>>
>> 2013/10/4 Dr. José A. Betancourt B. <jbetancourt en iscmc.cmw.sld.cu>
>>
>>> En el libro EPICALC (pagina 229-230) en el que está el siguiente script,
>>> todo nos funciona bien, pero cuando vamos a life table, ya allí no avanza,
>>> lo señalamos en el script, por favor quizá se nos haya ido algún detalle,
>>> pero fuimos siguiéndolo por el libro paso a paso y no no hemos percatado
>>>
>>> Todos los de el paquete survival de la ayuda del R funcionan perfectamente
>>>
>>>
>>>
>>> ###########################################
>>>
>>>
>>>
>>> library(survival)
>>>
>>>   library(epicalc)
>>>
>>> data(Marryage)
>>>
>>>   use(Marryage)
>>>
>>> attach(Marryage)
>>>
>>> summ(Marryage)
>>>
>>>
>>>
>>>
>>>
>>> age <- endyr - birthyr
>>>
>>> label.var(age, "Age")
>>>
>>>   summ(age, by = marital)
>>>
>>>
>>>
>>> age.marr <- maryr - birthyr
>>>
>>>   label.var(age.marr, "Age at marriage")
>>>
>>>   summ(.data[,8:9])
>>>
>>>
>>>
>>> married <- marital == "Married"
>>>
>>>   time <- ifelse(married, age.marr, age)
>>>
>>>
>>>
>>> (surv.marr <- Surv(time, married))
>>>
>>> head(data.frame(age, age.marr, married, surv.marr))
>>>
>>> .
>>>
>>>   cbind(age, sex, age.marr, married, surv.marr)[order(time),]
>>>
>>>
>>>
>>> #NO PODEMOS AVANZAR AQUI IT!!!!
>>>
>>> #Life table (there is no way we can make work the following line and the
>>> others depend on it, please help use
>>>
>>>
>>>
>>> fit <- survfit(surv.marr)  #Error in UseMethod("survfit", formula) :   no
>>> applicable method for 'survfit' applied to an object of class "Surv"
>>>
>>>
>>>
>>>   summary(fit, censor=TRUE)
>>>
>>>
>>>
>>>
>>>
>>> km1 <- summary(fit, censor=T)
>>>
>>>   attributes(km1)
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> --
>>>
>>> Este mensaje le ha llegado mediante el servicio de correo electronico que
>>> ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema
>>> Nacional de Salud. La persona que envia este correo asume el compromiso de
>>> usar el servicio a tales fines y cumplir con las regulaciones establecidas
>>>
>>> Infomed: http://www.sld.cu/
>>>
>>>
>>>
>>>
>>>          [[alternative HTML version deleted]]
>>>
>>>
>>> _______________________________________________
>>> R-help-es mailing list
>>> R-help-es en r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>>
>>>
>>
>>
>>
>> _______________________________________________
>> R-help-es mailing list
>> R-help-es en r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es