大家好，我是一个刚使用stata不久的小白，最近再使用bivariate probit模型时，不知道怎么求边际效应。我具体的分析步骤是这样的，1、首先估计bivariate probit模型：biprobit (Y1=X1 X2 X3 X4 X5) (Y2=X1 X6 X7  X8 X9),r partial difficult


2、计算在方程1和方程2中边际效应，我在researchgate上，看到有人是这样用的：
margins, dydx(X1) atmeans predict(pmarg1)  //估计在方程1(Y1=X1 X2 X3 X4 X5)中X1的边际效应
margins, dydx(X1) atmeans predict(pmarg2)  //估计在方程2(Y2=X1 X6 X7  X8 X9)中X1的边际效应

然而这样得出的第一个方向边际效应达到了19.5%，但是如果我只是做probit Y1 X1 X2 X3 X4 X5，然后margins,dydx(X1) atmeans，这样得出的边际效应只有2.5%。

虽然我知道不能直接对比19.5%和2.5%，但是两者的差距那么大，我是不是在第2步的时候，那里做错了，望各位会的大神指正一下！小弟非常感谢！！！