[求助] SAS regression的程序編碼

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我的data是

gvkey	year	quarter	v1	v2	v3	v4	v5
1	1992	1	213	32	67	678	45
1	1992	2	546	87	87	123	3
1	1992	3	345	756	7567	235	123
1	1992	4	47	867	6	346	4
1	1993	1	6	6	546	547	3456
1	1993	2	546	234	346	679	546
1	1993	3	7547	5	5467	678	567
1	1993	4	54	346	567	56	679
1	1994	1	123	568	678	546	78078
1	1994	2	65	346	567	5687	978
1	1994	3	867	789	67	679	9
1	1994	4	6	346	6	78	967
1	1995	1	867	47	345	9677	856
1	1995	2	978	235	234	545	754
1	1995	3	967	1234	26	234	5
2	1992	2	8	56	8	1234	34
2	1992	3	78	34669	56	346	5
2	1992	4	0	54235	547	568	2345
2	1993	1	68	47	568	67	346
2	1993	2	54	346	89	978	547
2	1993	3	65	58	78	78	68
2	1993	4	7	57	56	967	67
2	1994	1	679	36	34	854	8567
2	1994	2	547	56	6457	65	342
3	1993	1	569	9	786	756	756
3	1993	2	80	7	7876	967	856
3	1993	3	689	579	56	56	78
3	1993	4	78	2352	76	754	568
3	1994	1	568	57	78	65	7
3	1994	2	7	54	767	8	867
3	1994	3	856	5465	435	7	867
3	1994	4	967	76	456	568	856
3	1995	1	856	7	546	678	54
3	1995	2	97	78	5678	978	634
3	1995	3	67	3	568	967	654
3	1995	4	568	434	7	867	654
3	1996	1	8967	6	856	867	7
3	1996	2	967	578	856	456745	568
3	1996	3	978	58	87	36	5687

我得run出兩個model:
model 1: v1~v2+v2
model 2: v1~v2+v3+v4+v5

重點是：
1. model是run repeadly by gvkey
2. 每一次的regression只能有12個數據，也就是說，
第一次的regression只有1-12的數據，第二次的regression是2-13的數據，以此類推。

而我需要的output就會像是：

gvkey	year	quarter	r2_model_1	r2_model_2
1	1994	4	xxx	xxx		這個是gvkey 1的第12個數據
1	1995	1	xxx	xxx		第13個數據
1	1995	2	xxx	xxx		以此類推
1	1995	3	xxx	xxx		gvkey 2因為數據不足12個，所以沒有output
3	1995	4	xxx	xxx
3	1996	1	xxx	xxx
3	1996	2	xxx	xxx
3	1996	3	xxx	xxx

其中，
r2_model_1和r2_model_2依序是model 1與model2的R-squared的數值。

求請各位幫忙解惑，拜託了！感謝！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：regression regressio regress REG ESS 程序

1. data test;
   
2. input gvkey year quarter v1 v2 v3 v4 v5 @@;
   
3. cards;
   
4. 1
   
5. 1992
   
6. 1
   
7. 213
   
8. 32
   
9. 67
   
10. 678
    
11. 45
    
12. 
    
13. 1
    
14. 1992
    
15. 2
    
16. 546
    
17. 87
    
18. 87
    
19. 123
    
20. 3
    
21. 
    
22. 1
    
23. 1992
    
24. 3
    
25. 345
    
26. 756
    
27. 7567
    
28. 235
    
29. 123
    
30. 
    
31. 1
    
32. 1992
    
33. 4
    
34. 47
    
35. 867
    
36. 6
    
37. 346
    
38. 4
    
39. 
    
40. 1
    
41. 1993
    
42. 1
    
43. 6
    
44. 6
    
45. 546
    
46. 547
    
47. 3456
    
48. 
    
49. 1
    
50. 1993
    
51. 2
    
52. 546
    
53. 234
    
54. 346
    
55. 679
    
56. 546
    
57. 
    
58. 1
    
59. 1993
    
60. 3
    
61. 7547
    
62. 5
    
63. 5467
    
64. 678
    
65. 567
    
66. 
    
67. 1
    
68. 1993
    
69. 4
    
70. 54
    
71. 346
    
72. 567
    
73. 56
    
74. 679
    
75. 
    
76. 1
    
77. 1994
    
78. 1
    
79. 123
    
80. 568
    
81. 678
    
82. 546
    
83. 78078
    
84. 
    
85. 1
    
86. 1994
    
87. 2
    
88. 65
    
89. 346
    
90. 567
    
91. 5687
    
92. 978
    
93. 
    
94. 1
    
95. 1994
    
96. 3
    
97. 867
    
98. 789
    
99. 67
    
100. 679
     
101. 9
     
102. 
     
103. 1
     
104. 1994
     
105. 4
     
106. 6
     
107. 346
     
108. 6
     
109. 78
     
110. 967
     
111. 
     
112. 1
     
113. 1995
     
114. 1
     
115. 867
     
116. 47
     
117. 345
     
118. 9677
     
119. 856
     
120. 
     
121. 1
     
122. 1995
     
123. 2
     
124. 978
     
125. 235
     
