解决实际问题之填补确实月份末日数据

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

1. data actual;
   
2. informat ytime yymmdd10. stime yymmdd10.;
   
3. input ytime stime yvalue svalue;
   
4. format ytime yymmdd10.  stime yymmdd10.;
   
5. month=month(stime);
   
6. cards;
   
7. 2016-01-28         2016-05-13          10000        1000
   
8. 2016-01-28         2016-06-21          9000        1000
   
9. 2016-01-28         2016-08-03          8000        5000
   
10. 2016-01-28         2016-08-20          3000        3000
    
11. ;
    
12. 
    
13. /*生成一月份的末日*/
    
14. data first_mori;
    
15. informat ytime yymmdd10. etime yymmdd10.;
    
16. input ytime etime;
    
17. format ytime yymmdd10.  etime yymmdd10.;
    
18. cards;
    
19. 2016-01-28         2016-01-31
    
20. ;
    
21. run;
    
22. 
    
23. /*利用数组生成二月到十月的末日*/
    
24. data mori(keep=ytime etime);
    
25. ytime='28jan2016'd;
    
26. ytime1='2016-01-28';/*假如提供的是字符型日期*/
    
27. mori0=intnx('month',input(ytime1, anydtdte.),0,'e');
    
28. array mori(*) mori1-mori11;
    
29. do i=1 to dim(mori);/*数组下标不能从0开始*/
    
30. mori(i)=intnx('month',input(ytime1, anydtdte.),i,'e');/*间隔i个月最后一天*/
    
31. etime=mori(i);
    
32. format ytime mori0 mori1-mori11 etime yymmdd10.;
    
33. output;
    
34. end;
    
35. run;
    
36. 
    
37. /*生成所有月份末日*/
    
38. data allmori;
    
39. set first_mori mori;
    
40. month=month(etime);
    
41. run;
    
42. 
    
43. /*生成需要填补的每月末日*/
    
44. proc sql;
    
45. create table supplement as
    
46. select ytime,etime
    
47. from allmori where not exists
    
48. (select *
    
49. from actual
    
50. where allmori.month=
    
51. actual.month);
    
52. run;
    
53. 
    
54. /*月份相同的交集查询*/
    
55. proc sql;
    
56. create table month_intersect as
    
57. select *
    
58. from actual as a,
    
59. allmori as b
    
60. where a.month=b.month;
    
61. quit;
    
62. 
    
63. /*重命名为了后面的合并变量*/
    
64. data supplement_1;
    
65. set supplement(rename=(etime=stime));
    
66. svalue=0;
    
67. run;
    
68. 
    
69. /*set语句纵向拼接数据集*/
    
70. data youneed(drop=month);
    
71. set actual supplement_1;
    
72. run;
    
73. proc sort data=youneed;by descending stime;run;
    
74. 
    
75. /*retain函数复制上一条观测值，因为是之前是降序排序，相当于取下一条观测*/
    
76. data youneed1(drop=yvalue);
    
77. retain ytime stime yvalue1 svalue;
    
78. set youneed;
    
79. if yvalue^=. then yvalue1=yvalue;
    
80. run;
    
81. 
    
82. /*重命名换原你的变量名*/
    
83. data youneed1;
    
84. set youneed1(rename=(yvalue=yvalue1));
    
85. run;
    
86. /*最后升序排序一下，大功告成*/
    
87. proc sort data=youneed1;by stime;run;

复制代码

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：实际问题 日数据 Supplement intersect informat actual 9000