Question

对于喜欢在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 和侧翼之间的相互作用的估计（当 cue 具有级别“center”时）侧翼的级别为“全等”）为-3.9511,95％CI为[-12.997,5.535]

现在说我想使用ezANOVA按副科目和副作品运行anova，我希望得到95％的CI，用于按主题估算。这是我的模特：

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

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

model.f1$ANOVA

而且我不知道如何计算与这些估计相对应的95％CI。我想我应该可以使用ezBoot（），但我尝试了，我不知道如何实现它。

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

Answer 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"

如何使用ezANOVA获得95％的CI（）

1 个答案: