SAS 新手求教问题-经管之家官网！

免费学术公开课,扫码加入

DATA MY;
array x(1000) x1-x1000;
array y(1000) y1-y1000;
do i=1 to 1000;
do j=1 to 1000;
x(i)=rannor(0); /*生成1000个正态分布的随机数x*/
y(j)=rannor(0); /*生成1000个正态分布的随机数y*/
end;
end;
a=mean(of x1-x1000);/*计算1000个数x的均数*/
b=std(of x1- x1000);/*计算1000个数x的标准差*/
c=mean(of y1-y1000);/*计算1000个数y的均数*/
d=std(of y1-y1000);/*计算1000个数y的标准差*/
output;
%let m1= a;%let s1= b;%let m2=c;%let s2=d;%let l=-c/5;%let u=c/5;/*定义宏变量*/
array nt(2)nt1-nt2;
array nr(2)nr1-nr2;
do i=1to 2;
do j=1to 2;
nt(i)=&m1+&s1*rannor(0);/*生成2个正态分布的随机数，其均数是a，标准差是b*/
nr(j)=&m2+&s2*rannor(0);/*生成2个正态分布的随机数，其均数是c，标准差是d*/
end;
end;
Mt=mean(of nt1-nt2);/*计算前2个数的均数*/
St=std(of nt1-nt2);/*计算前2个数的标准差*/
Mr=mean(of nr1-nr2);/*计算后2个数的均数*/
Sr=std(of nr1-nr2);/*计算后2个数的标准差*/
SSt=St**2*(2-1);/*计算前一组的离均差平方和*/
SSr=Sr**2*(2-1);/*计算后一组的离均差平方和*/
Sc=(SST+SSr)/(4-2);/*计算2组的合并方差*/
Se=sqrt(Sc*(1/2+1/2));/*计算2组的标准误*/
T1=((Mt-Mr)-&l)/Se;/*对下限l进行假设检验，计算T1*/
T2=(&u-(Mt-Mr))/Se;/*对上线u进行假设检验，计算T2*/
P1=probt(-T1,4-2);/*计算-T1对应的T分布左侧累计概率*/
P2=probt(-T2,4-2);/*计算－T2对应的T分布右侧累计概率*/
PROC PRINT;
VAR P1 T1 P2 T2; /*P1 T1 P2 T2的值*/
run;
怎么将上面的程序重复1000次,求出1000个P1P2T1T1T2?

请求高人指导下 谢谢

「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP