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

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请大家帮忙看一下，假设现在有1个变量（变量名 var1）包含3条观测（a、b、c），具体如下：

1. var1
   
2. a
   
3. b
   
4. c

复制代码

现在想用 proc sql 生成一个新的表，表中的行、列（列名 new_var）如下：

1. new_var
   
2. a
   
3. ab
   
4. abc
   
5. ac
   
6. acb
   
7. b
   
8. ba
   
9. bac
   
10. bc
    
11. bca
    
12. c
    
13. ca
    
14. cab
    
15. cb
    
16. cba

复制代码

亦即由原来的3条观测通过排列的方式生成 \[P_{3}^{1}+P_{3}^{2}+P_{3}^{3}=15\] 条新的行（观测）。想了很长时间也没有解决，希望大家帮忙看看，万分感谢。

PS: 当然不用 proc sql 也行，只要最后能生成上述列就可以了，谢谢大家啦^^

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：proc sql sql ROC LaTeX atex

我来膜拜高手。

完美实现。若把SQL部分改成宏，就可以通用。

1. data a;
   
2.         input x $;
   
3. datalines;
   
4. a
   
5. b
   
6. c
   
7. ;
   
8. 
   
9. proc sql noprint;
   
10. 
    
11.         create table b as
    
12.                 select * from a
    
13.                 outer union
    
14.                 select * from (
    
15.                         select *
    
16.                                 from a(rename=(x=x1)) as l, a(rename=(x=y1)) as r
    
17.                                 where l.x1<> r.y1
    
18.                 )
    
19. 
    
20.                 outer union
    
21.                 select * from (
    
22.                         select * from a(rename=(x=z)) as l, (
    
23.                                 select * from a(rename=(x=x2)) as l,a(rename=(x=y2)) as r
    
24.                                 where l.x2<> r.y2
    
25.                         ) as r
    
26.                         where (l.z<> r.x2 and  l.z<> r.y2)
    
27.                 );
    
28. 
    
29. quit;
    
30. 
    
31. data b(keep=x);
    
32.         set b;
    
33.         length x $ 3;
    
34.         x=catt(x,x1,y1,z,x2,y2);
    
35. run;
    
36. 
    
37. proc sort data=b;by x;run;

复制代码

1. /* datastep solution just for the case, not for all cases */
   
2. 
   
3. data one;
   
4.    input var1 $;
   
5. datalines;
   
6. a
   
7. b
   
8. c
   
9. ; run;
   
10. 
    
11. data _null_;
    
12. length new $30;
    
13.    set one;
    
14.        new = var1;
    
15.        put new=;
    
16.       do p= 1 to nobs;
    
17.        set one(rename=var1=var2) nobs=nobs point=p;
    
18.          if trim(var1) ^= trim(var2) then do; new = cats(var1, var2);
    
19.          put new=; end;
    
20.         do p1 = 1 to nobs;
    
21.          set one(rename=var1=var3) nobs=nobs point = p1;
    
22.            if  trim(var1) ^= trim(var2)
    
23.            and trim(var1) ^= trim(var3)
    
24.            and trim(var2) ^= trim(var3) then do; new = cats(var1, var2, var3);
    
25.            put new=; end;
    
26.          end;
    
27.      end;
    
28. run;

复制代码

data have;
input var $;
cards;
a
b
c
;
run;

data have1;
set have have;
run;


data _null_;
if 0 then set have(rename=(var=_var));
if _n_=1 then do;
  declare hash h(dataset:'have(rename=(var=_var))',ordered:'a');
  h.definekey('_var');
  h.definedata('_var');
  h.definedone();
  declare hiter hi ('h');
end;
set have1 end=last;
h.find(key:var);
do rc=hi.next() by 0 while (rc=0);
  if findc(compress(_var),compress(var))=0 then do;
   _var=cats(var,_var);
    h.ref();
  end;
rc=hi.next();
end;
if last then h.output(dataset:'want(rename=_var=var)');
run;
   

zhengbo8 发表于 2014-4-6 01:17
完美实现。若把SQL部分改成宏，就可以通用。

嗯嗯~ 非常感谢，正在尝试宏的改写 ^^

yongyitian 发表于 2014-4-6 01:34

非常感谢!

farmman60 发表于 2014-4-6 06:55
data have;
input var $;
cards;

非常感谢!

you can leverage the proc plan to do it.


proc format ;
value trmt
   1='a'
   2='b'
   3='c'
   other=' '
   ;
run;

proc plan;
   factors Block=3 ordered
           Treat= 1 of 3 perm;
   ods output Plan=Combinations1;
run;

proc plan;
   factors Block=6 ordered
           Treat= 2 of 3 perm;
   ods output Plan=Combinations2;
run;

proc plan;
   factors Block=6 ordered
           Treat= 3 of 3 perm;
   ods output Plan=Combinations3;
run;

data all;
  length new_var $3;
  set Combinations1
  Combinations2
  Combinations3
  ;
  new_var=catt(put(Treat1,trmt.),put(Treat2,trmt.),put(Treat3,trmt.));

  keep new_var;
  run;

  proc print;run;

zhengbo8 发表于 2014-4-6 01:17
完美实现。若把SQL部分改成宏，就可以通用。

...表示 macro 的技术不到家，虽然知道要做什么，但完全不知道怎么改（因为 select * from 子块的数量是递增的）==...

能否给点提示和帮助，谢谢！