egen sic_year = group(粗分行业代码num year), label lname(sic_year)
         g 投资效率=.
         forvalues i = 1/218{
           qui reg 后一期新增投资额占均值资产比率 当期新增投资额占均值资产比率  Q 现金持有量t listyears 资产规模t 资产负债率 个股年收益率 if sic_year==`i', nocons
           qui predict e if e(sample), res
           qui replace 投资效率 = e if e(sample)
           drop e
         }