在logit模型中使用R中边距函数的"at“参数

Question

我希望能够在logit模型中分析连续变量和二元变量的边际效应。我希望R能提供hp的独立边际效应的平均值(在本例中为200)，同时还能找到等于1的vs变量的边际效应。我希望输出表还包括SE、p值和z得分。我在使用这个表时遇到了问题，当我让它运行时，它不会独立地计算这两个变量。下面是一个MRE。谢谢!

mod2 <- glm(am ~ hp + factor(vs), data=mtcars, family=binomial)

margins(mod2)
#> Average marginal effects
#> glm(formula = am ~ hp + factor(vs), family = binomial, data = mtcars)
#>        hp      vs1
#>  -0.00203 -0.03154

#code where I am trying to evaluate at the desired values. 
margins(mod2, at=list(hp=200, vs=1))

Answer 1

这是因为您已经将vs更改为一个因子。考虑以下几点

library(margins)
mod3 <- glm(am ~ hp + vs, data=mtcars, family=binomial)
margins(mod3, at=list(hp=200, vs=1))

# Average marginal effects at specified values
# glm(formula = am ~ hp + vs, family = binomial, data = mtcars)
# 
#  at(hp) at(vs)        hp       vs
#     200      1 -0.001783 -0.02803

这里没有真正的理由将vs变成一个因素；它是二分法的。