Question

假设我有一个因变量和4个独立变量。我怀疑只有3个自变量是显着的，所以我使用glm（y~x1 + x2 + x3 ...）函数。然后我得到这些变量的一些系数。现在我想运行glm（y~x1 + x2 + x3 + x4），但我想指定x1，x2，x3系数保持不变。我怎么能做到这一点？

谢谢！

Answer 1

我认为你不能适应某些自变量具有固定参数的模型。你可以做的是创建一个新变量y2，它等于你的第一个模型x1+x2+x3的预测值。然后，您可以使用第二个模型y~y2+x4将其作为一个独立变量与x4一起包含。

所以基本上是这样的：

m1 <- glm(y~x1+x2+x3...)
data$y2 <- predict(glm, newdata=data)
m2 <- glm(y~y2+x4...)