如何查看 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 和比较。您可以通过使用带有 by 参数的 pairs 有条件地比较一个因素的水平与其他两个因素的水平：

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)

本次通话中一些因素的交换可能会更令人满意。