[求助] SAS regression的程序編碼

h52101019

2717

收藏 2014-08-18

我的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的數值。

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

扫码加我拉你入群

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

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

全部回复

mingfeng07

2014-8-18 18:06:30

复制代码

扫码加我拉你入群

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

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

h52101019

2014-8-19 02:17:19

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

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

扫码加我拉你入群

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

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

jingju11

2014-8-19 05:07:31

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

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

复制代码

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;

扫码加我拉你入群

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

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

mingfeng07

2014-8-19 09:24:10

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

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

扫码加我拉你入群

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

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

jingju11

2014-8-20 05:26:43

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

方法不同而已.京剧

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群