有没有更好的方法在 R 中进行三向方差分析？

如何解决有没有更好的方法在 R 中进行三向方差分析？

现在我正在使用我自己的数据集，这与 https://www.datanovia.com/en/lessons/anova-in-r/ 中的示例非常相似，特别是三向方差分析部分。代码布局合理，但是当我进入使用箱线图进行可视化的部分时，遇到了意外错误。

headache %>%
  group_by(gender,risk,treatment) %>%
  get_summary_stats(pain_score,type = "mean_sd")

在他们的代码中，他们将 y 指定为疼痛评分，但对我来说，我得到了错误：

FUN (x[[i]],...) : object "pain_score" not found

因此，我无法理解他们制作的情节或进一步进行分析。

  bxp <- ggBoxplot(
  headache,x = "treatment",y = "pain_score",color = "risk",palette = "jco",facet.by = "gender"
  )
bxp

我拥有的所有软件包都是最新的，我的代码中没有看到任何错误。我已经在我的数据集中尝试了其他变量和同样的问题，当我将其更改为平均值时，我只得到了这条线，而不是箱线图。如果有人提供一些意见，将不胜感激！

解决方法

它工作得很好，你也没有在你的代码中提到库。如果您不提供所有完整信息，则很难提供帮助。无论如何，有人想重现这个。

请运行下面的，它应该在你的最后工作。您收到的错误表明您可能错误地删除了该列，或者会话中的数据已损坏。

这是完整的工作代码：

library(tidyverse)
library(rstatix)
library(ggpubr)

data("headache",package = "datarium")

headache %>%
  group_by(gender,risk,treatment) %>%
  get_summary_stats(pain_score,type = "mean_sd")


bxp <- ggboxplot(
  headache,x = "treatment",y = "pain_score",color = "risk",palette = "jco",facet.by = "gender"
)
bxp