* SAS code to find Multivariate t critical points;
 * Final PROBMC params. are the square roots of rho's;
 options linesize=70;

 data tcrit;
     rcon1 = sqrt(.4765);
      rcon2 = rcon1;
       rcon3 = rcon1;
 t331 = probmc('dunnett2', . , .95, 31, 3, OF rcon1-rcon3);

     lambd1 = sqrt(.4714);
      lambd2 = sqrt(.4714);
       lambd3 = sqrt(.4867);
 t331r = probmc('dunnett2',.,.95,31,3,OF lambd1-lambd3);

 cards;

 proc print;
 run;