Порядок в логистической регрессии

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