使用R进行ANOVA分析的总结数据

Question

是否可能使用R中的平均值、标准差和n值进行方差分析？以下是我的数据框：

q2data.mean <- c(90,85,92,100,102,106)
q2data.sd <- c(9.035613,11.479667,9.760268,7.662572,9.830258,9.111457)
q2data.n <- c(9,9,9,9,9,9)
q2data.frame <- data.frame(q2data.mean,q2data.sq,q2data.n)

我正在尝试寻找均方残差，因此我想查看ANOVA表。

非常感谢您的帮助！ :)

- y3trgfhsfgr

---

此外，通过谷歌搜索您的问题确切标题，会给出一个链接：http://stats.stackexchange.com/questions/95949/one-way-anova-from-summary-data-in-r，该链接指向：http://www.inside-r.org/packages/cran/rpsychi/docs/ind.oneway.second。 - thelatemail

Answer 1

使用 rspychi 包中的 ind.oneway.second：

library(rpsychi)
with(q2data.frame, ind.oneway.second(q2data.mean,q2data.sd,q2data.n) )

#$anova.table
#                SS df     MS     F
#Between (A) 2923.5  5 584.70 6.413
#Within      4376.4 48  91.18      
#Total       7299.9 53   
# etc etc

更新：rpsychi包已于2022年3月归档，但该函数仍可在此处使用：http://github.com/cran/rpsychi/blob/master/R/ind.oneway.second.R（感谢评论中的@jrcalabrese提供提示）

---

另外，与此无关的一点是，您的数据可能需要进行一些重命名。 q2data.frame是一个data.frame，无需将其放入标题中。此外，在q2data.frame内部指定q2data.mean是没有必要的，只需写mean即可。这将使您的代码变得复杂：

q2data.frame$q2data.mean

什么时候：

q2$mean

would give you all the info you need.

- thelatemail

---

你真是太棒了，谢谢！另外，对于那些没有包含rpsychi库的人，请在上述脚本之前添加以下代码：install.packages("rpsychi")install.packages("gtools") - y3trgfhsfgr

@BrettCochrane - install.packages 应该处理获取所有依赖项 - 不需要单独安装 gtools。 - thelatemail

不幸的是，rpsychi包在2022年3月被归档，但该函数仍然可在此处使用：https://github.com/cran/rpsychi/blob/master/R/ind.oneway.second.R - jrcalabrese