options ls=80;
data a;
p_value1=1-probf(2.5,2,20);  /* F-distribution */
p_value2=1-probchi(2.5,3);   /* Chisquare distribution */
p_value3=1-probt(2.5,16);    /*  t-distribution */


proc print data=a;
run;