求高手帮助

计算重合

sorry, i would like to help you i couldn't.

给你个sql的比较傻的方法，好的方法我自己留着，呵呵。
proc sql noprint;
        create table table_n as
                select time,id,code,count(distinct sd) as n from a group by time,id,code;
        create table table_b as
                select t1.time,t1.id,t1.code as code_s,t2.code as code_e
                from table_n t1,table_n t2
                where t1.time=t2.time and t1.id=t2.id and t1.code<t2.code;
        create table table_c as
                select t1.time,t1.id,t1.code_s,t1.code_e,t2.sd
                from table_b t1 left join a t2
                on t1.time=t2.time and t1.id=t2.id and (t2.code=t1.code_s or t2.code=t1.code_e);
        create table table_m as
                select time,id,code_s,code_e,count(distinct sd) as m label='总个数' from table_c group by time,id,code_s,code_e;
        create table table_m as
                select t1.*,t2.n as n1
                from table_m t1 left join table_n t2
                on t1.time=t2.time and t1.id=t2.id and t1.code_s=t2.code;       
        create table table_m as
                select t1.*,t2.n as n2,(t2.n+t1.n1-t1.m)/t1.m as chl label='重合率'
                from table_m t1 left join table_n t2
                on t1.time=t2.time and t1.id=t2.id and t1.code_e=t2.code;       
        create table table_clp as
                select time,id,code_s,code_e,avg(chl) as clp
                from table_m
                group by time,id,code_s,code_e;
quit;

learsaas 发表于 2015-5-6 12:41
给你个sql的比较傻的方法，好的方法我自己留着，呵呵。
proc sql noprint;
        create table table_n as

不好意思，id,code  可能都是很多的 ，对于每个id，要计算出code所有两两组合的重合率，然后除以这些组合数目球的对于每个id ，在某个时期的平均重合率

在线等

1. /* this is partial code that count the number of repeated sd */
   
2. /* you need to find the number of distinct sd in each time-id group
   
3.    then calculate the wanted ratio */
   
4. 
   
5. data have;
   
6. input time id code  sd;
   
7. datalines;
   
8. 199806        1        1        1
   
9. 199806        1        1        2
   
10. 199806        1        1        4
    
11. 199806        1        2        4
    
12. 199806        1        2        2
    
13. 199806        1        2        5
    
14. 199806        1        3        1
    
15. 199806        1        3        8
    
16. 199806        2        1        5
    
17. 199806        2        1        4
    
18. 199806        2        1        7
    
19. 199806        2        1        2
    
20. 199806        2        2        5
    
21. 199806        2        2        1
    
22. 199806        2        3        4
    
23. 199806        2        3        7
    
24. 199806        2        3        1
    
25. 199806        2        3        5
    
26. 199812        1        3        8
    
27. 199812        1        1        3
    
28. 199812        1        1        7
    
29. 199812        1        1        4
    
30. 199812        1        1        5
    
31. 199812        1        1        1
    
32. 199812        1        2        7
    
33. 199812        1        2        8
    
34. 199812        1        2        9
    
35. 199812        1        3        1
    
36. 199812        1        3        5
    
37. ; run;
    
38. 
    
39. 
    
40. proc sql noprint;
    
41.     select max(maxcode), max(maxsd) into: max_code,
    
42.                                         : max_sd
    
43.      from (select max(code) as maxcode, max(sd) as maxsd
    
44.            from have
    
45.            group by time, id);
    
46. quit;
    
47. %put &max_code;
    
48. %put &max_sd;
    
49. 
    
50. data want;
    
51.    array code_sd{&max_code, &max_sd} _temporary_ ;
    
52.    array c{&max_code} _temporary_ ;
    
53. 
    
54.    do i = 1 to &max_code;
    
55.      do j = 1 to &max_sd;
    
56.        code_sd[i,j]=0;
    
57.    end; end;
    
58. 
    
59.    do until (last.id);
    
60.      set have;
    
61.      by time id;
    
62.      code_sd[code, sd] = sd;
    
63.    end;
    
64. 
    
65.     do i = 1 to dim(c);
    
66.      c[i]= 0; end;
    
67.     do i = 1 to &max_code - 1;      
    
68.      do j = i+1 to &max_code;
    
69.       do k = 1 to &max_sd;
    
70.        if code_sd[i,k]=code_sd[j,k] and code_sd[i,k] > 0 then c[i]+1;
    
71.     end; end; end;
    
72. 
    
73.     total_count = 0;
    
74.     do i = 1 to dim(c);
    
75.       total_count = total_count + c[i];  
    
76.     end;
    
77.     output; keep time id total_count;
    
78. run;

复制代码