r - 如何查看 Ismeans 的所有输出？

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)

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