Hi
I ran a gologit2 model and then computed predicted probabilities for the ordinal dependent variable using the margins command. However, I found a ‘implausible” negative probability value for the outcome(2). How can I address this negative probability value? Any advice would be appreciated.

. margins, at(susr=3 offr=1 susw=3 vio=55.731)

Predictive margins Number of obs = 3,342
Model VCE : OIM

1._predict : Pr(forc==1), predict(pr outcome(1))
2._predict : Pr(forc==2), predict(pr outcome(2))
3._predict : Pr(forc==3), predict(pr outcome(3))
at : susr = 3
offr = 1
susw = 3
vio = 55.731

------------------------------------------------------------------------------
| Delta-method
| Margin Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
_predict |
1 | .4373939 .0940638 4.65 0.000 .2530323 .6217555
2 | -.2668683 .063431 -4.21 0.000 -.3911907 -.1425458
3 | .8294743 .0555001 14.95 0.000 .7206961 .9382526
------------------------------------------------------------------------------