ロジスティック回帰における順序付け

Question

ロジスティック回帰では、SASには、「降順」オプションを使用して0ではなく1をモデル化するオプションがあります。 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