求救贴：如何提高下面代码的效率？

data final;
input _name_ $;
cards;
test
;
run;
%macro samp(n,m,t);
data temp(drop=i);
do i= 1 to &n;
a=i;
output;
end;
run;
%macro repeat;
proc surveyselect data=temp noprint out=rand method=srs sampsize=&m seed=0;
run;
proc transpose data=rand out=tran;
run;
data trantemp(drop=i);
set tran;
run;
data final(drop=_name_);
set final trantemp;
run;
%mend repeat;
%do i=1 %to &t;
%repeat;
%end;
%mend samp;
%samp(33,7,1000000);
这段代码的大概意思是，参数宏SAMP（n,m,t）从N个号码中无放回抽取M个号码，然后抽取t次，并把t次的结果放在FINAL的DATASET里；我打算运行100万次，但我连运行10000次都很吃力，请问各位大虾如何提高效率呢？

1. %MACRO samp(n,m,t);
   
2.         data final;
   
3.         retain t;
   
4.         array x x1-x&m;
   
5.         do t=1 to &t by 1;
   
6.                 m=&m;
   
7.                 do n=1 to &n by 1;
   
8.                         r=ranuni(123);
   
9.                         if r<=m/(&n-n+1) then do;
   
10.                                 x[&m-m+1]=n;
    
11.                                 m=m-1;
    
12.                         end;
    
13.                 end;
    
14.         output;
    
15.         end;
    
16.         drop m n r;
    
17.         run;
    
18. %MEND samp;
    
19. %samp(30,17,1000000);

复制代码

不知道是不是准确的无放回
效率4秒多一点
希望对你有帮助

我的目的只需要从N个号码中无放回抽取M个号码，然后抽取t次，并把t次的结果放在FINAL的DATASET里
，目标的T为100W即可。各位大可不按我的代码的思路。

6203479170 发表于 2010-8-1 20:45
data final;
input _name_ $;
cards;
test
;
run;
%macro samp(n,m,t);
data temp(drop=i);
do i= 1 to &n;
a=i;
output;
end;
run;
%macro repeat;
proc surveyselect data=temp noprint out=rand method=srs sampsize=&m seed=0;
run;
proc transpose data=rand out=tran;
run;
data trantemp(drop=i);
set tran;
run;
data final(drop=_name_);
set final trantemp;
run;
%mend repeat;
%do i=1 %to &t;
%repeat;
%end;
%mend samp;
%samp(33,7,1000000);
这段代码的大概意思是，参数宏SAMP（n,m,t）从N个号码中无放回抽取M个号码，然后抽取t次，并把t次的结果放在FINAL的DATASET里；我打算运行100万次，但我连运行10000次都很吃力，请问各位大虾如何提高效率呢？

First there is no need/point to use macro.
Second why you need a million samplings. It sounds way too much.

Changing reps=4 to 1000000 will do it. Good luck!



   data temp;
      do i = 1 to 33;
         output;
      end;
    run;

     proc surveyselect  data=temp out=rand method=srs sampsize=7 seed=1
         reps=4;
   
   run;

   proc print; run;