我现在需要根据上图的(5)式进行回归，但现在有一些问题想要请教各位坛友：
首先，我的数据如下图：

也就是，有四个维度：年份-企业-出口目的国-贸易方式，即在某一年，a企业以某种贸易方式出口产品到某个国家（其中，上图的hs6代表产品的6位海关编码；company指的是企业；country指的是出口目的国；shipment指的是贸易方式）。
那么对于这种类型的数据，我要用什么命令回归(5)式呢（一定要保留贸易方式）？我最后的目的是要通过回归（5）式得到这个式子的残差。
ps:如果有坛友知道的话，希望能够帮助一下，万分感谢。最后，附上一部分数据。