在R中使用glm函数为零建模

Question

我正在尝试使用R中的glm函数构建逻辑回归模型。我的因变量是仅包含0和1的二项式。 此处0-不归，1-归。

我想建模不归（0），但R的glm函数默认构建为1。 就像在SAS中默认情况下会构建较低的值，并且我们可以在proc logistic中使用降序属性来更改顺序一样，glm中是否也有类似的东西？

我可以通过将原始数据中的0更改为1（反之亦然）来实现此目的，但是不想更改原始数据。

请帮助我或指导我如何在R中执行类似的操作。

提前致谢。

Answer 1

只需指定1 - y作为DV：

set.seed(42)
y <- sample(c(0, 1), 10, TRUE)
#[1] 1 1 0 1 1 1 1 0 1 1

fit <- glm(y ~ 1, family = binomial)
coef(fit)
# (Intercept) 
# 1.386294 
log(mean(y) / (1 - mean(y)))
#[1] 1.386294

1 - y
#[1] 0 0 1 0 0 0 0 1 0 0

fit1 <- glm(1 - y ~ 1, family = binomial)
coef(fit1)
#(Intercept) 
#-1.386294 
log(mean(1 - y) / (1 - mean(1 - y)))
#[1] -1.386294

Answer 2

另外，您可以使用... transform临时转换数据：

glm( data = transform( data.frame(y=0), y=y+1 ), ... )