data oc;
    keep beta nsample t plot2 arl;
    plot2=.;
    do nsample=1, 2, 3, 4, 5, 10, 15, 20;
       do j=0 to 400;
          t=j/100;
          beta=probnorm( 3-t*sqrt(nsample)) -
               probnorm(-3-t*sqrt(nsample));
          arl = 1/(1-beta);
          output;
          end;
       end;
    label t   ='Shift in Population Mean (Unit=Std Dev)'
         beta='Probability of Not Detecting Shift';

 symbol1 v=none i=join l=1  w=3 c=black;
 symbol2 v=none i=join l=2  w=3 c=black;
 symbol3 v=none i=join l=8  w=3 c=black;
 symbol4 v=none i=join l=20 w=3 c=black;
 symbol5 v=none i=join l=1  w=3 c=black;
 symbol6 v=none i=join l=2  w=3 c=black;
 symbol7 v=none i=join l=8  w=3 c=black;
 symbol8 v=none i=join l=20 w=3 c=black;

    axis1 order=(0.0 to 1.0 by 0.05)
          offset=(0,0);
    axis2 order= (0 to 4 by 0.5)
          minor= (n=3)
          offset=(0,0);
    legend1 label=('Sample Size n:');

 proc gplot data=oc;
    plot beta*t=nsample / frame
       vminor = 1
       hminor = 1
       lvref  = 2
       lhref  = 2
       vaxis  = axis1
       haxis  = axis2
       autovref
       autohref
       vzero
       hzero
       legend = legend1;
 title1 'OC Curves for Xbar Charts -- mu, sigma known';

    axis3 order=(0.0 to 380 by 20)
          offset=(0,0);
    axis4 order= (0 to 1.5 by 0.25)
          minor= (n=3)
          offset=(0,0);
    legend1 label=('Sample Size n:');

 proc gplot data=oc;
    plot arl*t=nsample / frame
       vminor = 3
       hminor = 1
       lvref  = 2
       lhref  = 2
       vaxis  = axis3
       haxis  = axis4
       autovref
       autohref
       vzero
       hzero
       legend = legend1;
 title1 'ARL Curves for Xbar Charts -- mu, sigma known';
 run;
 goptions reset=all;