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.