筛选观测

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

   a                     b             c
173.2465        5.641        28375
173.2908        5.645        7161
173.1294        5.642        767098
197.2151        9.998        51203
197.2644        9.999        13350
197.0907        9.999        1193382
197.3505        10.001        6035
173.8311        6.39        6832
197.899        10.34        7456

上述数据中，任意两行观测的变量a和b相减，若-0.5=<变量a的差值=<0.5, -0.3=<变量b的差值=<0.3，则将这些观测输出。
谢谢各位了！
结果是：

   a                                b                          c
173.2465           5.641           28375
173.2908           5.645           7161
173.1294           5.642          767098
197.2151          9.998            51203
197.2644          9.999            13350
197.0907          9.999            1193382
197.3505         10.001           6035

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：筛选观测

样本数据有问题吧？
只是小数部分相减的话，才可以得到期望的结果。
？？？？

1. data a;
   
2. input a b c;
   
3. cards;
   
4. 173.2465        5.641        28375
   
5. 173.2908        5.645        7161
   
6. 173.1294        5.642        767098
   
7. 197.2151        9.998        51203
   
8. 197.2644        9.999        13350
   
9. 197.0907        9.999        1193382
   
10. 197.3505        10.001        6035
    
11. 173.8311        6.39        6832
    
12. 197.899        10.34        7456
    
13. ;
    
14. run;
    
15. data a;
    
16. set a;
    
17. n=_n_;
    
18. run;
    
19. data b;
    
20. set a;
    
21. rename n=n1 a=a1 b=b1 c=c1;
    
22. run;
    
23. 
    
24. proc sql;
    
25. create table merge as
    
26. select a.*,b.* from a,b
    
27. where abs(a.a-b.a1)<=0.5 and (a.b-b.b1)<=0.3 and n^=n1;
    
28. quit;
    
29. 
    
30. proc sort data=merge(keep=n) nodupkey out=c;
    
31. by n;
    
32. run;
    
33. data final;
    
34. merge c(in=x) a;
    
35. by n;
    
36. if x;
    
37. drop n;
    
38. run;

复制代码

完美解决

   将上一帖结果中的重复值和 a, b 差为零的观测去掉，就是想要的结果。

proc sql noprint;
   select count(*) into :Nobs    from aaa;
   select a into :a separated by ' ' from aaa;
   select b into :b separated by ' ' from aaa;
quit;

data bbb;
   array aa{&nobs} _temporary_ (&a);
   array bb{&nobs} _temporary_ (&b);
   set aaa end=last;
   n=_N_;
     do i = 1  to &nobs;
       a_dif = abs(a - aa(i));
       b_dif = abs(b - bb(i));
       if  a_dif <= 0.5 or b_dif <= 0.3  then output;
     end;
run;

proc sort data=bbb out=ccc (keep=a b c) nodupkey;
   by n a b;
   where a_dif^=0 and b_dif^=0;
run;