/*
  tPvalue.sas
*/
options ls=80;
data myci;
  b1 = 0.12486;                   * slope estimate;
  SE = 0.01290;                   * Std. Error of b1;
  tcalc = 9.68;                 * t statistic value;
  df = 12;
  P_lower = probt(tcalc, df);  
  P_upper = 1-probt(tcalc, df);
  P_two_tail = 2*(1-probt(abs(tcalc),df));

* Note:  SAS version 9 uses 'CDF' as a generalization
         of 'probt';

proc print;
  run;