r - margins.plot：使用 'which' 参数来选择要包含在绘图中的边距

Question

我试图根据逻辑回归在 r中绘制边际效应。例如：

data <- mtcars
mod <- glm(am ~ cyl + hp + wt + mpg, family = binomial, data = data)

library(margins)
marg <- margins(mod, atmeans = TRUE)
summary(marg)

我可以运行边距绘图命令：

plot(marg)

它绘制了所有 IV 的边际效应和置信区间。我只想在情节中包含cyl我hp感兴趣的解释变量。根据 r 文档，这可以使用 'which' 参数来完成，该参数采用字符向量。但是，文档没有说明如何使用此参数。有谁知道如何使用“哪个”参数来要求 margins.plot 仅绘制选择的边际效应？不幸的是，上面链接的边距绘图帮助页面没有任何示例。

绘图图像

Answer 1

在绘图之前，我们可以使用函数中的variables选项指定感兴趣的变量margins()。

mod <- glm(am ~ cyl + hp + wt + mpg, family=binomial, data=mtcars)

library(margins)
marg <- margins(mod, variables=c("cyl", "hp"))

plot(marg)

给出：