我正在进行有序的概率回归分析。我现在使用R中的erer包计算每个解释变量对我的因变量（DV = 0,1,2）的边际效应。
现在，我的问题是如何绘制边际效应？我想得到的情节类似于这个here。
我使用命令polr来估计有序的概率回归。我在模型中有连续和二分解释变量。
我的模型看起来像这样：model<-polr(y~x1,x2,x3,x4,x5,data=mydata,Hess=TRUE, method="probit")
对于边际效应，我使用erer包中的命令：MEs<–ocME(model)
非常欢迎所有评论。请记住，我在问一个与有序概率回归分析有关的问题。