/* pdfcf.sas                        */
 /*                                  */
 /* compute central F probability    */
  options ls=80 ;
  data a ;
    input ndf ddf calph p0 ; 
    p1 = cdf('F',calph,ndf,ddf) ;    * generic cdf function ;
    p2 = probf(calph,ndf,ddf) ;      * cdf for central F ;
 cards ;
   4  4  4.11  .90
   3  5  3.62  .90
   3  5  7.76  .975
   4  5  3.52  .90
   4  5  7.39  .975
   3  6  3.29  .90
   3  6  6.60  .975
   4  6  3.18  .90
   4  6  6.23  .975 
 ;
  run ;
   proc print data=a ;
    title 'some central F probabilities' ;
 run ;