/* calculate p-values from = t distributions */

data Fstat;
  input df t;
  pt =  2*probt(-abs(t),df);
  
cards;
100 1.962
5 3.15
6 10.3
;

proc print;
  title 'P values for t statistics';
run;
        
proc print;
  format pt pvalue.;
  
  title 'P values for t statistics, nicer format';
run;