我正在做一个三元方差分析,我想对它进行事后分析。
为了按组进行比较,我发现emmeans包中的函数emmeans_test
是最佳选项。但是当我运行代码时,R找不到emmeans_test
。我不明白为什么。我已经安装了emmeans并完成了library(emmeans)
。
library(emmeans)
donnees_tot_g_J4 %>%
group_by(quality, temperature) %>%
emmeans_test(growth_rate ~ quantity, p.adjust.method = "bonferroni")
Error in emmeans_test(., growth_rate ~ quantity, p.adjust.method = "bonferroni"):
could not find function "emmeans_test"
函数emmeans_test
不是来自emmeans,而是来自rstatix。所以你需要:
if (!require(rstatix)) install.packages("rstatix")
library(rstatix)
谢谢你,就是这样!抱歉问了这个愚蠢的问题。
别担心。这个名字确实有误导性。事实上,这两个方案都不依赖或进口另一个rstatix仅表明它可以增强emmeans。
感谢dipetkov提供的有用反馈。
实际上,rstatix调用emmeans来进行实际分析;它并没有增强任何东西。例如,它不能处理省略三方交互的模型。
很抱歉造成混淆。我的声明是基于我在rstatix的CRAN页面上看到的内容。令人惊讶的是,emmeans既不在";取决于";也不是";"进口";;仅在";建议";。不知道如果包一开始就不导入emmeans,emmeans_test
将如何工作。
我强烈建议OP学习如何用
lm()
和emmeans()
进行分析,因为他们会有更多的灵活性(以及对结果的信心(。
谢谢。我更新了我的答案,将这个建议传递给OP,以及未来的读者。