Henrique, Funcionou. Valeu mesmo... eu sabia que era uma coisa que me faria me sentir um pateta, mas as vezes é preciso uma ajudinha. Abraço, Pedro Brasil Em 30 de novembro de 2016 15:04, Henrique Dallazuanna <wwwhsd@gmail.com> escreveu:
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