Log-rank a Wilcoxon
umístění v nabídce: Analysis_Survival_Log-rank a Wilcoxon.
Tato funkce poskytuje metody pro porovnání dvou nebo více přežití křivky, kde některé připomínky mohou být cenzurovány a kde celkový sdružení může být rozvrstvené. Metody jsou neparametrické v tom, že nevytvářejí předpoklady o distribuci odhadů přežití.
při absenci cenzury (např. ztráta sledovat, naživu na konci studie) metody zde prezentované snížit na Mann-Whitney (dva ukázkové Wilcoxonův) test pro dvě skupiny z doby přežití a Kruskal-Wallis test pro více než dvě skupiny, z doby přežití. StatsDirect dává komplexní sadu testů pro srovnání přežití data, která mohou být cenzurovány (Tarone a Keramiky, 1977; Kalbfleisch a Prentice, 1980; Cox a Oakes, 1984; Le, 1997).
zde testovaná nulová hypotéza je, že riziko úmrtí / události je ve všech skupinách stejné.
Peto je log-rank test je obecně nejvhodnější metodu, ale Prentice modifikovaný Wilcoxonův test je více citlivý, když poměr rizika je vyšší na počátku doby přežití než na konci (Peto a Peto, 1972; Kalbfleisch a Prentice, 1980). Test log-rank je podobný testu Mantel-Haenszel a někteří autoři jej označují jako test Cox-Mantel (Mantel a Haenszel, 1959; Cox, 1972).
Vrstev,
nepovinná proměnná, vrstev, umožňuje sub-zařadit do skupiny uvedené ve skupině identifikátor proměnné a testovat význam tohoto sub-klasifikace (Armitage a Berry, 1994; Lawless, 1982; Kalbfleisch a Prentice, 1980).
Wilcoxonův závaží
StatsDirect vám dává na výběr ze tří různých váhových metody pro všeobecné Wilcoxonův test, tyto jsou Peto-Prentice, Gehan-Breslow a Tarone-Ware. Ty Peto-Prentice metoda je obecně robustnější než ostatní, ale Gehan statistika je vypočtena běžně mnoho statistických softwarových balíků (Breslow, 1974; Tarone a Keramiky, 1977; Kalbfleisch a Prentice, 1980; Miller, 1981; Hosmer a Lemeshow 1999). Pokud plánujete použít jinou metodu vážení než Peto-Prentice, měli byste vyhledat statistické pokyny.
poměry nebezpečnosti
přibližný interval spolehlivosti pro log hazard-ratio se vypočítá pomocí následujícího odhadu standardní chyby (se):
– kde eij je rozsah expozice na riziko smrti (někdy se nazývá očekávané úmrtí), pro skupiny i k u j-tého odlišné pozorovány čas (někdy nazývané očekávané úmrtí), pro skupiny i k (Armitage a Berry, 1994).
je volitelně uveden přesný podmíněný odhad maximální věrohodnosti poměru rizik. Přesný odhad a jeho interval spolehlivosti (Fisher nebo mid-P) by měly být běžně používány přednostně před výše uvedenou aproximací. Exponenty Cox regresní parametry jsou také přesné odhady nebezpečnosti poměr, ale prosím, na vědomí, že nejsou přesné, pokud Breslow metoda byla použita ke korekci vazby v regresi. Pokud uvažujete o použití Cox regrese, poraďte se se statistikem.
Trend test
Pokud máte více než dvě skupiny, pak StatsDirect vypočítá variantu testu log-rank pro trend. Pokud se rozhodnete nezadávat skóre skupiny, jsou přiděleny jako 1,2,3 … n v pořadí skupin (Armitage and Berry, 1994; Lawless, 1982; Kalbfleisch and Prentice, 1980).
Technické validace
obecné testovací statistika je vypočtena kolem hypergeometrické rozdělení počet událostí v odlišné časy událostí:
– kde váhy wj pro log-rank testu je rovna 1, a wj pro všeobecné Wilcoxonův test je ni (Gehan-Breslow metoda); pro Tarone-Ware metoda wj je odmocnina z ni; a pro Peto-Prentice metoda wj je Kaplan-Meierova survivor funkce vynásobí (ni děleno ni +1). eij je očekávání smrti ve skupině i v JTH zřetelné pozorované době, kdy došlo k dj událostem/úmrtím. nij je číslo ohrožené ve skupině i těsně před J. odlišným pozorovaným časem. Statistika testu pro rovnost přežití napříč skupinami k (populace vzorkované) je přibližně chí-kvadrát distribuován na stupních volnosti k-1. Statistika testu pro monotónní trend je přibližně chí-kvadrát distribuován na 1 stupeň volnosti. c je vektor skóre, které jsou buď definovány uživatelem nebo přiděleny jako 1 až k.
rozptyl se odhaduje metodou, kterou Peto (1977) označuje jako „přesný“.
stratifikovaný test statistika je vyjádřena jako (Kalbfleisch a Prentice, 1980):
– kde statistiky definované výše jsou vypočteny do vrstev, které pak shrnul do vrstev před zobecněné inverzní a provést maticové operace.
příklad
od Armitage and Berry (1994, s. 479).
testovací sešit (pracovní list přežití: skupina fází, čas, cenzor).
následující údaje představují přežití ve dnech od vstupu do studie pacientů s difuzním histiocytárním lymfomem. Porovnávají se dvě různé skupiny pacientů, pacienti se stadiem III a pacienti s onemocněním stadia IV.
Stage 3: 6, 19, 32, 42, 42, 43*, 94, 126*, 169*, 207, 211*, 227*, 253, 255*, 270*, 310*, 316*, 335*, 346*
Stage 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*
* = cenzorovaných dat (pacient stále naživu, nebo zemřela z nesouvisejících příčina)
analyzovat tyto údaje v StatsDirect musíte nejprve připravit je na tři sloupce sešitu, jak je uvedeno níže:
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 pomocí funkce otevření souboru v nabídce Soubor. Poté vyberte Log-rank a Wilcoxon v části Analýza přežití v nabídce analýza. Vyberte sloupec označený „skupina fází“, když budete požádáni o identifikátor skupiny, vyberte „čas“, když budete požádáni o časy, a „cenzor“ pro cenzuru. Na otázku o vrstvách klikněte na tlačítko Zrušit.
příklad:
Logrank a Wilcoxonův testů,
Log Rank (Peto):
Pro skupinu 1 (Fáze skupina = 1)
Pozorované úmrtí = 8
Míry expozice k riziku úmrtí = 16.687031
Relativní rychlost = 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)
Podmíněné maximální věrohodnosti odhady:
Poměr rizik = 0.381485
Fisher Přesný 95% interval spolehlivosti = 0.154582 na 0.822411
Přesné Fisher jednostranné P = 0.0051, oboustranné P = 0.0104
Přesné mid-P 95% interval spolehlivosti = 0.167398 na 0.783785
Přesné mid-P jednostranné P = 0.0034, dva oboustranné P = 0.0068
Všeobecné Wilcoxonův (Peto-Prentice):
testovací statistika:
-5.19836, 5.19836
variance-kovarianční matice:
0.201506 | -4.962627 |
-4.962627 | 4.962627 |
Chi-square pro rovnocennost úmrtnost = 5.44529 P = 0.0196
Oba log-rank a Wilcoxonův testy prokázaly statisticky významný rozdíl v přežití mezi stage 3 a stage 4 pacientů v této studii.
stratifikovaný příklad
z 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 stratifikovaný log-rank test v StatsDirect. Tento příklad je zpracován ve druhém ze dvou klasických článků Richarda Petra a jeho kolegů (Peto et al., 1977, 1976). Vezměte prosím na vědomí, že StatsDirect používá přesnější vzorce rozptylu uvedené v části statistické poznámky na konci Peto et al. (1977).