/*
    TEC02_01.SAS
    Ｆ分布表を作成する
    by K.T. Ver1.0  2002/07/08（月）15:50:18
*/

DATA d1;
        file 'f_table.lst';
        prob=0.05; * 計算する確率; 
        pprob=1-prob;
        head1=' f2＼'; tate='f1';
        i=0; j=0 ; 
        PUT 80*' '@; 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(0.95 , i , j);
               PUT F5 8.2 @;
            END;PUT;
            * １％ ;
            PUT 5*' ' @; 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(0.99 , i , j);
               PUT F5 8.2 @;
            END;PUT ;
            * 空白行 ;
            PUT 5*' ' @; PUT '| ' @; PUT 168*' ';
        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;
            * １％ ;
            PUT 5*' ' @; 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(0.99 , i , j);
               PUT F5 8.2 @;
            END;PUT ;
            * 空白行 ;
            IF j<100000 THEN DO ;
               PUT 5*' ' @; PUT '| ' @; PUT 168*' ';
            END;
        END;
        
        PUT 175*'-';
        PUT ' ';
        PUT '自由度f1,f2から上側確率５％,および１％の対するＦの値を求める表 ';
        PUT '上段は５％,下段は１％水準 ';
RUN;