请看看这个问题

有   a1 a2 a3  b c d e1 e2 几个，现在求他们的宏变量组合 ，其中  a1 a2 a3 a3 其中必须选两个  e1 e2 其中必须选一个 ，b c d 是任意个数任意形式的组合 （ 包括三个中任意选一个  选两个 选三个 ），请生成在一个数据集中吧，变量取名为zuhe 。

小鳄鱼a 发表于 2015-12-14 19:27
非常感谢   请问   假如有这样一个数据集   我想提取他们的先前预定的各种组合   那么该怎样提取  一种组 ...

哦，如果已经有了组合，那么确实可以利用keep来做一下，方法分为两种：1）只用data部来实现，需要用到call execute语句。 2）利用宏。

1. /* 1.利用call execute 在一个data部中生成 */
   
2. 
   
3. data _null_;
   
4.   set all_combinations;  * 上一次生成的;
   
5. 
   
6.   call execute("data comb" || strip(_N_) || ";" ||
   
7.                "  set test;                   " ||
   
8.                "  keep " || zuhe || ";"         ||
   
9.                "run;"
   
10.               );
    
11. run;
    
12. 
    
13. /* 2.借助宏语句生成 */
    
14. 
    
15. %macro create_all_comb;
    
16.   proc sql
    
17.     noprint;
    
18.     select count(*) into :comb_count
    
19.     from all_combinations
    
20.     ;
    
21.   quit;
    
22. 
    
23.   %do i = 1 %to &comb_count.;
    
24.     data _null_;
    
25.       set all_combinations(firstobs = &i. obs = &i.);
    
26.       call symput("keep_var", strip(zuhe));
    
27.     run;
    
28. 
    
29.     data comb&i.;
    
30.       set test;
    
31.       keep &keep_var.;
    
32.     run;
    
33.   %end;
    
34. %mend;
    
35. 
    
36. %create_all_comb;

复制代码

用哪种，看楼主的喜好了。

请确保他们每一个都是宏变量   

1. %macro x(seed = 456,s2 = a1 a2 a3,s1 = e1 e2,r = b c d,rn = 2);
   
2.   data temp;
   
3.     do _n_ = 1 to countw("&s2. &s1. &r.");
   
4.       value = scan("&s2. &s1. &r.", _n_);
   
5.           order = ranuni(&seed.);
   
6.           select;
   
7.         when(_n_ <= countw("&s2.")) name = "s2";
   
8.         when(_n_ <= countw("&s2. &s1")) name = "s1";
   
9.                 otherwise name = "r&rn.";
   
10.           end;
    
11.           output;
    
12.         end;
    
13.   run;
    
14. 
    
15.   proc sort;
    
16.     by name order;
    
17.   run;
    
18. 
    
19.   data _null_;
    
20.     set temp;
    
21.     by name;
    
22. 
    
23.         length want $20.;
    
24.         retain want;
    
25.         if first.name then do;
    
26.       want = '';
    
27.           count = 0;
    
28.         end;
    
29. 
    
30.         want = catx(' ',want,value);
    
31.         count + 1;
    
32. 
    
33.         if cats(count) = compress(name,,'kd') then call symputx(name,want);
    
34.   run;
    
35. 
    
36.   %put s2 = &s2. s1 = &s1. r2 = &r2. ;
    
37. %mend;
    
38. 
    
39. %x;

复制代码

方法比较笨拙，如有更好的应该采纳更好的。

Tigflanker 发表于 2015-12-7 12:40
方法比较笨拙，如有更好的应该采纳更好的。

这个好像没形成组合吧   另外   做这个好像一般是用 ncomb  就是有些复杂    所以不会做了  。  

小鳄鱼a 发表于 2015-12-7 15:31
这个好像没形成组合吧   另外   做这个好像一般是用 ncomb  就是有些复杂    所以不会做了  。

奥，你是打算把所有组合罗列到一个数据集中啊

可是为什么要叫做‘宏变量组合’啊？

