data t_value;
file print;
k=8;
put  'v\单侧' @(k*1) '0.25' @(k*2) '0.20' @(k*3) '0.10' @(k*4) '0.05' @(k*5) '0.025'
     @(k*6) '0.01' @(k*7) '0.005' @(k*8) '0.0025' @(k*9) '0.001' @(k*10) '0.0005';
put  '  双侧' @(k*1) '0.05' @(k*2) '0.4' @(k*3) '0.20' @(k*4) '0.10' @(k*5) '0.05'
     @(k*6) '0.02' @(k*7) '0.01' @(k*8) '0.005' @(k*9) '0.002' @(k*10) '0.001';
do i=1 to 40;
   p=0;
   do j=0.25,0.2,0.1,0.05,0.025,0.01,0.005,0.0025,0.001,0.0005;
   p+1;
   t=round(abs(tinv(j,i)),0.001);
   t=trim(left(t));
   put  #(i+2) i @ (k*p) t @;
   end;
end;
q=i+2;
do i=50 to 90 by 10;
   p=0;
   q+1;
   do j=0.25,0.2,0.1,0.05,0.025,0.01,0.005,0.0025,0.001,0.0005;
   p+1;
   t=round(abs(tinv(j,i)),0.001);
   t=trim(left(t));
   put  #(q) i @ (k*p) t @;
   end;
end;
s=q+1;
do i=100,200,500,1000,10000;
   p=0;
   s+1;
   do j=0.25,0.2,0.1,0.05,0.025,0.01,0.005,0.0025,0.001,0.0005;
   p+1;
   t=round(abs(tinv(j,i)),0.001);
   t=trim(left(t));
   put  #(s) i @ (k*p) t @;
   end;
end;
title 't分布界值表';
run;