Overlevelsesanalyse
Overlevelsesanalyse
Overlevelsesanalyse
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Ventetidsanalyse Styrkeberegning 17<br />
Styrkeberegning for ventetidsdata<br />
Summen af varianserne, V , brugt ved beregningen af log rank<br />
teststørrelsen, bruges ogs˚a ved styrkeberegninger:<br />
Hvis man ønsker at kunne p˚avise en rate ratio p˚a RR for aktiv versus<br />
kontrol med et tosidet signifikansniveau p˚a α og en styrke p˚a 1−β, s˚a skal<br />
man inkludere tilstrækkeligt mange personer til, at V ≥<br />
( Uα/2 +U β<br />
ln(RR)<br />
(Ux=(1−x)-fraktilen i en normal fordeling, eksempelvis U 2.5%=1.96,<br />
U 5%=1.645, U 10%=1.282, U 20%=0.842)<br />
Hvis randomiseringsratioen (aktiv:kontrol) er a:1, og hændelsen er sjælden,<br />
er V ≈ a<br />
(a+1) 2 · M, hvor M=total antal hændelser, dvs. der skal inkluderes<br />
tilstrækkeligt mange personer til, at<br />
M ≥<br />
(a + 1)2<br />
a<br />
·<br />
( ) 2<br />
Uα/2 + Uβ<br />
ln(RR)<br />
Ventetidsanalyse Log rank 19<br />
Dele af output fra PROC PHREG<br />
The PHREG Procedure<br />
Model Information<br />
Data Set WORK.SKL<br />
Dependent Variable Day<br />
Censoring Variable Bld<br />
Censoring Value(s) 0<br />
Ties Handling DISCRETE<br />
Summary of the Number of Event and Censored Values<br />
Percent<br />
Total Event Censored Censored<br />
187 91 96 51.34<br />
Model Fit Statistics<br />
Without With<br />
Criterion Covariates Covariates<br />
-2 LOG L 738.406 737.488<br />
: : :<br />
Testing Global Null Hypothesis: BETA=0<br />
Test Chi-Square DF Pr > ChiSq<br />
Likelihood Ratio 0.9175 1 0.3381<br />
Score 0.9174 1 0.3382<br />
Wald 0.9124 1 0.3395<br />
) 2<br />
.<br />
Ventetidsanalyse Log rank 18<br />
Log rank test ved hjælp af PROC PHREG<br />
Log rank teststørrelsen kan beregnes som score-teststørrelsen i<br />
PROC PHREG med brug af option TIES=DISCRETE:<br />
PROC PHREG DATA=skl;<br />
MODEL day*bld(0) = sklero / TIES=DISCRETE;<br />
RUN;<br />
Hvis der er mere end 2 grupper, er det nødvendigt at bruge et CLASS statement,<br />
eksempelvis<br />
PROC PHREG DATA=skl;<br />
CLASS sklero;<br />
MODEL day*bld(0) = sklero / TIES=DISCRETE;<br />
RUN;<br />
Ventetidsanalyse Log rank 20<br />
Log rank teststørrelse=0.9174, p=0.3382, findes ud for “Score”<br />
under “Testing Global Null Hypothesis: BETA=0”<br />
Output fra PROC PHREG fortsat<br />
Analysis of Maximum Likelihood Estimates<br />
Parameter Standard Hazard<br />
Variable DF Estimate Error Chi-Square Pr>ChiSq Ratio<br />
Sklero 1 -0.20261 0.21212 0.9124 0.3395 0.817<br />
Log rank testet er d˚arligt til at detektere forskelle, der<br />
ændres markant med tiden, som en d˚arligere<br />
korttidsprognose samtidigt med en bedre langtidsprognose!