r - 使用方差分析在 R 语言中复制拉丁方分析

Question

我正在执行一个包含 2 个处理和 2 个块变量（参与者和过程）的 1 个因素（工具）实验。换句话说，一个 2x2 的拉丁方格。

块变量之一是参与者。由于我的样本有 8 个参与者，因此我有 4 个 2x2 拉丁方格。

为了对 1 个拉丁方格进行分析，我通过 R 中的以下命令使用 ANOVA：

result = aov(time~participant+process+tool);
print(summary(result));

我的问题是：如何在 R 中使用复制的拉丁方格执行 ANOVA 测试？

Answer 1

library(nlme) # Avoid lme4 for starter, documentation is not that good yet
result = lme(time~process*tool, data=...., random = ~1|participant)
summary(result)

Answer 2

我找到了一个很好的参考来了解如何分析复制的拉丁方格：

http://halweb.uc3m.es/esp/Personal/personas/jmmarin/esp/Disenno/CursoDisExpSAS.pdf

在第 42 页可以观察到我的情况是：“新行和相同的列”

因此，ANOVA 表应根据第 46 页中的表。特别注意来源和自由度 (df)。

为了在 R 中编码表达的概念，我做了以下工作：

结果=aov（时间~工具+参与者：副本+进程+副本）；

粗体部分是为了表示第 46 页中的表格而进行的更改。