data(ptitanic)
summary(ptitanic)
# survival rate was greater for females
rpart.rules(rpart(survived ~ sex, data = ptitanic))
# survival rate was greater for higher classes
rpart.rules(rpart(survived ~ pclass, data = ptitanic))
# survival rate was greater for children
rpart.rules(rpart(survived ~ age, data = ptitanic))
# main indicator of missing data is 3rd class esp. with many children
obs.with.nas <- rowSums(is.na(ptitanic)) > 0
rpart.rules(rpart(obs.with.nas ~ ., data = ptitanic, method = "class"))
Run the code above in your browser using DataLab