3路anova嵌套在r

时间:2014-01-17 00:33:39

标签: r lm anova

我正在试图找出完全因子实验的模型。

我有以下因素

Treatment   Day   Hour  Subject   ResponseVariable

10天的测量,每天4个不同的时间点,测量2个不同的治疗,12个受试者)治疗1中的6个受试者和治疗中的6个不同的受试者2)

我测量的每一天:治疗1的6名受试者,治疗2的其他6名,4个不同的时间点 对于受试者,我有12个不同的受试者,但受试者1-6在治疗-1中,受试者7-12在治疗-2中。受试者没有改变治疗方法,因此我在10天的每一次治疗中测量了相同的一组受试者

所以绊​​倒我的是指定正确的错误术语 我以为我有一般模型,但R给了我"Error() model is singular"

aov(ResponseVariable  ~ T + R + S + TR + TS + RS + Error(T/S)

任何想法都有帮助?

1 个答案:

答案 0 :(得分:1)

我遇到了同样的错误,我认为我的问题是缺少观察结果。你错过了任何观察吗?我相信它们对于线性混合效应来说不是一个问题,而且我已经读过有些人使用lme而不是重复测量方差分析。

您的错误术语可以解释为“每个T中的S效应”。这听起来好像是你想要的,所以我不认为这是导致你的错误信息的原因。

一个注意:我看到你有一个名为“T”的变量。 R让你这样做? T通常保留为“TRUE”。这可能是你问题的一部分。