[R-br] Extraindo parte de uma formula

Pedro Emmanuel Alvarenga Americano do Brasil emmanuel.brasil em gmail.com
Quarta Novembro 30 15:21:28 BRST 2016


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 em 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 em 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 em 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
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161130/2bd2a080/attachment.html>


Mais detalhes sobre a lista de discussão R-br