分行业回归。因为所用数据库的两位数行业代码不连续，所以分了两段：forvalues i=13/15{} 和forvalues i=17/37{}，大括号里面是xi:opreg lny if ind2==`i',exit(exit) state(age lnk) proxy(lnI) cvars(d_state d_ex)  vce(bootstrap, seed(1357) reps(5))

请问各位大神，这样的话，自变量系数该如何提取？

查看了其它帖子里，说提取系数用 statsby _b _se, by(id): reg Y X。可是把statsby放在opreg前面会报错。


还有帖子说用bcoeff或bcoeffs。因为我是分两段循环的，第一段的循环语句里加入bcoeff以后，运行时会显示_b_lnk和_b_lnl已经生成，下面就没法再运行下去了。也就是，得到的系数还是i = 13 时的，没有实现“分行业回归，得到不同自变量系数”的目的。


那OP法，如果是分段循环，自变量系数该用什么命令提取出来呢？