Das Skript beginnt mit der Erstellung von drei separaten Datensätzen ('control', 'rad', 'radbpa') unter Verwendung von DATA-Blöcken und Datalines-Anweisungen, um die Überlebenszeit- und Zensurdaten zu integrieren. Jeder Datensatz repräsentiert eine andere Behandlungsgruppe im Rahmen einer Studie. Diese drei Datensätze werden anschließend zu einem einzigen Datensatz namens 'prob7_7' zusammengeführt. Zwei Indikatorvariablen ('grp1', 'grp2') werden erstellt, um die Identifizierung der Gruppen 'rad' und 'radbpa' in späteren Analysen zu erleichtern. Nach einer kurzen Anzeige der kombinierten Daten über PROC PRINT zur Überprüfung wird der Kern der Analyse von PROC LIFETEST durchgeführt. Diese Prozedur ist so konfiguriert, dass sie die Überlebensfunktionen (Kaplan-Meier standardmäßig) und die Risikofunktionen zwischen den drei Gruppen ('trt') schätzt und vergleicht, mit einer Bonferroni-Korrektur für multiple Vergleiche. Sie generiert auch Überlebens- und Risikografiken zur Visualisierung der Ergebnisse. Symbolisierungsoptionen werden verwendet, um das Erscheinungsbild der Kurven in den Grafiken zu differenzieren.
Datenanalyse
Type : INTERNE_ERSTELLUNG
Die Daten (`time`, `censor`, `trt`) werden vollständig intern innerhalb des SAS-Skripts über DATA-Blöcke und `datalines`-Anweisungen erstellt. Es gibt keine Abhängigkeit von externen Dateien oder bereits vorhandenen SAS-Bibliotheken für die Erfassung der Rohdaten, außer den vom Skript selbst generierten Daten.
1 Codeblock
DATA STEP Data
Erklärung : Dieser DATA-Block erstellt einen Datensatz namens `control`. Er verwendet die `datalines`-Deklaration, um die Werte der Variablen `time` und `censor` zu lesen. Der Variablen `trt` wird der Wert 'control' zugewiesen, um diese Gruppe als die unbehandelte Gruppe zu kennzeichnen. Die im Originalcode in der Zeile `input` gefundene Syntax ` @code_sas_json/hsdua2304@gmail.com_SAS_Assignment_1.json` wurde entfernt, da sie ein Syntaxfehler ist und nicht dem Lesemechanismus der `datalines` entspricht. Der Wert 'controp' wurde zu 'control' korrigiert, um Konsistenz und Ausführbarkeit des Codes zu gewährleisten.
Kopiert!
data control; /*create data set for untreated group*/
input time censor;
trt='control'; /*Correction de 'controp' en 'control' pour cohérence avec le commentaire*/
datalines;
20 1 21 1 23 1 24 1 24 1 26 1 26 1 27 1 28 1 30
;
1
DATA control; /*create data set for untreated group*/
2
INPUT time censor;
3
trt='control'; /*Correction de 'controp' en 'control' pour cohérence avec le commentaire*/
4
DATALINES;
5
20121123124124126126127128130
6
;
7
2 Codeblock
DATA STEP Data
Erklärung : Dieser DATA-Block erstellt einen Datensatz namens `rad` für die bestrahlte Gruppe. Ähnlich wie der vorherige Block liest er die Daten für `time` und `censor` über `datalines` und weist der Variablen `trt` den Wert 'rad' zu. Die im Originalcode in der Zeile `input` gefundene Syntax ` @code_sas_json/hsdua2304@gmail.com_SAS_Assignment_1.json` wurde entfernt, da sie ein Syntaxfehler ist und nicht dem Lesemechanismus der `datalines` entspricht.
Kopiert!
data rad; /*create data set for the radiated group*/
input time censor;
trt='rad';
datalines;
26 1 28 1 29 1 29 1 30 1 30 1 31 1 31 1 32 1 35 0
;
1
DATA rad; /*create data set for the radiated group*/
2
INPUT time censor;
3
trt='rad';
4
DATALINES;
5
26128129129130130131131132135 0
6
;
7
3 Codeblock
DATA STEP Data
Erklärung : Dieser DATA-Block erstellt einen Datensatz namens `radbpa` für die mit Strahlung und BPA behandelte Gruppe. Er liest die Daten für `time` und `censor` über `datalines` und weist der Variablen `trt` den Wert 'radbpa' zu. Die im Originalcode in der Zeile `input` gefundene Syntax ` @code_sas_json/hsdua2304@gmail.com_SAS_Assignment_1.json` wurde entfernt, da sie ein Syntaxfehler ist und nicht dem Lesemechanismus der `datalines` entspricht.
Kopiert!
data radbpa; /*Create data set for the raduated + BPA group*/
input time censor;
trt='radbpa';
datalines;
31 1 32 1 34 1 35 1 36 1 38 1 38 1 39 1 42 0 42 0
;
1
DATA radbpa; /*Create data set for the raduated + BPA group*/
2
INPUT time censor;
3
trt='radbpa';
4
DATALINES;
5
31132134135136138138139142 0 42 0
6
;
7
4 Codeblock
DATA STEP Data
Erklärung : Dieser DATA-Block kombiniert die drei zuvor erstellten Datensätze (`control`, `rad`, `radbpa`) zu einem neuen, einzigen Datensatz namens `prob7_7`. Er erstellt außerdem zwei binäre Variablen `grp1` (ist 1, wenn `trt` 'rad' ist, sonst 0) und `grp2` (ist 1, wenn `trt` 'radbpa' ist, sonst 0), um Vergleiche zu erleichtern. Die Anweisungen `title` und `title3` definieren Titel für die nachfolgenden Ausgaben, insbesondere für die Prozeduren `PROC PRINT` und `PROC LIFETEST`.
Kopiert!
data prob7_7; set control rad radbpa; /*Combine all these three samples*/
grp1=(trt='rad');
grp2=(trt='radbpa');
title 'Problem 7_7, page 240 --Boron neutron capture therapy';
title3 'grp1=rad, grp2=radbpa';
;
1
DATA prob7_7; SET control rad radbpa; /*Combine all these three samples*/
2
grp1=(trt='rad');
3
grp2=(trt='radbpa');
4
title 'Problem 7_7, page 240 --Boron neutron capture therapy';
5
title3 'grp1=rad, grp2=radbpa';
6
;
5 Codeblock
PROC PRINT
Erklärung : Diese Prozedur zeigt den Inhalt des zuletzt erstellten oder aktiven Datensatzes an, in diesem Kontext `prob7_7`. Sie wird für eine schnelle Überprüfung der kombinierten Daten und zur Verifizierung der Korrektheit der Gruppenzusammenführung verwendet.
Kopiert!
proc print;
1
PROC PRINT;
6 Codeblock
PROC LIFETEST
Erklärung : Diese Prozedur führt eine Überlebenszeitanalyse auf dem Datensatz `prob7_7` durch. `plots=(s,h)` und `graphics` generieren Grafiken der Überlebensfunktionen (S, Kaplan-Meier standardmäßig) und der kumulativen Risikofunktionen (H). `time time*censor(0)` spezifiziert, dass `time` die Zeitvariable und `censor` die Ereignisvariable ist, wobei 0 eine zensierte Beobachtung anzeigt (das Ereignis ist nicht eingetreten). `strata trt/adjust=bon` fordert eine stratifizierte Analyse nach der Variablen `trt` (Gruppen 'control', 'rad', 'radbpa') an und wendet die Bonferroni-Korrektur für multiple Vergleiche zwischen den Strata an, um die p-Werte anzupassen. Die Anweisungen `symbol1`, `symbol2`, `symbol3` passen das Aussehen der Kurven in den generierten Grafiken an, um eine bessere visuelle Unterscheidung der Gruppen zu ermöglichen.
Kopiert!
/*Compare three groups, plot the survival functions and hazard functions*/
proc lifetest data=prob7_7 plots=(s,h) graphics;
time time*censor(0);
strata trt/adjust=bon;
symbol1 v=none color=black line=1;
symbol2 v=none color=black line=2;
symbol3 v=none color=black line=3;
run;
1
/*Compare three groups, plot the survival functions and hazard functions*/
2
PROC LIFETESTDATA=prob7_7 plots=(s,h) graphics;
3
time time*censor(0);
4
strata trt/adjust=bon;
5
symbol1 v=none color=black line=1;
6
symbol2 v=none color=black line=2;
7
symbol3 v=none color=black line=3;
8
RUN;
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.
Urheberrechtsinformationen : SAS-Code für Problem 7.7 aus dem Handbuch 'Survival Analysis Techniques for Censored and Truncated Data'.
SAS und alle anderen Produkt- oder Dienstleistungsnamen von SAS Institute Inc. sind eingetragene Marken oder Marken von SAS Institute Inc. in den USA und anderen Ländern. ® zeigt die Registrierung in den USA an. WeAreCAS ist eine unabhängige Community-Site und nicht mit SAS Institute Inc. verbunden.
Diese Website verwendet technische und analytische Cookies, um Ihre Erfahrung zu verbessern.
Mehr erfahren.