我正在执行一个 1 因子(工具)实验,其中包含 2 个处理和 2 个块变量(参与者和过程)。换句话说,一个 2x2 的拉丁方块。
块变量之一是参与者。由于我的样本有 8 名参与者,所以我有 4 个 2x2 拉丁方块。
为了仅对 1 个拉丁方块进行分析,我通过 R 中的以下命令使用 ANOVA:
result = aov(time~participant+process+tool);
print(summary(result));
我的问题是:如何在 R 中使用复制的拉丁方块执行方差分析测试?
library(nlme) # Avoid lme4 for starter, documentation is not that good yet
result = lme(time~process*tool, data=...., random = ~1|participant)
summary(result)
我找到了一个很好的参考来了解如何分析复制的拉丁方块:
http://halweb.uc3m.es/esp/Personal/personas/jmmarin/esp/Disenno/CursoDisExpSAS.pdf
在第 42 页可以观察到我的情况是:"新行和相同的列"
因此,方差分析表应根据第 46 页中的表。特别注意自由度(df)的来源和自由度。
为了在 R 中对表达的概念进行编码,我做了以下操作:
结果 = AOV(时间~工具+参与者:重复+过程+重复);
粗体部分是更改,以表示第 46 页中提供的表格。