Survival analyse

Auteur	dr. ir. N van Geloven
Co-Auteur	dr. R.B. Geskus
auteurschap op deze site

Survival analyse wordt gebruikt voor data die informatie geeft over de tijd tot het optreden van een bepaald event. Met tijd wordt in deze bedoeld het aantal jaren, maanden of weken vanaf de start van de follow-up van een patient tot aan het optreden van een event. Het event kan overlijden zijn (vandaar de naam survival analyse), maar ook een relapse, herstel of een ander helder gedefinieerd eindpunt. In het algemeen wordt er bij survival analyse vanuit gegaan dat er maar één type risico is en dat een event per individu maar hoogstens één keer kan optreden. Bij meerdere risico's of meerdere herhalingen van een event zijn meer geavanceerde analyse technieken nodig zoals recurrent event analyse of competing risk analyse.

Bij het meten van tijd-tot-event data is vaak sprake van gecensureerde waarnemingen. Hierbij heb je iets aan informatie over de 'survival time' van een patient, maar weet je niet precies hoe lang de tijd-tot-event is. Een voorbeeld. Stel dat in een studie patienten over een periode van twaalf weken gevolgd worden. Censurering (Eng.: censoring) kan op drie manieren optreden:

• Een patient is gevolgd vanaf het begin van de studie tot aan het eind van de 12 weekse studie periode. Op dat moment heeft de patient het event nog niet gekregen. We weten nu alleen dat de tijd-tot-event van deze persoon langer is dan 12 weken.
• Een patient trekt zich terug uit de studie na 6 weken, dus voor het eind van de studieperiode. Van deze persoon weten we alleen dat de tijd-tot event langer is dan 6 weken. Of een patient start aan de studie maar raakt na het laatste contact op 8 weken 'lost to follow up'. Ook van deze patient weten we alleen dat de eventvrije periode langer is dan 8 weken.
• Een patient kan niet tot het einde gevolgd worden doordat hij een ander type event krijgt. Bijvoorbeeld: we onderzoeken tijd tot relapse na operatie, maar een patient overlijdt tijdens de studie zonder een relapse gehad te hebben.

Kaplan Meier analyse

Wanneer gebruik ik een Kaplan Meier analyse?

Om de tijd tot het optreden van het event voor de gehele steekproef te kwantificeren, wordt er vaak een Kaplan Meier analyse gebruikt. De Kaplan Meier analyse is een niet-parametrische aanpak. Dat wil zeggen dat er geen aanname gedaan wordt ten aanzien van het soort verdeling (normaal, gamma etc.) van de event tijden. Toch zitten er wel een paar aannames achter deze methode:

• De prognose van patienten die vroeg of laat zijn geïncludeerd in de studie is gelijk.
• Patienten die gecensureerd zijn op een bepaald tijdstip hebben dezelfde prognose als patienten die op datzelfde tijdstip wel in de studie zijn gebleven (en op dat moment nog geen event hebben gehad). Er wordt niet aan deze aanname voldaan als bijvoorbeeld personen met een grote kans op een event vlak voordat het event optreedt de studie verlaten.
• De tijdstippen waarop de events optreden zijn exact bekend.
• Patienten worden vanaf het begin van de onderzochte periode gevolgd. Als er bijvoorbeeld pas 10 dagen na interventie gestart wordt met het volgen van survival stutus van patienten, dan zullen patienten die in de eerste 10 dagen reeds overleden zijn niet vertegenwoordigd zijn in de studie. Een Kaplan-Meier analyse geeft dan geen goede inschatting van de echte survival tijd. Survival modellen die rekening houden met dit verschijnsel worden late entry survival modellen genoemd.

De basis gedachte achter de Kaplan Meier procedure is het opdelen van de tijdsschaal. Op ieder moment dat er bij één of meerdere patienten een event optreedt, wordt de kans daarop geschat door het aantal events op dat moment te delen door het aantal personen dat op dat moment nog in de studie zit (i.e. niet eerder al een event heeft gehad of gecensureerd is). De survival functie (of cumulatieve proportie van 'survivors') combineert deze schattingen over alle tijdstippen. Op deze manier wordt het resultaat minimaal beïnvloed door patienten die halverwege de rit gecensureerd raken. Een Kaplan-Meier analyse resulteert in een 'life table' en een survival curve.

De volgende life table en survival curve laten bijvoorbeeld zien wat de verwachte overlevingskans is van patienten die een dotterbehandeling hebben ondergaan. Je ziet in de tweede kolom ("Time") van de life table dat er op dag 4 en dag 6 patienten gecensureerd raakten en dat er op dag 0 vijf events waren, op dag 1 een event, op dag 5 twee events en dag 7 een event optrad. Op de dagen dat er een event optrad, is de cumulatieve survival berekend. In de survival curve zie je de geschatte survival functie over de eerste 200 dagen na interventie. Na 200 dagen wordt de kans op overleven geschat op ongeveer 92%.

Hoe test ik of Kaplan Meier survival curves van elkaar verschillen?

We kunnen de curves ook apart schatten en vergelijken voor twee (of meer) subgroepen in onze data. Om te testen of de Kaplan Meier curves van 2 (of meer) groepen significant van elkaar verschillen, wordt vaak de log-rank test gebruikt. Hierbij is de nulhypothese dat de curves niet van elkaar verschillen. De test vergelijkt de geobserveerde verschillen tussen de survival curves met verschillen die zouden kunnen ontstaan als er geen verschil is tussen de groepen. In het volgende voorbeeld concluderen we bijvoorbeeld dat het verschil tussen treatment A en treatment B significant is:

Mijn Kaplan Meier curves kruisen, mag ik dan nog de log-rank test gebruiken?

Bij mijn Kaplan-Meier analyses voor het vergelijken van een factor, heb ik curves die elkaar kruisen. In statistische stukken wordt soms gezegd dat een log-rank test dan niet geschikt is om een potentieel verschil tussen curves te analyseren. Desondanks, als ik met SPSS een log-rank test uitvoer, krijg ik soms alsnog een significante p-value. Houdt de log-rank test (als uitgevoerd door SPSS) rekening met het feit dat de curves elkaar kruisen? Of zou ik eigenlijk geen log-rank mogen gebruiken dan?

De log-rank test houdt geen rekening met het kruisen van curves. De test bekijkt het verschil tussen de curves op ieder (event-)tijdspunt en telt de gevonden verschillen (in detail het verschil tussen observed en expected events) op over alle tijdspunten. Als 2 curves bijvoorbeeld in het midden kruisen, dan zal in de vroege periode de ene groep relatief meer events hebben en in de late periode de andere groep. Bij het optellen van de verschillen zullen deze twee periodes elkaar 'uitmiddelen'. Grote kans dat de log rank statistic dan niet significant is, terwijl er toch interessante verschillen kunnen zijn in de twee periodes apart. Om deze reden wordt de log rank test als minder geschikt (eigenlijk minder krachtig, minder power) gezien bij kruisende curves.

In jouw geval waarbij je kruisende curves hebt en toch een significante log rank test, zal het verschil in een van de twee periodes (voor of na het omslagpunt van de hazard) sterker zijn dan in de andere periode. Er is dan toch 'netto' een verschil 1 kant op gevonden. Dit kan informatief zijn, al moet je goed oppassen hoe zo'n verschil te interpreteren. Zie ook deze ingezonden brief nav van een verkeerde interpretatie in zo'n situatie.

Ik raad je aan iig de plot van de survival curves mede te rapporteren (en niet alleen de log rank statistic), zodat de lezer zelf kan beoordelen waar het verschil te verwachten is. Statistische alternatieven zijn het gebruik van testen die meer gewicht leggen op vroege of juist late periodes. Bijvoorbeeld de Wilcoxon test (in SPSS de Breslow test genoemd), deze test weegt verschillen in het begin van de survival curves zwaarder dan verschillen aan het einde van de curves. De test is bijvoorbeeld geschikt wanneer het effect van een behandeling op survival het sterkst is in de eerste fase en later minder effectief wordt en je wilt inzoomen op de eerste fase. Of de Tarone-Ware test. Deze test geeft net als de Wilcoxon test meer relatief meer gewicht aan verschillen op vroege momenten (maar minder uitgesproken dan de Wilcoxon test). Ook kun je overwegen de analyse apart te doen over twee tijdsintervallen, al is het dan heel lastig om te bepalen waar je de grens trekt. Het liefst zou je dergelijke alternatieve analyses al vooraf hebben gespecificeerd.

Ik heb veel gecensureerde patienten, kan ik de Breslow of Tarone-Ware test gebruiken in plaats van de log-rank test?

De test die ik nu gebruik is logrank. Ik zag dat er ook 2 andere testen zijn (breslow en Tarone-Ware). Het is me niet helemaal duidelijk wanneer je deze gebruikt behalve de vernoeming dat deze toetsen vooral het begin van de curve beoordelen. De reden dat ik mogelijk een andere test wil proberen is het feit dat ik veel censored zie in de curve. De patienten zijn niet allemaal even lang gevolgd en dus neemt het aantal patienten in follow up af naarmate de tijd verstrijkt. Kan/mag de breslow of tarone-ware in dit geval?

De Breslow en Tarone-Ware 'luisteren' inderdaad meer naar de beginperiode wanneer er nog meer personen in follow up zijn dan naar de eindperiode. Echter, het feit dat er later in de follow up veel gecensureerd wordt, is niet een reden om voor deze testen te kiezen. Ook de log-rank test zal namelijk het meest gevoed worden door de periodes waar er de meeste mensen nog in follow up zijn. De Breslow / Tarone-Ware zijn aan te raden wanneer er op basis van een biologische hypothese de verwachting is dat er in de beginperiode de grootste verschillen zullen zijn (bijvoorbeeld bij een eenmalige behandeling waarbij het effect uitvaagt over de tijd). Het zijn niet de standaard testen, dus het gebruik ervan zal duidelijk gemotiveerd moeten worden. Ik raad dus aan hier de 'gewone' log rank te gebruiken, de latere momenten met minder personen zullen vanzelf minder zwaar meetellen in de analyse.

Kan ik meerdere aparte Kaplan Meier curves bij elkaar in 1 grafiek plotten?

Voor een paper heb ik 3 kaplan meier curves die op de x-as allemaal even ver lopen (allemaal een tijdsduur van 27 maanden). Ik wil deze 3 curves in 1 plot zetten. Kan dat met SPSS en zo ja, hoe?

We zijn er zelf uitgekomen en wel op de volgende manier: In ons geval wilden wij 3 verschillende curves (survival, grote amputatie en kleine amputatie) in 1 KM grafiek zetten. We hebben de data opnieuw gerangschikt. 1 variable met time to event, 1 variable met event status (alledrie de events hadden status 1), en een nieuwe onderscheidende variable die duidelijk maakte dat de eerste time to event en event zelf sloeg op survival, tweede time to event en event zelf op grote ampuatie en drede time to event en event zelf op kleine amputatie. Invoeren gewoon time to event en event status en de nieuwe onderscheidende variable wordt de 'factor'. Mocht iemand anders ooit nog met hetzelfde probleem zitten, bij deze de oplossing.

Hoe kan ik mijn logrank analyse doen over alleen het eerste jaar van de follow up?

