*//power and sample size//;

data one;
  u0=8;
  std=sqrt(2.5);
  alpha=.05;
  n=10;
  *ua=16.5;
  zcl=probit(alpha/2);
  zcu=probit(1-(alpha/2));

 do n=10,50,100;
 *do std=1 to 5 by 1;
  do ua=4 to 12 by 0.5;
   xcl=u0 + zcl*(std/sqrt(n));
   xcu=u0 + zcu*(std/sqrt(n));
   zal= (xcl-ua) / (std/sqrt(n));
   zau= (xcu-ua) / (std/sqrt(n));
   effect=abs((ua-u0)/std);
   beta=probnorm(zau) - probnorm (zal);
   power=1-beta;
   output;
  end;
 end;

*proc print;
  *var n alpha beta u0 ua std effect;
*run;

proc gplot;
  symbol i=join v=star w=3;
  plot power*ua=n / vaxis=0 to 1 by .1 vminor=0;
run;