126. 234
     
127. 545
     
128. 754
     
129. 
     
130. 1
     
131. 1995
     
132. 3
     
133. 967
     
134. 1234
     
135. 26
     
136. 234
     
137. 5
     
138. 
     
139. 2
     
140. 1992
     
141. 2
     
142. 8
     
143. 56
     
144. 8
     
145. 1234
     
146. 34
     
147. 
     
148. 2
     
149. 1992
     
150. 3
     
151. 78
     
152. 34669
     
153. 56
     
154. 346
     
155. 5
     
156. 
     
157. 2
     
158. 1992
     
159. 4
     
160. 0
     
161. 54235
     
162. 547
     
163. 568
     
164. 2345
     
165. 
     
166. 2
     
167. 1993
     
168. 1
     
169. 68
     
170. 47
     
171. 568
     
172. 67
     
173. 346
     
174. 
     
175. 2
     
176. 1993
     
177. 2
     
178. 54
     
179. 346
     
180. 89
     
181. 978
     
182. 547
     
183. 
     
184. 2
     
185. 1993
     
186. 3
     
187. 65
     
188. 58
     
189. 78
     
190. 78
     
191. 68
     
192. 
     
193. 2
     
194. 1993
     
195. 4
     
196. 7
     
197. 57
     
198. 56
     
199. 967
     
200. 67
     
201. 
     
202. 2
     
203. 1994
     
204. 1
     
205. 679
     
206. 36
     
207. 34
     
208. 854
     
209. 8567
     
210. 
     
211. 2
     
212. 1994
     
213. 2
     
214. 547
     
215. 56
     
216. 6457
     
217. 65
     
218. 342
     
219. 
     
220. 3
     
221. 1993
     
222. 1
     
223. 569
     
224. 9
     
225. 786
     
226. 756
     
227. 756
     
228. 
     
229. 3
     
230. 1993
     
231. 2
     
232. 80
     
233. 7
     
234. 7876
     
235. 967
     
236. 856
     
237. 
     
238. 3
     
239. 1993
     
240. 3
     
241. 689
     
242. 579
     
243. 56
     
244. 56
     
245. 78
     
246. 
     
247. 3
     
248. 1993
     
249. 4
     
250. 78
     
251. 2352
     
252. 76
     
253. 754
     
254. 568
     
255. 
     
256. 3
     
257. 1994
     
258. 1
     
259. 568
     
260. 57
     
261. 78
     
262. 65
     
263. 7
     
264. 
     
265. 3
     
266. 1994
     
267. 2
     
268. 7
     
269. 54
     
270. 767
     
271. 8
     
272. 867
     
273. 
     
274. 3
     
275. 1994
     
276. 3
     
277. 856
     
278. 5465
     
279. 435
     
280. 7
     
281. 867
     
282. 
     
283. 3
     
284. 1994
     
285. 4
     
286. 967
     
287. 76
     
288. 456
     
289. 568
     
290. 856
     
291. 
     
292. 3
     
293. 1995
     
294. 1
     
295. 856
     
296. 7
     
297. 546
     
298. 678
     
299. 54
     
300. 
     
301. 3
     
302. 1995
     
303. 2
     
304. 97
     
305. 78
     
306. 5678
     
307. 978
     
308. 634
     
309. 
     
310. 3
     
311. 1995
     
312. 3
     
313. 67
     
314. 3
     
315. 568
     
316. 967
     
317. 654
     
318. 
     
319. 3
     
320. 1995
     
321. 4
     
322. 568
     
323. 434
     
324. 7
     
325. 867
     
326. 654
     
327. 
     
328. 3
     
329. 1996
     
330. 1
     
331. 8967
     
332. 6
     
333. 856
     
334. 867
     
335. 7
     
336. 
     
337. 3
     
338. 1996
     
339. 2
     
340. 967
     
341. 578
     
342. 856
     
343. 456745
     
344. 568
     
345. 
     
346. 3
     
347. 1996
     
348. 3
     
349. 978
     
350. 58
     
351. 87
     
352. 36
     
353. 5687
     
354. ;
     
355. run;
     
356. 
     
357. proc sort data=test;by gvkey year;run;
     
358. proc sql;
     
359. create table test1 as select *,count(gvkey) as count from test group by gvkey having count(gvkey)>=12;
     
360. select count(distinct gvkey) into:count_var from test1;
     
361. %let count_var=&count_var;
     
362. select  count into:count1-:count&count_var from (select distinct gvkey,count from test1);
     
363. quit;
     
364. %macro test;
     
365. %do i=1 %to &count_var;/*按组*/
     
366.    %let k=%eval(&i-1);
     
367.    %let count0=0;
     
368.    %do j=1 %to &&count&i-11;/*组内*/
     
369.       %let var_=%eval(&j+&&count&k);
     
370.           %let var_1=%eval(&var_+11);
     
371.       data test2;
     
372.             set test1(firstobs=&var_ obs=&var_1);
     
373.       run;
     
374. 
     
375.      proc reg data=test2 outest=test3(keep=_RSQ_) edf;
     
376.       model v1=v2 v3;
     