Op een polikliniek voor HIV-geinfecteerde kinderen in Malawi hebben we in de periode van 2004-2008 een cohort kinderen gevolgd vanaf het moment dat zij startten met antiretrivirale therapie. We zijn geinteresseerd in hun uitkomsten en risicofactoren voor mortaliteit en uitval uit behandeling (LTFU). 1) Ik probeer uitkomsten na 12 maanden op een rij te zetten dmv analyze_survival_lifetables. Dit lijk te lukken. 2) Ik wil naar voorspellende factoren kijken voor overlijden en dood op 12 maanden. Voor univariate analyses zou ik hiervoor willen gebuiken: analyze_survival -> kaplanmeier met compare factor log-rank test. Hiermee doe ik een analyse voor uitkomst na totale duur folllow up. Het lukt mij niet hierbij specifiek naar uitkomsten tot 12 maanden te kijken.

Als er een langere follow up is, en je wilt in een analyse alleen focussen op de periode tot 12 maanden, dan is het nodig om de informatie van de personen op 12 maanden te censureren. Als het goed is heb je nu in de dataset de volgende 2 variabelen:

-tijd tot laatste follow up of tijd tot optreden event (mortaliteit of uitval behandeling): FUtime

-event indicator (wel of niet het event opgetreden): FUstatus (1=event, 0= geen event)

Om alleen op de eerste 12 maanden te focussen, is het nodig om de volgende aanpassingen te maken:

indien FUtime>12 en FUstatus =0 -> FUtime=12 FUstatus=0

indien FUtime>12 en FUstatus =1 -> FUtime=12 FUstatus=0

indien FUtime=<12 en FUstatus =1 -> variabelen blijven hetzelfde

indien FUtime=<12 en FUstatus =0 -> variabelen blijven hetzelfde

Met de nieuwe FUtime (die nooit groter is dan 12) en de nieuwe status (die nu iets minder events bevat) kan een Kaplan-Meier analyse gedaan worden voor alleen het eerste jaar.

Hoe ga ik om met patienten die vanwege overlijden niet meer het event van interesse kunnen krijgen?

Ik heb een RCT gedaan met een follow-up van 1 jaar met als belangrijkste uitkomstmaat: het aantal dagen vanaf randomisatie tot terugkeer naar werk. Nu wil ik onderzoeken of de tijd tot terugkeer naar werk significant verschilt tussen beide groepen (controlegroep en interventiegroep). Hiervoor heb ik KM-survival analyse gedaan met Log rank toets. Nu is het zo dat een aantal patiënten gedurende de follow-up periode overleden zijn en een aantal patiënten palliatief worden behandeld. Deze patiënten kunnen de positieve uitkomstmaat (weer aan het werk) nooit meer bereiken. Hoe ga je om met deze patiënten in de survival analyse? Moet ik deze patiënten helemaal uit de analyse halen? Moet ik deze als censored beschouwen? Moet ik deze patiënten zien als negatieve uitkomst (=niet aan het werk). Kortom wat is het effect van censored data op een survival analyse en logrank toets en welke afweging is het meest juist?

De voor- en nadelen van de verschillende alternatieven:

1. de overleden/palliatieve patienten in de KM-analyse laten en censureren op het moment van dit event

Voordeel: de informatie van deze patienten voordat ze overleden / voordat ze palliatief waren doet mee in de analyse / geen verlies van data Nadeel: de KM-analyse veronderstelt dat de prognose van deze patienten, nadat ze gecensureerd worden, gelijk is aan patienten die nog wel in de analyse zitten (aanname van non-informative censoring). De schattingen (en vergelijkingen) op latere follow up momenten worden daardoor (onevenredig) veel bepaald door relatief goede patienten die wellicht niet representatief zijn voor de totale groep.

2. de overleden/palliatieve patienten in de KM-analyse uit de analyse laten

Voordeel: eenvoudig uit te leggen aanpak, je beperkt de studiegroep alsware het een exclusiecriterium is Nadeel: alleen als op t=0 al bekend is dat iemand tot de 'exclusiegroep' behoort is dit een logische aanpak, voor deze personen zou je kunnen zeggen dat de interventie toch geen invloed meer kan hebben op terugkeer naar werk omdat dat (ongeacht welke behandeling ook) geen resultaat hierop kan hebben. Voor personen die in het begin nog alle kanten op kunnen en pas later achteruit gaan zou het zonde zijn/kan het vertekening geven die info niet te gebruiken.

