Rでタイタニックの生存者分析

Program R

z <- data.frame(Titanic)

Titanic1 <- data.frame(Class = rep(z[,1],z[,5]), Sex = rep(z[,2], z[,5]), Age = rep(z[,3], z[,5]), Survived = rep(z[,4],z[,5]))

 

########補足：データの型の修正はこちらの方がスマート ##########

class(Titanic)

TN1 <- data.frame(Titanic)

TN2 <- data.frame(lapply(TN1, function(i) rep(i, TN1[,5])))

TN3 <- TN2[,-5]

#検算

TN4 <- xtabs(~Class+Sex+Age+Survived,data=TN3)

TN4 == Titanic

########以上補足##########

(r <- glm(Survived ~ Class, data = TN3, family = "binomial"))

table(r$fitted.values)

(q <- glm(Survived ~ ., data = TN3, family = "binomial"))

summary(q)

cbind(TN1, SurveRate = predict(q, Titanic, type = "response"))

exp(q$coefficient)

#exp(coefficient)で各変数についてのオッズ比を算出出来る。