¿Cómo especificar múltiples pendientes aleatorias en glmmPQL (paquete MASS)?

Question

Me gustaría especificar más de un efecto aleatorio en un modelo lineal mixto generalizado con glmmPQL (paquete MASS). Actualmente estoy utilizando los datos de Heck, Thomas y Tabata (2012).

Así que lo que me gustaría conseguir es probar, si las pendientes de los predictores de nivel 1 (femeninos, menores y bajos) varían entre las unidades de nivel 2.

¿Es este el código correcto para hacerlo? Si no - ¿Cómo puedo especificar múltiples efectos aleatorios simultáneamente?

PQL<-glmmPQL(readprof~female+minor+lowses+schcomp+smallsch, random=(~1|schcode/minor/female/lowses), family = binomial(link = logit), data = Heck)

Muchas gracias.

Answer 1

Así que lo que me gustaría conseguir es probar si las pendientes de los predictores del nivel 1 (femenino, menor y menor) varían entre las unidades del nivel 2.

Soy nuevo en glmmPQL y en los modelos de efectos mixtos, y no estoy del todo seguro de lo que quiere decir con "unidades de nivel 2" aquí, así que no puedo comprobar que su modelo está construido correctamente. Sin embargo, creo que puedo ayudarte con tu pregunta:

¿Cómo puedo especificar varios efectos aleatorios simultáneamente?

Según el glmmPQL documentación , es necesario crear una matriz de las fórmulas para introducir múltiples efectos aleatorios.

PQL<-glmmPQL(readprof~female+minor+lowses+schcomp+smallsch, random=list(schcode = ~1, minor = ~1, female = ~1, lowses = ~1), family = binomial(link = logit), data = Heck)

Encontré la información anterior en stackoverflow .

Las barras inclinadas se utilizan para covariables en un modelo, en lugar de permitirle simplemente introducir los múltiples efectos aleatorios. Puedes ver cómo funcionan en el enlace.

Espero que esto ayude.