3. de overleden/palliatieve patienten in de KM-analyse als negatieve uitkomst zien / Competing risk analyse

Voordeel: Je zou kunnen ageren dat deze patienten een ultiem 'slechte' uitkomst hebben en ook als zodanig meegerekend moeten worden. Je zegt dan eigenlijk dat ze het event 'terug naar werk' nooit meer zullen bereiken. Hun status staat dan op 0, en de follow up tijd zou je (gekunsteld) op oneindig lang kunnen zetten (of iig langer dan de langste follow up). Je kunt dit ook formeel in een competing risk analyse doen. Nadeel: de analyse is minder standaard / competing risk analyse wordt niet standaard in SPSS aangeboden.

Hoe bereken ik de 'events per person year'?

'Ik wil graag mijn uitkomsten (opnames in verband met verklevingen in de buik en het vóórkomen van littekenbreuken) weergeven als 'event per person-years at risk', om te corrigeren voor de wisselende follow up duur. Hoe kan ik dit het beste doen? Ik zou een Kaplan-Meier kunnen maken van time-to-event, maar ik wil ook graag deze maat weergeven. Moet ik dan een variable aanmaken waarin ik event ja/ nee deel door de tijd 'at risk'? Of is er een functie in SPSS waarmee ik dit makkelijker kan berekenen?

De uitkomst 'events per person-years at risk' is simpel te berekenen door alle events bij elkaar op te tellen (som_events), alle tijd at risk bij elkaar op te tellen (som_tijd) en die twee op elkaar te delen: som_events/som_tijd.

Omdat SPSS gewend is om te rekenen 'per rij' en niet zozeer per kolom, zijn deze optellingen makkelijker te doen als je de data even in een ander pakket (excel of R) zet.

Bij de maat 'events per person-years at risk' wil je waarschijnlijk ook een betrouwbaarheidsinterval berekenen? Die zou je kunnen berekenen mbv de calculator op www.openepi.com. Kies in de linkerbalk voor person time -> 1 rate. Geef het aantal events en de totale tijd in. Je krijgt dan de rate en een betrouwbaarheidsinterval terug. Let goed op, de rate wordt per 10, 100, 1000 person years gegeven (afhankelijk van het aantal informatieve decimalen). Het is dus goed om de getallen naast je eigen berekening (zoals boven) te leggen.

Waar vind ik de Kaplan-Meier-analyse in SPSS?

Je vindt de Kaplan-Meier-analyse in SPSS 28 onder Analyze -> Survival -> Kaplan-Meier.

Cox regressie

Wanneer gebruik ik een Cox proportional hazards model?

Wanneer er bij een survival analyse slechts één voorspellende factor onderzocht wordt die slechts een paar mogelijke waardes kan aannemen, kan een Kaplan Meier analyse volstaan. Zijn er echter meerdere factoren die de overleving kunnen beinvloeden, dan is een uitgebreider model nodig. Een veel gebruikt model is het Cox proportional hazards model. Hiermee kun je net als met andere multivariabele regressie modellen meerdere voorspellers tegelijk bestuderen. Voordeel van het Cox model boven een Kaplan Meier analyse is bovendien dat er naast de vraag "is er een effect" ook antwoord gegeven wordt op de vraag "hoe groot is het effect" van een bepaalde voorspeller. Het Cox model geeft de grootte van een effect weer met behulp van hazard ratio's (HR), waarbij een HR van 1 aangeeft dat er geen verschil is. De hazard is de instantane kans om een event te krijgen, en een hazard ratio is vergelijkbaar met een relatief risico. Bijvoorbeeld: "During the 20-year follow-up, 83 women developed lung cancer. There was a moderate association between hormone use and lung cancer (hazard ratio, 1.13; 95% CI, 1.08-1.73)."

