0

现在我正在使用我自己的数据集，该数据集与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

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

4

1 回答 1

2

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

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

这是完整的工作代码：

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
于 2021-03-20T06:42:57.393 回答