如何根据混合模型结果计算均值差异

编程问答 2022-08-11

问题描述

我正在使用glmmTMB软件包运行混合模型，并使用预测代码使用以下代码从该模型计算预测的均值和置信区间：

运行模型：

mod.1 <- glmmTMB(step.rate ~ Treatment*Week + (1|Animal.ID),data = data.df,family = poisson)

计算预测均值

new.data <- data.frame(Week = data.df$Week,Treatment = data.df$Treatment,Animal.ID = 
            data.df$Animal.ID) 

new.data <- cbind(new.data,predict(mod.1,new.data = new.data,type = "response",re.form = NA,se.fit = TRUE))

但是，我现在想绘制均值和置信区间的差异，并想知道是否存在一种方法来计算预测均值（和置信区间）的差异。

解决方法

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

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

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

glmmtmbpredictrr