Tigflanker 发表于 2015-12-7 15:35
奥，你是打算把所有组合罗列到一个数据集中啊

可是为什么要叫做‘宏变量组合’啊？

我想着到时候可以引用   表达可能不清楚   比如如果能够排列成 88个组合的话   假如用data a_i ; keep + 宏变量组合。 利用i ,这样可以循环出 88个组合

看了一下，是好像可以通过comb做出，但对于allcomb这个函数在PDV中的运行机制，还是不太理解

小鳄鱼a 发表于 2015-12-7 20:27
我想着到时候可以引用   表达可能不清楚   比如如果能够排列成 88个组合的话   假如用data a_i ; keep +  ...

楼主你好，我写了个代码，但是根据我的理解，只有42中组合，请楼主确认下我说的对不对：
1） a1 a2 a3中必须选两个，也就是说只能选择2个，那么它的组合数是3
2） e1 e2中必须选一个，组合数是2
3）b c d 任意组合形式那么可能的状况是7
4）3 x 2 x 7 = 42

不太确定我理解的是否正确，我觉得其中有歧义的地方在于1）和2） “必选两个”这个说法，是只能选择两个，还是说2个及以上？如果我的理解正确，下面的代码应该就能满足楼主的需求了。

1. data set_a(keep = zuhe_a);
   
2.   array x(3) $2. ('a1' 'a2' 'a3');
   
3. 
   
4.   n = dim(x);
   
5. 
   
6.   k = 2;
   
7. 
   
8.   ncomb = comb(n, k);
   
9. 
   
10.   do i = 1 to ncomb;
    
11.     rc = allcomb(i, k, of x[*]);
    
12.     zuhe_a = strip(x1) || " " || strip(x2);
    
13.         output;
    
14.   end;
    
15. run;
    
16. 
    
17. %macro create_set_bcd;
    
18.   %do i = 1 %to 3;
    
19.     data set_bcd_&i.(keep = zuhe_bcd);
    
20.       array x(3) $1. ('b' 'c' 'd');
    
21.       retain ncomb n;
    
22. 
    
23.       n = dim(x);
    
24.       ncomb = comb(n, &i.);
    
25.       do j = 1 to ncomb;
    
26.         rc = allcomb(j, &i., of x[*]);
    
27.         zuhe_bcd = substr(strip(x1) || " " || strip(x2) || " " || strip(x3), 1, &i. * 2 - 1);
    
28.         output;
    
29.       end;
    
30.     run;
    
31. 
    
32.     proc sort
    
33.       data = set_bcd_&i.;
    
34.       by zuhe_bcd;
    
35.     run;
    
36.   %end;
    
37. 
    
38.   data set_bcd;
    
39.     set set_bcd_1 - set_bcd_3;
    
40.   run;
    
41. %mend;
    
42. 
    
43. %create_set_bcd;
    
44. 
    
45. data set_e(keep = zuhe_e);
    
46.   array x(2) $2. ('e1' 'e2');
    
47. 
    
48.   n = dim(x);
    
49. 
    
50.   k = 1;
    
51. 
    
52.   ncomb = comb(n, k);
    
53. 
    
54.   do i = 1 to ncomb;
    
55.     rc = allcomb(i, k, of x[*]);
    
56.     zuhe_e = strip(x1);
    
57.         output;
    
58.   end;
    
59. run;
    
60. 
    
61. proc sql;
    
62.   create table all_combinations as
    
63.   select strip(zuhe_a) || " " || strip(zuhe_bcd) || " " || strip(zuhe_e) as zuhe
    
64.   from set_a a, set_bcd b, set_e c
    
65.   ;
    
66. quit;

复制代码

sniperhgy 发表于 2015-12-11 14:12
楼主你好，我写了个代码，但是根据我的理解，只有42中组合，请楼主确认下我说的对不对：
1） a1 a2 a3中 ...

非常感谢   请问   假如有这样一个数据集   我想提取他们的先前预定的各种组合   那么该怎样提取  一种组合生成一个数据集    类似用 keep之类的功能