我正在使用一项来自随机对照试验的数据集，该数据旨在评估巴基斯坦女子乡村中学的存在对提供私立教育的影响。作为分析的一部分，我试图预测变量的边际效应，在本例中为变量GSSexist存在的1981年人口（pop81）使用不同的百分比值。执行此操作的Stata命令如下：
margins, dydx(i.GSSexist) at((p25)pop81) at((median)pop81) at((p75)pop81) at((p90)pop81)
其中px代表特定的百分位数。如何在R中使用“边距”命令执行此操作？
我尝试了以下代码，但没有得到任何有意义的结果。为了明确起见，Probit_1指的是我用来提取分别在第25、50、75和90％百分数值处的边际效应的概率模型。
summary(margins(probit_1, probs = c(0.25,0.50,0.75,0.90), variables = "pop81"))
关于如何改善代码以获得所需结果的任何想法吗？
非常感谢！