Het Cox model maakt een belangrijke aanname: de proportional hazards assumptie. Deze zegt dat de HR's constant zijn in de tijd. Dit betekent dat het relatieve verschil in hazard veroorzaakt door verschillende waardes van een voorspeller gelijk blijft in de tijd. Er wordt in het bovenstaande voorbeeld dus vanuit gegaan dat het wel of niet gebruiken van hormonen een zelfde invloed heeft na 1 jaar als na 20 jaar follow up: op ieder tijdstip wordt de kans op het optreden van longkanker bij hormoongebruik 1.13 keer hoger geschat dan dezelfde kans als er geen hormonen werden gebruikt. Bij een Cox regressie, zijn er verder geen specifieke aannames over de verdeling van de hazard functie of voorspellende variabelen. Ze hoeven dus niet normaal verdeeld te zijn .

Hoe test ik de proportional hazards aanname binnen een Cox model?

Er zijn meerdere manieren om de PH aanname van een Cox model te testen:

• Grafisch

Als een variabele in het model voldoet aan de proportional hazards (PH) aanname, dan zullen de geschatte overlevingscurves, op de zogenaamde log-schaal, voor de waardes van deze variabele parallel lopen. Je kunt hier naar kijken door in het cox model de variabele voor welke je de PH aanname wilt onderzoeken niet meer als covariaat, maar als stratum mee te geven en te vragen om log minus log survival plots. De resulterende lijnen moeten op het oog parallel lopen (zeker elkaar niet kruisen) om te voldoen aan de PH aanname voor deze variabele. Voor iedere variabele in het model is een dergelijke plot nodig. In deze voorbeeldsyntax kun je zien hoe het maken van een log minus log plot er in SPSS syntax uit komt te zien.

• Testen van (Schoenfeld) residuen

Er kan een formele test gedaan worden mbv Schoenfeld residuen. Er is geen standaard functie hiervoor in SPSS. In Stata (functie estat phtest) of R (functie cox.zph) kan het wel direct. Download hier een script voor het testen van de PH aanname in R.

• Toevoegen van een tijdsafhankelijke variabele

Anticiperend op het mogelijk niet voldoen aan de PH aanname kun je ook een time dependent Cox model draaien met daarin een tijdsafhankelijke variabele opgenomen (bijvoorbeeld GENDER*T). Als deze tijdsafhankelijke variabele significant bijdraagt aan het model, is er blijkbaar niet aan de PH aanname voldaan.

Mijn variabelen lijken niet te voldoen aan de proportional hazards-aanname, mag ik ze nu niet meenemen in de Cox-regressie?

Ik ben aan de slag gegaan met het testen van de proportional hazards-aanname. In de log-log-plots kruisen de lijnen van bijna alle covariaten. Betekent dit dat ik al deze covariaten niet mee mag nemen in de analyse?

Het is niet zo dat als de lijnen kruisen dat je de variabele dan sowieso niet mee mag nemen. Je moet per variabele bekijken wat er aan de hand is. Mogelijk opties:

1. Een variabele is niet gerelateerd aan de uitkomstmaat. De log-minus-log survival curves lopen dan naar verwachting over elkaar heen. Door random fluctuatie kunnen de curves ook 'om elkaar heen wiebelen' en elkaar dus meerdere keren kruisen. Conclusie is hier dat er geen effect is (ook geen proportioneel effect).

