options ls=75;

  
data powercomp;
   drop mu1-mu5;
   t=5;
do sigma=1 to 2 by .5;
do n=2 to 10;
   df1=t-1;   /* numerator df */
   df2=(n-1)*t;
   array muvec{5} mu1-mu5 (12,12,11,10,9);
   sumtau2=css(of mu1-mu5);  /* css is a built-in SAS function for corrected sum of squares */

   ncp=sumtau2*n/(sigma*sigma);     /* noncentrality ("shift") parameter */

   critval=finv(.95,df1,df2);     /* critical value */
   beta = probf(critval,df1,df2,ncp);
   power=1-beta;
   output;
end;
end;
run;
proc print  data=two;run;
symbol value=dot i=join;
proc gplot;
   *plot power*n;
   title "HW 3, Problem 2";
   plot power*n=sigma;
run;