如何找出每个个体中x值小于2并离2最近的观测

1. data aaa;
   
2. input no x;
   
3. cards;
   
4. 1  1
   
5. 1  1.5
   
6. 1  1.8
   
7. 1  2
   
8. 1  3
   
9. 2  1.5
   
10. 2  2
    
11. 2  2.5
    
12. 3  1.5
    
13. 3  1.9
    
14. 3  2.5
    
15. 3  4
    
16. ;
    
17. run;

复制代码

如何找出每个个体中x值小于2并离2最近的观测，用一个新变量标记。
谢谢！

问题转化为每个ID的x<2的min(2-x)OBS.
proc sql;
create table obs as select min(2-x) as min,no,x from aaa where x<2 group by no having min=(2-x);
quit;

wang1839 发表于 2016-10-31 23:31
问题转化为每个ID的x

数据集aaa中多一个变量，比如用0表示每个个体中x值小于2并离2最近的观测。

proc sql;
create table wanted as
select a.*, minx
from aaa a left join
(select *, abs(x-2) as y, 0 as minx
from aaa
where x<2
group by 1
having y=min(y)) b
on a.no=b.no and a.x=b.x;
quit;

merge回去 by no x，然后min 不为空FLAG=0;ELSE flag=1

我的答案

1. data bbb;
   
2.         set aaa(where=(x<2));
   
3.         by no x;
   
4.         if last.no then flag='Y';
   
5. run;
   
6. 
   
7. data ccc;
   
8.         merge aaa(in=a) bbb(in=b);
   
9.         by no x;
   
10. run;

复制代码

笨方法，供参考！

1. /*原始数据input*/
   
2. data aaa;
   
3. input no x;
   
4. cards;
   
5. 1  1
   
6. 1  1.5
   
7. 1  1.8
   
8. 1  2
   
9. 1  3
   
10. 2  1.5
    
11. 2  2
    
12. 2  2.5
    
13. 3  1.5
    
14. 3  1.9
    
15. 3  2.5
    
16. 3  4
    
17. ;
    
18. run;
    
19. /*对x小于2的记录标记flag=1，计算距离值2-x，大于等于2的记录标记flag=2，距离值默认999*/
    
20. data bbb;
    
21.         set aaa;
    
22.         if x<2 then do;
    
23.                 flag=1;
    
24.                 distance=2-x;
    
25.                 end;
    
26.         else do;
    
27.                 flag=2;
    
28.                 distance=999;
    
29.                 end;
    
30. run;
    
31. /*排序*/
    
32. proc sort data=bbb;
    
33.         by no flag distance;
    
34. run;
    
35. /*赋值最小距离值标签*/
    
36. data ddd;
    
37.         set bbb;
    
38.         by no flag distance;
    
39.         if first.no & first.flag & first.distance then group=0;
    
40.         keep no x group;
    
41. run;
    
42. proc print data=ddd;
    
43. run;

复制代码