comparison - 使用对比度作为效果大小的度量 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，并在 med_Hobo=30 时增加到约 2 倍，ci 相应增加

基于 glmmTMB 模型的权重预测

提前致谢！

标签： comparisonmixed-modelscompare-contrast