Especificación del término Error () en medidas repetidas ANOVA en R

Tengo problemas para definir los términos de error para un ANOVA de medidas repetidas bidireccionales en R. Mis datos consisten en estimaciones de densidad de madera para tres posiciones radiales (interior, medio y exterior) a lo largo de un núcleo extraído de un árbol. Hay un total de 20 especies de árboles, 6 individuos de cada especie y dos núcleos de cada árbol.

Para probar el efecto de la posición radial en la densidad de la madera, utilizo el siguiente modelo ANOVA de dos vías con un término de error que explica la variabilidad entre individuos:

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

Sin embargo, no estoy seguro de si mi especificación del término de error es adecuada. ¿También debo dar cuenta de la variabilidad dentro de un núcleo? Para mí, esta variabilidad es la misma debido a la posición radial, que es el factor principal que me interesa.

Aunque he dedicado algo de tiempo a leer acerca de cómo especificar el término de error en ANOVA de medidas repetidas, todavía tengo problemas para especificar prácticamente el término de error. Apreciaré algo de ayuda con esto.

Podría ser

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

Eso explica el error dentro del sujeto de Radialposition. Si tiene otros factores dentro del sujeto, tírelos (en una interacción) con Radialpositionel denominador de error, como + Error(Individual/(Radialpostion*wiFactorA)). Esa es mi comprensión de eso. Eso coincide con las medidas repetidas de SPSS GLMsi no le faltan datos.