如何用sas实现排列组合

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

从0-9这10个数字中选取任意的4个数字组成一个四位数，问这些四位数是哪些？哪位大侠帮忙写下程序，谢谢

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：如何用 如何

Example
The following is an example of the ALLCOMB function.
data _null_;
   array x[5] $3 ('ant' 'bee' 'cat' 'dog' 'ewe');
   n=dim(x);
   k=3;
   ncomb=comb(n,k);
   do j=1 to ncomb+1;
      rc=allcomb(j, k, of x[*]);
      put j 5. +3 x1-x3 +3 rc=;
   end;
run;
SAS writes the following output to the log:
    1   ant bee cat    rc=0
    2   ant bee ewe    rc=3
    3   ant bee dog    rc=3
    4   ant cat dog    rc=2
    5   ant cat ewe    rc=3
    6   ant dog ewe    rc=2
    7   bee dog ewe    rc=1
    8   bee dog cat    rc=3
    9   bee ewe cat    rc=2
   10   dog ewe cat    rc=1
   11   dog ewe cat    rc=-1

1. data test(keep=aa);
   
2.         array x[10] $1 ('0' '1' '2' '3' '4' '5' '6' '7' '8' '9');
   
3.         n=dim(x);
   
4.         k=4;
   
5.         ncomb=comb(n,k);
   
6. 
   
7.         do j=1 to ncomb;
   
8.                 rc=allcomb(j, k, of x[*]);
   
9. 
   
10.                 if x1>0 then
    
11.                         do;
    
12.                                 aa=input(cats(x1,x2,x3,x4),4.);
    
13.                                 output;
    
14.                         end;
    
15.         end;
    
16. run;

复制代码

上面写得不对，改一下

1. %let s=4;
   
2. 
   
3. data test(keep=aa);
   
4.         array x[10] $1 ('0' '1' '2' '3' '4' '5' '6' '7' '8' '9');
   
5.         array y[&s] $1;
   
6.         n=dim(x);
   
7.         k=&s;
   
8.         ncomb=comb(n,k);
   
9. 
   
10.         do i=1 to ncomb;
    
11.                 rc=allcomb(i, k, of x[*]);
    
12. 
    
13.                 do j=1 to k;
    
14.                         y[j]=x[j];
    
15.                 end;
    
16. 
    
17.                 nfact=fact(k);
    
18. 
    
19.                 do j=1 to nfact;
    
20.                         call allperm(j, of y[*]);
    
21. 
    
22.                         if (y1>0) then
    
23.                                 do;
    
24.                                         aa=input(cats(of y1-y&s),10.);
    
25.                                         output;
    
26.                                 end;
    
27.                 end;
    
28.         end;
    
29. run;
    
30. 
    
31. proc sort data=test;
    
32.         by aa;
    
33. run;

复制代码

1. data a(keep=value);
   
2. do i=1 to 9;
   
3.   do j=0 to 9;
   
4.    do k=0 to 9;
   
5.     do l=0 to 9;
   
6.          if i^=j and i^=k and i^=l and j^=k and j^=l and k^=l then do;
   
7.             value=i*1000+j*100+k*10+l;
   
8.             output;
   
9.          end;
   
10.     end;
    
11.    end;
    
12.   end;
    
13. end;
    
14. run;

复制代码

/* note: 10选4的组合应该有210种, ncomb(10,4) = 210 */

下面是不用 comb(), allcomb() functions 而生成 n 选４的全部组合的程序，其中n可以是大于４的任意数。

1. data comb_n_4;
   
2.   n=10; r=4;  
   
3.   do k1 = 0 to (n-1) - (r-1);              
   
4.    do k2 = k1+1 to (n-1) - (r-2);         
   
5.      do k3 = k2+1 to (n-1) - (r-3);      
   
6.        do k4 = k3+1 to (n-1) - (r-4);         
   
7.        index=catx(',', k1, k2, k3, k4);
   
8.    put index=;
   
9.    output;  
   
10.   end; end; end; end;
    
11. run;

复制代码

mingfeng07 发表于 2014-12-20 19:17

要是可以重复取值
data a(keep=value);
do i=1000 to 9999;
            value=i;
            output;
         end;
run;

%macro test1();
data test2(keep=aa);
    array x[5] $1 ('age','wealth','sex','x1','x2');
    n=dim(x);
        %do k=1 %to 3;
    ncomb=comb(n,&k.);
       
    do j=1 to ncomb;
            
          rc=allcomb(j, &k., of x[*]);
          if x1>0 then
           do;
                aa=putn(cats(of x1-x&k.),10.);
                output;
            end;
     end;
         %end;
run;

%mend;
%test1();
连接字符有问题~~~~~~~~~~~~~~·