margins.plot:使用' which‘参数选择要包含在绘图中的页边距

Question

我正在尝试基于逻辑回归来绘制r中的marginal effects。例如：

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

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

我可以运行margins plot命令：

plot(marg)

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

plot image

Answer 1

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

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

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

plot(marg)

提供：

​

​