377.           model v1=v2 v3 v4 v5;
     
378.      run;
     
379.      proc transpose data=test3 out=test3(keep=col1 col2  rename=(col1=r2_model_1 col2=r2_model_2));
     
380.      run;
     
381.          data test4;
     
382.          set test3;
     
383.          set test1(firstobs=&var_1 obs=&var_1);
     
384.          run;
     
385.          %if &i=1 and &j=1 %then %do;
     
386.          data real_test;
     
387.          set test4;
     
388.          run;
     
389.          %end;
     
390.          %else %do;
     
391.          proc append data=test4 base=real_test;
     
392.          run;
     
393.          %end;
     
394.    %end;
     
395. %end;
     
396. %mend;
     
397. %test;
     
398. proc sql;
     
399. create table test_ as select gvkey,year,quarter,r2_model_1,r2_model_2 from real_test;
     
400. quit;

复制代码

mingfeng07 发表于 2014-8-18 18:06

感謝！非常感謝！感謝這好高級的編程！真的非常感謝！

h52101019 发表于 2014-8-19 02:17
感謝！非常感謝！感謝這好高級的編程！真的非常感謝！

indeed.非常好的程序.
这是我的程序.为什么我选择这样的程序呢? 当然不时因为和以上的程序不一样. 原因是它给出一个非常清晰的LOG 窗口展示.程序良好的解析使得运行过程一目了然. 当然前提是保持好的效率.但是,很显然,使用宏和程序运行效率并不是一个概念.更多的宏程序的运行效率非常低下.因为较好的效率往往来源于对SAS运行的良好理解从而导致的对某些高效率机制的正确使用.
京剧

1. ods select FitStatistics; ods graphics off;
   
2. %macro mdl(fobs);
   
3. proc reg data =test(firstobs =%eval(&fobs-12+1) obs =&fobs);
   
4. g_&fobs._1: model v1 =v2 v3;
   
5. g_&fobs._2: model v1 =v2 v3 v4 v5;
   
6. %mend mdl;
   
7. data t_run;
   
8. set test end =Eof;
   
9. by gvkey notsorted;
   
10. if _n_ =1 then call execute('options nonotes ls =150; ods output FitStatistics(persist=proc) =fs(keep =model Label2 nValue2 where =(Label2="R-Square"));');
    
11. if first.gvkey then n =0;
    
12. n ++1;
    
13. if n >= 12 then do;
    
14. _obs_ =_n_; output;
    
15. call execute('%mdl('||cats(_obs_)||')');
    
16. end;
    
17. if Eof then call execute('run;quit; ods output clear; options notes ls =120;');
    
18. run;
    
19. proc sql;
    
20. create table final as
    
21. select t1.gvkey, t1.year, t1.quarter, t2.nvalue2 as r2_model_1, t3.nvalue2 as r2_model_2
    
22. from t_run as t1
    
23. left join fs(where=(model like '%1') ) as t2 on t1._obs_ =input(scan(t2.model, 2, '_'), best.)
    
24. left join fs(where=(model like '%2') ) as t3 on t1._obs_ =input(scan(t3.model, 2, '_'), best.)
    
25. ;
    
26. quit;

复制代码

LOG Window:
NOTE: CALL EXECUTE generated line.
1 + options nonotes ls =150; ods output FitStatistics(persist=proc) =fs(keep =model Label2 nValue2 where =(Label2="R-Square"));
2 + proc reg data =test(firstobs =1 obs =12); g_12_1: model v1 =v2 v3; g_12_2: model v1 =v2 v3 v4 v5;
3 + proc reg data =test(firstobs =2 obs =13); g_13_1: model v1 =v2 v3; g_13_2: model v1 =v2 v3 v4 v5;
4 + proc reg data =test(firstobs =3 obs =14); g_14_1: model v1 =v2 v3; g_14_2: model v1 =v2 v3 v4 v5;
5 + proc reg data =test(firstobs =4 obs =15); g_15_1: model v1 =v2 v3; g_15_2: model v1 =v2 v3 v4 v5;
6 + proc reg data =test(firstobs =25 obs =36); g_36_1: model v1 =v2 v3; g_36_2: model v1 =v2 v3 v4 v5;
7 + proc reg data =test(firstobs =26 obs =37); g_37_1: model v1 =v2 v3; g_37_2: model v1 =v2 v3 v4 v5;
8 + proc reg data =test(firstobs =27 obs =38); g_38_1: model v1 =v2 v3; g_38_2: model v1 =v2 v3 v4 v5;
9 + proc reg data =test(firstobs =28 obs =39); g_39_1: model v1 =v2 v3; g_39_2: model v1 =v2 v3 v4 v5;
10 + run;
10 + quit;
10 + ods output clear; options notes ls =120;

jingju11 发表于 2014-8-19 05:07
indeed.非常好的程序.
这是我的程序.为什么我选择这样的程序呢? 当然不时因为和以上的程序不一样. 原因 ...

比我写的有条理、清晰多了，还是京剧大神厉害！

mingfeng07 发表于 2014-8-19 09:24
比我写的有条理、清晰多了，还是京剧大神厉害！

方法不同而已.京剧