为什么我会为“plot.margins”和“margins”获得不同的 AME？答案

【问题标题】：Why am I obtaining different AME's for "plot.margins" and "margins"?为什么我会为“plot.margins”和“margins”获得不同的 AME？
【发布时间】：2020-12-21 04:18:55
【问题描述】：

我正在使用“margins”包来获得逻辑回归的平均边际效应。当我尝试绘制这些边际效应时，绘图的结果与我从“边际”函数获得的结果不同。我不知道为什么会这样 D：

我用于回归的公式如下：

mod_7 <- glm(ID_izq_WV7 ~ sex_7 + edad_7 + edusup_7 + ingresos_7 + indiv_eco_7 +
             interv_gob_7 + lib_ind_WV7 + prefdem_7 + postmat_WV7 + confgob_WV7 +
             agnostic_7, family = binomial)

然后，我使用“margins”函数来获得平均边际效应：

summary(margins(mod_7, type = "link"))

从这段代码中，我得到以下结果： Results: Average Marginal Effects

然后，我使用以下代码为相同的回归公式创建一个图：

par(mar = c(8, 5, 4, 3))
plot(margins(mod_7), las = 2)

从这段代码中，我得到以下情节： Plot: Average Marginal Effects

图中显示的结果与我从“summary.margins”得到的结果不匹配。有谁知道为什么会这样？非常感谢您的帮助！

【问题讨论】：

如果您包含一个简单的reproducible example 以及可用于测试和验证可能的解决方案的示例输入，则更容易为您提供帮助。

标签： r plot statistics

【解决方案1】：

我仍然不知道为什么会发生这种情况，但我意识到图中显示的值与表中显示的值的顺序相同。变量的名称只是以不同的顺序打乱了。为了解决这个问题，我只是按正确的顺序标记了变量：

plot(margins_mod_europa_4,
 labels = c("Secularismo", "Confianza en instituciones gubernamentales",
           "Edad", "Educación Superior", "Individualismo económico",
           "Ingresos altos", "Intervención gubernamental", 
           "Valoración de libertades individuales", "Valores postmaterialistas",
           "Preferencia por la democracia", "Sexo = Mujer"),
 las = 2)

【讨论】：