由于没有数据，只好给个例子供参考：

*- Reproduce the results in Roychowdhury(2006, JAE)
*  Roychowdhury, S., 2006,
*    Earnings management through real activities manipulation,
*    Journal of Accounting and Economics, 42 (3): 335-370.

   *-行业重新编码 1,2,3 ……
         cap drop sic123
         egen sic123 = group(sic2 year), label lname(sic_year)
            qui sum sic123
         global N = r(max)
         
   *-分行业, 分年度回归分析  
         dropvars DACC e
         gen DACC = .
         forvalues i = 1/$N{
           qui reg acc invA DS_DAR PPE if (sic123==`i')
           qui predict e if e(sample), res
           qui replace DACC = e if e(sample)
           drop e
         }

连老师，您写的我看的不是太懂。不过，我按照您之前教我的方法goolgle stata gini save.找到如下命令：
levels city, local(reg)
ge gini_reg = .
foreach r of local reg {
  ineqdeco income if city == `r'
  replace gini_reg = $S_gini if city == `r'
}
成功将各个city基尼系数保存，且与ineqdeco income, by(city)屏幕显示的一样。现在贴出来，请连老师多批评指正，也供广大学友借鉴

基本上是同一个思路，不错。