Rでロジスティック回帰分析

[:ja]library(MASS)
#「低体重出生とそのリスク因子の関連」を調べたデータ
head(birthwt)
str(birthwt)
#列の選択
sample_logi<-birthwt[,c(1,2,3,5,6,7,8)]
#ポンドをキログラムに変換
sample_logi$lwt<-sample_logi$lwt * 0.454
#データフレームの項目を直接使えるようにする。{sample_logi$lwt}→{lwt}
attach(sample_logi)

train.glm <- glm(order_of_finish~., family=binomial,data= train)
summary(output.glm)

対数odds{log(2.5kg未満の新生児が生まれる確率/生まれない確率)=
1.38186+(-0.04223)*年齢+(-0.03154)*体重+0.55076*喫煙有無+0.59316*早産経験有無+1.86364*高血圧有無+0.73675*子宮神経過敏有無}

#対数oddsから確率へ変換 自然対数表示exp
確率=exp(対数odds)/(1+exp(対数odds))

# 高血圧以外の条件が同じ場合に、要素の有無が、確率を~倍引き上げる
exp(output.glm$coefficients)[:]

コメントを残す

メールアドレスが公開されることはありません。 * が付いている欄は必須項目です