1 ATIVIDADE 2 - ANÁLISE DE SOBREVIVÊNCIA

1.1 Dados do Exercício

1.1.1 Os dados da tabela abaixo referem-se aos tempos de sobrevivência (em dias) de pacientes com câncer submetidos à radioterapia (o numero 0 indica censura).

2 Exercício - Análise de Sobrevivência

2.1 Obtenha: O ajuste dos modelos exponencial e Weibull para o conjunto de dados de radioterapia, compare os modelos pelas estatísticas AIC, BIC e CAIC e pelo teste TRV. Compare o ajuste dos modelos com o Kaplan-Meier.

2.1.1 Obtendo ajustes dos modelos exponencial

## $par
## [1] 432.9823
## 
## $value
## [1] 283.4069
## 
## $counts
## function gradient 
##       11       10 
## 
## $convergence
## [1] 0
## 
## $message
## NULL
## 
## $hessian
##              [,1]
## [1,] 0.0002195435
##             l      AIC     CAIC      BIC
## [1,] 566.8138 568.8138 569.3471 570.7056
##      parametros       EP   tvalue       valorp       LI       LS
## [1,]   432.9823 67.49004 6.415498 2.669135e-08 297.3559 568.6086

2.1.2 Obtendo ajustes dos modelos Weibull

## $par
## [1] 359.271369   1.086988
## 
## $value
## [1] 265.8835
## 
## $counts
## function gradient 
##      160      100 
## 
## $convergence
## [1] 1
## 
## $message
## NULL
## 
## $hessian
##              [,1]         [,2]
## [1,] 0.0002152802  0.001074071
## [2,] 0.0010740706 60.984123223
##             l      AIC    CAIC      BIC
## [1,] 531.7669 535.7669 536.676 539.5506
##      parametros         EP   tvalue       valorp          LI         LS
## [1,] 359.271369 68.1580253 5.271153 1.516755e-06 222.3026894 496.240049
## [2,]   1.086988  0.1280592 8.488173 1.729388e-11   0.8296439   1.344333

3 Sobrevivência estimada versus o Kaplan-Meier

require(survival)
ekm=survfit(Surv(t,censur)~1)
t=seq(0,1500,0.01)
plot(ekm,conf.int=F, xlab="t",ylab="S(t)")
lines(t,s,col='blue', lty=1, lwd=2)
lines(t,s1,col='red', lty=1, lwd=2)

legend(600,1,c('Kaplan-Meier','Weibull','Exponencial'),bty='o',cex=1,
col=c('black','blue','red'), lty=c(1,1,1), lwd=c(2,2,2))

4 Créditos

Material criado por MATEUS ELIAS.