28.07.2013 Views

Overlevelsesanalyse

Overlevelsesanalyse

Overlevelsesanalyse

SHOW MORE
SHOW LESS

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!

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!