使用多因子预测器从 GLM 中删除截距

Question

我正在 R 中运行带有 logit 链接函数的二项式逻辑回归。我的响应是阶乘 [0/1] 并且我有两个多级阶乘预测变量 - 我们称它们为 a 和 b，其中 a 有 4 个因子水平(a1、a2、a3 ,a4) 和 b 有 9 个因子水平 (b1,b2...b9)。所以:

mod <- glm(y~a+b, family=binomial(logit),data=pretend) summary(mod)

然后模型输出将显示有关模型的所有信息以及系数。

摘要输出中缺少 a 和 b(a1 和 b1)的因子水平。我知道它固定在模型的“拦截”中。我已经读过，如果我想删除截距项并查看这些因子水平的估计值，我可以在模型公式中添加 -1 或 +0，即:

mod2 <- glm(y~a+b-1, family=binomial(logit),data=pretend) 

...或者...
mod2 <- glm(y~a+b+0, family=binomial(logit),data=pretend)
总结(mod2)

在新模型 (mod2) 中，截距项消失了，变量 a 的因子水平 a1 在系数列表中给出。但是，变量 b 的因子水平 b1 仍然缺失，并且鉴于不再有截距项，那么我如何解释该因子水平的优势比？

有人可以向我解释如何获得 b1 的系数以及为什么会发生这种情况吗？

谢谢你。

Answer 1

有趣的是a1给出。人们会期望一个因子水平作为“引用”，因此在输出中没有任何 OR(因为它是 1.0)。

我想 b1是您的引用，因此是隐藏的，因此是 1.0。