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


10

У мене виникають проблеми з визначенням помилок для двосторонніх повторних заходів ANOVA в Р. Мої дані складаються з оцінок щільності деревини для трьох радіальних положень (внутрішнього, середнього та зовнішнього) уздовж ядра, витягнутого з дерева. Всього існує 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
amoeba
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.