Tengo datos recopilados de un experimento organizado de la siguiente manera:
Dos sitios, cada uno con 30 árboles. 15 son tratados, 15 son control en cada sitio. De cada árbol, tomamos muestras de tres piezas del tallo y tres piezas de las raíces, por lo que 6 muestras de nivel 1 por árbol, que está representado por uno de los dos niveles de factores (raíz, tallo). Luego, de esas muestras de tallo / raíz, tomamos dos muestras diseccionando diferentes tejidos dentro de la muestra, que está representado por uno de los dos niveles de factores para el tipo de tejido (tipo de tejido A, tipo de tejido B). Estas muestras se miden como una variable continua. El número total de observaciones es 720; 2 sitios * 30 árboles * (tres muestras de tallo + tres muestras de raíz) * (una muestra de tejido A + una muestra de tejido B). Los datos se ven así ...
ï..Site Tree Treatment Organ Sample Tissue Total_Length
1 L LT1 T R 1 Phloem 30
2 L LT1 T R 1 Xylem 28
3 L LT1 T R 2 Phloem 46
4 L LT1 T R 2 Xylem 38
5 L LT1 T R 3 Phloem 103
6 L LT1 T R 3 Xylem 53
7 L LT1 T S 1 Phloem 29
8 L LT1 T S 1 Xylem 21
9 L LT1 T S 2 Phloem 56
10 L LT1 T S 2 Xylem 49
11 L LT1 T S 3 Phloem 41
12 L LT1 T S 3 Xylem 30
Estoy intentando ajustar un modelo de efectos mixtos usando R y lme4, pero soy nuevo en los modelos mixtos. Me gustaría modelar la respuesta como el Tratamiento + Factor de Nivel 1 (tallo, raíz) + Factor de Nivel 2 (tejido A, tejido B), con efectos aleatorios para las muestras específicas anidadas dentro de los dos niveles.
En R, estoy haciendo esto usando lmer, de la siguiente manera
fit <- lmer(Response ~ Treatment + Organ + Tissue + (1|Tree/Organ/Sample))
Según tengo entendido (... lo cual no es seguro, y por qué estoy publicando) el término:
(1|Tree/Organ/Sample)
Especifica que 'Muestra' está anidada dentro de las muestras de órganos, que está anidada dentro del árbol. ¿Es este tipo de anidamiento relevante / válido? Lo siento si esta pregunta no está clara, si es así, especifique dónde puedo dar más detalles.