J'ai des données recueillies à partir d'une expérience organisée comme suit:
Deux sites de 30 arbres chacun. 15 sont traités, 15 sont contrôlés sur chaque site. De chaque arbre, nous échantillonnons trois morceaux de la tige et trois morceaux des racines, soit 6 échantillons de niveau 1 par arbre, représentés par l'un des deux niveaux de facteur (racine, tige). Ensuite, à partir de ces échantillons de tige / racine, nous prélevons deux échantillons en disséquant différents tissus dans l'échantillon, ce qui est représenté par l'un des deux niveaux de facteur pour le type de tissu (type de tissu A, type de tissu B). Ces échantillons sont mesurés en tant que variable continue. Le nombre total d'observations est 720; 2 sites * 30 arbres * (trois échantillons de tige + trois échantillons de racines) * (un échantillon de tissu A + un échantillon de tissu B). Les données ressemblent à ceci ...
ï..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
J'essaie d'adapter un modèle à effets mixtes en utilisant R et lme4, mais je suis novice dans les modèles mixtes. J'aimerais modéliser la réponse sous la forme traitement + facteur de niveau 1 (tige, racine) + facteur de niveau 2 (tissu A, tissu B), avec des effets aléatoires pour les échantillons spécifiques imbriqués dans les deux niveaux.
En R, je le fais en utilisant lmer, comme suit
fit <- lmer(Response ~ Treatment + Organ + Tissue + (1|Tree/Organ/Sample))
D'après ma compréhension (... ce qui n'est pas certain, et pourquoi je poste!), Le terme:
(1|Tree/Organ/Sample)
Spécifie que 'Sample' est imbriqué dans les échantillons d'organes, ce qui est imbriqué dans l'arborescence. Ce type d'imbrication est-il pertinent / valide? Désolé si cette question n'est pas claire, dans l'affirmative, veuillez préciser où je peux élaborer.