Análise de Sobrevivência
Life Tables
&
Kaplan-Meier
Prof. Ivan Balducci
FOSJC / Unesp
Análise de Sobrevivência
Todo organismo vivo acaba morrendo, e a morte é função
do tempo
O tempo faz sempre o mesmo, quer rapidamente quer devagar,
com tudo o que tem vida, por meio de uma ferida, por meio de
uma doença, pelo fogo, pela fome ou por qualquer outra coisa:
em algum momento, toda essa vida florescente será cinza.
Caducidade é o que a cinza exprime
Os sinais sagrados. Romano Guardini. Ed. Quadrante, SP, 1993: A cinza (p.34-35)
Análise de Sobrevivência
Todo organismo vivo acaba morrendo, e a morte é função
do tempo
Tudo se converterá em cinza: a minha casa, a minha roupa, o meu
dinheiro: campo, prados e bosques; a mão com que agora escrevo, os
olhos com que leio, o meu corpo inteiro; as pessoas que amei, que
odiei, que temi. O que sobre a face da terra me pareceu grande e o que
me pareceu pequeno e desprezível – tudo cinza..., tudo...
Os sinais sagrados. Romano Guardini. Ed. Quadrante, SP, 1993: A cinza (p.34-35)
Não há nada que resista ao tempo
Análise de Sobrevivência
•Dados de tempo de sobrevivência
•tempo de sobrevivência é definido como
o tempo para a ocorrência de um evento específico,
que pode ser o desenvolvimento de uma doença,
resposta a um tratamento, reincidência, ou morte.
•Análises de tempo de sobrevivência usados em estudos
biomédicos, Engenharia, Sociologia, Companhia de
Seguros, Marketing etc...
•Por exº: para a Sociologia, o tempo de sobrevivência
pode ser a duração do 1º casamento
Análise de Sobrevivência
A Life Table dá a curva de sobrevivência para um grupo
de indivíduos. Ela pode ser usada para estimar a
sobrevivência ao longo de qualquer tempo da escala.
O método da Life Table dá a probabilidade de morte (ou
outro evento) para cada intervalo de tempo designado.
Há dois métodos para estimar a sobrevivência:
(1) o atuarial, ou Life Table, e (2) o Kaplan-Meier
Análise de Sobrevivência
•Para um grupo de indivíduos, há uma distribuição de
tempos para a morte, isto é, a idade não é a mesma para
todos os indivíduos no grupo.
•A distribuição de mortes pode ser usada para derivar a
probabilidade de mortes antes de uma idade específica
para o grupo.
•A técnica mais comumente usada para determinar a
probabilidade da morte (ou de algum end-point) é a Life
Table.
Análise de Sobrevivência
A palavra sobrevivência é usada, mas outros end-points
além da morte podem ser usados.
A curva de sobrevivência será denotada por S(t), a probabilidade
de sobrevivência t ou mais anos após a entrada no estudo.
Há 4 critérios para calcular uma curva de sobrevivência:
1. Um evento que indica o ponto inicial (exº: início de administrar a
droga;
2. End point (exº: morte, reaparição de uma doença);
3. Entrada de pessoas num estudo pode ser em qualquer tempo
durante o período de follow-up após a 1ª pessoa entrar;
4. Nem todos os pacientes podem ter sido acompanhados durante o
período de tempo.
Análise de Sobrevivência
Os pacientes são acompanhados até algum end-point
Os pacientes que não alcançam o end-point caem
em duas categorias:
1. Perde-se a pista deles por algum motivo;
2. Abandonam o estudo (desistem)
O termo para perdido ou desistente é censurado
Life Tables: Observações Censuradas
•Uma aplicação comum de dados de sobrevivência
são as observações censuradas. Nem todos os
elementos selecionados para o estudo atendem ao
seguimento ao longo do período previamente
estabelecido, podendo, a qualquer momento,
desistir por motivos diversos.
•
observações censuradas ocorrem quando o
dado de interesse não foi registrado:
1. seja devido ao fato de desistência do sujeito em
estudo, ou porque
2. o estudo terminou antes de que o evento ocorresse.
Modelos de sobrevivvência
3 patterns de mortalidade
1. Baixa early in life, Alta later in life
 2. Mortalidade Constante
 3. Alta early in life, Baixa later in life

Análise de Sobrevivência
Objetivo: calcular a probabilidade de
ocorrência de eventos ao longo de certo
período, o qual é dividido em vários intervalos
de tempo
Os eventos podem ser de diversas naturezas, tais como:
a) períodos de remissão de uma doença;
b) tempo de sobrevivência de indivíduos acometidos
por neoplasia;
c) duração de matrimônios;
d) efeitos de drogas e outros.
Life Tables
Life Tables: é uma tabela de distribuição de freqüência aumentada.
É o método mais direto para descrever a sobrevivência em uma amostra.
A distribuição dos tempos de sobrevivência é dividido em um certo nº de
intervalos.
Para cada intervalo podemos calcular:
• nº (%) de casos que entram no respectivo intervalo “vivos”,
• nº (%) ............................................................................“mortos”,
• nº (%) ..........................................................“perdidos” ou “censurados”.
Exemplo de uma Life Table
nºs vivos no
início do
intervalo
“at risk”
ni
nº de
mortos
durante o
intervalo
di
nº de
censurados
durante o
intervalo
ci
Probab.
estimada de
sobreviver
1º Intervalo
0 – 11
15
2
0
0.87
2º Intervalo
12 – 23
13
3º Intervalo
24-35
10
4º Intervalo
36-47
6
Meses após
cirurgia
(i)
Probab. (Si)
estimada de
sobreviver
desde o
início até o
fim do
intervalo i
pi
=
0.87
=
0.73
=
0.73
x
2
1
0.84
x
0
4
1.00
x
3
3
0.33
=
0.24
Fórmulas para Life Tables
pi = (ni-di-ci/2) / ni – ci/2
ni = nº de pacientes vivos no início do intervalo i
pi = probabilidade de sobreviver o intervalo i
di= nº de mortos no intervalo i
ci = nº de censurados no intervalo i
Si = probabilidade de sobrevivência desde o início até o
fim do intervalo(i)
Si = p1 p2p3... pi
Exemplo: dados de sobrevivência p/ Life Table
Survival of Patients Folowing Radical Mastectomy for Breast Cancer
Exemplo
alfa
Tempo (meses)
Tempo (meses)
Paciente
Morte
Censurado
1
6 meses
2
8
3
20
4
20
5
20 meses
6
24
7
25
8
30
9
31
10
37
11
37
12
38
13
40
14
15
42
45
Survival of Patients Folowing Radical Mastectomy for Breast Cancer
Exemplo alfa
Tempo (meses)
Tempo (meses)
Paciente
Morte
Censurado
1
6 meses
2
8
1º intervalo:
3
20
0 - 11 meses
4
20
5
20 meses
•Intervalos escolhidos: 12 meses
•Após 12 meses da cirurgia (0-11 meses):
2 pacientes morreram;
Zero foram censurados, conseqüentemente:
n1 = 15, d1 = 2, c1 = 0.
“d = death”
Survival of Patients Folowing Radical Mastectomy for Breast Cancer
Exemplo alfa
Tempo (meses)
Tempo (meses)
Paciente
Morte
Censurado
1
6 meses
2
8
1º intervalo:
3
20
0 - 11 meses
4
20
5
20 meses
Fórmula para cada intervalo(i):
pi = (ni-di-ci/2) / ni – ci/2
Si = p1
A probabilidade estimada de sobreviver no final do 1º
ano (S1) é: p1 = [(15-2-0/2) / (15-0/2) ] = 0.87
S1= 12meses = 0.87 = 87% é a chance de sobreviver
Survival of Patients Folowing Radical Mastectomy for Breast Cancer
Exemplo alfa
Tempo (meses)
Tempo (meses)
Paciente
Morte
Censurado
1
6 meses
2
8
3
20
4
20
1º intervalo
2º intervalo:
5
20 meses
6
24
12-23 mes
Para o 2º intervalo: sobraram 13 = n2= “at risk”
•Durante o 24º mês da cirurgia (12-23 meses):
2 pacientes morreram;
1 paciente foi censurado, conseqüentemente:
n2 = 13, d2 = 2, c2 = 1.
Survival of Patients Folowing Radical Mastectomy for Breast Cancer
Exemplo alfa
Tempo (meses)
Tempo (meses)
Paciente
Morte
Censurado
1
6 meses
2
8
3
20
4
20
2º intervalo:
5
20 meses
6
24
12-23 mes
pi = (ni-di-ci/2) / ni – ci/2
S2 = p1 x p2
A probabilidade estimada de sobreviver desde o começo
até o final do 2º ano (S2) é: p2 = [(13-2-1/2) / (13-1/2) ] = 0.84
S2= 24 meses = 0.87 x 0.84 = 0.73 = 73%é a chance de viver
Survival of Patients Folowing Radical Mastectomy for Breast Cancer
Exemplo alfa
Tempo (meses)
Tempo (meses)
Paciente
Morte
Censurado
1
6 meses
2
8
3
20
4
20
5
20 meses
6
24
7
25
8
30
9
31
3º
intervalo:
24-35
•Durante o 36º mês da cirurgia (24 – 35 meses):
• Zero pacientes morreu;
• 4 pacientes foram censurados,
• conseqüentemente: n3 = 10, d3 = 0, c3 = 4.
Survival of Patients Folowing Radical Mastectomy for Breast Cancer
Exemplo alfa
Tempo (meses)
Tempo (meses)
Paciente
Morte
Censurado
1
6 meses
2
8
3
20
4
20
5
20 meses
6
24
7
25
8
30
9
31
3º
intervalo:
24-35
p3 = [(10-0-4/2) / (10-4/2) ] = 1.0
A probab. estimada de sobreviver desde o começo até o final do 3º ano S3 é:
S3 = p1 x p2 x p3
S3= 36 meses = 0.87 x 0.84 x 0.73 x 1 = 0.73
Life Table final
Exemplo alfa
nºs vivos no
início do
intervalo
“at risk”
ni
nº de
mortos
durante o
intervalo
di
nº de
censurados
durante o
intervalo
ci
Probab.
estimada de
sobreviver
1º Intervalo
0 – 11
15
2
0
0.87
2º Intervalo
12 – 23
13
3º Intervalo
24-35
10
4º Intervalo
36-47
6
Meses após
cirurgia
(i)
Probab. (Si)
estimada de
sobreviver
desde o
início até o
fim do
intervalo i
pi
=
0.87
=
0.73
=
0.73
x
2
1
0.84
x
0
4
1.00
x
3
3
0.33
=
0.24
Survival of Patients Folowing Radical Mastectomy for Breast
Cancer
Exemplo alfa
Paciente
Morte
1
Censurado
Tempo
Evento
6 meses
6
1
2
8
8
1
3
20
20
1
4
20
20
1
5
20
20
0
6
24
24
0
7
25
25
0
8
30
30
0
9
31
31
0
10
37
37
1
11
37
37
1
12
38
38
0
13
40
40
0
42
1
45
0
14
15
42
45
para
entrada
com os
dados em
um
programa,
por exº:
MINITAB
Exemplo alfa
Resolução via MINITAB for Windows
STAT>> RELIABILITY / SURVIVAL > Nonparametric Distribution Analysis Right Censoring
Specify time intervals: 0 12 24 36 48
Actuarial Table
Interval
Lower
Upper
0.000000
12.0000
12.0000
24.0000
24.0000
36.0000
36.0000
48.0000
Number
Entering
15
13
10
6
Survival
Time Probability
12.0000
0.8667
24.0000
0.7280
36.0000
0.7280
48.0000
0.2427
Number
Failed
2
2
0
3
Number
Censored
0
1
4
3
Standard
Error
0.0878
0.1162
0.1162
0.1664
Conditional
Probability Standard
of Failure
Error
0.1333
0.0878
0.1600
0.1037
0.0000
0.0000
0.6667
0.2222
95.0% Normal CI
Lower
Upper
0.6946
1.0000
0.5002
0.9558
0.5002
0.9558
0.0000
0.5687
Curva de sobrevivência do MINITAB for Windows
Exemplo alfa
Nonparametric Survival Plot for DATA
Actuarial Method
Censoring Column in CENSUR
1.0
Median
0.9
Probability
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0
10
20
30
months
40
41.637
FIM DO MÉTODO LIFE TABLE
INÍCIO DO MÉTODO KAPLAN-MEIER
MÉTODO DE KAPLAN-MEIER
•O método de estimar a sobrevivência é similar à análise atuarial
diferindo apenas em que o tempo que entra no estudo não é
dividido em intervalos para a análise
•É especialmente apropriado em estudos
que envolvem um nº pequeno de pacientes
•A sobrevivência é estimada cada vez que o paciente morre,
assim as desistências são ignoradas na análise
MÉTODO DE KAPLAN-MEIER
•Método de estimar a curva de sobrevivência
que usa o tempo exato da morte
Exemplo
beta
MÉTODO DE KAPLAN-MEIER
Survival of Patients Folowing Radical Mastectomy for Breast Cancer
Tempo (meses)
Tempo (meses)
Paciente
Morte
Censurado
1
6 meses
2
8
3
20
4
20
5
20 meses
6
24
7
25
8
30
9
31
10
37
11
37
12
38
13
40
14
15
42
45
Exemplo beta
MÉTODO DE KAPLAN-MEIER
Survival of Patients Folowing Radical Mastectomy for Breast Cancer
Tempo (meses)
Tempo (meses)
Paciente
Morte
Censurado
1
6 meses
2
8
3
20
4
20
5
20 meses
A primeira morte ocorre 6 meses após a cirurgia, e 14 pacientes estão
ainda vivos: p1 = 14/(14+1) = 0.93
A segunda morte ocorre 8 meses após a cirurgia, e 13 pacientes estão ainda
vivos: p2 = 13/(13+1) = 0.93
Duas mortes (nº s 3 e 4) ocorrem 20 meses após a cirurgia, p3 = 11/(11+2) =
0.85. Observe que o paciente censurado em 20 meses é considerado vivo
para o cálculo a 20 meses, ou seja, ele não é considerado no cálculo de p3.
Exemplo beta
MÉTODO DE KAPLAN-MEIER
St ... Curva de sobrevivência no instante t
A sobrevivência cumulativa é dada pelo produto:
St = p1 p2p3... pi
p1 = 14/(14+1) = 0.93  S1= 6 meses = 0.93
p2 = 13/(13+1) = 0.93  S2 = 8 meses = 0.93 x 0.93 = 0.86
p3 = 11/(11+2) = 0.85.  S3 = 20 meses = 0.93 x 0.93 x 0.85 = 0.72
Exemplo beta
KAPLAN-MEIER Final
Meses
após
cirurgia
nºs vivos
6
14
1
0.93
8
13
1
0.93
ni = at risk
nº de
mortos
di
Si
pi
[1 – (di/ni)]
=
x
=
0.93
0.86
x
20
11
2
0.85
37
4
2
0.67
=
x
=
42
1
1
0.50
=
0.73
0.49
x
0.24
Survival of Patients Folowing Radical Mastectomy for Breast
Cancer
Exemplo beta
Paciente
Morte
1
Censurado
Tempo
Evento
6 meses
6
1
2
8
8
1
3
20
20
1
4
20
20
1
5
20
20
0
6
24
24
0
7
25
25
0
8
30
30
0
9
31
31
0
10
37
37
1
11
37
37
1
12
38
38
0
13
40
40
0
42
1
45
0
14
15
42
45
para entrada
com os
dados em
um
programa,
por exº:
MedCalc;
Statistix;
Statistica;
Minitab, etc...
“statistical’s
software”
STATISTIX 7.0
KAPLAN-MEIER - EXEMP..., 9/12/2003, 5:02:00 PM
KAPLAN-MEIER PRODUCT-LIMIT SURVIVAL DISTRIBUTION
Exemplo beta
TIME VARIABLE:
TIME
EVENT VARIABLE: EVENT
CENTIME
DIED
SORED
6
1
0
8
1
20
AT
RISK
LOWER
UPPER
95% C.I.
S(t)
95% C.I.
SE S(t)
H(t)
15
0.7119
0.9333
0.9875
0.0644
0.0690
0
14
0.6314
0.8667
0.9610
0.0878
0.1431
2
1
13
0.5006
0.7333
0.8830
0.1142
0.3102
24
0
1
10
25
0
1
9
30
0
1
8
31
0
1
7
37
2
0
6
0.2475
0.4889
0.7356
0.1603
0.7156
38
0
1
4
40
0
1
3
42
1
0
2
0.0698
0.2444
0.5824
0.1905
1.4088
45
0
1
1
Exemplo beta
Curva de sobrevivência do STATISTICA for windows
Life Tables vs Kaplan-Meier
Método Life Tables para dados agrupados
Método Kaplan-Meier para dados agrupados e não agrupados
Os dois métodos podem dar diferentes resultados mesmo se
os mesmos dados forem considerados.
Os resultados diferirão se observações censuradas e mortes
ocorrerem no mesmo intervalo. No exemplo beta, duas mortes
e uma observação censurada ocorreu 20 meses após a
cirurgia.
Os resultados irão diferir, portanto, para os dois métodos
após 20 meses.
Life Tables vs Kaplan-Meier
Kaplan-Meier:
- a vantagem sobre o Life Tables é que o resultado da
estimativa não depende do agrupamento dos dados (em um
certo nº de intervalos).
“Os dois métodos são idênticos se os intervalos da
Life Table contém no máximo uma observação”
Exemplo gama
Método de Kaplan-Meier
Patient 1
TIME
EVENT
1
4
1:
DIED
2
4
0:
CENS.
3
10
1:
DIED
4
14
1:
DIED
5
14
0:
CENS.
6
24
1:
DIED
died
Lost to follow-up
Patient 2
Patient
Patient
Patient
Patient
Nº
3
4
5
6
died
died
Lost to follow-up
died
4
10
14
Months Since Enrollment
24
Exemplo gama
Método de Kaplan-Meier
(1)
ni
di
qi
pi = 1- qi
tratamento
(meses)
(2)
Nº alive e
followed
durante
cada time
intervalo
(3)
Nº que
died at
each time
point
Proporção
que died at
that time:
Proporção
que
survived at
that time:
(3)/(2)
1.00-(4)
Cumulative
Survival
4
6
1
.167
.833
.833
10
4
1
.250
.750
?
14
3
1
.333
.667
?
24
1
1
1.00
.000
?
Times to
death
desde o
início do
at risk
(5)
(4)
Si = p1 p2..pI
(6)
Método de Kaplan-Meier
(1)
Times to
death
(3)
desde o
(2)
Nº who
início do
died at
tratamento Nº alive at
each time each time
(meses)
(4)
Proporção
que died
at that
time:
(3)/(2)
4
6
1
.167
10
4
1
.250
14
3
1
.333
(5)
Proporção
que
survived at
that time:
Cumulative
1.00-(4)
Survival
=
x
=
.750
x
.667 =
.833
x
Exemplo gama
(6)
.833
.625
.417
=
24
1
1
1.00
.000
.000
Exemplo gama
Kaplan-Meier Plot (N=6)
% Cumulative Surviving
100
.833
80
.625
60
.417
40
20
0
.0
0
4
10
14
24
Months After Enrollment
Exemplo gama
Nº
TIME
EVENT
1
4
1:
DIED
2
4
0:
CENS.
3
10
1:
DIED
4
14
1:
DIED
5
14
0:
CENS.
6
24
1:
DIED
Forma de
entrada dos
dados nos
Statistical´s
softwares:
Statistix
Minitab
MedCalc, etc...
Exemplo gama
STATISTIX 7.0
KAPLAN-MEIER PRODUCT-LIMIT SURVIVAL DISTRIBUTION
TIME VARIABLE: TIME
EVENT VARIABLE: EVENT
TIME
4
10
14
24
DIED
1
1
1
1
CENSORED
1
0
1
0
AT
RISK
6
4
3
1
LOWER
95% C.I.
0.4688
0.2779
0.1534
0.0000
S(t)
0.8333
0.6250
0.4167
0.0000
UPPER
95% C.I.
0.9659
0.8783
0.7379
0.0000
Termos que devem ser familiares
Life Tables
Kaplan-Meier
Dados censurados
Curvas de sobrevivência
Download

4 - FOSJC UNESP