Spécification du terme Error () dans les mesures répétées ANOVA dans R

J'ai des problèmes avec la définition des termes d'erreur pour une mesure répétée bidirectionnelle ANOVA dans R. Mes données consistent en des estimations de la densité du bois pour trois positions radiales (intérieure, moyenne et extérieure) le long d'un noyau extrait d'un arbre. Il y a un total de 20 espèces d'arbres, 6 individus de chaque espèce et deux noyaux de chaque arbre.

Pour tester l'effet de la position radiale sur la densité du bois, j'utilise le modèle ANOVA bidirectionnel suivant avec un terme d'erreur qui tient compte de la variabilité entre les individus:

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

Cependant, je ne sais pas si ma spécification du terme d'erreur est adéquate. Suis-je également censé tenir compte de la variabilité au sein d'un noyau? Pour moi, cette variabilité est la même que celle due à la position radiale qui est le principal facteur qui m'intéresse.

Bien que j'aie consacré du temps à lire sur la spécification du terme d'erreur dans ANOVA à mesures répétées, j'ai toujours des problèmes avec la spécification pratique du terme d'erreur. J'apprécierai de l'aide avec ceci.

Ce serait

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

Cela explique l'erreur dans le sujet de Radialposition. Si vous avez d'autres facteurs intra-sujet, jetez-les (dans une interaction) avec Radialpositiondans le dénominateur d'erreur, comme + Error(Individual/(Radialpostion*wiFactorA)). C'est ma compréhension de cela. Cela correspond aux mesures répétées de SPSS GLMsi vous n'avez aucune donnée manquante.