编程forvalues循环问题求助，内外循环，年份i循环上加不了，循环中断

在做环境不确定性的回归时，是要对满足YEAR为1 ~5的代码进行回归，并且展开循环，分为内外两个循环。但是当内部循环遇到不满足条件的reg的时候就会自动终止，提示no observation,或者insufficient observation.猜测由此终止了循环，更不用提外部的年份叠加的循环了，论坛中有老师指出可以用levelsof，但不知具体如何处理，是不是需要加一个内循环，进行continue或者break.附上数据和代码。谢谢各位！
bysort stkcd:egen n=count( year)
drop if n<5
encode(stkcd),gen(id)
xtset id year
egen idnew = group(id)
sum idnew
local N = r(max)
gen msale=.
gen dev=.
forvalues i = 2000(1)2011{
forvalues j = 1/876{
gen YEAR= year-`i'
qui reg sale YEAR if YEAR>0& YEAR<6 & idnew==`j'
qui predict p if YEAR<6 & YEAR>0& idnew==`j',r
qui sum p if idnew==`j',de
qui dis "i=`i' j=`j'"
qui replace dev=r(sd) if idnew==`j' & YEAR==5
qui sum sale if idnew==`j' & YEAR<6 & YEAR>0 ,de
qui replace  msale=r(mean) if idnew==`j' & YEAR==5
drop p YEAR
}
}

2001-2016中小板创业板sale.xls (531.5 KB, 需要: 5000 个论坛币)

2018-6-14 19:07:24 上传

需要: 5000 个论坛币  [购买]

2001-2016中小板创业板sale.dta (144.25 KB, 需要: 5000 个论坛币)

2018-6-14 19:07:24 上传

需要: 5000 个论坛币  [购买]

1. bysort stkcd:egen n=count( year)
   
2. drop if n<5
   
3. encode(stkcd),gen(id)
   
4. 
   
5. xtset id year
   
6. egen idnew = group(id)
   
7. sum  idnew
   
8. local N = r(max)
   
9. gen msale=.
   
10. gen dev=.
    
11. 
    
12. 
    
13. forvalues i = 2000(1)2011{
    
14. forvalues j = 1/876{
    
15. cap drop p YEAR
    
16. gen YEAR= year-`i'
    
17. cap qui reg sale YEAR if YEAR>0& YEAR<6 & idnew==`j'
    
18. cap qui predict p     if YEAR<6 & YEAR>0& idnew==`j',r
    
19. cap qui sum p if idnew==`j',de
    
20. cap qui dis "i=`i' j=`j'"
    
21. cap qui replace dev=r(sd) if idnew==`j' & YEAR==5
    
22. cap qui sum sale if idnew==`j' & YEAR<6 & YEAR>0 ,de
    
23. cap qui replace  msale=r(mean) if idnew==`j' & YEAR==5
    
24. cap drop p YEAR
    
25. }
    
26. }

复制代码

在reg前加capture noisily也试过，但是在YEAR为6时也就直接终止了，循环进行不下去。请教各位！

于果 发表于 2018-6-14 19:09
在reg前加capture noisily也试过，但是在YEAR为6时也就直接终止了，循环进行不下去。请教各位！

1. 你若要问程序，永远附上资料；若要附上资料，永远用 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。2. 類似這種情況，用 loop 基本上都是不好的作法，總是有其他更簡潔與迅速之作法！

蓝色 发表于 2018-6-15 10:34

非常感谢您的回复！谢谢！刚刚按照您的程序运行了一下，出现 variable YEAR already defined 的提示，但是命令中已有在代码前进行cap drop YEAR p的步骤，不知为何仍然会出现YEAR已经存在的错误提示，还麻烦您指点~谢谢！

1. bysort stkcd:egen n=count( year)
   
2. drop if n<5
   
3. encode(stkcd),gen(id)
   
4. 
   
5. xtset id year
   
6. egen idnew = group(id)
   
7. sum  idnew
   
8. local N = r(max)
   
9. gen msale=.
   
10. gen dev=.
    
11. 
    
12. 
    
13. forvalues i = 2000(1)2011{
    
14. forvalues j = 1/876{
    
15. cap drop p
    
16. cap drop YEAR
    
17. gen YEAR= year-`i'
    
18. cap qui reg sale YEAR if YEAR>0& YEAR<6 & idnew==`j'
    
19. cap qui predict p     if YEAR<6 & YEAR>0& idnew==`j',r
    
20. cap qui sum p if idnew==`j',de
    
21. cap qui dis "i=`i' j=`j'"
    
22. cap qui replace dev=r(sd) if idnew==`j' & YEAR==5
    
23. cap qui sum sale if idnew==`j' & YEAR<6 & YEAR>0 ,de
    
24. cap qui replace  msale=r(mean) if idnew==`j' & YEAR==5
    
25. cap drop p
    
26. cap drop YEAR
    
27. }
    
28. }

复制代码

cap drop p
cap drop YEAR

拆成2行应该就没有问题

蓝色 发表于 2018-6-15 12:10
cap drop p
cap drop YEAR

谢谢蓝色大大！感恩~

谢谢大大！