我在 R 中定义了一个 lmer 模型,它具有 2 个固定的相互作用效应和三个随机效应。第一个固定效应“A”是分类的,而第二个固定效应“B”是连续的:
library(lme4)
model1<- lmer(log(response)~A*B+(1|random1)+(1|random2)+(1|random3), data=Data)
我曾尝试使用 lsmeans(现在称为 emmeans):
lsmeans(model1, pairwise~A*B, adjust="tukey")
我也尝试过使用 glht:
summary(glht(model1, linfct = mcp(A*B ="Tukey")))
但这些选项都不起作用。
我是R新手。如何对我的 LMM 执行事后测试?