各位前辈：以下是我编的程序，其中暂元里面包括了6个变量，与主题无关。
我的编程目的是，有13对变量，也就是26个变量，任取每三对（也就是aj ak aL bj bk bL）与暂元里的6个变量一起组成自变量，对因变量缩尾ROIC做分组回归，然后怎么把这暂元里的6个变量的系数和这三对变量的6个系数以及拟合优度提出来。
如果用 _b 命令的话，aj ak aL bj bk bL的系数所在列是不固定的，我想把这6个系数固定在6个列里。

local 暂元 缩尾上一年ξ-CRDA
      forvalues j=1/11{  
      local j1=`j'+1   
              forvalues k=`j1'/12{   
                   local k1=`k'+1
                       forvalues L=`k1'/13{   
                       statsby  拟合优度=e(r2)  自变量j=`j' 自变量k=`k' 自变量L=`L'  ,saving(D:\result\r`j'_`k'_`L'.dta,replace): reghdfe 缩尾ROIC `暂元' a`j' a`k' a`L' b`j' b`k' b`L', absorb(公司 年份) vce(cluster 公司)
                  
                  }
                  }
                  }


这是在statsby语句里使用_b 之后，可见生成了26列，不便于最后比较结果                  



这是我的目标样式，即在自变量L列之后，生成暂元里的6个变量系数和固定的6个aj ak aL bj bk bL系数