使用已知系数和未知截距拟合glm

Question

我正在尝试使用glm拟合逻辑回归模型，其中我只对截距感兴趣 - 但我仍然希望模型与已知系数拟合。 例如：

或

beta <- c(24.5,3.6,2.87,7.32)

所以我想用

model <- glm(y~x_1+x_2+x_3+x_4, family=binomial(link="logit"), data=dt)

并以某种方式包含已知的beta，因此glm函数仅适合alpha。我怎样才能做到这一点？

Answer 1

使用偏移量，将已知术语添加到GLM的线性预测变量（公式的RHS，logit标度）。

beta <- c(24.5, 3.6, 2.87, 7.32)
dt <- transform(dt,
         pred=beta[1]*x_1+beta[2]*x_2+beta[3]*x_3+beta[4]*x_4)
model <- glm(y~1+offset(pred), family=binomial(link="logit"), data=dt)