Amigos de R, Estou tentando fazer uma função com alguns loops que comparam modelos em diferentes subsets de dados. No caminho, eu precisaria extrair a parte dependente das formulas dos modelos, mas não estou sabendo como. O que me interessa é repetir a parte dependente da fomrula nos diferente subsets.
library(survival) f <- coxph(Surv(time, status) ~ ph.ecog + age, data=lung) f$formula Surv(time, status) ~ ph.ecog + age
# Eu gostaria de extrair o "Surv(time, status)" para um objeto de tal forma que eu consiga seguir um modelo igual ao abaixo. objeto <- Surv(time, status) # Esse passo que eu não estou conseguindo. E deveria vir do primeiro fit "f" f2 <- survfit(objeto ~ 1, data=lung, subset = lung$age >= 63 , type = "kaplan-meier") f2b <- survfit(objeto~ 1, data=lung, subset = lung$age < 63 , type = "kaplan-meier") Abraço forte, Pedro Brasil