如何查看lsmeans的所有输出？

Question

我使用lme4包运行了线性混合模型。然后，要查看所有对比，我必须使用lsmeans包运行lsmeans函数。这是我运行的功能：

library('lme4')
library('lsmeans')
lsmeans(lmer1, pairwise ~ vowel * experience * lang_sound, adjust="tukey")

然而，我看不到整个输出，因为有太多对。任何人都可以告诉我应该怎么做才能从这个函数中获得输出？

我尝试了'sink()'，但它不起作用。也许是因为我没有运行lsmeans命令的名称。我在Windows上使用RStudio。

Answer 1

sink函数执行，但您必须提供文件名以保存输出。例如

sink(file = "lsm-output.txt")
lsmeans(...)
sink()

最后一次sink()调用会将输出恢复为控制台。

那就是说，我认为你不需要看到所有那些成对的比较。我建议在单独的调用中进行lsmeans和比较。您可以使用带有pairs参数的by来有条件地对其他两个因素的级别进行比较：

library("lsmeans")
lmer1.lsm <- lsmeans(lmer1, ~ vowel * experience * lang_sound)
lmer1.lsm     # display the means
pairs(lmer1.lsm, by = c("experience", "lang_sound"))
pairs(lmer1.lsm, by = c("vowel", "lang_sound"))
pairs(lmer1.lsm, by = c("vowel", "experience"))

我还认为你应该想象你的结果，例如：通过构建交互图：

lsmip(lmer1.lsm, experience ~ vowel | lang_sound)

在这次电话会议中，一些因素的互换可能会更令人满意。