prediction_data <- explanatory_data %>% mutate( has_churned = predict(mdl_churn_vs_relationship, explanatory_data, type = "response"), most_likely_outcome = round(has_churned) )