Określanie terminu Błąd () w ANOVA w powtarzanych pomiarach w R.

Mam problem z określeniem terminów błędów dla ANOVA z powtarzanymi dwukierunkowymi pomiarami w R. Moje dane obejmują szacunki gęstości drewna dla trzech położeń promieniowych (wewnętrzna, środkowa i zewnętrzna) wzdłuż rdzenia wydobytego z drzewa. Istnieje w sumie 20 gatunków drzew, 6 osobników każdego gatunku i dwa rdzenie z każdego drzewa.

Aby przetestować wpływ położenia promieniowego na gęstość drewna, używam następującego dwukierunkowego modelu ANOVA ze składnikiem błędu, który uwzględnia zmienność między osobnikami:

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

Nie jestem jednak pewien, czy moja specyfikacja terminu błędu jest odpowiednia. Czy mam także uwzględniać zmienność w rdzeniu? Dla mnie ta zmienność jest taka sama jak ze względu na położenie promieniowe, które jest głównym czynnikiem, który mnie interesuje.

Chociaż poświęciłem trochę czasu na przeczytanie o określaniu terminu błędu w ANOVA z powtarzanymi pomiarami, wciąż mam problemy z praktycznym określeniem terminu błędu. Doceniam pomoc w tym zakresie.

To byłby

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

Odpowiada to za błąd wewnętrzny w wysokości Radialposition. Jeśli masz inne czynniki wewnątrz podmiotu, wrzuć je (w interakcji) z Radialpositionmianownikiem błędu, np + Error(Individual/(Radialpostion*wiFactorA)). Tak rozumiem. To pasuje do powtarzanych pomiarów SPSS, GLMjeśli nie brakuje danych.