Log-rank and Wilcoxon
localização do Menu: Analysis_Survival_Log-rank and Wilcoxon.
Esta função fornece métodos para comparar duas ou mais curvas de sobrevivência onde algumas das observações podem ser censuradas e onde o agrupamento global pode ser estratificado. Os métodos são não-paramétricos na medida em que não fazem suposições sobre as distribuições de estimativas de sobrevivência.
In the absence of censorship (e.g. perda de seguimento, viva no final do estudo) os métodos aqui apresentados reduzem a um teste de Mann-Whitney (duas amostras de Wilcoxon) para dois grupos de tempos de Sobrevivência e um teste de Kruskal-Wallis Para mais de dois grupos de tempos de sobrevivência. StatsDirect fornece um conjunto abrangente de testes para a comparação de dados de sobrevivência que pode ser censurado (Tarone e Ware, 1977; Kalbfleisch e Prentice, 1980; Cox e Oakes, 1984; Le, 1997).
a hipótese nula testada aqui é que o risco de morte/evento é o mesmo em todos os grupos.
o teste log-rank de Peto é geralmente o método mais apropriado, mas o teste de Wilcoxon modificado pelo Prentice é mais sensível quando a razão dos perigos é maior nos tempos de sobrevivência precoces do que nos finais (Peto e Peto, 1972; Kalbfleisch e Prentice, 1980). O teste log-rank é semelhante ao teste Mantel-Haenszel e alguns autores referem-se a ele como o teste Cox-Mantel (Mantel e Haenszel, 1959; Cox, 1972).
Estratos
Uma variável opcional, estratos, permite sub-classificar os grupos especificados no grupo, identificador de variável e para testar a significância desta sub-classificação (Armitage e Berry, 1994; Proença, 1982; Kalbfleisch e Prentice, 1980).os pesos do Wilcoxon dão-lhe uma escolha de três métodos diferentes para o teste generalizado do Wilcoxon, estes são o Peto-Prentice, Gehan-Breslow e Tarone-Ware. O método Peto-Prentice é geralmente mais robusto do que os outros, mas a estatística Gehan é calculada rotineiramente por muitos pacotes de software estatístico (Breslow, 1974; Tarone e Ware, 1977; Kalbfleisch e Prentice, 1980; Miller, 1981; Hosmer e Lemeshow 1999). Você deve procurar orientação estatística Se você planeja usar qualquer método de ponderação que não o Peto-Prentice.
razões de perigo
calcula-se um intervalo de confiança aproximado para a razão de perigo log utilizando a seguinte estimativa de erro-padrão (SE)::
– onde eij é a extensão da exposição ao risco de morte (às vezes chamado esperado de mortes) para o grupo i de k no jth distintos de tempo observado (às vezes chamado esperado de mortes) para o grupo i de k (Armitage e Berry, 1994).é facultativamente indicada uma estimativa condicional da probabilidade máxima exacta da razão de perigo. A estimativa exata e o seu intervalo de confiança (Fisher ou mid-P) devem ser usados rotineiramente em preferência à aproximação acima. Os expoentes dos parâmetros de regressão de Cox também são estimadores exatos da taxa de risco, mas por favor note que eles não são exatos se o método de Breslow foi usado para corrigir laços na regressão. Por favor, consulte um estatístico se você está considerando usar regressão Cox.
Teste de tendência
Se você tem mais de dois grupos, então StatsDirect irá calcular uma variante do teste log-rank para a tendência. Se você optar por não inserir as pontuações do grupo, então elas são alocadas como 1,2,3 … n in group order (Armitage and Berry, 1994; Lawless, 1982; Kalbfleisch and Prentice, 1980).
validação Técnica
A geral estatística de teste é calculada em torno de um hypergeometric distribuição do número de eventos em distintos tempos de evento:
– onde o peso wj para o teste log-rank é igual a 1, e wj para a generalização do teste de Wilcoxon é ni (Gehan-Breslow método); para o Tarone-Ware método wj é a raiz quadrada de ni; e para o Peto-Prentice método wj é o de Kaplan-Meier sobrevivente de função multiplicado por (ni dividido por ni +1). eij é a expectativa de morte no grupo i na jth distinct observed time where dj events/deaths occurred. nij é o número em risco no grupo I pouco antes do jth tempo observado distinto. A estatística de teste para a igualdade de sobrevivência entre os grupos k (populações amostradas) é aproximadamente qui-quadrado distribuído em K-1 graus de liberdade. A estatística do teste para a tendência da monótona é aproximadamente qui-quadrado distribuído em 1 grau de liberdade. c é um vetor de pontuações que são definidas pelo usuário ou atribuídas como 1 A K.
variância é estimada pelo método que Peto (1977) se refere como “exato”.
a estatística do ensaio estratificada é expressa como (Kalbfleisch e Prentice, 1980):
– onde as estatísticas acima definidas são calculadas dentro dos estratos e então somadas entre estratos antes das operações de inversa generalizada e transposição da matriz.exemplo de Armitage e Berry (1994, p. 479).
Test workbook (Survival worksheet: Stage Group, Time, Censor).
os dados seguintes representam a sobrevivência nos dias desde o início do ensaio de doentes com linfoma histiocítico difuso. Dois grupos diferentes de pacientes, aqueles com Fase III e aqueles com Fase IV da doença, são comparados.
Fase 3: 6, 19, 32, 42, 42, 43*, 94, 126*, 169*, 207, 211*, 227*, 253, 255*, 270*, 310*, 316*, 335*, 346*
Fase 4: 4, 6, 10, 11, 11, 11, 13, 17, 20, 20, 21, 22, 24, 24, 29, 30, 30, 31, 33, 34, 35, 39, 40, 41*, 43*, 45, 46, 50, 56, 61*, 61*, 63, 68, 82, 85, 88, 89, 90, 93, 104, 110, 134, 137, 160*, 169, 171, 173, 175, 184, 201, 222, 235*, 247*, 260*, 284*, 290*, 291*, 302*, 304*, 341*, 345*
* = dados censurados (paciente ainda vivo ou morreu de outra causa)
Para analisar estes dados em StatsDirect primeiro você deve prepará-los em livro três colunas, como mostrado abaixo:
Stage group | Time | Censor |
1 | 6 | 1 |
1 | 19 | 1 |
1 | 32 | 1 |
1 | 42 | 1 |
1 | 42 | 1 |
1 | 43 | 0 |
1 | 94 | 1 |
1 | 126 | 0 |
1 | 169 | 0 |
1 | 207 | 1 |
1 | 211 | 0 |
1 | 227 | 0 |
1 | 253 | 1 |
1 | 255 | 0 |
1 | 270 | 0 |
1 | 310 | 0 |
1 | 316 | 0 |
1 | 335 | 0 |
1 | 346 | 0 |
2 | 4 | 1 |
2 | 6 | 1 |
2 | 10 | 1 |
2 | 11 | 1 |
2 | 11 | 1 |
2 | 11 | 1 |
2 | 13 | 1 |
2 | 17 | 1 |
2 | 20 | 1 |
2 | 20 | 1 |
2 | 21 | 1 |
2 | 22 | 1 |
2 | 24 | 1 |
2 | 24 | 1 |
2 | 29 | 1 |
2 | 30 | 1 |
2 | 30 | 1 |
2 | 31 | 1 |
2 | 33 | 1 |
2 | 34 | 1 |
2 | 35 | 1 |
2 | 39 | 1 |
2 | 40 | 1 |
2 | 41 | 0 |
2 | 43 | 0 |
2 | 45 | 1 |
2 | 46 | 1 |
2 | 50 | 1 |
2 | 56 | 1 |
2 | 61 | 0 |
2 | 61 | 0 |
2 | 63 | 1 |
2 | 68 | 1 |
2 | 82 | 1 |
2 | 85 | 1 |
2 | 88 | 1 |
2 | 89 | 1 |
2 | 90 | 1 |
2 | 93 | 1 |
2 | 104 | 1 |
2 | 110 | 1 |
2 | 134 | 1 |
2 | 137 | 1 |
2 | 160 | 0 |
2 | 169 | 1 |
2 | 171 | 1 |
2 | 173 | 1 |
2 | 175 | 1 |
2 | 184 | 1 |
2 | 201 | 1 |
2 | 222 | 1 |
2 | 235 | 0 |
2 | 247 | 0 |
2 | 260 | 0 |
2 | 284 | 0 |
2 | 290 | 0 |
2 | 291 | 0 |
2 | 302 | 0 |
2 | 304 | 0 |
2 | 341 | 0 |
2 | 345 | 0 |
Alternatively, open the test workbook usando a função Ficheiro Abrir do menu ficheiro. Em seguida, selecione Log-rank e Wilcoxon da seção de Análise de sobrevivência do menu de análise. Selecione a coluna marcada ” Stage group “quando for solicitado o identificador de grupo, selecione” Time “quando for pedido o times e” Censor ” para censura. Carregue no botão Cancelar quando perguntar sobre os estratos.
Por exemplo:
Logrank e testes de Wilcoxon
Log Rank (Peto):
Para o grupo 1 (Fase de grupo = 1)
Observadas mortes = 8
Extensão da exposição ao risco de morte = 16.687031
taxa Relativa = 0.479414
For group 2 (Stage group = 2)
Observed deaths = 46
Extent of exposure to risk of death = 37.312969
Relative rate = 1.232815
test statistics:
-8.687031, 8.687031
variance-covariance matrix:
0.088912 | -11.24706 |
-11.24706 | 11.24706 |
Chi-square for equivalence of death rates = 6.70971 P = 0.0096
Hazard Ratio, (approximate 95% confidence interval)
Group 1 vs. Group 2 = 0.388878, (0.218343 to 0.692607)
Condicional estimativas de máxima verossimilhança:
Hazard Ratio = 0.381485
Exato de Fisher 95% intervalo de confiança = 0.154582 para 0.822411
Exato de Fisher, de um lado P = 0.0051, dois lados P = 0.0104
meio Exato-P 95% intervalo de confiança = 0.167398 para 0.783785
meio Exato-P um lado P = 0.0034, dois lados P = 0.0068
Generalizada de Wilcoxon (Peto-Prentice):
teste estatístico:
-5.19836, 5.19836
de variância-covariância da matriz:
0.201506 | -4.962627 |
-4.962627 | 4.962627 |
Qui-quadrado para a equivalência de taxas de mortalidade = 5.44529 P = 0,0196
log-rank e Wilcoxon testes demonstraram uma diferença estatisticamente significativa na sobrevivência de experiências entre a fase 3 e fase 4 pacientes deste estudo.
exemplo estratificado
de Peto et al. (1977):
Group | Trial Time | Censorship | Stratum |
1 | 8 | 1 | 1 |
1 | 8 | 1 | 2 |
2 | 13 | 1 | 1 |
2 | 18 | 1 | 1 |
2 | 23 | 1 | 1 |
1 | 52 | 1 | 1 |
1 | 63 | 1 | 1 |
1 | 63 | 1 | 1 |
2 | 70 | 1 | 2 |
2 | 70 | 1 | 2 |
2 | 180 | 1 | 2 |
2 | 195 | 1 | 2 |
2 | 210 | 1 | 2 |
1 | 220 | 1 | 2 |
1 | 365 | 0 | 2 |
2 | 632 | 1 | 2 |
2 | 700 | 1 | 2 |
1 | 852 | 0 | 2 |
2 | 1296 | 1 | 2 |
1 | 1296 | 0 | 2 |
1 | 1328 | 0 | 2 |
1 | 1460 | 0 | 2 |
1 | 1976 | 0 | 2 |
2 | 1990 | 0 | 2 |
2 | 2240 | 0 | 2 |
Censorship 1 = death event
Censorship 0 = lost to follow-up
Stratum 1 = renal impairment
Stratum 2 = no renal impairment
The table above shows you how to prepare data for a teste de log-rank estratificado no StatsDirect. Este exemplo é trabalhado no segundo de dois trabalhos clássicos de Richard Peto e colegas (Peto et al., 1977, 1976). Note – se que StatsDirect utiliza as fórmulas de variância mais precisas mencionadas na secção notas estatísticas no final de Peto et al. (1977).