ฉันมีข้อมูลที่รวบรวมจากการทดสอบที่จัดระเบียบดังนี้:
สองไซต์แต่ละแห่งมีต้นไม้ 30 ต้น ได้รับการปฏิบัติ 15, 15 คือการควบคุมในแต่ละเว็บไซต์ จากต้นไม้แต่ละต้นเราสุ่มตัวอย่างก้านสามชิ้นและรากสามชิ้นดังนั้น 6 ระดับ 1 ตัวอย่างต่อต้นซึ่งแสดงโดยหนึ่งในสองระดับปัจจัย (รากลำต้น) จากนั้นตัวอย่างต้นกำเนิด / รากเราใช้สองตัวอย่างโดยการผ่าเนื้อเยื่อต่าง ๆ ภายในตัวอย่างซึ่งแสดงโดยหนึ่งในสองระดับปัจจัยสำหรับประเภทเนื้อเยื่อ (ประเภทเนื้อเยื่อ A, ประเภทเนื้อเยื่อ B) ตัวอย่างเหล่านี้วัดเป็นตัวแปรต่อเนื่อง จำนวนการสังเกตทั้งหมดคือ 720; 2 ไซต์ * ต้นไม้ 30 ต้น * (ตัวอย่างลำต้นสามชิ้น + ตัวอย่างรากสามต้น) * (เนื้อเยื่อหนึ่งตัวอย่าง + เนื้อเยื่อหนึ่งตัวอย่าง B) ข้อมูลมีลักษณะเช่นนี้ ...
ï..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 และ lme4 แต่เป็นของใหม่สำหรับโมเดลผสม ฉันต้องการสร้างแบบจำลองการตอบสนองเช่นการรักษา + ระดับ 1 ปัจจัย (ลำต้น, ราก) + ปัจจัยระดับ 2 (เนื้อเยื่อ A, เนื้อเยื่อ B) โดยมีเอฟเฟกต์แบบสุ่มสำหรับตัวอย่างเฉพาะที่ซ้อนกันภายในสองระดับ
ใน R ฉันกำลังทำสิ่งนี้โดยใช้ lmer ดังนี้
fit <- lmer(Response ~ Treatment + Organ + Tissue + (1|Tree/Organ/Sample))
จากความเข้าใจของฉัน (... ซึ่งไม่แน่ใจและทำไมฉันโพสต์!) คำ:
(1|Tree/Organ/Sample)
ระบุว่า 'ตัวอย่าง' ซ้อนอยู่ภายในตัวอย่างอวัยวะซึ่งซ้อนอยู่ภายในต้นไม้ การทำรังแบบนี้เกี่ยวข้องหรือไม่ถูกต้อง? ขออภัยหากคำถามนี้ไม่ชัดเจนถ้าเป็นเช่นนั้นโปรดระบุว่าฉันสามารถทำอย่างละเอียดได้ที่ไหน