2. Een variabele is gerelateerd aan de uitkomstmaat, maar de curves kruisen op een gegeven ogenblik. Hier is het lastiger. Als de curves kruisen, dan is het effect van een variabele niet proportioneel. Dat wil zeggen, op het ene tijdsinterval is de ene groep beter af en op het andere tijdsinterval de andere groep. Als de Cox-analyse toch een proportioneel effect van een variabele aanwijst, dan betekent dit zoveel als dat er gemiddeld over de hele tijdslijn toch 1 van de 2 groepen beter uit de bus komt. 'Uitgemiddeld' doet deze groep het dan beter. Dit kan in bepaalde situaties een valide conclusie zijn, maar wanneer alleen de Cox-resultaten (HRs) getoond worden vertelt dit niet het hele verhaal. Bijvoorbeeld als een groep in het begin een groot voordeel laat zien, maar richting het einde van de follow-up het slechter doet dan zijn tegenhanger, dan is het te verwachten dat in een nog langere follow-up-periode het effect om kan gaan slaan. Deze informatie is goed om erbij te vermelden; als een lezer alleen de HR ziet, denkt ie dat de groep het overal in de tijd beter zal doen. Wanneer de variabele duidelijk (en verklaarbaar) een differentieel effect in de tijd laat zijn, dan kan het verstandig te zijn om 2 HRs te schatten, 1 voor de vroege periode en 1 voor de late. Dit kan m.b.v. het aanmaken van een tijdsafhankelijk covariaat.

Hoe kan ik de output lezen van de analyse van een categorische variabele in een Cox model?

Voor een cohort HIV geinfecteerde kinderen, gevolgd vanaf het moment dat ze met antiretrovirale medicatie starten, ben ik op zoek naar voorspellers van overlijden. Ik ben nu bezig met de Cox regressie (vooralsnog bezig met univariate analyse). Hieronder een voorbeeld van output van een categoriale variable (leeftijdscatgorieën) met meer dan 2 categorieën (referentie = kinderen>12 van wie we verwachten dat ze de beste prognose hebben) (1) (2) (3) betreffen 0-1jaar, 1-5jaar en 5-12jaar oude kinderen. Wat betekent de p-waarde die genoemd staat op de regel van de referentiewaarde?

Variables in the Equation
		B	S.E.	Wald	df	Sig.	Exp(B)	95.0% C.I.for Exp(B)
								Lower	Upper
Leeftijd	leeftijd			23.539	3	.000
	leeftijd(1)	1.094	.261	17.613	1	.000	2.985	1.791	4.795
	leeftijd(2)	.773	.211	13.379	1	.000	2.166	1.431	3.276
	leeftijd(3)	.351	.236	2.206	1	.137	1.420	.894	2.256

De eerste regel laat het resultaat zien van de overall test voor de variabele. Dit is de test die als nulhypothese heeft dat alle 4 de categorieën dezelfde kans (eigenlijk hazard) hebben op overlijden. De toets kijkt of er ergens een verschil is in categorieën, je ziet ook dat de test 3 'degrees of freedom' heeft. Je kunt het vergelijken met een ANOVA test bij het vergelijken van gemiddelde waarden in 4 groepen. De test hoort dus niet specifiek bij de referentiecategorie, dat vindt SPSS gewoon een handige locatie om 'm neer te zetten. In principe ga je als volgt te werk: je kijkt eerst of de overall test significant is. Als dat zo is, dan ben je geoorloofd om door te gaan en te kijken naar de verschillen tussen de categorieën onderling.

Hoe rapporteer ik de resultaten van een categorische variabele in een Cox model?

Ik heb een coxregressie uitgevoerd met een dichotome uitkomstmaat (wel/geen artritis ontwikkeling) en de variabele is een categorische maat (0-4). Normaal gesproken ben ik gewend om de uitslag met een continue variabele als volgt te noteren: HR x.x (95% CI x.x to x.x); p=0.xxx. Nu las ik echter dat je vanwege de categorische maat moet kijken naar de Wald uitslag; hier komt wel een p-waarde uit, maar geen CI e.d. Hoe moet je deze uitslag dan noteren in een manuscript?

