3 way anova inested in r

Question

Я пытаюсь выяснить модель для полностью факториального эксперимента.

У меня следующие факторы

Treatment Day Hour Subject ResponseVariable 

10 дней измерений, 4 различные моменты времени в течение каждого дня, 2 различные виды лечения, измеренные, 12 субъектов) 6 субъектов в лечении 1 и 6 различных предметов в лечении 2)

за каждый день, который я измерил: 6 субъектов в лечении 1, остальные 6 в лечении 2, в 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) 

любые мысли бы помочь?

Answer 1

У меня такая же ошибка, и я думаю, что в моей проблеме отсутствовали наблюдения. Вам не хватает каких-либо наблюдений? Я считаю, что они не являются проблемой для линейных смешанных эффектов, и я читал, что некоторые люди используют lme вместо повторных измерений ANOVA для этих случаев.

Ваш ошибочный термин может быть интерпретирован как «S-эффект внутри каждого T». Это звучит из вашего описания, как будто это то, что вы хотите, поэтому я не думаю, что это вызывает ваше сообщение об ошибке.

Одно замечание: я вижу, у вас есть переменная с именем «T». R позволит вам это сделать? T обычно зарезервирован для значения «ИСТИНА». Это может быть частью вашей проблемы.