各位前辈和老师好，我在看一篇新冠肺炎疫情对国际贸易的影响的论文中，发现作者的模型中包含了虚拟变量向量，如下图，作者使用的是PPML方法，那么这个在用stata进行回归的时候，该如何实现呢？

我试着用以下代码，但是发现结果并不理想，并且还有一个变量cdi8被omited掉了。请问各位大佬，该怎么写代码实现呢？
forvalues j = 1/8{
           gen d`j' = 1 if month == `j'
                   replace d`j' = 0 if d`j' == .
         }
forvalues j = 1/8{
           gen cdi`j' = d`j'*casei
                   gen cdj`j' = d`j'*casej
         }
ppmlhdfe exp case cdi2 cdi3 cdi4 cdi5 cdi6 cdi7 cdi8 cdj1 cdj2 cdj3 cdj4 cdj5 cdj6 cdj7 cdj8,absorb(code#year code#month) nolog