请教这个程序应该如何编写

JINGJU学习到了
谢谢

1. data class;
   
2.    input class $ @@;
   
3.    datalines;
   
4.    A A A B B C
   
5.    ;
   
6. run;
   
7. 
   
8. %macro  class;
   
9. proc sql;   
   
10.    select count(distinct class) into :class_n from class;
    
11.    select count(*) into :number from class ;
    
12. quit;
    
13. 
    
14. %let class_n=&class_n;
    
15. %let number=&number;
    
16. 
    
17. proc sql;
    
18.    select distinct class into :class1-:class&class_n  from class;
    
19.    %do i=1 %to &class_n;
    
20.    select count(*) into :&&class&i from class where class="&&class&i";
    
21.    %end;
    
22. quit;
    
23. 
    
24. proc sql;
    
25.     create table test1 as
    
26.         select distinct
    
27.         %do i=1 %to &number-1;
    
28.             a&i..class as a&i,
    
29.         %end;
    
30.         a&number..class as a&number
    
31.         from
    
32.         %do i=1 %to &number-1;
    
33.             class as a&i,
    
34.         %end;
    
35.            class as a&number;
    
36. quit;
    
37. 
    
38. data test2;
    
39.       set test1;
    
40. 
    
41.           %do i=1 %to &class_n;
    
42.       n_&&class&i=0;
    
43.           %end;
    
44. 
    
45.       %do i=1 %to &number;
    
46.             %do j=1 %to &class_n;
    
47.         if a&i="&&class&j" then n_&&class&j=n_&&class&j+1;
    
48.                 %end;
    
49.       %end;
    
50. 
    
51.       %do i=1 %to &class_n;
    
52.             if n_&&class&i=&&&&&&class&i;
    
53.         %end;
    
54. run;
    
55. %mend class;
    
56. 
    
57. %class;

复制代码

harlon1976 发表于 2010-3-12 09:05
感谢楼上兄弟的无私奉献，我要仔细学习，我再问一下，如果还是2类，但观测个数不是五个，例如共有n个，其中A有x个，B有y个，则总共的结果有N=n!/(x!*y!)，那么这个程序中的一些参数应该如何调整，最好能将这些参数放在宏中，能够任意附值。请兄弟再赐教啊。

1. 
   
2. DATA RAW;
   
3. INPUT X  class  $;
   
4. DATALINES;
   
5. 1      A
   
6. 5      B
   
7. 4      A
   
8. 2      B
   
9. 8      A
   
10. 9     A
    
11. ;
    
12. PROC SQL NOPRINT;
    
13. SELECT DISTINCT CLASS, COUNT(*) INTO :C1-:C2, :N1-:N2 FROM RAW GROUP BY CLASS;
    
14. QUIT;
    
15. %LET N = %EVAL(&N1+&N2);
    
16. %LET Ncombinations = %SYSFUNC(COMB(&N, &N1));
    
17. DATA RAW1;
    
18. SET RAW;
    
19. IF _N_ = 1 THEN DO;
    
20.   NA = 0; NB = 0;
    
21. END;
    
22.   NA +(class = "&C1"); NB +(class = "&C2");
    
23.   N = NA*(class = "&C1")+NB*(class = "&C2");
    
24. RUN;
    
25. PROC PLAN;
    
26.       factors combination = &Ncombinations ordered ZZ = &N1 of &N comb;
    
27.       ods output Plan = Combinations;
    
28. RUN;
    
29. DATA  Combinations_1;
    
30. SET Combinations;
    
31. ARRAY XX{&N} $; ARRAY  ZZ{&N1};
    
32. DO I = 1 TO DIM(XX);
    
33.   DO J = 1 TO DIM(ZZ);
    
34.    IF ZZ[J] = I THEN DO;
    
35.     XX[I] = "&C1";
    
36.     LEAVE;
    
37.    END;   
    
38.   END;
    
39.   IF  XX[I] NE "&C1" THEN XX[I] = "&C2";;
    
40.   class = XX[I];
    
41.   NA = SUM(NA, (class = "&C1")); NB = SUM(NB, (class = "&C2"));
    
42.   N = NA*(class = "&C1")+NB*(class = "&C2");
    
43.       OUTPUT;
    
44. END;
    
45. KEEP combination class N;
    
46. RUN;
    
47. PROC SQL;
    
48. CREATE TABLE Total AS
    
49.   SELECT a.combination, a.class, b.X FROM Combinations_1 AS a, RAW1 AS b
    
