OPTION LS=76 NOCENTER;

/* Model y=metal + ingot +e  */
/* Levels    3       7       */

%LET SEED      = 1234567;
%LET n         = 200; 
%LET c         = 112;   
%LET p0        = 0.5;
%LET p1        = 0.6; 
%LET alpha     = 0.05;
%LET power     = 0.99;
%let nsim      = 10000;

data key;
   sig=0.0;
   do sim = 1 to ≁
      y=0;
      do i=1 to &n;
         a=ranuni(&seed);
         IF (a<&p0) THEN y=y+1;
      end;     
         IF (y>&c) THEN sig=sig+1;
    end;
    sig=sig/&nsim; OUTPUT;
         
   power=0.0;
   do sim = 1 to &nsim;
      y=0;
      do i=1 to &n;
         a=ranuni(&seed);
         IF (a<&p1) THEN y=y+1;
      end;     
         IF (y>&c) THEN power=power+1;
    end;
    power=power/&nsim;OUTPUT;
run;
       
PROC PRINT DATA=key;