我正在审查R包OpenMx进行遗传流行病学分析,以了解如何指定和拟合SEM模型。我对此很陌生,所以请多多包涵。我正在遵循《OpenMx用户指南》第59页上的示例。他们在这里绘制以下概念模型:
在指定路径时,他们将潜在的“一个”节点对显示的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
(严格的单合子双胞胎)。我有两个问题:
当他们说路径的“开始”值为0.6时,是否像在估计GLM时那样设置具有初始值的数值积分例程?
为什么严格根据单卵双胞胎估算这个值?