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]]

