理学療法士がまとめたノート

統計学備忘録 since2016

Rを使って統計学を勉強するブログです

Rで簡単 多重ロジスティック回帰分析

サンプルはRのmtcarsを使用します.車に詳しくないので回帰分析に相応しいサンプルか分かりませんが、学習のため数値のみ使用させていただきます.

サンプルから直接取り込み

dat <- mtcars[,c(8,2,5,6)]    #ここがポイント8列目を応答変数として、2列目・5列目・6列目を説明変数として設定
ans <- glm(dat$vs~.,data=dat,family=binomial) 
ans <- glm(dat$vs~dat$cyl+dat$drat+dat$wt,family=binomial)  

これで終了です.

summary(ans)

Call:
glm(formula = dat$vs ~ dat$cyl + dat$drat + dat$wt, family = binomial)

Deviance Residuals: 
     Min        1Q    Median        3Q       Max  
-2.43462  -0.15612  -0.06808   0.22595   1.35668  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)  
(Intercept)   22.302     11.577   1.926   0.0541 .
dat$cyl       -3.375      1.544  -2.186   0.0288 *
dat$drat      -2.189      1.874  -1.168   0.2428  
dat$wt         1.777      1.660   1.070   0.2845  
---
Signif. codes:  0***0.001**0.01*0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 43.86  on 31  degrees of freedom
Residual deviance: 13.92  on 28  degrees of freedom
AIC: 21.92

Number of Fisher Scoring iterations: 7

オッズ比

exp(ans$coefficients)

 (Intercept)      dat$cyl     dat$drat       dat$wt 
4.849125e+09 3.423400e-02 1.120633e-01 5.911782e+00