Я подгоняю свои данные к общей линейной смешанной модели с фиксированным эффектом Обработка и случайным эффектом Сцепление.Вот мой код:
model<-glmer(cbind(Successes,Failures)~Treatment+(1|Clutch),
data = cont, family = "binomial")
Моя работа связана с соотношениями полов, и я определяю женщину как успех, а мужчину - как неудачу в каждом наблюдении.У меня есть 4 разных лечения.Я хочу построить (предпочтительно используя ggplot) прогнозируемое соотношение полов из модели для каждой процедуры с учетом сцепления (с 95% доверительными интервалами).Я понимаю, что это, вероятно, большой вопрос, но может ли кто-нибудь помочь мне с кодом, который мне понадобится для этого?Я искал онлайн последние несколько дней.Спасибо!