ggplot GLMM二元结果

Question

我正在寻找一种方法来绘制具有二元结果变量的GLMM，其中它们正确或错误地响应（响应）和1个因子，8个级别（级别），每个重复24次，20个人（ID）。

我希望制作一个图表，其中正确响应水平的概率在y轴上，并且每个人在x轴上绘制为包含水平的不同颜色的点。这有点让我陷入困境，但它并不是我想要的，我很确定它不是正确的方法。有20个人，颜色变成了渐变，无法分辨出谁在哪里。任何帮助将不胜感激。

ggplot(df, aes(x = level, y = response, group = id, colour = id)) +
       geom_point() + stat_smooth(method = "glm", family= "binomial", se = FALSE) 

这是一个可重复的小例子，包含3个人和2个级别。

df <- data.frame("id"=c(rep(1,10), rep(2,10), rep(3,10)), 
      "level"=c(1,1,1,1,1,2,2,2,2,2,1,1,1,1,1,2,2,2,2,2,1,1,1,1,1,2,2,2,2,2),
      "response"=c(1,0,0,0,0,1,1,1,1,0,0,0,1,1,1,1,1,0,0,0,1,1,1,0,0,0,1,1,0,1))
       df$level <- factor(df$level)
       df$id <- factor(df$id)