Modelo de sobrevivência Weibull com covariáveis variáveis no tempo em R

7

Estou tentando executar um modelo de sobrevivência usando a abordagem Weibull, mas o problema é que tenho covariáveis que variam no tempo. Estou usando o pacote de sobrevivência em R. Minha ligação é:

output <- survreg(Surv(start, stop, fail) ~ gdppc + [...] + cluster(name), data = mydata, dist="weibull")

que gera o seguinte erro:

Error in survreg(Surv(start, stop, fail) ~ gdppclag + : 
  Invalid survival type

O coxphprocedimento funciona bem, mas quero usar o weibull.

Minha primeira pergunta é: a abordagem Weibull pode explicar covariáveis variantes no tempo? Examinei alguns textos e vejo que a abordagem Cox PH pode ser estendida a covariáveis variantes no tempo. É menos claro se a abordagem Weibull pode fazer isso.

Segundo, se de fato o Weibull pode funcionar, quais são os pacotes em R que podem processá-lo?

r survival cox-model weibull

— George
fonte

Obrigado, reformulei / esclarei a questão para colocá-la no tópico.

— George

Eu não ficaria apressado que o coxph funcionasse ou que o sobrevivente não pudesse lidar com isso. A mensagem de erro implica que seu código de falha possui valores estranhos. Você poderia elaborar como codificar iniciar, parar e falhar? Além disso, você leu a vinheta de sobrevivência em covariáveis variáveis no tempo? As vinhetas de sobrevivência são muito boas.

— Wayne

2

O flexurvpacote pode fazer isso: https://cran.r-project.org/web/packages/flexsurv/index.html

Basta chamar a flexsurvregfunção em vez de survreg.

2

Você pode fazer isso com o Survivalpacote da seguinte maneira:

Srv <- Surv(start, stop, fail, type="interval" )

e então você pode usar Srvno seu modelo como:

output <- survreg(Surv(start, stop, fail) ~ gdppc + [...] + cluster(name), data = mydata, dist="weibull")

Espero que isso possa ajudar :)

— Abdal
fonte

1

Isso também pode ser alcançado com o aftregcomando do pacote eha . Para obter a mesma parametrização normalmente usada pelo survregcomando do pacote de sobrevivência e também pelo flexsurvregcomando do pacote flexsurv , a paramopção deve ser configurada "LifeExp"como, conforme explicado na documentação do pacote . Uma adaptação do seu código seria, portanto, a seguinte:

output <- aftreg(Surv(start, stop, fail) ~ gdppc + [...] + cluster(name), data = mydata, dist="weibull", param="lifeExp")

Uma vantagem do eha::aftregcomando é que ele é compatível com o stargazer e, portanto, sua saída pode ser facilmente exportada para o formato LaTeX.

— George Costanza
fonte

Modelo de sobrevivência Weibull com covariáveis ​​variáveis ​​no tempo em R

Modelo de sobrevivência Weibull com covariáveis variáveis no tempo em R