Hi! I ran a multivariate probit regression on 3 binary dependent variables using Capellari and Jenkins' -mvprobit- package and calculated the predicted marginal probability for each of the 3 outcome variables using the postestimation command -mvppred pmargm, pmarg-. How can I compute the marginal effects of a regressor(s) on each predicted marginal probability? I tried using -margins- as follows but the error message indicates that "option pmargm1 not allowed."

margins, dydx(*) atmeans predict(pmargm1) post