Veröffentlicht am :
Statistik SASHELP

Überlebensanalyse mit Kurven und Risikotabellen

Dieser Code ist auch verfügbar auf: English Español Français
Wartet auf Validierung
Das Skript beginnt mit der Konfiguration der ODS-Umgebung (Output Delivery System) für die grafische Ausgabe. Anschließend wird die LIFETEST-Prozedur für den Datensatz 'sashelp.BMT' ausgeführt, um eine nach der Variablen 'Group' stratifizierte Überlebensanalyse durchzuführen. Die Daten des Überlebensdiagramms werden in einem Datensatz namens 'SurvivalPlotData' erfasst. Danach wird die SGPLOT-Prozedur dreimal verwendet, um verschiedene Versionen eines Überlebensdiagramms zu generieren, wobei die Position und der Stil der Tabelle der Probanden im Risiko (außerhalb oder innerhalb des Diagramms) und der allgemeine ODS-Stil variiert werden.
Datenanalyse

Type : SASHELP


Das Skript verwendet ausschließlich den Datensatz BMT aus der Standard-SASHELP-Bibliothek, der Daten zur Knochenmarktransplantation enthält.

1 Codeblock
PROC LIFETEST Data
Erklärung :
Dieser Block initialisiert die grafische Ausgabeumgebung und führt die LIFETEST-Prozedur aus. Die Überlebensanalyse wird auf den Daten sashelp.BMT unter Verwendung der Variablen T für die Zeit und Status für das Ereignis durchgeführt. Die Analyse ist nach 'Group' stratifiziert. Die grafischen Ergebnisse, einschließlich der Punkte der Überlebenskurve, werden dank der 'ods output'-Anweisung in der Tabelle 'SurvivalPlotData' gespeichert.
Kopiert!
1%let gpath='.';
2%let dpi=300;
3ods html close;
4ods listing style=htmlblue image_dpi=&dpi gpath=&gpath;
5 
6ods graphics on;
7ods OUTPUT Survivalplot=SurvivalPlotData;
8PROC LIFETEST DATA=sashelp.BMT plots=survival(atrisk=0 to 2500 BY 500);
9 time T * STATUS(0);
10 strata Group / test=logrank adjust=sidak;
11RUN;
2 Codeblock
PROC SGPLOT
Erklärung :
Dieser Block verwendet PROC SGPLOT, um eine Überlebenskurve (STEP-Anweisung) mit Markierungen für zensierte Daten (SCATTER-Anweisung) zu generieren. Das Hauptmerkmal dieses Diagramms ist das Hinzufügen einer Tabelle der Probanden im Risiko (xaxistable), die außerhalb, unterhalb der X-Achse positioniert ist.
Kopiert!
1ods graphics / reset width=5in height=3in imagename='4_5_1_Survival Plot_SG_V94';
2title 'Product-Limit Survival Estimates';
3title2 h=0.8 'With Number of AML Subjects at Risk';
4PROC SGPLOT DATA=SurvivalPlotData;
5 step x=time y=survival / group=stratum lineattrs=(pattern=solid) name='s';
6 scatter x=time y=censored / markerattrs=(symbol=plus) name='c';
7 scatter x=time y=censored / markerattrs=(symbol=plus) GROUP=stratum;
8 xaxistable atrisk / x=tatrisk location=outside class=stratum colorgroup=stratum;
9 keylegend 'c' / location=inside position=topright;
10 keylegend 's';
11RUN;
3 Codeblock
PROC SGPLOT
Erklärung :
Ähnlich wie der vorherige Block generiert dieser Code dieselbe Überlebenskurve. Der Hauptunterschied besteht darin, dass die Tabelle der Probanden im Risiko (xaxistable) mit 'location=inside' konfiguriert ist, wodurch sie unten im Diagramm, über dem Plotbereich, platziert wird.
Kopiert!
1ods listing style=htmlBlue;
2ods graphics / reset width=5in height=3in imagename='4_5_2_Survival Plot_Inner_SG_V94';
3title 'Product-Limit Survival Estimates';
4title2 h=0.8 'With Number of AML Subjects at Risk';
5PROC SGPLOT DATA=SurvivalPlotData;
6 step x=time y=survival / group=stratum lineattrs=(pattern=solid) name='s';
7 scatter x=time y=censored / markerattrs=(symbol=plus) name='c';
8 scatter x=time y=censored / markerattrs=(symbol=plus) GROUP=stratum;
9 xaxistable atrisk / x=tatrisk location=inside class=stratum colorgroup=stratum
10 separator valueattrs=(size=7 weight=bold) labelattrs=(size=8);
11 keylegend 'c' / location=inside position=topright;
12 keylegend 's';
13RUN;
4 Codeblock
PROC SGPLOT
Erklärung :
Dieser letzte Block ändert den ODS-Stil auf 'journal' und generiert eine dritte Version des Diagramms. Die Tabelle der Probanden im Risiko befindet sich immer noch innen. Ästhetische Optionen werden hinzugefügt, wie z.B. Beschriftungen direkt auf den Kurven (curvelabel) und ein anderer Markierungsstil für zensierte Daten (vollständiger Kreis).
Kopiert!
1ods listing style=journal;
2ods graphics / reset width=5in height=3in imagename='4_5_3_Survival Plot_Inner_Journal_SG_V94';
3title 'Product-Limit Survival Estimates';
4title2 h=0.8 'With Number of AML Subjects at Risk';
5PROC SGPLOT DATA=SurvivalPlotData;
6 step x=time y=survival / group=stratum lineattrs=(pattern=solid) name='s'
7 curvelabel curvelabelattrs=(size=6) splitchar='-';
8 scatter x=time y=censored / name='c' markerattrs=(symbol=circlefilled size=4);
9 xaxistable atrisk / x=tatrisk location=inside class=stratum colorgroup=stratum
10 separator valueattrs=(size=7 weight=bold) labelattrs=(size=8);
11 keylegend 'c' / location=inside position=topright;
12RUN;
Dieses Material wird von We Are Cas "wie besehen" zur Verfügung gestellt. Es gibt keine ausdrücklichen oder stillschweigenden Garantien hinsichtlich der Marktgängigkeit oder Eignung für einen bestimmten Zweck in Bezug auf die hierin enthaltenen Materialien oder Codes. We Are Cas ist nicht verantwortlich für Fehler in diesem Material, wie es jetzt existiert oder existieren wird, noch bietet We Are Cas technischen Support dafür an.