Задание члена Error () в повторных измерениях ANOVA в R


10

У меня проблемы с определением условий ошибки для двухсторонних повторяющихся измерений ANOVA в R. Мои данные состоят из оценок плотности древесины для трех радиальных положений (внутреннего, среднего и внешнего) вдоль керна, извлеченного из дерева. Всего существует 20 видов деревьев, 6 особей каждого вида и два ядра от каждого дерева.

Чтобы проверить влияние радиального положения на плотность древесины, я использую следующую двухстороннюю модель ANOVA с погрешностью, которая учитывает изменчивость среди людей:

radpos.aov <- aov(WD ~ Species*Radialposition + Error(Individual), data=Radpos)

Однако я не уверен, что моя спецификация термина ошибки адекватна. Должен ли я также учитывать изменчивость в ядре? Для меня эта изменчивость такая же, как и из-за радиального положения, которое является основным фактором, который меня интересует.

Хотя я посвятил некоторое время чтению об указании термина ошибки в ANOVA с повторными измерениями, у меня все еще есть проблемы с практическим определением термина ошибки. Я буду признателен за помощь с этим.

Ответы:


8

Это было бы

 radpos.aov <- aov(WD ~ Species*Radialposition + Error(Individual/(Radialposition)), data=Radpos)
 summary(radpos.aov, type=3)

Что составляет внутри предметную ошибку Radialposition. Если у вас есть другие внутри-субъектные факторы, добавьте их (во взаимодействии) со Radialpositionзнаменателем ошибок, например + Error(Individual/(Radialpostion*wiFactorA)). Это мое понимание этого. Это соответствует повторным измерениям SPSS, GLMесли у вас нет недостающих данных.


1
+1. Смотрите подробную шпаргалку здесь: dwoll.de/rexrepos/posts/anovaMixed.html . А также здесь: personal-project.org/r/r.guide/r.anova.html
амеба
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.