DM 'LOG;CLEAR;OUT;CLEAR;';
OPTIONS NODATE NONUMBER LS=96 PS=60;
/*
FILENAME TEMP 'EWMA2.EPS';
GOPTIONS DEVICE=PSEPSF GSFNAME=TEMP GSFLEN=132
VSIZE=4.75 GSFMODE=REPLACE;
*/
data piston;
     do sample=1 to 25;
     do item=1 to 5; n= item + 5*(sample-1);
        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
 ;
symbol v=dot width=3;

proc macontrol data=piston;
      ewmachart diameter*sample / weight  = 0.25
                          mu0 = 74
                          outtable=ewma
                          xsymbol = mu0  reset
                          table tableoutlim
                          outlimits = limset
                          haxis=1 to 25;
     label diameter='Diameter Ewma'
           sample = 'Piston Ring Sample';
TITLE 'EWMA Chart for Piston-Ring Diameters with Resets';
TITLE2 'lambda weight=0.25 (sigma unknown)';

proc print data=limset;
run;