2014-10-02 120 views
3

是否可以在r中运行方差分析,仅使用均值,标准差和n值?这里是我的数据框:R中的方差分析使用汇总数据

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) 

我想找到平均方差残差,所以我想看看方差分析表。

任何帮助将非常感谢! :)

+2

这是第一个谷歌结果为“平均值和标准差的R anova in”:http://r.789695.n4.nabble.com/Perform-one-way-ANOVA-using-standard-deviation-and -mean-td3265548.html – thelatemail 2014-10-02 22:25:53

+0

我在网上也找不到任何东西,这就是我发布这个问题的原因。 :) – y3trgfhsfgr 2014-10-02 23:11:23

+1

另外,在谷歌搜索这个问题的确切标题给出了一个链接到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 2014-10-02 23:28:48

回答

4

在这里你去,用ind.oneway.secondrspychi包:

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 

到无关的方面说明,您的数据可能会与一些做重命名。 q2data.framedata.frame,不需要放在标题中。此外,无需在q2data.frame内指定q2data.mean - 当然mean就足够了。它只是意味着你最终像复杂的代码:

q2data.frame$q2data.mean 

时:

q2$mean 

会给你你需要的所有信息。

+0

你是天赐之物,谢谢!此外:对于那些没有rpsychi库的人,在上面的脚本之前包含以下代码:install.packages(“rpsychi”) install.packages(“gtools”) – y3trgfhsfgr 2014-10-03 00:16:38

+1

@BrettCochrane - 'install.packages'应该处理获取所有依赖关系 - 无需单独安装gtools。 – thelatemail 2014-10-03 01:12:55

+0

+1新信息给我 – akrun 2014-10-03 03:10:30