求用HASH做闭区间求和

求助一个问题，总目标是做分类合计的，但是他要求的分类全是闭区间，就有点麻烦了，想了一会感觉用HASH表做比较好，但是本人HASH水平也能做个inner join，故还请求助。。。

数据如下：（subjid是观测的ID，在这里只列出了 ‘第一个人’ 的情况）
data a;
input value x subjid;
cards;
1 1 001
2 1 001
3 2 001
4 2 001
5 2 001
6 3 001
7 3 001
8 5 001
;
run;

想做成的结果：

Subjid	X	SUM	说明
001	[1,2]	15	1+2+3+4+5
001	[2,3]	25	3+…+7
001	[3,4]	13	6+7+0
001	[4,5]	8	0+8
001	[1,5]	36	1+…+8，觉得麻烦就不做哈
002	[1,2]	…
…	…	…

请朋友们帮助哈，谢谢！！

更新：
大概是这个意思，
data b;
if 0 then set a;
if _n_=1 then do;
  dcl hash h(dataset:'a');
  h.definekey('x');
  h.definedata(all:'Y');
  h.definedone();
end;

do i = 1 to 4;
  j = i + 1;

  if not h.find(key:i) or not h.find(key:j) then sum+value;
  output;
end;
run;

不过值不对，而且也没考虑到subjid。。

在线等呵，ths~~

Tigflanker 发表于 2013-12-20 12:54
不好意思，我没选你的答案，因为我想通过这个学习HASH，这个我已经用数组（count1 到 count5） + Transpo ...

SAS HASH Object 还有一个名字叫 associativearray，　所以原则上说用array可以做的事情，用hash object 应该也可以做到。

下面程序是此问题用 hash object 的做法, 其中包括了用hash做分组求和，填充缺失值，和指定闭区间求和。

1. data a;
   
2. input value x subjid $;
   
3. cards;
   
4. 1 1 001
   
5. 2 1 001
   
6. 3 2 001
   
7. 4 2 001
   
8. 5 2 001
   
9. 6 3 001
   
10. 7 3 001
    
11. 8 5 001
    
12. 1 1 002
    
13. 2 1 002
    
14. 3 2 002
    
15. 4 2 002
    
16. 5 2 002
    
17. 6 3 002
    
18. 7 3 002
    
19. 8 5 002
    
20. ; run;
    
21. 
    
22. data b(keep = subjid xx sum_x);
    
23. length subjid xx $3.;
    
24. 
    
25.   declare hash h(ordered: 'Y');
    
26.   h.definekey('subjid', 'x');
    
27.   h.definedata('subjid', 'x', 'vsum');
    
28.   h.definedone();
    
29.   declare hiter hi ('h');
    
30. 
    
31.     do until (eof);        /* load the sum of values for each subjid_x */
    
32.       set a end = eof;
    
33.          rc = h.check();
    
34.          if rc ne 0 then       do;
    
35.             vsum = value;
    
36.             rc1 = h.add();    end;
    
37.     else if rc eq 0 then       do;
    
38.             vsum ++ value;
    
39.           rc1 = h.replace();  end;
    
40.     end;
    
41. 
    
42.    do subjid = '001', '002';     /* fill the missing x */
    
43.        do x = 1 to 5;
    
44.           rc = h.find();
    
45.           if rc ne 0 then    do;
    
46.           vsum = 0;
    
47.           rc = h.replace(); end;
    
48.        end;
    
49.    end;
    
50.   
    
51.     do subjid = "001", "002";   /* calculating expected sum */
    
52.        do i = 1 to 4;
    
53.             x = i;
    
54.            rc = h.find();
    
55.         sum_x = vsum;
    
56.             x = x + 1;
    
57.            rc = h.find();
    
58.         sum_x = sum_x + vsum;
    
59.            xx = catx('_', i, x);
    
60.            output;
    
61.        end;
    
62.          sum_x = 0;
    
63.        do x = 1 to 5;
    
64.          rc = h.find();
    
65.          sum_x ++ vsum;
    
66.        end;
    
67.           xx = catx('_', '1', '5');
    
68.           output;
    
69.     end;  
    
70. run;

复制代码

考虑到subjid是什么意思？subjid也是有不同值？

1. /* 是不是这样, 没有用hash */
   
2. data a;
   
3. input value x id;
   
4. cards;
   
5. 1 1 001
   
6. 2 1 001
   
7. 3 2 001
   
8. 4 2 001
   
9. 5 2 001
   
10. 6 3 001
    
11. 7 3 001
    
12. 8 5 001
    
13. 1 1 002
    
14. 2 1 002
    
15. 3 2 002
    
16. 4 2 002
    
17. 5 2 002
    
18. 6 3 002
    
19. 7 3 002
    
20. 8 5 002
    
21. ; run;
    
22. 
    
23. proc sql;
    
24. create table aaa as
    
25.     select distinct id, x, sum(value) as svalue
    
26.     from a
    
27.     group by id, x
    
28.     order by id, x;
    
29. run;
    
30. 
    
31. data bbb;
    
32.     set aaa;
    
33.        by id;
    
34.        value_lag = lag(svalue);
    
35.        x_lag = lag(x);
    
36.        retain first_x first_value total_x;
    
37.     if first.id then do;
    
38.        first_value = svalue;
    
39.        first_x = x;  Total_x = svalue;
    
40.        end;
    
41.     else do;
    
42.            x_x = catx('_',  put(x_lag, 1.), put(x, 1.) );
    
43.          sum_x = svalue + value_lag;
    
44.        total_x = Total_x + svalue;
    
45.        output;
    
46.     end;
    
47.     if last.id then do;
    
48.        x_x = catx('_', put(first_x, 1.), put(x, 1.));
    
49.        sum_x = total_x;
    
50.        output;
    
51.     end;
    
52.    keep id x_x sum_x;
    
53. run;

复制代码

yongyitian 发表于 2013-12-20 10:46

不好意思，我没选你的答案，因为我想通过这个学习HASH，这个我已经用数组（count1 到 count5） + Transpose做好了

你写的我只能给你加这么多，不好意思，谢谢啊。