2014-10-02 54 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