data _null_;
   input scale $ k1 obs_a1 n1 obs_a2 n2;
   /* set up the width of the desired confidence interval */
   ci_alph = .95;
   retain ci_alph;
   k2     = k1;

   if (obs_a1 > obs_a2) then do;
      /* necessary to reverse everything */
      tmp    = obs_a1;
      obs_a1 = obs_a2;
      obs_a2 = tmp;

      tmp = k1;
      k1  = k2;
      k2  = tmp;

      tmp = n1;
      n1  = n2;
      n2  = tmp;
      tmp = 1;
   end;
   else tmp = 0;


   w     = (1 - obs_a2)/(1 - obs_a1);
   ci    = 100*ci_alph;
   format ci 2.;

   df1 = n1 - 1;
   df2 = n2 - 1;

   p_null = probf(w,df2,df1);
   format scale $12.;
   format n1 n2 5.;
   format k1 k2 4.;
   format p_null 4.2;

   if (tmp = 0) then do;
      put scale n1 n2 p_null;
/*
      put "Given:";
      put "   Observed alpha 1 = " obs_a1;
      put "   n1               = " n1;
      put "   Observed alpha 2 = " obs_a2;
      put "   n2               = " n2;
      put "   N(items) 1       = " k1;
      put "   N(items) 2       = " k2;
      put "Probability of difference greater than observed on null";
      put "  hypothesis of no difference = " p_null;
*/
   end;
   else do;
      put scale n1 n2 p_null;
/*
      put "Given:";
      put "   Observed alpha 1 = " obs_a2;
      put "   n1               = " n2;
      put "   Observed alpha 2 = " obs_a1;
      put "   n2               = " n1;
      put "   N(items) 1       = " k2;
      put "   N(items) 2       = " k1;
      put "Probability of difference greater than observed on null";
      put "  hypothesis of no difference = " p_null;
*/
   end;

   cards;
EAT_bul 6 .19 48 .36 103
EAT_bul 6 .19 48 .16 33
EAT_bul 6 .19 48 .70 213
EAT_bul 6 .36 103 .16 33
EAT_bul 6 .36 103 .70 213
EAT_bul 6 .16 33 .70 213
EAT_diet 13 .30 48 .58 103
EAT_diet 13 .30 48 .82 33
EAT_diet 13 .30 48 .85 213
EAT_diet 13 .58 103 .82 33
EAT_diet 13 .58 103 .85 213
EAT_diet 13 .82 33 .85 213
EAT_oral 7 -.01 48 .38 103
EAT_oral 7 -.01 48 .48 33
EAT_oral 7 -.01 48 .61 213
EAT_oral 7 .38 103 .48 33
EAT_oral 7 .38 103 .61 213
EAT_oral 7 .48 33 .61 213
EAT_tot 26 .54 48 .70 103
EAT_tot 26 .54 48 .74 33
EAT_tot 26 .54 48 .87 213
EAT_tot 26 .70 103 .74 33
EAT_tot 26 .70 103 .87 213
EAT_tot 26 .74 33 .87 213
HAD_anx 7 .53 48 .35 103
HAD_anx 7 .53 48 .79 33
HAD_anx 7 .53 48 .82 213
HAD_anx 7 .35 103 .79 33
HAD_anx 7 .35 103 .82 213
HAD_anx 7 .79 33 .82 213
HAD_dep 7 .53 48 .41 103
HAD_dep 7 .53 48 .71 33
HAD_dep 7 .53 48 .75 213
HAD_dep 7 .41 103 .71 33
HAD_dep 7 .41 103 .75 213
HAD_dep 7 .71 33 .75 213
;

run;