在下有数千家公司的数据(附件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值?
我已经把我的数据截了一部分放在附件里,希望高手能够帮我解答一下!?万分感激!