在Logistic回归中排序

Question

在逻辑回归中,SAS可以选择使用"降序"选项来模拟1而不是0.R中有什么方法,我们可以做同样的事情吗？

我正在使用的代码如下:

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

此致,Ari

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