我明白你想要在margins包生成的边际效应图里只保留cyl和hp这两个变量，确实可以通过plot()函数的which参数轻松实现这个需求，具体操作方法如下：

核心操作逻辑 #

plot.margins()（也就是调用plot(marg)时实际运行的函数）的which参数接受一个字符向量，你只需要把想要展示的自变量名称按顺序传入即可。变量名称要和模型中使用的名称完全一致（注意大小写敏感）。

完整示例代码 #

# 加载数据并拟合logistic回归模型
data <- mtcars
mod <- glm(am ~ cyl + hp + wt + mpg, family = binomial, data = data)

# 加载margins包并计算均值处的边际效应
library(margins)
marg <- margins(mod, atmeans = TRUE)

# 仅绘制cyl和hp的边际效应及置信区间
plot(marg, which = c("cyl", "hp"))

小提示 #

如果不确定变量的准确名称，可以先运行summary(marg)查看输出结果里的变量列，确保which参数中的名称和输出里的完全匹配，避免因名称拼写错误导致的问题。

内容的提问来源于stack exchange，提问作者dc_ps_stat