再问! 如何计算newey-west 的t值?

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

在下有数千家公司的数据（附件original中只摘了其中五家），对每家公司CO，进行time-series regression，按以下面这模型：Y=X1;

通过这个code，我可以得到，每家公司在每个模型下的截距intercept(我定义为alpha)以及t-statsitics of intercept; 还有其他一些统计值：

proc reg data=origianl;

ods output nobs=a anova=b fitstatistics=c parameterestimates=d;

model y=x1;

by co;

quit;

data final;

set a(where=(Label="Number of Observations Read") keep=co Label N);

set b(where=(Source="Model") keep=source FValue);

set c(where=(Label2="Adj R-Sq") keep=Label2 nValue2 rename=(nValue2=adj_R_sq));

set d(where=(Variable="Intercept") keep=Variable Estimate tValue Probt rename=(Estimate=alpha tValue=t_alpha Probt=p_alpha));

set d(where=(Variable="X1") keep=Variable Estimate tValue Probt rename=(Estimate=beta1 tValue=t_beta1 Probt=p_beta1));

drop Label Source Label2 Variable;

run;

但是发现很多人用的是Newey-West t值来看intercept (alpha)显著性。在下想请教一下：(1)为什么需要用Newey-West Heteroskedasticity and Autocorrealtion consistent standard error 来对t进行调整？我算出来的t值和这个Newey-West的t值有什么区别吗？ (2)如何用SAS CODE 来计算每家公司的newey-west t值？

我已经把我的数据截了一部分放在附件里，希望高手能够帮我解答一下！？万分感激！

扫码加我 拉你入群

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

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

分享0 收藏3 回帖

关键词：newey-west newey west Est NEW original source where 如何 模型

你如果用新方法或者才出来不久的方法就不要指望SAS自带的程序,你自己需要编写

Newey-west是1987年就出来的方法，在spss里有自带程序。用sas是可以做到的，我大概了解到可能用这个语句：
Proc model data=original;
  Param b0 b1;
  y=b0+b1*x1;
  fit y/gmm kernel=(bart,1, lnM/lnT);
  instruments x1;
by co;
run;
主要对kernel后面三个参数是什么意思不太了解，如何规定bart, 1, lnM/lnT (M=bandwidth,T=no.of the observations) , 好像还有用kernel=(bart, laglength+1,0) 的.
就我的例子，应该如何规定kernel?望高手指教一下！

真的没有人能够指点一下吗？版主能否解答一下？

期待高手来解决

1# funwin LZ，这个问题我知道，bart就这样打上去。 laglength要看你的数据多少。一般400个设为3。你可以参见eviews的newey and west操作。最后一项，你已经把公式列上去了。你就根据你的数据去算就可以了。

另外，我有一个问题，想向LZ请教一下，如何让SAS自动将一千家公司的回归都做完？如果一个一个的点，貌似要很久。

lisutong1987 发表于 2010-7-18 19:12
另外，我有一个问题，想向LZ请教一下，如何让SAS自动将一千家公司的回归都做完？如果一个一个的点，貌似要很久。

貌似把数据集按公司名排序，回归时用个by，就可以了，当然1000多个回归，你可能要指定个输出数据集，输出你要的值，

要是全在output窗口，会溢出。