Pedro, 

Você pode usar a função update:

update(f$formula, ~ 1)

2016-11-30 14:26 GMT-02:00 Pedro Emmanuel Alvarenga Americano do Brasil via R-br <r-br@listas.c3sl.ufpr.br>:

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

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

--

Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40" S 49° 16' 22" O

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.