Aşağıdaki şekilde düzenlenen bir deneyden toplanan veriler var:
Her biri 30 ağaçlı iki alan. 15 işlenir, 15 her yerde kontrol edilir. Her bir ağaçtan, gövdenin üç parçasını ve üç parçanın köklerini örnekliyoruz; bu nedenle, iki faktör seviyesinden (kök, gövde) biriyle temsil edilen ağaç başına 6 seviye 1 örnek. Daha sonra, bu kök / kök örneklerinden, doku tipi için iki faktör seviyesinden biri ile temsil edilen, örnek içindeki farklı dokuları disekte ederek iki örnek alırız (doku tipi A, doku tipi B). Bu örnekler sürekli değişken olarak ölçülür. Toplam gözlem sayısı 720'dir; 2 alan * 30 ağaç * (üç kök örnek + üç kök örnek) * (bir doku A örnek + bir doku B örnek). Veriler şuna benziyor ...
ï..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
R ve lme4 kullanarak karışık efektler modeline uymaya çalışıyorum, ancak karışık modeller için yeni. Yanıtı, iki seviyeye yerleştirilmiş spesifik örnekler için rastgele etkilerle Tedavi + Seviye 1 Faktörü (kök, kök) + Seviye 2 Faktörü (doku A, doku B) olarak modellemek istiyorum.
R, bunu aşağıdaki gibi lmer kullanarak yapıyorum
fit <- lmer(Response ~ Treatment + Organ + Tissue + (1|Tree/Organ/Sample))
Anladığım kadarıyla (... kesin değil ve neden gönderiyorum!) Terimi:
(1|Tree/Organ/Sample)
'Örnek' öğesinin, ağacın içine yerleştirilen organ örneklerinde iç içe olduğunu belirtir. Bu tür yuvalama uygun mu / geçerli mi? Üzgünüm, eğer bu soru net değilse, eğer öyleyse, lütfen nerede çalışabileceğimi belirtin.