ฉันกำลังทบทวนแพคเกจ R OpenMx สำหรับการวิเคราะห์ทางระบาดวิทยาทางพันธุกรรมเพื่อเรียนรู้วิธีการระบุและเหมาะสมกับแบบจำลอง SEM ฉันยังใหม่กับสิ่งนี้ดังนั้นทนกับฉัน ฉันกำลังตัวอย่างต่อไปนี้ในหน้า 59 ของคู่มือการใช้งาน OpenMx ที่นี่พวกเขาวาดโมเดลแนวคิดต่อไปนี้:
และในการระบุเส้นทางพวกเขาตั้งค่าน้ำหนักของโหนด "หนึ่ง" แฝงไปยังโหนด bmi "T1" และ "T2" ที่ประจักษ์เป็น 0.6 เพราะ:
เส้นทางหลักที่น่าสนใจคือจากตัวแปรแฝงแต่ละตัวไปยังตัวแปรที่สังเกตได้ สิ่งเหล่านี้ได้รับการประเมิน (ซึ่งทั้งหมดถูกตั้งค่าไว้ฟรี) รับค่าเริ่มต้น 0.6 และป้ายกำกับที่เหมาะสม
# path coefficients for twin 1
mxPath(
from=c("A1","C1","E1"),
to="bmi1",
arrows=1,
free=TRUE,
values=0.6,
label=c("a","c","e")
),
# path coefficients for twin 2
mxPath(
from=c("A2","C2","E2"),
to="bmi2",
arrows=1,
free=TRUE,
values=0.6,
label=c("a","c","e")
),
ค่าของ 0.6 นี้ได้มาจากความแปรปรวนโดยประมาณbmi1
และbmi2
(จากเคร่งครัดโมโนคู่แฝด Zygotic) ฉันมีสองคำถาม:
เมื่อพวกเขาบอกว่าเส้นทางได้รับค่า "เริ่มต้น" เป็น 0.6 เช่นนี้การตั้งค่ารูทีนการรวมตัวเลขที่มีค่าเริ่มต้นเช่นในการประมาณค่าของ GLMs?
เหตุใดค่านี้จึงถูกประเมินอย่างเข้มงวดจากแฝดแฝด