Annotierte Statistische Visualisierung mit PROC GPLOT

Dieses Skript generiert zunächst einen Testdatensatz mit verschiedenen physikalischen Messungen, klassifiziert nach Typ (TEST). Anschließend wird die PROC MEANS Prozedur verwendet, um den Mittelwert, das Minimum und das Maximum für jede Kategorie zu berechnen. Ein komplexer DATA-Schritt erstellt einen Annotationsdatensatz (Annotate Facility), um manuell horizontale Fehlerbalken zu zeichnen, die den Min-Max-Bereich darstellen, und den Mittelwert zu markieren. Schließlich zeigt die PROC GPLOT Prozedur (SAS^©/GRAPH) das Ergebnis an. Hinweis: Dieser Code verwendet ältere Grafikprozeduren (SAS^©/GRAPH) anstelle der modernen ODS Graphics Prozeduren (SGPLOT).

Datenanalyse

Type : INTERNE_ERSTELLUNG

Die Daten werden direkt im Skript über die CARDS-Anweisung im DATA-Schritt 'a' definiert.

1 Codeblock

DATA STEP Data

Erklärung :
Initialisierung globaler Grafikoptionen und Erstellung des Quelldatensatzes 'a', der die Variablen TEST (Kategorie) und BREAKS (Wert) enthält.

Kopiert!

1	/* Set graphics options */
2	goptions reset=all cback=white border;
3
4	/* Set system options */
5	DATA a;
6	INPUT TEST $ BREAKS ;
7	CARDS;
8	Cold 5
9	Cold 12
10	Cold 14
11	Cold 22
12	Cold 52
13	Heat 20
14	Heat 25
15	Heat 10
16	Heat 22
17	Heat 47
18	Gases 12
19	Gases 25
20	Gases 33
21	Gases 48
22	Gases 24
23	Pressure 10
24	Pressure 12
25	Pressure 14
26	Pressure 22
27	Pressure 60
28	Xrays 20
29	Xrays 25
30	Xrays 14
31	Xrays 22
32	Xrays 29
33	Humidity 20
34	Humidity 25
35	Humidity 33
36	Humidity 40
37	Humidity 24
38	;
39	RUN;

2 Codeblock

PROC SORT

Erklärung :
Sortierung der Daten nach der Gruppenvariable 'TEST' zur Vorbereitung der Aggregation.

Kopiert!

1	/* Sort
2	data by variable TEST */
3	PROC SORT;
4	BY TEST;
5
6	RUN;
7

3 Codeblock

PROC MEANS Data

Erklärung :
Berechnung deskriptiver Statistiken (Mittelwert, Standardabweichung, Min, Max), gruppiert nach 'TEST'. Die Ergebnisse werden in der Ausgabetabelle 'b' gespeichert.

Kopiert!

1	/**************************************************/
2	/* Create an output data set, B using PROC MEANS */
3	/* that contain new variables, MEAN, STD, STDERR, */
4	/* MIN, and MAX. */
5	/**************************************************/
6
7	PROC MEANS mean std stderr min max DATA=a;
8	BY TEST;
9	OUTPUT out=b mean=mean min=min max=max;
10	RUN;

4 Codeblock

DATA STEP Data

Erklärung :
Erstellung eines speziellen 'Annotate'-Datensatzes. Dieser definiert grafische Befehle (move, draw), um blaue Linien vom Minimum zum Maximum zu zeichnen und rote Markierungen für den Mittelwert, das Minimum und das Maximum zu setzen. Verwendet das Datenkoordinatensystem (xsys='2', ysys='2').

Kopiert!

1	/****************************************************************/
2	/* Create an annotate data set, ANNO to draw the bars at +/- 1, */
3	/* 2, or 3 Standard Deviation or Standard Error of the mean. */
4	/****************************************************************/
5
6	DATA anno;
7	retain xsys ysys '2' when 'a';
8	LENGTH color function $8 ;
9	SET b;
10
11	/* Draw the horizontal line from min to max */
12	function='move'; xsys='2'; ysys='2'; yc=TEST; x=min; color='blue'; OUTPUT;
13	function='draw'; x=max; color='blue'; size=2; OUTPUT;
14
15	/* Draw the MEAN horizontal line making the SIZE bigger */
16	function='move'; xsys='2';ysys='2';yc=TEST;x=mean; color='red'; OUTPUT;
17	function='draw'; x=mean; ysys='9'; y=+2; size=4; OUTPUT;
18	function='draw'; x=mean; y=-4; size=4; OUTPUT;
19
20	/* Draw the line for the MIN value */
21	function='move';xsys='2';ysys='2';yc=TEST;x=min;color='red';OUTPUT;
22	function='draw';x=min;ysys='9';y=+2;size=2;OUTPUT;
23	function='draw';x=min;y=-4;size=2;OUTPUT;
24
25	/* Draw the line for the MAX value */
26	function='move';xsys='2';ysys='2';yc=TEST;x=max;color='red';OUTPUT;
27	function='draw';x=max;ysys='9';y=+2;size=2;OUTPUT;
28	function='draw';x=max;y=-4;size=2;OUTPUT;
29	RUN;

5 Codeblock

PROC GPLOT

Erklärung :
Generierung der finalen Grafik. Die PROC GPLOT Prozedur verwendet den Datensatz 'anno', um die benutzerdefinierten Zeichnungen über die definierten Achsen zu legen. Die horizontale Achse ist von 0 bis 100 festgelegt.

Kopiert!

1	axis1 order=(0 to 100 BY 10);
2	symbol1 i=none v=none c=black;
3
4	PROC GPLOT DATA=b ;
5	plot testmean / anno=anno haxis=axis1 href=30 60 90; / The HREF= option draws reference lines */
6	RUN;
7	QUIT;

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.

Retour à la liste