[求助] PROC SQL 生成关于同一变量各观测的排列组合

if you use SAS 9.4,  call lexperk routine will make it easy.

data a;
   length new_var $3;
   array x[3] $3 ('a' 'b' 'c' );
   n=dim(x);

   k=1;
   nperm=perm(n, k);
   do j=1 to nperm;
      call lexperk(j, k, of x[*]);
      new_var=x[1];
          output;
   end;

   k=2;
   nperm=perm(n, k);
   do j=1 to nperm;
      call lexperk(j, k, of x[*]);
      new_var=catt(x[1],x[2]);
          output;
   end;

   k=3;
   nperm=perm(n, k);
   do j=1 to nperm;
      call lexperk(j, k, of x[*]);
      new_var=catt(x[1],x[2],x[3]);
          output;
   end;
   
   keep new_var;
run;

proc print;run;

bobguy 发表于 2014-4-6 10:11
if you use SAS 9.4,  call lexperk routine will make it easy.

data a;

Sorry, but I'm currently using SAS 9.3 ^^

我重新了写了一下，改成通用版本了。例如：计算a、b、c、d。
你仔细看看递归，找到规律。
看明白了，就好写宏了。

1. data a;
   
2.         input x $;
   
3. datalines;
   
4. a
   
5. b
   
6. c
   
7. d
   
8. ;
   
9. 
   
10. proc sql noprint;
    
11. 
    
12.         create table b as
    
13. 
    
14.           select catt(x1,x2,x3,x4) as x  from (
    
15. 
    
16.                 select x as x1 from a
    
17. 
    
18.                 outer union
    
19.                                 
    
20.                 /* 第一次递归*/
    
21.                 select catt(l.x,r.x) as x2 from (
    
22.                         select *  
    
23.                                 from a as l, a as r
    
24.                                 where r.x not like catt("%",l.x,"%")
    
25.                 )
    
26. 
    
27.                 outer union
    
28. 
    
29.                 /* 第二次递归*/         
    
30.                 select catt(l.x,r.x) as x3 from (
    
31.                         select * from a as l, (
    
32.                                 select catt(l.x,r.x) as x from a as l,a as r
    
33.                                 where r.x not like catt("%",l.x,"%")
    
34.                         ) as r
    
35.                         where r.x not like catt("%",l.x,"%")
    
36.                  )
    
37. 
    
38.                 outer union
    
39. 
    
40.                 /* 第三次递归*/         
    
41.                 select catt(l.x,r.x) as x4 from (
    
42.                         select * from a as l, (
    
43. 
    
44.                                         select catt(l.x,r.x) as x from (
    
45.                                                 select * from a as l, (
    
46.                                                         select catt(l.x,r.x) as x from a as l,a as r
    
47.                                                         where r.x not like catt("%",l.x,"%")
    
48.                                                 ) as r
    
49.                                                 where r.x not like catt("%",l.x,"%")
    
50.                                           )
    
51. 
    
52.                         ) as r
    
53.                         where r.x not like catt("%",l.x,"%")
    
54.                    )
    
55. 
    
56.                );
    
57. 
    
58. quit;

复制代码