我正在研究响应变量为体重 (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,并在 med_Hobo=30 时增加到约 2 倍,ci 相应增加
提前致谢!