Of course, probabilities greater 1 do not make sense. That’s the reason why we prefer a “bended” graph, such as the s-type ogive in logistic regression. Let’s plot that instead.

First, we need to get the survival probabilities:

`d %>% mutate(pred_prob = predict(glm1, type = "response")) -> d`

Notice that

`type = "response`

gives you the probabilities of survival (ie., of the modeled event).

Read the whole thing.