如何使用ezANOVA獲得95％CI（）

這是一個編程問題，適用於喜歡在R中使用ez包的人。我習慣於在lmer（）中使用線性混合效果模型。在lmer（）的有用輸出中，我得到了每個實驗因素的係數值，並且使用pvals.fnc（）我可以很容易地得到95％的置信區間（CI）與模型係數一起報告。如何使用ezANOVA獲得95％CI（）

我最近開始使用ezANOVA，我想知道：有沒有一種主流的方式來獲得相同的輸出？也就是說，我想獲得一個實驗因子的係數值和一個CI值。下面是示例代碼，以使該混凝土：

library(languageR) #necessary to use pvals.fnc() 
library(lme4)  #necessary for lmer() 
library(ez)  #necessary for ezANOVA 
data(ANT)   #load sample data

如果我使用11聚物，我估計我的模型，然後得到的95％CI爲係數：

model_lmer = lmer(formula = rt ~ cue*flank + (1|subnum), data = ANT) 
pvals.fnc(model_lmer, withMCMC=T)$fixed

因此，舉例來說，我知道的互動提示和之間側翼的估計值（當線索具有電平「中心」和側翼具有電平「全等」）是-3.9511並且95％CI是[-12.997， 5.535 ]

現在說我想用ezANOVA來運行anova by-subjects和by-items，並且我想獲得95％的可信區間來進行按主題的估計。這是我的模型：

model.f1 = ezANOVA(data=ANT, dv=rt,wid=subnum,within=.(cue,flank),return_aov=T)

但在輸出中，我沒有看到模型估算，當我做：

model.f1$ANOVA

而且我不知道如何將95％可信區間對應的計算到那些估計。我想我應該能夠使用ezBoot（），但我試過了，我不知道如何實現它。

有什麼建議嗎？謝謝你的幫助！

來源

2013-07-30 Sol

此答案由其他論壇上的「ez」軟件包作者提供。我在這裏複製它以防別人認爲它有用：

"One somewhat hacky way to get CIs for effects is to use ezStats() to get the means 
and FLSD, compute the difference between the means to get the effect, 
and divide the FLSD by sqrt(2) to get the CI"

來源

2013-08-01 17:21:11 Sol

如何使用ezANOVA獲得95％CI（）

回答

相關問題