DM 'LOG;CLEAR;OUT;CLEAR;';
OPTIONS NODATE NONUMBER LS=72 PS=66;
/*
FILENAME TEMP 'CUSUM3.PS';
GOPTIONS DEVICE=PSEPSF GSFNAME=TEMP GSFLEN=132
VSIZE=4.75 HSIZE=6 GSFMODE=APPEND;
*/
DATA IN;
  DO SAMPLE = 1 TO 48;
     INPUT SOLIDS @@; OUTPUT;
  END;
 CARDS;
 43.7  44.4  45.0  44.1  46.4  43.6  46.2  43.5  44.5  46.3
 45.9  45.3  44.2  44.4  46.8  44.2  45.6  44.9  46.1  46.4
 43.8  44.3  44.5  46.0  47.2  46.1  45.9  45.3  46.8  45.1
 46.1  43.8  47.8  43.4  46.1  45.9  44.7  44.2  45.9  46.9
 45.8  47.1  44.6  47.6  44.6  46.1  45.8  44.9

SYMBOL1 VALUE=DOT WIDTH=3;

proc cusum data=in;
     xchart solids*sample
           / mu0=45 sigma0=1   h=4.0 k=0.5
             delta=1.0 dataunits tablesummary;
     label solids='Percent Solids Cusum'
           sample = 'Sample';
TITLE 'CUSUM for Percent Solids (h=4 k=.5)';

proc cusum data=in;
     xchart solids*sample
           / mu0=45 sigma0=1   h=8.0 k=0.25
             delta=1.0 dataunits tablesummary;
     label solids='Percent Solids Cusum'
           sample = 'Sample';
TITLE 'CUSUM for Percent Solids (h=8 k=.25)';

proc cusum data=in;
     xchart solids*sample
           / mu0=45 sigma0=1  alpha=.05 beta=.10
             delta=1.0 dataunits tablesummary;
     label solids='Percent Solids Cusum'
           sample = 'Sample';
TITLE 'CUSUM for Percent Solids (alpha=.05 beta=.10)';

proc cusum data=in;
     xchart solids*sample
           / mu0=45 sigma0=1  alpha=.01 beta=.10
             delta=1 dataunits tablesummary;
     label solids='Percent Solids Cusum'
           sample = 'Sample';
TITLE 'CUSUM for Percent Solids (alpha=.01 beta=.10)';

run;