********************************************************;
*** GENERATE NORMAL, EXPONENTIAL, AND GAMMA VARIATES ***;
***          FROM DISTRIBUTION WITH MEAN = 20        ***;
********************************************************;
DATA IN;
  DO N = 1 TO 250;
   _NORMAL = 20 + RANNOR(5510);     ** NORMAL(20,1)   **;
   _GAMMA  = .5*RANGAM(20921,40);   ** GAMMA(.5,40)   **;
    OUTPUT;
  END;

SYMBOL1 VALUE=DOT WIDTH=3 L=1;

TITLE 'PROCESS CAPABILITY COMPARISON';

PROC CAPABILITY DATA=IN GRAPHICS;
     VAR _NORMAL _GAMMA;
     SPEC LSL=17 USL=23 TARGET=20 ALPHA=.05;

     HISTOGRAM  _NORMAL _GAMMA / NORMAL(INDICES);
     INSET MEAN (5.3) STD='Std Dev' (5.3)
           SKEWNESS (5.3) KURTOSIS (5.3) /
           HEADER = 'Summary Statistics'  POS = NE;
     INSET N CP (4.2) CPK (4.2) CPM (4.2) / POS = NW;

     HISTOGRAM  _NORMAL _GAMMA / GAMMA(THETA=0 INDICES);
     INSET MEAN (5.3) STD='Std Dev' (5.3)
           SKEWNESS (5.3) KURTOSIS (5.3) /
           HEADER = 'Summary Statistics'  POS = NE;
     INSET N CP (4.2) CPK (4.2) CPM (4.2) / POS = NW;

     CDFPLOT    _NORMAL _GAMMA / NORMAL;
     QQPLOT     _NORMAL / NORMAL;
     PPPLOT     _NORMAL / NORMAL;
RUN;