50.    WHERE a.class = b.class AND a.N = b.N;
    
51. QUIT;
    
52. %MACRO OUTPUTmcr;
    
53. %DO I = 1 %TO &Ncombinations;
    
54.   DATA Dset&I;
    
55.    SET Total (WHERE = (combination = &I));
    
56.   RUN;
    
57. %END;
    
58. %MEND OUTPUTmcr;
    
59. %OUTPUTmcr
    
60. ;
    
61. /*Assuming there is no missing class variable (2 levels) and no replicated records, that is X by Class unique*/

复制代码

12# frackdeng

Hi, I spent too much coding on connecting variable x and class. Why you just ignore X?

sorry. this is a very good code. Thanks, I learn a lot.

我没理解错的话，应该就是把class排序而已，之后只要把X对应上去就可以了。在这里A对应1，4，8；B对应，5，2
14# jingju11

frackdeng 发表于 2010-3-12 13:46
我没理解错的话，应该就是把class排序而已，之后只要把X对应上去就可以了。在这里A对应1，4，8；B对应，5，2
14# jingju11

You are right. But as far as I can image, it is not a work of a couple of lines of coding. Why don't you just have a try.

1. data class;
   
2. input x class $ @@;
   
3. datalines;
   
4. 1 A 4 A 8 A 5 B 2 B
   
5. ;
   
6. run;
   
7. 
   
8. proc sort data=class;
   
9. by class;
   
10. run;
    
11. 
    
12. %macro class;
    
13. proc sql;
    
14. select count(distinct class) into :class_n from class;
    
15. select count(*) into :number from class ;
    
16. quit;
    
17. 
    
18. %let class_n=&class_n;
    
19. %let number=&number;
    
20. 
    
21. proc sql;
    
22. select distinct class into :class1-:class&class_n from class;
    
23. %do i=1 %to &class_n;
    
24. select count(*) into :&&class&i from class where class="&&class&i";
    
25. %end;
    
26. quit;
    
27. 
    
28. proc sql;
    
29. create table test1 as
    
30. select distinct
    
31. %do i=1 %to &number-1;
    
32. a&i..class as a&i,
    
33. %end;
    
34. a&number..class as a&number
    
35. from
    
36. %do i=1 %to &number-1;
    
37. class as a&i,
    
38. %end;
    
39. class as a&number;
    
40. quit;
    
41. 
    
42. data test2;
    
43. set test1;
    
44. %do i=1 %to &class_n;
    
45. n_&&class&i=0;
    
46. %end;
    
47. 
    
48. %do i=1 %to &number;
    
49. %do j=1 %to &class_n;
    
50. if a&i="&&class&j" then n_&&class&j=n_&&class&j+1;
    
51. %end;
    
52. %end;
    
53. 
    
54. %do i=1 %to &class_n;
    
55. if n_&&class&i=&&&&&&class&i;
    
56. %end;
    
57. run;
    
58. 
    
59. data _null_;
    
60. set test2 end=eof;
    
61. if eof then call symput('total',_n_);
    
62. run;
    
63. 
    
64. proc transpose data=test2 out=test3;
    
65. var
    
66. %do i=1 %to &number;
    
67. a&i
    
68. %end;;
    
69. run;
    
70. 
    
71. %do i=1 %to &total;
    
72. data class&i;
    
73. set test3;
    
74. obs=_n_;
    
75. rename col&i=class;
    
76. run;
    
77. 
    
78. proc sort data=class&i;
    
79. by class;
    
80. run;
    
81. 
    
82. data class&i;
    
83. merge class class&i;
    
84. by class;
    
85. run;
    
86. 
    
87. proc sort data=class&i out=class_&i(keep=x class);
    
88. by obs;
    
89. run;
    
90. %end;
    
91. 
    
92. %mend class;
    
93. 
    
94. %class;

复制代码

16# jingju11

jingju11 的程序可行，只是我的数据是16个，是7+9类型的，运算的时间比较长，不知是否还有更简洁的程序，同时也谢谢frackdeng 的程序。

还有个问题请教，上面程序生成了10个数据集dest1-dest10，如何编写这样一个程序：随机生成1-10中的任意一个整数，比如是2则就调用数据集dest2。谢谢！

harlon1976 发表于 2010-3-12 19:42
还有个问题请教，上面程序生成了10个数据集dest1-dest10，如何编写这样一个程序：随机生成1-10中的任意一个整数，比如是2则就调用数据集dest2。谢谢！

调用是什么意思？