Angabe des Termes Error () in ANOVA mit wiederholten Messungen in R.

Ich habe Probleme beim Definieren der Fehlerterme für eine ANOVA mit zwei wiederholten Messungen in R. Meine Daten bestehen aus Holzdichteschätzungen für drei radiale Positionen (innen, mittel und außen) entlang eines aus einem Baum extrahierten Kerns. Es gibt insgesamt 20 Baumarten, 6 Individuen jeder Art und zwei Kerne von jedem Baum.

Um die Auswirkung der radialen Position auf die Holzdichte zu testen, verwende ich das folgende Zwei-Wege-ANOVA-Modell mit einem Fehlerterm, der die Variabilität zwischen Individuen berücksichtigt:

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

Ich bin mir jedoch nicht sicher, ob meine Angabe des Fehlerterms angemessen ist. Soll ich auch die Variabilität innerhalb eines Kerns berücksichtigen? Für mich ist diese Variabilität dieselbe wie die aufgrund der radialen Position, die der Hauptfaktor ist, an dem ich interessiert bin.

Obwohl ich einige Zeit damit verbracht habe, über die Angabe des Fehlerausdrucks in ANOVA mit wiederholten Messungen zu lesen, habe ich immer noch Probleme mit der praktischen Angabe des Fehlerausdrucks. Ich werde mich über Hilfe dabei freuen.

Es wäre

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

Das erklärt den innerhalb des Subjekts liegenden Fehler von Radialposition. Wenn Sie andere subjektinterne Faktoren haben, werfen Sie diese (in einer Interaktion) Radialpositionin den Fehler-Nenner, wie z + Error(Individual/(Radialpostion*wiFactorA)). So verstehe ich das. Dies entspricht den wiederholten Messungen von SPSS, GLMwenn keine Daten fehlen.