请教一下，通过retain变量累加，如何设置累加次数

lovexialulu 发表于 2017-3-9 15:54
这种应该是你要的吧，组内累积衰减50%

还不是，可能我之前没说清楚，有点不好描述，又上传了图片

还有，想问一下，logistic回归模型是不是主要用于横截面数据，可不可以用在时间序列数据上

1. data zip;
   
2. input zip$ q_1 q_2 q_3;
   
3. cards;
   
4. 50003 0 0 0
   
5. 50003 0 1 0
   
6. 50003 1 1 0
   
7. 50003 1 1 0
   
8. 50003 0 1 0
   
9. 50003 0 0 0
   
10. 50003 1 0 0
    
11. 50003 0 1 0
    
12. 50003 0 1 0
    
13. 50003 1 0 0
    
14. 50003 0 0 0
    
15. 50003 2 1 0
    
16. 50003 0 0 0
    
17. 50003 1 0 0
    
18. 50010 0 0.8 4
    
19. 50010 1 1.8 0
    
20. 50010 0 2.8 0
    
21. 50010 1 2 4
    
22. 50010 0 3 4
    
23. 50010 4 3 6
    
24. 50010 0 2.8 0
    
25. 50010 1 2 4
    
26. 50010 0 3 4
    
27. 50010 3 3 6
    
28. 50010 0 2.8 0
    
29. 50010 1 2 4
    
30. 50010 0 3 4
    
31. 50010 2 3 6
    
32. 50010 0 2.8 0
    
33. 50011 0 3 6
    
34. 50011 2 2.8 0
    
35. 50011 0 2 4
    
36. 50011 0 3 4
    
37. 50011 1 3 6
    
38. 50010 1 2.8 0
    
39. ;
    
40. proc sort data=zip out=b_;by zip;run;
    
41. data b;
    
42. set b_;
    
43. order=_n_;
    
44. run;
    
45. 
    
46. %macro jj(indata,outdata,case,od,out,ct);
    
47. data &outdata.;
    
48. set &indata.(where=(&case.));
    
49. by zip;
    
50. 
    
51. retain num&od. ;
    
52. if q_1 ne 0  then  num&od.+1;
    
53. if first.zip then do;
    
54.    if q_1 ne 0 then num&od.=1;else
    
55.    num&od.=0;
    
56. end;
    
57. run;
    
58. 
    
59. data del_&out.;
    
60. set &outdata. ;
    
61. if num&od.=0;
    
62. run;
    
63. 
    
64. data cumpter_c&out.;
    
65. set &outdata.(where=(num&od. ne 0));
    
66. by zip;
    
67. retain num&ct.;
    
68. 
    
69. if first.zip then  num&ct.=1;else num&ct.+1;
    
70. 
    
71. if num&ct. le 8 then do;
    
72. rate=1/(2**num&ct.);
    
73. if num&ct.=1 then want=q_1;else want=sum(q_1,rate);
    
74. end;
    
75. 
    
76. run;
    
77. %mend jj;
    
78. %jj(b,c,%str(1=1),0,1,1);
    
79. %jj(cumpter_c1,d,%str(num1 gt 8),2,2,3);
    
80. %jj(cumpter_c2,e,%str(num3 gt 8),4,3,5);
    
81. 
    
82. data final;
    
83. set del_1  cumpter_c1(where=(num1 le 8))
    
84.     del_2  cumpter_c2(where=(num3 le 8))
    
85.     del_3  cumpter_c3(where=(num5 le 8));
    
86.         if want=. then want=0;
    
87.         keep want q_1 zip order;
    
88. proc sort;by order;
    
89. run;

复制代码

data a;
x=100;
run;
data b;
set a;
sum=0;
retain sum;
do i=1 to 8;
sum+x;
x=x/2;
output;
end;
run;

data a;
x=100;
run;
data b;
set a;
sum=0;
retain sum;
do i=1 to 9;
sum+x;
x=x/2;
output;
end;
run;