/* calculate p-values from F distributions */

data Fstat;
  input df1 df2 F;
  pf =  1-probf(f,df1,df2);
  
cards;
2 6 5.14
5 8 5.14
5 343 57.10
;

proc print;
  title 'P values for F statistics';
run;
        
proc print;
  format pf pvalue.;
  
  title 'P values for F statistics, nicer format';
run;