This script first generates a test dataset containing various physical measurements classified by type (TEST). It then uses the MEANS procedure to calculate the mean, minimum, and maximum for each category. A complex DATA step constructs an annotation dataset (Annotate Facility) to manually draw horizontal error bars representing the Min-Max interval and mark the mean. Finally, the GPLOT procedure (SAS©/GRAPH) displays the result. Note: This code uses older graphical procedures (SAS©/GRAPH) instead of modern ODS Graphics procedures (SGPLOT).
Data Analysis
Type : CREATION_INTERNE
The data is defined directly in the script via the CARDS statement in the DATA step 'a'.
1 Code Block
DATA STEP Data
Explanation :
Initialization of global graphics options and creation of the source dataset 'a' containing the variables TEST (category) and BREAKS (value).
Initialization of global graphics options and creation of the source dataset 'a' containing the variables TEST (category) and BREAKS (value).
Copied!
| 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 Code Block
PROC SORT
Explanation :
Sorting data by the group variable 'TEST' to prepare for aggregation.
Sorting data by the group variable 'TEST' to prepare for aggregation.
Copied!
| 1 | /* Sort |
| 2 | data by variable TEST */ |
| 3 | PROC SORT; |
| 4 | BY TEST; |
| 5 | |
| 6 | RUN; |
| 7 |
3 Code Block
PROC MEANS Data
Explanation :
Calculation of descriptive statistics (mean, standard deviation, min, max) grouped by 'TEST'. Results are stored in output table 'b'.
Calculation of descriptive statistics (mean, standard deviation, min, max) grouped by 'TEST'. Results are stored in output table 'b'.
Copied!
| 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 Code Block
DATA STEP Data
Explanation :
Creation of a special 'Annotate' dataset. It defines graphical commands (move, draw) to draw blue lines from min to max, and red markers for the mean, min, and max. Uses the data coordinate system (xsys='2', ysys='2').
Creation of a special 'Annotate' dataset. It defines graphical commands (move, draw) to draw blue lines from min to max, and red markers for the mean, min, and max. Uses the data coordinate system (xsys='2', ysys='2').
Copied!
| 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 Code Block
PROC GPLOT
Explanation :
Generation of the final graph. The GPLOT procedure uses the 'anno' dataset to superimpose custom drawings on the defined axes. The horizontal axis is set from 0 to 100.
Generation of the final graph. The GPLOT procedure uses the 'anno' dataset to superimpose custom drawings on the defined axes. The horizontal axis is set from 0 to 100.
Copied!
| 1 | axis1 order=(0 to 100 BY 10); |
| 2 | symbol1 i=none v=none c=black; |
| 3 | |
| 4 | PROC GPLOT DATA=b ; |
| 5 | plot test*mean / anno=anno haxis=axis1 href=30 60 90; /* The HREF= option draws reference lines */ |
| 6 | RUN; |
| 7 | QUIT; |
This material is provided "as is" by We Are Cas. There are no warranties, expressed or implied, as to merchantability or fitness for a particular purpose regarding the materials or code contained herein. We Are Cas is not responsible for errors in this material as it now exists or will exist, nor does We Are Cas provide technical support for it.
