2012-06-27 273 views
1

我正在执行2个处理和2个块变量(参与者和处理)的1因素(工具)实验。换句话说,一个2x2拉丁方格。使用ANOVA在R语言中复制拉丁方分析

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

做的分析只是1拉丁方,我通过在R中的以下命令使用ANOVA:

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

我的问题是:我如何执行与R中复制的拉丁方ANOVA检验?

+0

看那包'lme4'和'nlme'在页面呈现表中的变化。 – mnel

回答

0

我找到了一个很好的参考,以了解如何分析复制拉丁方:

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

在42页是可以观察到我的情况是:“新行和同列“

因此,方差分析表应按照第46页的表格。特别注意来源和自由度(df)。

为了译码R中的表达的概念,我做了以下:

结果= AOV(时间〜工具+ 参与者:replic + +处理replic);

中的加粗部分是为了表示46

+1

严格地说,这没问题,但我认为拉丁广场是一个重复随机化的特例,我对该专栏不感兴趣。所以我更喜欢把重复放入重复中(见下文)。 –

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