多元Logit模型回归中因变量预测概率问题

         在多元Logit模型回归中，在回归结果分析后，想预测多分类因变量在某一自变量条件下的预测概率(包括置信区间)，
请问需要什么样的命令？本人用了“predict”和“adjust”命令没有正确算出来，此问题源自于一篇文章“农村独生子女生育选择影响因素计测度研究（作者：杜本峰）”，先截图如下，并上传该文章，请各位高手指点，谢谢！

margins, at(婚姻状态 = (1(1) 3)) predict(outcome(1)) vsquish
margins, at(婚姻状态 = (1(1) 3)) predict(outcome(2)) vsquish
margins, at(婚姻状态 = (1(1) 3)) predict(outcome(3)) vsquish