r - Logistic 回归中的排序

Question

在逻辑回归中，SAS 可以选择使用“降序”选项对 1 而不是 0 进行建模。 R 中有什么方法可以让我们做同样的事情吗？

我正在使用的代码如下:

glm(y~x1+x2+x3, family=binomial(link="logit"), na.action=na.pass)

问候，
阿里

Answer 1

该选项与建模完全相同1-y , 并将返回相同的系数，但具有不同的符号。所以要么你把1-y在模型中，或者你只是反转你的系数:

Data <- data.frame(
    y = rbinom(100,1,0.5),
    x1 = rnorm(100),
    x2 = rnorm(100),
    x3 = rnorm(100)
)

mod1 <- glm(y~x1+x2+x3, family=binomial(link="logit"), 
           na.action=na.pass,data=Data)

mod2 <- glm((1-y)~x1+x2+x3, family=binomial(link="logit"),
           na.action=na.pass,data=Data)

> all.equal(coef(mod2),-coef(mod1))
[1] TRUE