Especificando o termo Error () em medidas repetidas ANOVA em R

Estou tendo problemas para definir os termos de erro para uma ANOVA de medidas repetidas em dois sentidos em R. Meus dados consistem em estimativas da densidade da madeira para três posições radiais (interna, média e externa) ao longo de um núcleo extraído de uma árvore. Há um total de 20 espécies de árvores, 6 indivíduos de cada espécie e dois núcleos de cada árvore.

Para testar o efeito da posição radial na densidade da madeira, utilizo o seguinte modelo ANOVA de duas vias com um termo de erro que explica a variabilidade entre indivíduos:

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

No entanto, não tenho certeza se minha especificação do termo de erro é adequada. Também devo explicar a variabilidade dentro de um núcleo? Para mim, essa variabilidade é a mesma devido à posição radial, que é o principal fator no qual estou interessado.

Embora tenha dedicado algum tempo a ler sobre a especificação do termo de erro na ANOVA de medidas repetidas, ainda tenho problemas em praticamente especificar o termo de erro. Agradeço alguma ajuda com isso.

Seria

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

Isso explica o erro dentro do assunto de Radialposition. Se você tiver outros fatores dentro do assunto, jogue-os (em uma interação) com Radialpositionno denominador de Erro, como + Error(Individual/(Radialpostion*wiFactorA)). Essa é a minha compreensão disso. Isso corresponde às medidas repetidas do SPSS GLMse você não tiver dados ausentes.