如何解决在Logit模型中使用R中margins函数的“ at”参数
我希望能够分析logit模型中连续变量和二进制变量的边际效应。我希望R提供hp的独立边际效应的平均值(在此示例中为200),同时还找到vs变量的边际效应等于1。我希望输出表还包括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))
解决方法
这是因为您已将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
变成一个因素;这是二分法。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。