请教：用stata做事件研究法用到循环语句，因为样本量大，回归结果很慢怎么改进？

请教：我用stata做事件研究法计算CAR，用到循环语句时大概要执行16000+次，速度非常慢，我等了有三个多小时都还没回归完。希望大家能帮我看看，该如何改进？只有10个论坛币，别嫌少啊

gen predicted_return=.
forvalues i=1(1)16307 {
qui reg ret market_ret if ID==`i' & estimation_window==1
predict p if ID==`i'
replace predicted_return=p if ID==`i' & event_window==1
drop p
}

尔后建议请用 dataex (先 ssc install dataex 并见说明) 将原始 Stata 资料中具有”代表性”的一部分资料列出，以供有意回答者实验之用，并能提供具体操作指令。并请参考 http://www.jianshu.com/p/9870080fe769,  https://bbs.pinggu.org/thread-5048204-1-1.html, 与 https://bbs.pinggu.org/thread-5917273-1-1.html。

顶一下

大神在哪里？在线等~求助

黃河泉 发表于 2018-1-11 07:09
尔后建议请用 dataex (先 ssc install dataex 并见说明) 将原始 Stata 资料中具有”代表性”的一部分资料列 ...

赞！其实作为提问者，也是苦于不知道如何把自己的问题说清楚，以便于更好地得到大家的帮助。谢谢啦~我现在就来试试，哈哈！

好吧，失败了。说我的字数超过15000，心累~

试一下

1. * Example generated by -dataex-. To install: ssc install dataex
   
2. clear
   
3. input float ID str6 stkcd str7 id float(date event_date dif ret) double market_ret float(event_window estimation_window)
   
4. 1 "000001" "0000011" 17170 17611 -266 -.024879               .008512 0 0
   
5. 1 "000001" "0000011" 17171 17611 -265 -.074415                .00978 0 0
   
6. 1 "000001" "0000011" 17174 17611 -264  .023737                .03352 0 0
   
7. 1 "000001" "0000011" 17175 17611 -263  .032162               .029485 0 0
   
8. 1 "000001" "0000011" 17176 17611 -262   .02029               .025921 0 0
   
9. 1 "000001" "0000011" 17177 17611 -261  .049716 -.0057540000000000004 0 0
   
10. 1 "000001" "0000011" 17178 17611 -260  .050068              -.028983 0 0
    
11. 1 "000001" "0000011" 17181 17611 -259  .050258               .055619 0 0
    
12. 1 "000001" "0000011" 17182 17611 -258  .044172               .030431 0 0
    
13. 1 "000001" "0000011" 17183 17611 -257        0              -.015147 0 0
    
14. 1 "000001" "0000011" 17184 17611 -256 -.049354               .009873 0 0
    
15. 1 "000001" "0000011" 17185 17611 -255   .04759               .035119 0 0
    
16. 1 "000001" "0000011" 17188 17611 -254  .047788                .04398 0 0
    
17. 1 "000001" "0000011" 17189 17611 -253  .019144               .006285 0 0
    
18. 1 "000001" "0000011" 17190 17611 -252 -.021547                 .0108 0 0
    
19. 1 "000001" "0000011" 17191 17611 -251 -.017504               -.03783 0 0
    
20. 1 "000001" "0000011" 17192 17611 -250      .05               .018095 0 0
    
21. 1 "000001" "0000011" 17195 17611 -249  .049808               .028906 0 0
    
22. 1 "000001" "0000011" 17196 17611 -248  .050052              -.007395 0 0
    
23. 1 "000001" "0000011" 17197 17611 -247 -.050149              -.059873 0 0
    
24. end
    
25. format %dCY-N-D date
    
26. format %dCY-N-D event_date

复制代码

循环语句

1. gen predicted_return=.
   
2. forvalues i=1(1)16307 {
   
3.         qui reg ret market_ret if ID==`i' & estimation_window==1
   
4.         predict p if ID==`i'
   
5.         replace predicted_return=p if ID==`i' & event_window==1
   
6.         drop p
   
7.         }

复制代码

lihao_ 发表于 2018-1-11 09:56
试一下

循环语句

1. 这个资料不具代表性。请寄一个 ID=1＆2 的档给我 (river@mail.tku.edu.tw)，我来试试(或问问)看。2. 这种类似问题用 loop 来处理大概都不是好 idea 。

黃河泉 发表于 2018-1-11 10:39
1. 这个资料不具代表性。请寄一个 ID=1＆2 的档给我 (river@mail.tku.edu.tw)，我来试试(或问问)看。2. 这 ...

数据已经发到您邮箱，特别感谢！

lihao_ 发表于 2018-1-11 11:24
数据已经发到您邮箱，特别感谢！

我将你的问题 post 美国 Stata 论坛去寻求帮忙，请摒息以待。