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

在Logistic回归中排序

1 个答案: