0

在 R 中更改 margins-package 中的类型参数时,我有一个关于不同输出的问题。这是一个可重现的示例:

library(margins)
x <- glm(am ~ cyl + hp * wt, data = mtcars, family = binomial)
me <- margins(x, type = "response") 
me2 <- margins(x, type = "link")
summary(me)
summary(me2)

两个摘要在水平和重要性上都产生了不同的输出。所以我的问题是,为什么两种类型的 p 值不同?这是各自的输出:

> summary(me)
 factor     AME     SE       z      p   lower  upper
    cyl  0.0216 0.0493  0.4377 0.6616 -0.0750 0.1181
     hp  0.0027 0.0023  1.1596 0.2462 -0.0018 0.0072
     wt -0.5158 0.2685 -1.9209 0.0547 -1.0421 0.0105

> summary(me2)
 factor      AME     SE       z      p    lower  upper
    cyl   0.5156 1.1695  0.4409 0.6593  -1.7765 2.8077
     hp   0.0515 0.0357  1.4430 0.1490  -0.0185 0.1215
     wt -12.2426 7.6784 -1.5944 0.1108 -27.2920 2.8067

这是包的手册:

https://cran.r-project.org/web/packages/margins/vignettes/Introduction.html#Average_Marginal_Effects_and_Average_Partial_Effects

他们解释说“响应”类型返回部分效应,即每个变量对结果量表的贡献,以线性预测变量的链接函数转换中涉及的其他变量为条件,并且“链接”类型返回真正的边际效应,即每个变量在线性预测器尺度上的边际贡献。那么,如果只是尺度变换,为什么两种类型之间的显着性水平会发生变化呢?我想我不太了解真正的边际效应和部分效应之间的区别。任何帮助,将不胜感激。谢谢!

4

0 回答 0