无法使用lsmeans获取glmer的调整平均值

我有一个glm，我想得到使用lsmeans的调整平均值。下面的代码创建了模型，并且似乎做得很正确：

library(lmerTest)
data$group <- as.factor(data$grp)
data$site <- as.factor(data$site)
data$stimulus <- as.factor(data$stimulus)

data.acc1 = glmer(accuracy ~ site + grp*stimulus + (1|ID), data=data, family=binomial)

或

---

有什么建议吗？

问题在于，您已经使用glmer创建了一个广义线性混合模型。在这种情况下，是一个混合逻辑回归模型，而不是使用lmer的线性混合模型。lsmeans函数不接受glmer创建的对象，因为它们不是线性混合模型

这篇文章中的答案可能会有所帮助：

---

如果您想了解/计算混合GLM的边际效应，这篇文章可能会很有用：

错误消息的可能副本来自lmerTest包。试试lsmeans包。

lsmeans(data.acc1, "stimulus")

data.lsm <- lsmeans(data.acc1, accuracy ~ stimulus ~ grp)
pairs(data.lsm)