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]]
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@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@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > >-- Daniel [[alternative HTML version deleted]]
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