如何使用 ezANOVA() 获得 95% CI答案

【问题标题】：How to get 95% CIs using ezANOVA()如何使用 ezANOVA() 获得 95% CI
【发布时间】：2013-07-30 17:14:05
【问题描述】：

这是一个编程问题，适合喜欢在 R 中使用 ez 包的人。我习惯于使用带有 lmer() 的线性混合效果模型。在 lmer () 的有用输出中，我得到了每个实验因子的系数值，并且使用 pvals.fnc() 我可以轻松获得 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

如果我使用 lmer，我会估计我的模型，然后得到系数的 95% CI：

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

因此，例如，我知道 cue 和 flank 之间的交互作用的估计（当 cue 具有“中心”级别时） flank 的级别为“全等”）为 -3.9511，95% CI 为 [-12.997, 5.535]

现在假设我想使用 ezANOVA 运行 anova by-subjects 和 by-items，并且我想获得 95% CI 的 by-subject 估计值。这是我的模型：

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

但在输出中，我没有看到模型估计值：

model.f1$ANOVA

而且我不知道如何计算与这些估计值对应的 95% CI。我想我应该可以使用 ezBoot() 但我试过了，但我不确定如何实现它。

有什么建议吗？感谢您的帮助！

【问题讨论】：

标签： r anova confidence-interval

【解决方案1】：

这个答案是由另一个论坛的“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"

【讨论】：