如何使用emmeans软件包绘制均值和置信区间的差异？

问题描述

我正在运行一个包含2向交互的混合模型。模型看起来像这样

model.1 <- glmmTMB(Weight ~ Treatment*Week + (1|Animal.ID),data=df,family=nbinom1) 

治疗是一个因素，有2个级别（感染和未感染），而一周是一个因素，有9个级别（9周）。

模型的输出表明，治疗与3-7周之间的星期之间存在显着的交互作用。

但是，当我绘制预测的均值和置信区间时，条形图重叠并且该图看起来不太有说服力。

我想知道是否有更好的方法来展示此模型的结果。我最初的想法是使用emmeans软件包绘制均值和置信区间的差异，但我不知道该怎么做。

我尝试使用以下代码（均不正确）计算均值差，但其值与预测的均值相同，并且两个图看起来相同。这是正确的吗？

plot.1 <- emmip(model.1,Treatment~Week,type = "response",CIs = TRUE,plotit = FALSE) 

我想知道如何使用emmeans软件包或其他软件包（如果可能）计算均值和置信区间的差异。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）