问题描述
现在我正在使用我自己的数据集,这与 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