利用工业企业数据库面板数据，做了一个回归，想实现：
xi：reg l lnsize capital growth i.cic i.year 其中，cic为四位数行业
然后调用各解释变量的系数constant,a,b,c,d,e....,f..，生成：
l_expect = constant+a*lnsize + b*capital +c*growth +e*cic(每个企业的四位数行业代码对应的上述回归的系数e)+f*year(如：上述回归2005年系数为0.2 ，此处则为0.2*2004)

请问应当如何实现？

回归部分结果如图：

其他的变量好说，主要是想把各行业的系数保存并与企业自身的行业代码匹配，然后生成l_expect。
自己能想到的方法是先把各行业代码和对应系数保存出来，然后与原数据基于行业代码merge，然后再算。但是第一步的保存不知道如何实现。求解答，跪谢！