我将数据拟合为一般线性混合模型，其中“治疗”为固定效应，而“离合器”为随机效应。这是我的代码：
model<-glmer(cbind(Successes,Failures)~Treatment+(1|Clutch),
data = cont, family = "binomial")
我的工作涉及性别比例，每次观察我都将女性定义为成功，将男性定义为失败。我有4种不同的治疗方法。我想绘制模型（最好使用ggplot），并考虑离合器（95％置信区间）下每种治疗模型的预测性别比。我意识到这可能是一个很大的问题，但是有人可以帮助我提供执行此操作所需的代码吗？我过去几天一直在网上搜索。谢谢！