DM 'LOG;CLEAR;OUT;CLEAR;';
OPTIONS NODATE NONUMBER PS=66 LS=108;

DATA PISTON;
     DO SAMPLE = 1 TO 25;
     DO ITEM = 1 TO 5;
        INPUT DIAMETER @@;
        DIAMETER = DIAMETER + 70;
        OUTPUT;
     END; END;
 CARDS;
 4.030 4.002 4.019 3.992 4.008   3.995 3.992 4.001 4.011 4.004
 3.988 4.024 4.021 4.005 4.002   4.002 3.996 3.993 4.015 4.009
 3.992 4.007 4.015 3.989 4.014   4.009 3.994 3.997 3.985 3.993
 3.995 4.006 3.994 4.000 4.005   3.985 4.003 3.993 4.015 3.988
 4.008 3.995 4.009 4.005 4.004   3.998 4.000 3.990 4.007 3.995
 3.994 3.998 3.994 3.995 3.990   4.004 4.000 4.007 4.000 3.996
 3.983 4.002 3.998 3.997 4.012   4.006 3.967 3.994 4.000 3.984
 4.012 4.014 3.998 3.999 4.007   4.000 3.984 4.005 3.998 3.996
 3.994 4.012 3.986 4.005 4.007   4.006 4.010 4.018 4.003 4.000
 3.984 4.002 4.003 4.005 3.997   4.000 4.010 4.013 4.020 4.003
 3.988 4.001 4.009 4.005 3.996   4.004 3.999 3.990 4.006 4.009
 4.010 3.989 3.990 4.009 4.014   4.015 4.008 3.993 4.000 4.010
 3.982 3.984 3.995 4.017 4.013
 ;

PROC MEANS DATA=PISTON NOPRINT; BY SAMPLE;
     VAR DIAMETER;
     OUTPUT OUT=MSET MEAN=MEAN STD=STD RANGE=RANGE;

PROC PRINT DATA=MSET;
     ID SAMPLE; VAR MEAN STD RANGE;
TITLE 'PISTON RING DATA SUMMARY: MEANS, STANDARD DEVIATIONS AND RANGES';

SYMBOL V=DOT WIDTH=3;

PROC SHEWHART DATA=PISTON GRAPHICS;
     XCHART DIAMETER*SAMPLE / MU0=74 SIGMA0=.01 SIGMAS=3;
     RCHART DIAMETER*SAMPLE / SIGMA0=.01 SIGMAS=3;
     SCHART DIAMETER*SAMPLE / SIGMA0=.01 SIGMAS=3;
TITLE 'SHEWHART CHART for Mean Piston-Ring Diameters across Samples';
RUN;