Calculadora de sobrevivência de Kaplan–Meier
Estime a curva de sobrevivência a partir de dados de tempo até o evento com censura. Informe o tempo de acompanhamento e o status de evento de cada indivíduo, e a calculadora retorna as probabilidades de sobrevivência de Kaplan–Meier, o valor de S(t) no tempo escolhido e o tempo mediano de sobrevivência.
Estimativa da sobrevivência ao longo do tempo com dados censurados
Fórmula de sobrevivência de Kaplan–Meier
Onde:
- \( t_i \) — os tempos distintos em que ocorreram eventos, até o tempo \( t \)
- \( d_i \) — número de eventos no tempo \( t_i \)
- \( n_i \) — número de indivíduos em risco (ainda sob observação e sem evento) imediatamente antes de \( t_i \)
- \( \hat{S}(t) \) — probabilidade estimada de permanecer sem evento além do tempo \( t \)
O estimador produto-limite multiplica as probabilidades condicionais de sobreviver a cada tempo de evento. Observações censuradas não são eventos; elas apenas saem do conjunto de risco \( n_i \) no seu tempo de censura — é assim que o método lida com o acompanhamento incompleto.
Sobrevivência de Kaplan–Meier – exemplo de cálculo
Exemplo: 10 indivíduos com tempos de acompanhamento 3; 4; 5; 6; 6; 8; 8; 9; 11; 14 e status de evento 1; 1; 0; 1; 0; 1; 1; 0; 1; 0 (6 eventos, 4 censuras), com a sobrevivência avaliada em \( t_0 = 6 \).
\[ \hat{S}(3) = 1 - \tfrac{1}{10} = 0.9 \quad\Rightarrow\quad \hat{S}(4) = 0.9 \left(1 - \tfrac{1}{9}\right) = 0.8 \]
\[ \hat{S}(6) = 0.8 \left(1 - \tfrac{1}{7}\right) = 0.6857 \]
\[ \hat{S}(8) = 0.6857 \left(1 - \tfrac{2}{5}\right) = 0.4114 \quad\Rightarrow\quad \hat{S}(11) = 0.4114 \left(1 - \tfrac{1}{2}\right) = 0.2057 \]
Resultado: \( \hat{S}(6) = 0.6857 \) (68,6% esperados sem evento após o tempo 6). A curva de sobrevivência cai para 0,5 ou menos pela primeira vez em \( t = 8 \); portanto, o tempo mediano de sobrevivência é 8.
Como usar a calculadora de sobrevivência de Kaplan–Meier
- Digite os tempos de acompanhamento de todos os indivíduos, um por linha ou separados por ponto e vírgula ou espaços. Use a mesma unidade (dias, semanas, meses) do início ao fim.
- Digite os valores de status de evento (1 = evento, 0 = censurado) na mesma ordem dos tempos: 1 se o indivíduo sofreu o evento (morte, recidiva, falha de equipamento…), 0 se o acompanhamento terminou sem o evento.
- Defina no campo avaliar sobrevivência no tempo o ponto \( t_0 \) em que você quer a probabilidade de sobrevivência \( S(t_0) \).
- Leia os resultados: \( S(t_0) \), o tempo mediano de sobrevivência, as contagens de indivíduos, eventos e censuras, a tabela completa das estimativas em cada tempo distinto e uma interpretação do resultado.
O que é o estimador de Kaplan–Meier?
O estimador de Kaplan–Meier (produto-limite) é o método não paramétrico padrão para estimar uma função de sobrevivência a partir de dados de tempo até o evento. Aqui, "sobrevivência" significa qualquer evento de interesse — não apenas a morte, mas também recidiva, recuperação, falha de equipamento, cancelamento de clientes e assim por diante. O resultado é uma função em degraus que começa em 1 e só cai nos tempos em que ocorrem eventos; entre um evento e outro, a estimativa permanece constante.
O que torna o método essencial é a censura. Na maioria dos estudos, alguns indivíduos saem antes de o evento acontecer: o estudo termina, um paciente desiste ou uma máquina é aposentada por motivos não relacionados. Você não pode simplesmente contá-los como livres de evento para sempre, porque não se sabe o seu futuro. O método de Kaplan–Meier resolve isso removendo os indivíduos censurados do conjunto de risco no seu tempo de censura, aproveitando a informação de que eles sobreviveram pelo menos até ali.
Sem nenhuma censura, o estimador se reduz à simples proporção de indivíduos ainda sem evento em cada tempo. Com censura, ele continua sem viés, desde que a censura seja independente do processo do evento. A curva também é a grandeza comparada entre grupos pelo teste de log-rank, e \( S(t) \) pode ser lida diretamente: por exemplo, \( \hat{S}(6) = 0.6857 \) significa que uma proporção estimada de 68,6% dos indivíduos permanece sem evento além do tempo 6.
Como interpretar os resultados
Sobrevivência no tempo escolhido \( t_0 \): a probabilidade estimada de estar sem evento após \( t_0 \). Se nenhum evento ocorreu em \( t_0 \) ou antes dele, a estimativa é 1.
Tempo mediano de sobrevivência: o menor tempo observado em que \( \hat{S}(t) \) cai para 0,5 ou menos — o tempo até o qual se espera que metade dos indivíduos tenha sofrido o evento. Se a curva nunca chega a 0,5 durante o acompanhamento, a mediana é "não alcançada"; isso é comum em estudos curtos e com poucos eventos, e não significa que a mediana seja infinita — apenas que o acompanhamento foi curto demais para estimá-la.
A tabela mostra, em cada tempo distinto, o número em risco, os eventos e as censuras, e a atualização de \( \hat{S}(t) \). Note que a estimativa muda apenas nos tempos de evento; tempos com apenas censura aparecem na tabela por transparência, mas não alteram a curva.
Premissas e dicas práticas
- Censura não informativa: os indivíduos censurados devem ter o mesmo risco futuro de evento que os que continuam sob observação. Uma censura ligada ao evento (por exemplo, pacientes que abandonam o estudo porque pioraram) gera viés na estimativa.
- Unidades consistentes e indivíduos independentes: todos os tempos em uma única unidade, uma linha por indivíduo, sem observações duplicadas ou correlacionadas.
- Empates no mesmo tempo: se eventos e censuras ocorrem no mesmo tempo, a convenção (adotada aqui) é que as censuras saem do conjunto de risco depois dos eventos daquele tempo — os indivíduos censurados ainda contam como em risco quando os eventos do mesmo instante são processados.
- Amostras pequenas: com poucos eventos, a função em degraus é grosseira e \( S(t) \) fica imprecisa; a mediana só pode ser estimada quando a curva cruza 0,5.
- O estimador é descritivo para um único grupo. Para comparar formalmente duas curvas de sobrevivência, use um teste de log-rank; para resumir a ocorrência de eventos entre dois grupos, veja a calculadora de razão de chances e risco relativo. Quando se prefere um modelo paramétrico, os tempos de sobrevivência costumam ser modelados com a distribuição de Weibull ou a distribuição exponencial; para comparações de resultados categóricos, a calculadora de teste qui-quadrado é a ferramenta padrão.
Perguntas frequentes
A ordem em que eu digito os indivíduos importa?
Não. A calculadora ordena todas as observações pelo tempo de acompanhamento antes de calcular qualquer coisa, então você pode listar os indivíduos em qualquer ordem — desde que cada tempo esteja pareado com o seu próprio status de evento na mesma posição dentro da lista. Só o pareamento entre as duas listas precisa estar correto.
Posso digitar dados agrupados ou resumidos em vez de uma linha por indivíduo?
Não, a calculadora precisa de um tempo de acompanhamento e um status por indivíduo. Contagens agregadas como "12 eventos no mês 6" não podem ser digitadas diretamente. Se os seus dados estiverem agrupados, seria preciso reconstruir as linhas individuais (por exemplo, digitar 12 vezes o mesmo tempo com status 1), mas saiba que o agrupamento costuma perder informação sobre os tempos exatos dos eventos e as censuras.
Por que a calculadora não mostra intervalos de confiança para S(t)?
Intervalos de confiança exigem métodos adicionais, mais comumente a fórmula de Greenwood para o erro padrão de \( \hat{S}(t) \), muitas vezes em escala log-menos-log para manter os limites entre 0 e 1. Sem o intervalo, lembre-se de que \( \hat{S}(t) \) é apenas uma estimativa pontual: com poucos eventos ela pode ser bem imprecisa, então um valor como 0,70 não deve ser lido como exato.
Qual é a diferença entre tempo mediano e tempo médio de sobrevivência?
A mediana é o tempo em que \( \hat{S}(t) \) cai pela primeira vez para 0,5 — espera-se que metade dos indivíduos tenha sofrido o evento até ali. O tempo médio de sobrevivência é a área sob a curva de sobrevivência, mas não pode ser estimado de forma confiável a menos que o acompanhamento continue até quase todos os indivíduos terem sofrido o evento. É por isso que a mediana é o resumo padrão na análise de sobrevivência.
De quantos indivíduos ou eventos eu preciso para uma estimativa confiável de Kaplan–Meier?
Não existe um mínimo fixo, mas a precisão depende principalmente do número de eventos, e não apenas do número de indivíduos. Com menos de aproximadamente 10–15 eventos, a curva é uma escadaria grosseira e a mediana pode ficar fora de alcance. Como regra prática, o erro padrão de \( \hat{S}(t) \) diminui aproximadamente com a raiz quadrada do número de eventos, então mais eventos significam uma curva visivelmente mais estável.
Posso converter os resultados para outra unidade de tempo?
Sim. Se todos os tempos digitados usarem uma unidade única e consistente, todas as saídas (o tempo escolhido \( t_0 \), a mediana e os tempos da tabela) estarão nessa mesma unidade. Para trocar de unidade, converta os tempos de entrada antes de digitá-los ou converta os tempos de saída depois — por exemplo, uma mediana de 540 dias equivale a 78,6 semanas (540 ÷ 7). As probabilidades de sobrevivência em si não dependem de unidade e nunca mudam.
Quais são os erros mais comuns ao digitar dados de sobrevivência?
Os erros típicos são: listar tempos e status em ordens diferentes (o pareamento precisa corresponder), usar valores diferentes de exatamente 0 ou 1 como status, misturar unidades como dias e meses no mesmo conjunto de dados e contar por engano um indivíduo censurado como evento. Verifique também linhas duplicadas — o mesmo indivíduo digitado duas vezes infla o conjunto de risco e eleva \( \hat{S}(t) \) com viés para cima.