01.03.2013 Views

Applied Statistics Using SPSS, STATISTICA, MATLAB and R

Applied Statistics Using SPSS, STATISTICA, MATLAB and R

Applied Statistics Using SPSS, STATISTICA, MATLAB and R

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

362 9 Survival Analysis<br />

A: The Heart Valve Survival datasheet contains the computed final date<br />

for the study (variable DATE_STOP). This is the date of the first occurring event,<br />

if it did occur, or otherwise, the last date the patient was known to be alive <strong>and</strong><br />

well. The survivor function estimate shown in Figure 9.4 is obtained by using<br />

<strong>STATISTICA</strong> with DATE_OP <strong>and</strong> DATE_STOP as initial <strong>and</strong> final dates, <strong>and</strong><br />

variable EVENT as censored data indicator. From this figure, one can estimate that<br />

about 85% of patients survive five years (1825 days) without any event occurring.<br />

1.2<br />

1.1<br />

1.0<br />

0.9<br />

0.8<br />

0.7<br />

0.6<br />

0.5<br />

0.4<br />

S(t)<br />

Complete Censored<br />

t (days)<br />

0 1000 2000 3000 4000 5000 6000 7000 8000<br />

Figure 9.4. Kaplan-Meier estimate of the survivor function for the event-free<br />

survival of patients with heart valve implant, obtained with <strong>STATISTICA</strong>.<br />

9.2.3 <strong>Statistics</strong> for Non-Parametric Analysis<br />

The following statistics are often needed when analysing survival data:<br />

1. Confidence intervals for S(t).<br />

For the Kaplan-Meier estimate, the confidence interval is computed assuming that<br />

the estimate Sˆ ( t)<br />

is normally distributed (say for a number of intervals above 30),<br />

with mean S(t) <strong>and</strong> st<strong>and</strong>ard error given by the Greenwood’s formula:<br />

[ Sˆ<br />

( t)<br />

]<br />

2<br />

ˆ ⎪⎧<br />

k d j ⎪⎫<br />

s ≈ S(<br />

t)<br />

⎨∑<br />

j=<br />

1 ⎬ , for tk ≤ t < tk+1. 9.13<br />

⎪⎩<br />

n j ( n j − d j ) ⎪⎭

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

Saved successfully!

Ooh no, something went wrong!