/* calculate quantiles for t distributions */

/* input: df and probability (e.g. 0.975 to compute 95% ci */
data Fstat;
  input df p;
  qt =  tinv(p,df);
cards;
100 0.975
5 0.995
128 0.975
;

proc print;
  title 'quantiles for t statistics';
run;