求助Mlogit模型交互项的边际效用的stata命令

大家好，最近正在写毕业论文遇到一个棘手的问题，希望大家能帮帮忙。
我的模型是多项logit模型，具体为
Y（0，1，2）=aX1+bX2+cX3+dX4+eOther, 其中X1，X2，X3，X4都是虚拟变量的交互项，比如X1=r*D1,X2=r*D2,而D1,D2是分位数虚拟变量。那么如何用stata计算出Y=1,Y=2时X1，X2等这些交互项的边际效应？

我用的Mlogit命令为  Mlogit X1 X2 X3 X4 Other, robust 关键是怎么求边际效应？

试试在最后加上rrr，我也正在准备做这个模型，我的问题更基础，我的原始数据是8000+微观截面数据，有15个选择，怎么整理成条件logit可以用的数据格式？