![[PukiWiki]](image/stat_logo.jpg "[PukiWiki]"){kind=logo}
ロジスティック回帰分析のサンプルデータ
タイタニック号
統計パッケージRの中にタイタニック号での生存者・非生存者のデータがあります。
データの内容
- Class（乗客のクラス／乗組員） 1st, 2nd, 3rd, Crew
- Sex（性別） Male, Female
- Age（成人／子供） Child, Adult
- Survived（生存） No, Yes
CSVデータにしたものを添付しておきます。
Titanic.csv低体重出生とリスク因子
library(MASS) data(birthwt)
- low ：低体重出生の有無を示す２値変数（児の出生時体重 2.5 kg 未満が 1 ）
- age ：年齢
- lwt ：最終月経時体重
- race ：人種（１＝白人，２＝黒人，３＝その他）
- smoke ：喫煙の有無（１＝あり）
- ptl ：非熟練労働経験数，
- ht ：高血圧の既往（１＝あり）
- ui ：子宮神経過敏の有無（１＝あり）
- ftv ：妊娠の最初の３ヶ月の受診回数
- bwt ：児の出生児体重 (kg)
library(MASS)
data(birthwt)
#
attach(birthwt)
low<-factor(low)
race<-factor(race, labels=c("white","black" ,"other"))
ptd<-factor(ptl>0)
smoke<-(smoke>0)
ht<-(ht>0)
ui<-(ui>0)
ftv<-factor(ftv)
levels(ftv)[-(1:2)]<-"2+"
# 必要な変数だけのデータフレーム bwt を定義する
bwt<-data.frame(low,age,lwt,race,smoke,ptd,ht,ui,ftv)
detach(birthwt)
# データの確認
str(bwt)
summary(res<-glm(low ~ ., family=binomial, data=bwt))
# 変数選択するなら
summary(res2<-step(res))