Veröffentlicht am :
Statistik CREATION_INTERNE

Freeman-Tukey- und t-Tests mit Bootstrap-Resampling (PROC MULTTEST)

Dieser Code ist auch verfügbar auf: English Español Français
Wartet auf Validierung
Dieses Skript analysiert Daten zur Karzinogenität bei Tieren (Vorhandensein von Tumoren S1 und S2, Todeszeitpunkt T) in Abhängigkeit von der verabreichten Dosis. Es verwendet PROC MULTTEST mit der BOOTSTRAP-Option, um angepasste p-Werte zu berechnen, wobei ein Freeman-Tukey-Test für binäre Variablen und ein t-Test für die kontinuierliche Variable angewendet werden.
Datenanalyse

Type : CREATION_INTERNE


Die Daten werden über einen DATA STEP generiert, der DATALINES mit dem Operator ' @@' verwendet, um mehrere Beobachtungen pro Zeile zu lesen.

1 Codeblock
DATA STEP Data
Erklärung :
Erstellung des Datensatzes 'a'. Jede Beobachtung repräsentiert ein Tier mit zwei Tumorindikatoren (S1, S2), der Todeszeit (T) und der erhaltenen Dosis. Der Operator ' @@' wird verwendet, um mehrere Einträge aus einer einzigen Datenzeile zu lesen.
Kopiert!
1title 'Freeman-Tukey and t Tests with Bootstrap Resampling';
2 
3DATA a;
4 INPUT S1 S2 T Dose @code_sas_json/8_SAS_Intro_ReadFile_MultiCol_@@.json;
5 DATALINES;
60 1 104 1 1 0 80 1 0 1 104 1
70 1 104 1 0 1 100 1 1 0 104 1
81 0 85 2 1 0 60 2 0 1 89 2
91 0 96 2 0 1 96 2 1 0 99 2
101 0 60 3 1 0 50 3 1 0 80 3
110 1 98 3 0 1 99 3 1 0 50 3
12;
2 Codeblock
PROC MULTTEST Data
Erklärung :
Ausführung mehrerer statistischer Tests. Die Option 'bootstrap' aktiviert das Resampling (10.000 Stichproben, fester Seed). Ein Freeman-Tukey-Test (ft) wird für die binären Variablen S1/S2 und ein Mittelwerttest (mean) für T angegeben. Ein linearer Kontrast wird definiert, um den Trend über die Dosisstufen hinweg zu testen.
Kopiert!
1PROC MULTTEST DATA=a bootstrap nsample=10000 seed=37081 outsamp=res;
2 test ft(S1 S2 / lowertailed) mean(T / lowertailed);
3 class Dose;
4 contrast 'Linear Trend' 0 1 2;
5RUN;
3 Codeblock
PROC PRINT
Erklärung :
Anzeige der ersten 36 Beobachtungen des Ausgabedatensatzes 'res', der die Resampling-Ergebnisse enthält.
Kopiert!
1PROC PRINT DATA=res(obs=36);
2RUN;
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 SAMPLE LIBRARY - Example provided by Keith Soper, Merck