R中有序probit/logit模型中每个变量的单一边际效应

我必须重新创建一个使用R在Stata中构建的模型。我想在使用有序probit/logit模型后计算边际效应。我有一个因变量（序数），有5个结果和几个自变量

我尝试了“MASS”和“erer”以及“oglmx”软件包，所有这些都为我的DV的每一个结果提供了边际效果。比如：

x