Meestal worden de significantietoetsen van HR's in Cox regressie met een Wald statistic gegenereerd. Dat geldt evengoed voor continue als voor categorische voorspellers. Wat hier wellicht het issue is, is dat er bij het evalueren van een categorical voorspeller eerst bekeken moet worden of de 'overal' toets significant is, dwz of de toets die kijkt of het opsplitsen in de 5 categorien (met 4 df) bijdraagt aan het model sign is. Als die uitslag positief is, kun je door gaan met inzoomen op de verschillende contrasten (HR's) apart en deze op gebruikelijke wijze met 95% CI rapporteren. De p-waardes van de losse HR's zijn doorgaans wat minder relevant omdat deze zeer afhankelijk zijn van de gekozen referentiecategorie.

Mijn Cox regressie convergeert niet, heb ik wel de juiste analyse gebruikt?

Ik wil graag een Cox regression doen op een patiënten cohort om te bepalen of sommige variabelen een independent predictive value hebben. Ik heb een cohort van 86 patiënten. Voor de cohort wil ik 8 variabelen vergelijken, 7 binaire en een continue (leeftijd). In de cohort waren er 6 events. Als ik het Cox regression model dan probeer te runnen krijg ik de volgende warning: ‘Since coefficients did not converge no further models will be fitted’. De tabel geeft dan voor alle variabelen een significantie van 1 en een hazard ratio van 1. Wat betekent dit? En is een Cox regression wel de geschikte test voor mijn data dan?

Ik denk dat je last hebt van een ernstige vorm van te weinig power. Je hebt maar 6 events in de data. Voor een Cox regressie zijn er volgens de vuistregels 10 a 20 events per parameter die je wilt testen nodig. Met 6 events kun je helaas niet zinvol naar effecten van variabelen kijken mbv een Cox model. Uiteraard kun je wel de ruwe data beschrijven.

Waar vind ik het Cox proportional hazards model in SPSS?

Je vindt de Kaplan Meier analyse in SPSS 16 onder Analyze -> Survival -> Cox Regression of Cox w./ Time Dep Cov.

Referenties

• Altman DG, Bland JM. Time to event (survival) data. BMJ 1998;317:468-9.
• Bewick V, Cheek L, Ball J. Statistics review 12: Survival analysis. Critical Care 2004, 8:389-394.
• Bland JM, Altman DG. Survival probabilities. The Kaplan-Meier method. BMJ 1998;317:1572.
• Bland JM and Altman DG. The log rank test. BMJ 2004;Vol 328 may 2004.
• An Introduction to Survival Analysis - by Mark Stevenson from EpiCentre, IVABS, Massey University. 31 pagina's cursusmateriaal (Engels) met voorbeeld syntax in R.
• Singh R and Mukhopadhyay K. Survival analysis in clinical trials: Basics and must know areas. Perspect Clin Res. 2011 Oct-Dec; 2(4): 145–148.
• Hosmer, D. W., and S. Lemeshow. 1999. Applied Survival Analysis. New York: John Wiley and Sons.
• Kleinbaum, D. G. 1996. Survival Analysis: A Self-Learning Text. New York: Springer-Verlag.
• Norusis, M. 2004. SPSS 13.0 Advanced Statistical Procedures Companion. Upper Saddle-River, N.J.: Prentice Hall, Inc..
• R tutorial survival analyse - youtube filmpje met uitleg over uitvoeren van een survival analyse in R.
• R tutorial Cox analyse - youtube filmpje met uitleg over uitvoeren van een Cox analyse in R.

Klaar met lezen? Je kunt naar het OVERZICHT van alle statistische onderwerpen op deze wiki gaan of naar de pagina KEUZE TOETS voor hulp bij het uitzoeken van een geschikte toets of analyse. Wil je meer leren over biostatistiek? Volg dan de AMC e-learning Practical Biostatistics. Vind je op deze pagina's iets dat niet klopt? Werkt een link niet? Of wil je bijdragen aan de wiki? Neem dan contact met ons op.

De wiki biostatistiek is een initiatief van de voormalige helpdesk statistiek van Amsterdam UMC, locatie AMC. Medewerkers van Amsterdam UMC kunnen via intranet ondersteuning aanvragen.