url <- “https://raw.githubusercontent.com/jbryer/CompStats/master/Data/titanic3.csv” df <- read.csv(url, stringsAsFactors = TRUE) head(df)
df <- df[complete.cases(df\(age, df\)sex, df\(pclass, df\)survived), ] df\(survived <- as.factor(df\)survived)
model <- glm(survived ~ age + sex + pclass + sibsp + parch, data = df, family = binomial)
df$predicted_prob <- predict(model, type = “response”)
library(pROC) roc_obj <- roc(df\(survived, df\)predicted_prob) plot(roc_obj, col = “blue”, main = “ROC Curve - Titanic”) auc(roc_obj)