   Hola a todos, mi consulta va relacionada con los paquete nlme o lmer. Mi
idea va relacionada a correr un modelo lineal mixto en el cual una de las
matrices de incidencia (Z), sea previamente formada por la suma de dos
matrices o por su puperposición (Z1+Z2=Z) dentro de esas funciones. Con un
ejemplo creo que me explicaré mejor

outputF<-lmer(formula=ALT~ (1|MADRE)+(1|PADRE)+(1|FAMILIA), data=datosF)

El modelo anterior me permite correr el efecto aleatorio de padre y madre
por separado, considerando que cada uno corresponde a una columna de mi
base de datos, pero mi intención es superponer estas matrices o sumarlas
(si es que no estoy errado en el término exacto) y correrlas como un efecto
aleatorio único.

La matriz en vez de asociar a un padre y una madre por separado
considerarlos conjuntamente, tras sumar las matrices.

1 1 0 0 0 0 0
0 0 1 0 0 1 0
0 1 0 1 0 0 0

Manualmente lo podría realizar, pero la idea es integrarlo en la función,
pre-ajuste de alguna forma, para obtener las varianzas. Saludos, atte

Felipe

	[[alternative HTML version deleted]]

