Cox proportional hazards regression
Schat hoe de hazard van een gebeurtenis verandert met een of meer factoren, bij gelijkblijvende overige factoren.
Klinisch voorbeeld
Verschilt de hazard op overlijden tussen twee behandelingen als leeftijd en tumorstadium worden meegenomen?
Wanneer te gebruiken
- De uitkomst is de tijd tot een gebeurtenis, met sommige patiënten gecensureerd.
- U moet corrigeren voor andere variabelen.
- Concurrerende gebeurtenissen zijn zeldzaam of komen niet voor.
Aannames om te controleren
Proportionele hazards
De hazard ratio wordt als constant in de tijd aangenomen. Controleer dat, bijvoorbeeld met Schoenfeld-residuen (cox.zph in R).
Niet-informatieve censurering
Patiënten die gecensureerd zijn, bijvoorbeeld wegens verlies uit follow-up, moeten hetzelfde toekomstige risico hebben als patiënten die onder observatie blijven.
Controleer eerst
- Aanname van proportionele hazards
Wat te rapporteren
- Hazard ratio met 95%-BI voor uw hoofdfactor
- De variabelen waarvoor u hebt gecorrigeerd, en waarom
- Het aantal events, en hoe u proportionele hazards hebt gecontroleerd
Veelgemaakte fouten
Controleer proportionele hazards
Cox regression gaat ervan uit dat de hazard ratio in de tijd constant is. Controleer dat, bijvoorbeeld met Schoenfeld-residuen.
Let op het aantal events
Bij ja/nee- en overlevingsuitkomsten bepaalt het aantal events hoeveel variabelen het model betrouwbaar kan dragen, niet het aantal patiënten.
Maak van tijd-tot-gebeurtenis geen ja/nee
Tellen wie de gebeurtenis aan het eind had, negeert wanneer die plaatsvond en hoe lang elke patiënt gevolgd is.
Alternatieven
| Als uw situatie is… | Gebruik |
|---|---|
| Corrigeren is niet nodig | Kaplan-Meier curves and log-rank test |
| Een concurrerende gebeurtenis komt vaak voor | Competing risks: cumulative incidence and Fine-Gray |
In uw software
# R
library(survival)
fit <- coxph(Surv(time_months, death) ~ treatment + age + stage, data = cohort)
summary(fit)
cox.zph(fit)
COXREG time_months
/STATUS=death(1)
/METHOD=ENTER treatment age stage
/PRINT=CI(95).