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