/* Saved as probability in SAS */
option nocenter ls=76;

/* Distribution functions */

    data dist;input name1 $ name2  $ df1 df2 x p;
    IF name1='normal' THEN p=PROBNORM(x);
    IF name1='inverse' AND name2='normal' THEN x=PROBIT(p);
    IF name1='t' THEN p=PROBT(x,df1);
    IF name1='inverse' AND name2='t' THEN x=TINV(p,df1);
    IF name1='chisqr' THEN p=PROBCHI(x,df1);
    IF name1='inverse' AND name2='chisqr' THEN x=CINV(p,df1);   
    IF name1='F' THEN p=PROBF(x,df1,df2);
    IF name1='inverse' AND name2='F' THEN x=FINV(p,df1,df2);   
    CARDS;

    normal . . .   1.96 .
    inverse normal . . .  0.95 
    t .  5 .  1.476 .
    inverse t 10 . .  0.95
    chisqr  . 7 .  12.0 .
    inverse chisqr 11 . . 0.975
    F     . 5 10 3.33 .
    inverse F 3 5 . 0.95
    inverse chisqr 1.56 . . 0.975
    inverse chisqr 1.56 . . 0.025
    inverse chisqr 2 . . 0.975
    inverse chisqr 2 . . 0.025
    proc print;VAR name1 name2  df1 df2 x p;
    title1 '  SAS for z, t, chi-square and F distributions';
    title2 ' The relation is Prob(Y < x) = p ';