我有一个 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)
但是,当我尝试使用以下任何代码来调整模型的方法时,我会收到错误消息
lsmeansLT(model, test.effs = test.effs, ddf = ddf) 中的错误:
该模型不是线性混合效应模型。
lsmeans(data.acc1, "stimulus")
或者
data.lsm <- lsmeans(data.acc1, accuracy ~ stimulus ~ grp)
pairs(data.lsm)
有什么建议吗?