/*
    TEC02_02.SAS
    Ｆ分布表を作成する2
    （５％と１％の表を別々に出力するバージョン）
    by K.T. Ver1.0  2002/07/08（月）15:50:18
*/

DATA d1;
        file 'f_table2.lst';
        prob=5; * 計算する確率 ここを 1 に変更すると1%の表になる  ;  
        pprob=1-prob/100;
        head1=' f2＼'; tate='f1';
        i=0; j=0 ; 
        PUT 78*' '@; PUT 'ＳＡＳによるＦ分布表（' @ ; PUT prob @; PUT  '%）'; PUT;
 
        * 見出し部分 ;
	PUT '     (分子の自由度) ';
        PUT head1 5.0  tate 2.0 @;
        DO i=1,2,3,4,5,6,7,8,9,10,11,12,14,16,20,24,30,40,60,120;
            PUT i 8.0 @; 
        END; PUT '     ∞';
        * 桁数  5+2+ 21*8=175 ;
        PUT 175*'-';

        * 数表本体 ;
        DO j=1 TO 30 BY 1;
            PUT j 5.0 @ ; PUT '| ' @;
            * ５％ ;
            DO i=1,2,3,4,5,6,7,8,9,10,11,12,14,16,20,24,30,40,60,120,1000000;
               F5=finv(pprob , i , j);
               PUT F5 8.2 @;
            END;PUT;
        END;

        DO j=40,60,120,1000000 ;
            IF j<1000000 THEN DO ;
                PUT j 5.0 @ ; PUT '| ' @; 
            END; ELSE DO;
                PUT '   ∞' @ ; PUT '| ' @; 
            END;
            * ５％ ;
            DO i=1,2,3,4,5,6,7,8,9,10,11,12,14,16,20,24,30,40,60,120,1000000;
               F5=finv(0.95 , i , j);
               PUT F5 8.2 @;
            END;PUT;
        END;
        
        PUT 175*'-';
        PUT ' ';
        PUT '自由度f1,f2から上側確率' @; PUT prob @ ; PUT '%に対するＦの値を求める表 ';
RUN;