[R-es] pregunta

[neo]

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