想要分年度分行业的回归，设想，每次回归返回系数 用的forvalue ,但是这执行了一次，程序如下

egen newvar = group(year nindcd1)
qui su newvar
global N = r(max)
forvalue i= 0/143{
   reg Y X1 X2 X3 ROA if newvar==i
   gen ncoeff=_b[ROA]
  replace ncoeff=_b[ROA] if newvar==i
   }