DM 'LOG;CLEAR;OUT;CLEAR;';
OPTIONS NODATE NONUMBER LS=96 PS=60;
/*
FILENAME TEMP 'EWMA1.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 sigma0=.005
                          outtable=ewma
                          xsymbol = mu0
                          table tableoutlim
                          haxis=1 to 25;
     label diameter='Diameter Ewma'
           sample = 'Piston Ring Sample';
TITLE 'EWMA Chart for Piston-Ring Diameters';
TITLE2 'lambda weight=0.25 (mu, sigma known)';

proc macontrol data=piston;
      ewmachart diameter*sample / weight  = 0.25
                          mu0 = 74 sigma0=.005
                          outtable=ewma
                          xsymbol = mu0
                          table tableoutlim reset
                          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 (mu, sigma known)';

run;