Plot KI Kaplan-Meier

[swimfan]

Hallo,
weiß jemand hier, wie man von einer Überlebenskurve nach Kaplan-Meier die 95%KI als Grafik bekommen kann? Bisher habe ich diese Syntax gefunden:

km tZ5JLoPro /status=Z5JLoPro(2)
/print = none
/save survival(surv_km) se(surv_se).

compute var2_6 = surv_se**2/(surv_km**2*(ln(surv_km))**2).
compute cl =ln(-ln(surv_km)) - 1.96*sqrt(var2_6).
compute cu =ln(-ln(surv_km)) + 1.96*sqrt(var2_6).
compute l = exp(-exp(cu)).
compute u = exp(-exp(cl)).
exe.

graph /scatterplot(overlay)=tZ5JLoPro tZ5JLoPro tZ5JLoPro with surv_km l u (pair).

tZ5JLoPro ist meine Zeitvariable, Z5JLoPro(2) ist mein Zensor.

Nun ist das Problem, dass die Grafik mir nur die Zensiert-Kreuzchen anzeigt und obwohl in der Legende die Überlebenskurven an sich auftauchen, ich es nicht schaffe, dass man die Linien auch tatsächlich im Diagramm sehen kann.

Frage 2: Weiß jemand, wie ich die Werte für die 95%KI zu beliebigem Zeitpunkt in einer Tabelle bekommen kann?

Danke für eure Hilfe, Brit

[drfg2008]

http://spssx-discussion.1045642.n5.nabble.com/

[swimfan]

Danke! Nun hab ich die Linien.

Neue Frage: Weiß zufällig jemand, wie ich die Farbe der Interpolationslinien separat verändern kann? Nun hab ich nämlich in meinem Graph eine Variable mit vier Ausprägungen, d.h. inclusive KI´s 12 schwarze Linien... sehr unübersichtlich

Viele Grüße und Dank, Brit