救急啊，谢谢各位

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

数据结构如下
data test;
input x y z@@;
cards;
5.88  173.4723  5094
5.88  173.2468  51751
5.881  173.1299  1203539
5.883  173.2951  13767
5.883  173.4073  6816
5.884  173.7063  3000
5.887  173.9003  4000
5.882  123.100   10000
5.933  173.1298   20000
run;
说明：要求先按变量X和Y排序，把任意两行变量X的差值落在0到0.003之间，同时Y差值落在范围0到0.3425之间的行输出，在满足这些条件的行中用一个新行输出，新行中Z取所有满足条件行中Z的最大值,X为最大Z值对应的值，Y为满足条件行的所有Y值的平均值。
例如
5.88 173.2468   51751
5.881 173.1299  1203539
这两行X的差值为0.001，Y的差值为0.1169，同时满足两个条件，这两行中Z的最大值是1203539.故X取值为5,881，输出新行为
X           Y                            Z
5.881    173.18835  1203539
输出数据集结果应为：
5.881   173.13      1203539
5.887   173.8033  4000
5.882   123.100   10000
5.933   173.1298  20000
这个问题纠结很久了，请各位大侠指教！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Input cards Data test card 平均值 最大值

无解啊

1. %let n=9; /*定义你原数据有多少观测*/
   
2. %let num=%eval(&n*(&n-1)/2); /*排列组合*/
   
3. %put #
   
4. 
   
5. /*生成排列组合数据集*/
   
6. data sequence;
   
7. informat c $8.;
   
8. do i=1 to &n-1;
   
9. do j=1 to &n-i;
   
10. a=i;b=i+j;
    
11. c=catx(",",a,b);
    
12. output;
    
13. end;
    
14. end;
    
15. drop i j;
    
16. run;
    
17. 
    
18. 
    
19. proc sql;
    
20. select c into:c1-:c%left(&num) from sequence;
    
21. quit; /*将排列组合的观测序号赋给多个宏变量 */
    
22. 
    
23. data test;
    
24. input x y z@@;
    
25. cards;
    
26. 5.88  173.4723  5094
    
27. 5.88  173.2468  51751
    
28. 5.881  173.1299  1203539
    
29. 5.883  173.2951  13767
    
30. 5.883  173.4073  6816
    
31. 5.884  173.7063  3000
    
32. 5.887  173.9003  4000
    
33. 5.882  123.100   10000
    
34. 5.933  173.1298   20000
    
35. ;
    
36. run;
    
37. 
    
38. proc sort data=test;
    
39. by z;
    
40. run;
    
41. 
    
42. %macro xunhuan;
    
43. proc delete data=null; run;
    
44. %do j=1 %to #
    
45. data tem;
    
46. do n=&&c&j;
    
47. set test point=n;
    
48. output;
    
49. end;
    
50. stop;
    
51. run;
    
52. data tem1;
    
53. set tem;
    
54. x1=x-lag(x);
    
55. y1=y-lag(y);
    
56. x1=round(x1,0.001);
    
57. y1=round(y1,0.0001);
    
58. drop x1 y1;
    
59. y=(lag(Y)+y)/2;
    
60. if x1^=. and abs(x1)<0.003 and y1^=. and abs(y1)<0.3425 then output;
    
61. run;
    
62. proc append base=null data=tem1 force;
    
63. run;
    
64. %end;
    
65. %mend;
    
66. %xunhuan;

复制代码

希望大家多多指教，谢谢