SAS PROGRAMMING
                                    HANDOUT #14

     This handout will illustrate some of the features of PROC GPLOT and also show
how do to a normal probablity plot.

LIBNAME JIM 'A:\STAT534';
DATA DGNP; SET JIM.GNP; ONE=1;

PROC UNIVARIATE DATA=DGNP NOPRINT; VAR INVEST;
 OUTPUT OUT=AA PCTLPTS=0 20 40 60 80 100 PCTLPRE=INV_;
PROC PRINT DATA=AA; RUN;

DATA BB; SET AA; ONE=1;

DATA CC; MERGE DGNP BB; BY ONE;
 CINV1=0; CINV2=0; CINV3=0; CINV4=0; CINV5=0;
 IF INV_0<=INVEST<INV_20 THEN CINV1=1;
 IF INV_20<=INVEST<INV_40 THEN CINV2=1;
 IF INV_40<=INVEST<INV_60 THEN CINV3=1;
 IF INV_60<=INVEST<INV_80 THEN CINV4=1;
 IF INV_80<=INVEST<=INV_100 THEN CINV5=1;

 INV=1*CINV1+2*CINV2+3*CINV3+4*CINV4+5*CINV5;

PROC REG DATA=CC;  ID INV;                            * This will do regression and;
 MODEL GNP=INVEST/P;                                   * put the predicted and residuals;
 OUTPUT OUT=JIM R=R P=P;  RUN;                         * in the data set  JIM;

PROC SORT DATA=JIM OUT=JIM2; BY R;
DATA BOB; SET JIM2;
 NN=(_N_-3/8)/(126+.25);                              *Here 126 is the numbers of obs;
 Z=PROBIT(NN);

PROC GPLOT DATA=BOB;
 TITLE J=L 'NORMAL PROBABILITY PLOT';
 PLOT R*Z/ HAXIS=AXIS1  VAXIS=AXIS2 VREF=0;
  AXIS1 LABEL=('STANDARD NORMAL') ORDER=-3 TO 3 BY 1  ;
  AXIS2 LABEL=(ANGLE=90 ROTATE=0 'RESIDUAL VALUES')  ;
   SYMBOL V=o I=RL L=1 C=RED;
   FOOTNOTE J=L 'SAS EMPLOYMENT DATA';   RUN;

PROC UNIVARIATE PLOT NORMAL DATA=JIM; VAR R;    RUN;

PROC GCHART DATA=BOB;
 TITLE 'HISTOGRAM OF RESIDUALS';                   *Produces histogram of residuals;
 VBAR R/MIDPOINTS=-500 TO 1500 BY 300;        RUN;

           * If your program has two graphs then you have to page-down in the graph;
           * window to see the second graph;

PROC GCHART DATA=JIM;
 TITLE 'HISTOGRAM OF LOG SALARY';
 VBAR GNP/SUBGROUP=INV;
RUN;

-------------------------------------------------------------------------------------------------------

PROC GPLOT DATA=DGNP;
 TITLE 'INVESTMENT AND GOVERNMENT SPENDING';
 PLOT INVEST*DATE GOVT*DATE/OVERLAY;           *OVERLAY causes both plots to be on the same page;
 SYMBOL1 V=I H=.5 I=JOIN L=1 C=RED;
 SYMBOL2 V=G H=.5 I=JOIN L=2 C=BLACK;
RUN;