2

我在 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 执行事后测试?

4

0 回答 0