面板数据做多层回归分析交互效应如何进行？

求各位老师帮忙！

非平衡面板数据
主效应包括三个自变量一个因变量
调节效应1加入一个公司层面的变量
以上选择了xtreg fe r进行了回归。

调节效应2想加入一个区域层面的变量，做一个调节的调节，三重交互的那种

请问是否可以用xtmixed进行回归分析呢？如果可以的话，stata语句这样写是否正确？
xtmixed dv iv1 iv2 iv3 ctrl mod1 iv1*mod1 iv2*mod1 iv3*mod1 mod2 iv1*mod2 iv2*mod2 iv3*mod2 mod1*mod2 iv1*mod1*mod2 iv2*mod1*mod2 iv3*mod1*mod2||cid:||id: ,cov(un)

有几个问题：
1.xtmixed的两层如何定义，上面指令中我的cid是区域分类 id是企业个人分类 是写成||cid:||id: 还是||id:||cid:
2.两个层次后是否需要写iv1 iv2 iv3？是一起放进模型还是三个自变量分别进行三重交互分析？
3.是否i.year还需要进入方程?
4. 是否还有更好的办法?

谢谢！！