The learner status should be removed first as it provides the greatest adjusted \(R^2\) without it being included in the model.

p <- c(51, 53, 55)

exp(11.6630-0.2162*p) / (1 + exp(11.6630-0.2162*p))
## [1] 0.6540297 0.5509228 0.4432456
data <- data.frame(prob=c(0.341, 0.251, 0.179, 0.124, 0.084, 0.056, 0.037, 0.024, 0.654, 0.551, 0.443), temp=c(57, 59, 61, 63, 65, 67, 69, 71, 51, 53, 55))

ggplot(data, aes(x=temp, y=prob)) +
  geom_point() +
  geom_smooth()
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'