data _null_;
   obs_a = .82;    /* observed alpha value in sample           */
   n     = 213;      /* n observations in sample                 */
   k     = 7;       /* n items in scale                         */
   p     = .95;     /* spread for confidence interval           */
   nul_a = .9;      /* null model population alpha to be tested */

   ci    = 100*p;
   format ci 2.;

   f     = (1 - obs_a) / (1 - nul_a);

   n1    = (n - 1) * (k - 1);
   n2    = n - 1;
   p_null = 1 - probf(f,n2,n1);
   p1    = (1 - p)/2;
   p2    = p + p1;
   f1    = finv(p1,n2,n1);
   f2    = finv(p2,n2,n1);

   lwr   = 1 - (1 - obs_a)*f2;
   upr   = 1 - (1 - obs_a)*f1;
   format lwr upr 5.2;

   put;
   put "Given:";
   put "   Observed alpha  = " obs_a;
   put "   n               = " n;
   put "   Number of items = " k;
   put "   and null hypothesis population alpha = " nul_a;
   put;
   put "P against one-tailed alternative populn. alpha of " nul_a " is " p_null;
   put ci "% confidence interval from " lwr " to " upr;
run;