每个公司同时分行业分国家回归，求每个公司的残差。代码怎么写？

数据分不同行业不同国家。现在要对每个公司同时按行业按国家回归，想要得到的是每家公司的残差是多少。这个代码怎么写？小白求指教。

egen group=group(country industry)
sum group
local n=r(max)
gen e=.
forvalue i=1/`n' {
reg  y x1 x2 x3 x4 if g==`i'   
est store m_`i'
predict r if e(sample), res
replace e=r if e(sample)
drop r
}
drop g

est store m_`i' 是我为了看回归系数加的，如果只需要求残差，可以去掉