这个怎么用sas实现？很实际的一个问题。。

数据如下：含有两个字段，第一个是手机号码，第二个是优惠代码
13900000000 10102345
13900000000 10102346
13900000000 10102347
13900000000 10102348
13900000001 10102345
13900000001 10102346
13900000001 10102349
13900000002 10102345
13900000002 10102347
13900000002 10102349
13900000003 10102345
13900000003 10102346
13900000003 10102348
13900000003 10102349
需求如下：算出同时使用两个优惠组合的客户数

prom_1	prom_2	客户数

如：10102345和10102346同时使用的客户数是3个。


sql解法很精彩，现在求data步实现方法，data步悬赏500币！

1. data a;
   
2.         input nbr code;
   
3.         datalines;
   
4.         13900000000 10102345
   
5. 13900000000 10102346
   
6. 13900000000 10102347
   
7. 13900000000 10102348
   
8. 13900000001 10102345
   
9. 13900000001 10102346
   
10. 13900000001 10102349
    
11. 13900000002 10102345
    
12. 13900000002 10102347
    
13. 13900000002 10102349
    
14. 13900000003 10102345
    
15. 13900000003 10102346
    
16. 13900000003 10102348
    
17. 13900000003 10102349
    
18. ;
    
19. 
    
20. run;
    
21. 
    
22. proc sql;
    
23.         create table b as
    
24.                 select
    
25.                         a.code as prom_1
    
26.                         ,b.code as prom_2
    
27.                         ,count(distinct a.nbr) as Csr_Cnt
    
28.                 from a a
    
29.                 inner join a b
    
30.                 on a.code<b.code
    
31.                         and a.nbr=b.nbr
    
32.                 group by
    
33.                 a.code
    
34.                 ,b.code
    
35.         ;
    
36. quit;

复制代码

用SAS中的描述性统计工具可以看到吧

-                                                             -

这回复速度也太快了吧？不知道谁的时钟不准了

{:3_42:}

{:3_53:}

500，这么快就被秒杀了

data raw;

  input nbr code;

  datalines;

13900000000 10102345

13900000000 10102346

13900000000 10102347

13900000000 10102348

13900000001 10102345

13900000001 10102346

13900000001 10102349

13900000002 10102345

13900000002 10102347

13900000002 10102349

13900000003 10102345

13900000003 10102346

13900000003 10102348

13900000003 10102349

;
run;
proc sort data=raw nodup;
  by code nbr;
run;
proc summary data=raw nway;
  var nbr;
  class code;
  output out=nbr_cnt(drop=_:) n=cnt;
run;
data cnt_2;
  set nbr_cnt;
  where cnt=2;
run;