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

Respostas:

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