Используйте контраст как меру размера эффекта glmmTMB

Question

Я изучаю эффект экспериментального лечения с переменной ответа Вес (Wt). У меня есть 2 контрольных лечения (A,B) и 2 экспериментальных лечения (C,D). Я хочу посмотреть, увеличивается ли эффект лечения D по сравнению с контрольным лечением с температурой (med_Hobo), и то же самое для лечения C по сравнению с лечением D. (Дополнительные фиксированные эффекты: Atl.Pac, NS и случайный эффект Сайт)

Матрица контрастности:

mat<-rbind(c(-0.5,-0.5,0,1),
           c(0,0,-1,1))
Cmat<-ginv(mat)

Модель:

WW<-glmmTMB(Wt~med_Hobo*Treatment*Atl.Pac*NS+(1|Site)+ (1|Site:Treatment), 
            data=WetWtsEnv,contrasts=list(Treatment=Cmat),family=Gamma(link="log"))

Используя emmeans и emtrends, я могу продемонстрировать, что средний контраст 1 является значимым, а средний контраст 2 - нет. Но я изо всех сил пытаюсь оценить какие-либо изменения в контрасте с температурой. На рисунке ниже, обработка D (синий) явно отличается от температуры, как мне продемонстрировать это цифрами / цифрами? Например, я бы ожидал, что контрасты при med_Hobo=10 будут равны 1 и увеличатся в ~2 раза при med_Hobo=30, с соответствующим увеличением ci.

Прогноз веса на основе модели glmmTMB

Заранее спасибо!