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

假设数据集中有变量两个，分别为数值型的x和字符型的class。假设有如下5个观测：
x     class
1      A
5      B
4      A
2      B
8      A
在这个数据集中，变量class有两个取值，分别为A、B，如果现在考察对变量class进行重新排序，根据概率论知识：这里有5条观测，分成两类，一类（A类）有3个，第二类（B类）有2个，则可以重新组合的结果为10种，例如下面的结果就是其中的两种：
x         class                x           class
1            A                   5              B
4            A                   2              B
8            A                   8              A
5            B                  1               A
2            B                  4               A
现在的问题是：如何根据原始的数据集，编写SAS程序生成所有的10种结果。
如果现在不只是A、B两类，比如3类，类似的程序又如何编写呢？请高手给予解决。

我的理解是10=5!/(3!*2!)

即只与A、B的顺序有关，和x的排列没关系，即A间无差别

LZ是不是这个意思

这里不考虑x的影响，当然x随着class的变动而同步变动，你的理解公式完全正确，请给出这个程序，好象你以前就帮我解决个问题，希望你这次还能帮帮我，谢谢！

data a;
input x class $@;
cards;
1      A
5      B
4      A
2      B
8      A
;
run;
proc sort;
by class x;
run;
proc append base=a data=a;
run;
data b;
set a;
class1=lag1(class);
class2=lag2(class);
class3=lag3(class);
class4=lag4(class);
class5=lag5(class);
xn=class;
x1=lag1(x);
x2=lag2(x);
x3=lag3(x);
x4=lag4(x);
x5=lag5(x);
run;
data c(drop=class5 xn x5);
set b;
where class5 ^=' ' and x5^=.;
run;
proc print;
run;

结果出来了，但是组合成数据集，还要努力，时间有限，期待下面的人写完

最好能写成一个宏文件出来
抛砖引玉。。。。

貌似不对吧......

x1-class1
x2-class2
x3-class3
x4-class4
x-class
可以优化一下
呵呵

欢迎高手继续努力，我想大家主要是学知识，钱多少无所谓的，6楼还在考虑吧，谢谢你啊

harlon1976 发表于 2010-3-11 17:22
欢迎高手继续努力，我想大家主要是学知识，钱多少无所谓的，6楼还在考虑吧，谢谢你啊

1. 
   
2. title 'All Combinations of (2 Choose 5) Integers';
   
3. PROC PLAN;
   
4.       factors combination=10 ordered
   
5.              group = 2 of 5 comb;
   
6.       ods output Plan = Combinations;
   
7. RUN;
   
8. DATA  Combinations_1;
   
9. SET Combinations;
   
10. ARRAY XX{5} $;
    
11. DO I = 1 TO 5;
    
12.   IF group1 = I THEN XX[I] = 'B';
    
13.    ELSE IF group2 = I THEN XX[I] = 'B';
    
14.     ELSE XX[I] = 'A';
    
15.   class = XX[I];
    
16.   NA = SUM(NA, (class = 'A')); NB = SUM(NB, (class = 'B'));
    
17.   N = NA*(class = 'A')+NB*(class = 'B');
    
18.       OUTPUT;
    
19. END;
    
20. KEEP combination class N;
    
21. RUN;
    
22. DATA RAW;
    
23. INPUT X  class  $;
    
24. IF _N_ = 1 THEN DO;
    
25.   NA = 0; NB = 0;
    
26. END;
    
27.   NA +(class = 'A'); NB +(class = 'B');
    
28.   N = NA*(class = 'A')+NB*(class = 'B');
    
29. DATALINES;
    
30. 1      A
    
31. 5      B
    
32. 4      A
    
33. 2      B
    
34. 8      A
    
35. ;
    
36. PROC SQL;
    
37. CREATE TABLE Total AS
    
38.   SELECT a.combination, a.class, b.X FROM Combinations_1 AS a, RAW AS b
    
39.    WHERE a.class = b.class AND a.N = b.N;
    
40. QUIT;
    
41. %MACRO OUTPUTmcr;
    
42. %DO I = 1 %TO 10;
    
43.   DATA D&I;
    
44.    SET Total (WHERE = (combination = &I));
    
45.   RUN;
    
46. %END;
    
47. %MEND OUTPUTmcr;
    
48. %OUTPUTmcr
    
49. ;

复制代码

I think here is binomial and otherwise multinomial if class has more than 2 levels.
If more than 2 levels for the class, probably I will do the combination each and then combine them together by Cartesian method provided by proc SQL.
JingJu

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