生成t双侧分位数表。

shiyiming

题目如下：
生成t分布的双侧分位数表。水平取0.001, 0.002, 0.005, 0.01, 0.02, 0.05,0.10, 0.20，自由度1-100。分位数精确到小数点后3位，表格应为行列对齐格式，并有列标题，写出生成这样的表格并放到一个文本文件中的sas程序。

各位有什么好的做法？

shiyiming

[code:d7d28]
data prob;
input p;
cards;
0.001
0.002
0.005
0.01
0.02
0.05
0.10
0.20
;
run;
data result;
set prob;
do df=1 to 100;
result=-tinv(p/2,df);
output;
end;
run;
[/code:d7d28]
生成文本文件，我用proc tabulate，不过太土了，就不贴上来现眼了。是不是应该用proc report 和ods ？我在这方面很弱。  <!-- s:oops: --><img src="{SMILIES_PATH}/icon_redface.gif" alt=":oops:" title="Embarassed" /><!-- s:oops: -->

shiyiming

Data example.a;
do p=0.001, 0.002, 0.005, 0.01, 0.02, 0.05, 0.1, 0.2;
df=1;
  do until (df&gt;100);
    q1=TINV(p/2,df);
    q2=TINV((1-p/2),df);
    output;
   df=df+1;
  end;
end;
proc export data=example.a outfile='d:\sas\exercise\example\c.txt' DBMS='Tab' replace;
run;

shiyiming

Data example.a;
do p=0.001, 0.002, 0.005, 0.01, 0.02, 0.05, 0.1, 0.2;
  df=1;
  do until (df&gt;100);
     left=TINV(p/2,df);
     right=TINV((1-p/2),df);
     output;
      file 'd:\sas\exercise\example\c.txt';
      put p 1-5 .3 df 10-15 left 20-30 .3 right 35-45 .3;
     df=df+1;
end;
end;
run;