Ungleichgewichtete Varianzanalyse mit EFFECTPLOT (ICOMEP2)

Schwierigkeitsgrad

Anfänger

Veröffentlicht am : 21/02/2021

Dieses Skript veranschaulicht, wie ein zweifaktorielles Versuchsdesign (Medikament und Krankheit) mit Interaktionen auf ungleichgewichteten Daten analysiert wird. Es verwendet die GENMOD-Prozedur, um das Modell anzupassen, und die EFFECTPLOT-Anweisung, um verschiedene Diagramme zu generieren: Effektkurven, Interaktionen mit Aufteilung (slice), Boxplots und Mosaikdiagramme.

Datenanalyse

Type : CREATION_INTERNE

Die Daten werden dynamisch in der DATA-Schritt 'a' über die DATALINES-Anweisung generiert. Sie simulieren eine klinische Studie mit den Variablen 'drug', 'disease' und der Antwort 'y'.

1 Codeblock

DATA STEP Data

Erklärung :
Erstellung der Datentabelle 'a'. Der DATA-Schritt liest die Medikamenten- und Krankheits-IDs und iteriert dann, um bis zu 6 'y'-Messungen pro Eingabezeile zu lesen.

Kopiert!

1	DATA a;
2	INPUT drug disease @;
3	DO i=1 to 6;
4	INPUT y @;
5	OUTPUT;
6	END;
7	DATALINES;
8	1 1 42 44 36 13 19 22
9	1 2 33 . 26 . 33 21
10	1 3 31 -3 . 25 25 24
11	2 1 28 . 23 34 42 13
12	2 2 . 34 33 31 . 36
13	2 3 3 26 28 32 4 16
14	3 1 . . 1 29 . 19
15	3 2 . 11 9 7 1 -6
16	3 3 21 1 . 9 3 .
17	4 1 24 . 9 22 -2 15
18	4 2 27 12 12 -5 16 15
19	4 3 22 7 25 5 12 .
20	;

2 Codeblock

PROC GENMOD

Erklärung :
ODS-Grafiken aktivieren. Ausführung von PROC GENMOD zur Modellierung der Antwortvariable 'y' (Standard: Normalverteilung). 'effectplot / obs' zeigt das Effekt-Diagramm mit den Beobachtungen. 'effectplot interaction(sliceby=disease)' zeigt Interaktionen, aufgeteilt nach der Variablen 'disease', mit Konfidenzintervallen (clm).

Kopiert!

1	ods graphics on;
2	PROC GENMOD DATA=a;
3	class drug disease;
4	model y=disease drug disease*drug / d=n;
5	effectplot / obs;
6	effectplot interaction(sliceby=disease) / clm;
7	RUN;

3 Codeblock

PROC GENMOD

Erklärung :
Zweite Analyse mit zusätzlichen Visualisierungen: Boxplots, gekreuzte Interaktionen mit Beobachtungen, Mosaikdiagramme und Interaktionen, facettiert nach 'disease' (plotby).

Kopiert!

1	PROC GENMOD DATA=a;
2	class drug disease;
3	model y=drug disease drug*disease / d=n;
4	effectplot box;
5	effectplot interaction(x=drug*disease) / obs;
6	effectplot mosaic;
7	effectplot interaction(plotby=disease);
8	RUN;
9	ods graphics off;

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.

Zugehörige Dokumentation

Aucune documentation spécifique pour cette catégorie.