悬赏500论坛币求 SAS MACRO 高手！ 第二部分！

%macro test(n);
    data test;
            %do i=1 %to &n;
                    do v&i =1 to &n;
                    %if &i=1 %then %do;
                           tmp&i=cats("*",v&i,"*");
                        %end;
                        %else %do;
                            if index(tmp%eval(&i-1), cats("*",v&i,"*")) then continue;
                                else tmp&i="*"||catx("*",of v1-v&i)||"*";
                        %end;
                        %if &i=&n %then output;;
                 %end;
                 %do i=1 %to &n;
                end;
          %end;
                  drop tmp:;
        run;
%mend;

%test(4)
%test(6)

1. * Note: call allperm 最多可以对18个变量进行排列. Don't try it!;
   
2. 
   
3. %let n=4;
   
4. data PermFact;
   
5.    array x[&n] (1:&n);
   
6.      nfact=fact(&n);
   
7.       do i=1 to nfact;
   
8.       call allperm(i, of x[*]);
   
9.      output;
   
10.     end;
    
11.         keep x1-x&n;
    
12. run;
    
13. proc sort data=PermFact;
    
14.   by x1-x&n;
    
15. run;

复制代码

%macro perm(r) / parmbuff;   
      %let i=2;               
      %let things=;                     
      %do %while (%scan(&syspbuff,&i) ne );     
        %let p&i="%scan(&syspbuff,&i)";         
        %if &i=2 %then %let things=&&p&i;   
        %else %let things=&things,&&p&i;        
        %let i=%eval(&i+1);
      %end;
      %let n=%eval(&i-2);
      data perm;
         drop i j copy;
         array check (*) $ 10 r1-r&r;        
          %do m=1 %to &r;                 
            do r&m = &things;
          %end;
          copy=0;
            do i=2 to &r;               
              do j=1 to i-1;            
                if check(j)=check(i) then copy+1;
              end;
            end;
          if copy = 0 then output;        
          if copy = 0 then put r1-r&r;      
          %do m=1 %to &r;
            end;                             
          %end;
        run;
   %mend perm;
%perm(5,1,2,3,4,5);

yongyitian 发表于 2015-4-9 20:00

感谢您的回答！看了您的代码查了一下才知道了allperm这个功能的存在。
不知道有没有什么办法可以不依靠allperm这个功能来实现。

jl60156 发表于 2015-4-10 03:12
%macro perm(r) / parmbuff;   
      %let i=2;               
      %let things=;                  ...

感谢您的回答！有些地方看的一知半解。
最大的问题在于呼出宏时只能给出一个数值而非给定一串数据。

silentzilch 发表于 2015-4-10 10:41
感谢您的回答！看了您的代码查了一下才知道了allperm这个功能的存在。
不知道有没有什么办法可以不依靠a ...

楼上的 macro 不用 call allperm routine。
也可用于其他字符串, 并指定需要排列得字符串得个数. 例如 %permute(3,aaa, bbb, ccc,ddd)
http://support.sas.com/documentation/cdl/en/lrdict/64316/HTML/default/viewer.htm#a002208648.htm

http://www.urz.uni-heidelberg.de ... /ts498-combperm.txt

yongyitian 发表于 2015-4-10 11:37
楼上的 macro 不用 call allperm routine。
也可用于其他字符串, 并指定需要排列得字符串得个数. 例如  ...

感谢您的回答！

pobel 发表于 2015-4-10 13:27
%macro test(n);
    data test;
            %do i=1 %to &n;

感谢您的回答！

%let n=5  ;   /*n should <=18*/
proc iml;
  n = &n;
  p = allperm(n);
  print p;  /* this is the result*/
create jiecheng var{p};   
   append;
   show contents;
   close jiecheng;
run;

/* the following step is just for output dataset ;*/
data result;
  set jiecheng;
  array xxx[&n] a1-a&n;
  retain xxx;
  count+1;
  xxx[count]= p;
  if count=&n then do;
    call missing (count);
    output;
  end;
  keep a: ;
run;