SAS编写 t 界值表

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

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;

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