全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学
2010-4-9 22:01:23
以下只是我个人的一些思路,希望高手能够指教一下!

Proc sql;


Create table newdata (drop=_1-_100) as


Select *, pre+_1 as y1,


pre+_2 as y2,



………


pre+_100 as y100,

From originaldata;

Quit;

怎么能把这个100步,用loop或宏语句表示呢?

Proc reg data=newdata;


Model y1=x1 x2 x3 x4


By co;

Run;

……….

Proc reg data=newdata;


Model y100=x1 x2 x3 x4;


By co;

Run;

这里也涉及到怎么用loop或宏语句表示.还有如何把我要的每家公司的每次sampleintercept, t value of intercept显示出来?用ods语句?

Data pvalue;


Count=0;


Numiters=100;


If t_new>t_obs then count=count+1;



End;





p_value=count/numiters;

Run;

用这个语句能够计算出我所要的p-value吗?
二维码

扫码加我 拉你入群

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

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

2010-4-10 00:34:57
31# funwin
什么是pre?

I see, the PRE is the old Y without including intercept
二维码

扫码加我 拉你入群

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

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

2010-4-10 02:51:54
yeah, pre is just a variable, equals to b*X1+c*X2+d*X3..., that's right, you got it,
二维码

扫码加我 拉你入群

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

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

2010-4-10 03:13:20
funwin 发表于 2010-4-9 22:01
以下只是我个人的一些思路,希望高手能够指教一下!

Proc sql;



Create table newdata (drop=_1-_100) as



Select *, pre+_1 as y1,



pre+_2 as y2,




………



pre+_100 as y100,


From originaldata;


Quit;


怎么能把这个100步,用loop或宏语句表示呢?

Proc reg data=newdata;



Model y1=x1 x2 x3 x4



By co;


Run;


……….


Proc reg data=newdata;



Model y100=x1 x2 x3 x4;



By co;


Run;


这里也涉及到怎么用loop或宏语句表示.还有如何把我要的每家公司的每次sampleintercept, t value of intercept显示出来?用ods语句?

Data pvalue;



Count=0;



Numiters=100;



If t_new>t_obs then count=count+1;




End;






p_value=count/numiters;


Run;


用这个语句能够计算出我所要的p-value吗?
复制代码
二维码

扫码加我 拉你入群

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

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

2010-4-10 03:25:41
For each Y, you have 100 corresponding bootstrap-Y’s and thus 100 bootstrap-intercepts. Among those 100 intercepts' T statistics, how many exceed the T of the original intercept (only 1 for each company)? If that is what you are asking, so
复制代码
二维码

扫码加我 拉你入群

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

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

2010-4-10 07:48:49
34# jingju11

So useful. It's been done. I really appreciate your kindly help.
Have a good weekend.
二维码

扫码加我 拉你入群

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

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

2010-4-10 07:53:28
35# jingju11

Yeah, that's what I need. You've totally grabbed my idea for data processing. Cheers.
二维码

扫码加我 拉你入群

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

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

2010-5-30 05:09:10
上次多谢你指点,帮我解决了bootstrap中的残值再抽样问题,非常感谢。
现在想用另一种bootstrap方法(这simulation方法,我称之为jointly sampling.):
1, 对每家公司,按公式:Yi,t=ai+bi*Bi,t+ci*Ci,t 做时间序列回归。得出每家公司的ai, 以及 Zi,t=Yi,t-ai;
2, 然后对每家公司的时间t进行有放回的再抽样(抽取300次,我的研究时间段一共300个月),不管每家公司原来有多少期,都抽出同样300期。比如有一家公司虽然只有70期,但也抽出300期来。
(上次进行的是残值再抽样,而且抽取的次数和公司原来的期数一致)
3;做 Zi,t =alphai+bi*Bi,t+ci*Ci,t 时间序列回归(每个公司都有抽取出来的300期了),得出一个新的alphai and t-statistic of alpha ;
5. 步骤2-3,重复1000次。也就是说每家公司最后有1000个 new alpha and t-statitics of alpha.
不知如何能够实现?望jingju11能否再次指教一下!非常感激
二维码

扫码加我 拉你入群

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

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

栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群