Estoy revisando el paquete R OpenMx para un análisis de epidemiología genética con el fin de aprender a especificar y ajustar modelos SEM. Soy nuevo en esto, así que tengan paciencia conmigo. Estoy siguiendo el ejemplo en la página 59 de la Guía del usuario de OpenMx . Aquí dibujan el siguiente modelo conceptual:
Y al especificar las rutas, establecen el peso del nodo "uno" latente en los nodos bmi manifestados "T1" y "T2" en 0.6 porque:
Las principales rutas de interés son las de cada una de las variables latentes a la respectiva variable observada. Estos también se estiman (por lo tanto, todos se liberan), obtienen un valor inicial de 0.6 y las etiquetas apropiadas.
# 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")
),
El valor de 0.6 proviene de la covarianza estimada de bmi1
y bmi2
(de pares gemelos estrictamente mono zigóticos). Tengo dos preguntas:
Cuando dicen que a la ruta se le da un valor "inicial" de 0.6, ¿es esto como establecer una rutina de integración numérica con valores iniciales, como en la estimación de GLM?
¿Por qué este valor se estima estrictamente a partir de los gemelos monocigóticos?