首页 文章

3路anova嵌套在r

提问于
浏览
0

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

我有以下因素

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 回答

  • 1

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

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

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

相关问题