如何将回归得到的变量系数保存

ywh19860616 发表于 2014-5-22 21:22

您好，我您给我的代码试验了一次，有一个错误提示conformability error，请问这个怎么办？可不可以分两次做，就是说一套代码专门获得系数的，一套代码专门用来获得P值得

1. mat b = J(1000,1,0)
   
2. mat se = J(1000,1,0)
   
3. mat p = J(1000,1,0)
   
4. 
   
5. preserve
   
6. forvalues i=1/1000{
   
7. sort year
   
8. sample 100, count by(year)
   
9. save `i'.dta,replace
   
10. restore,preserve
    
11. 
    
12. merge 1:1 num year using `i'.dta
    
13. sort year
    
14. gen part`i'=1 if _merge==3
    
15. replace part`i'=0 if part`i'==.
    
16. drop _merge
    
17. regress DA2 part`i',level(95)
    
18. mat b[`i',1] = _b[part`i']
    
19. mat se[`i',1] = _se[part`i']
    
20. mat p[`i',1] = 2*ttail(e(df_r),abs(_b[part`i']/_se[part`i']))
    
21. test part`i'
    
22. }
    
23. 
    
24. matrix list b

复制代码

这个是可以的吧，如果不行，把你的数据上传。

ywh19860616 发表于 2014-5-22 21:34

太好了，您真强，成功了。请问您是怎么学的STATA，感觉您逻辑太清晰了。

jinlong1835 发表于 2014-5-20 23:50
多谢高人，在你的帮助下，刚才我试出来了，谢谢您。

求问，您试出来最终是代码是怎么编的呢？还想请教下您后来二项检验的代码是怎么编的？以及文献中说对jones模型进行分年度分行业回归，最终得出的系数是一个，怎么理解您解决了吗？