mod1 <- glm(outcome == "Alive" ~ smoker + age, data = Whickham)
mod2 <- glm(outcome == "Alive" ~ smoker * age, data = Whickham)
mod3 <- glm(outcome == "Alive" ~ smoker * ns(age,3), data = Whickham)
cv_pred_error(mod1, mod2, mod3, output = "likelihood")