有沒有更好的方法在 R 中進行 3 路方差分析？

Question

現在我正在使用我自己的數據集，該數據集與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

因此，我無法獲得他們制造的 plot 或進一步進行分析。

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

我擁有的所有軟件包都是最新的，我的代碼中沒有看到任何錯誤。 我在我的數據集中嘗試了其他變量和相同的問題，當我將其更改為表示時，我只是得到了這條線，而不是框 plot。 如果有人有一些意見，將不勝感激！

Answer 1

它工作得很好，你也沒有在你的代碼中提到這些庫。 如果您不提供所有完整信息，則很難提供幫助。 無論如何，有人想重現這一點。

請運行以下命令，它應該在您的最后工作。 您收到的錯誤表明您可能錯誤地刪除了該列，或者 session 中的數據已損壞。

這是完整的工作代碼：

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