2010-04-26 52 views
2

我正在嘗試計算2x2析因設計所需的樣本量。我有兩個問題。雙向Anova功耗分析

1)我使用的包PWR和單向ANOVA功能用下面的代碼

pwr.anova.test(k = , n = , f = , sig.level = , power =) 

但是計算所需的樣本大小,我想看兩個因素ANOVA,因爲這是更有效地估計羣體意味着比單向anova更有效。我找不到雙向anova函數。 [R]中是否有包裝或常規?

2)此外,我是否安全地假設,因爲我使用單向anova權力計算,樣本量會更保守(即更大)?

+3

這不是編程有關,但你應該何去何從羅布海德門的網站提議投票:http://meta.stackexchange.com/questions/5547/proposal-for-統計的數據挖掘,stackexchange現場。 – Shane 2010-04-26 01:16:20

回答

0

搜索後 - 我無法找到這個網上的任何解決方案。

我建議你做什麼(如果你知道怎麼做的話)就是使用模擬來編程。 如果你不知道該怎麼做,那麼寫一個關於「我怎樣編寫雙向anova仿真,實現功耗分析」的SO問題,看看有什麼人可以幫助你:)

此外,您可以通過這裏檢查代碼開始:

http://www.rforge.net/doc/packages/NCStats/power.sim.html

對於通過模擬在功率計算一個開始。

請注意Jeromy寫的 - 這種功效分析適用於多種結果。

有趣的主題 - 我很想跟進它。

最佳,

塔爾

+0

謝謝你們兩位的輸入傑里米和塔爾。我會考慮你的建議以及如何融入它們。 – Thomas 2010-04-27 12:07:07

+0

@Thomas你是否爲你的問題開了另一篇文章? R中只需要一條線來計算F檢驗的功率。這不是一個好的地方,但是這裏有一個功能可以完成這項工作:功率< - 函數(alpha,eff,n,m,l)df1 < - m -l df2 < - n - m c < - qf (1-alpha,df1,df2) lambda-eff^2 * n pow <-pf(c,df1,df2,ncp = lambda,lower.tail = FALSE) return(pow) } – 2013-07-16 20:43:28