请问如何做笛卡尔乘积的dummy程序

例如我有数据集a：

1. data a;
   
2. input a s d;
   
3. cards;
   
4. 1 1 1
   
5. 1 2 2
   
6. 1 3 1
   
7. 2 4 2
   
8. 2 2 3
   
9. ;run;

复制代码

我想得到b数据集，包括24条，是a数据集中 a变量的所有情况 * s变量的所有情况 * d变量的所有情况（另外我需要考虑11组之多，希望方法能通用）

请问应该如何来做呢？

另外hard code不行喔，也就是说用do循环是不行的。

谢谢。

用cross join 即可，可探索效率更高的算法

1. proc sql noprint;
   
2.    create table B as
   
3.       select X.*,Y.*,Z.*
   
4.          from (select distinct a from a) as X
   
5.                 cross join
   
6.                (select distinct s from a) as Y
   
7.                 cross join
   
8.                (select distinct d from a) as Z
   
9.         order by a,s,d
   
10.     ;
    
11. quit;

复制代码

楼主你好，请试试看我的代码：

1. data a;
   
2. input a s d;
   
3. cards;
   
4. 1 1 1
   
5. 1 2 2
   
6. 1 3 1
   
7. 2 4 2
   
8. 2 2 3
   
9. ;run;
   
10. 
    
11. %macro CreateDescartes(in_dt, out_dt);
    
12.   proc contents
    
13.     noprint
    
14.     data = &in_dt.
    
15.     out = &in_dt._summary(keep = NAME VARNUM);
    
16.   run;
    
17. 
    
18.   proc sort
    
19.     data = &in_dt._summary;
    
20.     by VARNUM;
    
21.   run;
    
22. 
    
23.   data _null_;
    
24.     set &in_dt._summary end = eof;
    
25.     call symput(strip("var") || strip(_N_), strip(NAME));
    
26.         if eof then
    
27.       call symput("var_count", strip(_N_));
    
28.   run;
    
29. 
    
30.   proc sql;
    
31.     create table &out_dt. as
    
32.     select
    
33.       %do i = 1 %to &var_count.;
    
34.         &&var&i.
    
35.         %if &i. lt &var_count. %then %do;
    
36.          ,
    
37.         %end;
    
38.       %end;
    
39.     from
    
40.       %do i = 1 %to &var_count.;
    
41.         (select distinct &&var&i.
    
42.          from &in_dt.)
    
43.         %if &i. lt &var_count. %then %do;
    
44.          ,
    
45.         %end;
    
46.       %end;
    
47.     ;
    
48.   quit;
    
49. %mend;
    
50. 
    
51. %CreateDescartes(a, b);

复制代码

data temp;
     set a;
     dummy=1;
run;
proc means data=temp  completetypes nway noprint;
     class a s d;
     var dummy;
     output out=want(keep=a s d) n=n;
run;

1. proc sql;
   
2.    create table asd as
   
3.    select a.a, b.s, c.d
   
4.    from a(keep=a) as a,
   
5.         a(keep=s) as b,
   
6.         a(keep=d) as c;
   
7. quit;
   
8. proc sort data=asd out=asd_sort nodup;
   
9.      by a s d;
   
10. run;

复制代码

thanks ,  learned