Comprobación de supuestos para efectos aleatorios en modelos de efectos mixtos anidados en R / S-Plus

Question

Estoy elaboración de modelos anidados modelo de efectos mixtos con solo el intercepto en la parte aleatoria, de la forma:

fit4<-lme(fixed = Stdquad~factor(LayoutN)+factor(nCHIPS.fixed), random = ~1|Class.Ordered/student)

Cuando intento comprobar la hipótesis sobre los independientes, idénticamente distribuidas (iid) de efectos aleatorios en SPlus usando ranef(), sigo recibiendo el mensaje de error:

Problem in sort.list: ordering not defined for mode "list": sort.list(x, partial)

Sospecho que la anidación es el problema, ya que si me quito el estudiante, las parcelas venir hasta bien.

Me gustaría saber cómo puedo comprobar que los efectos aleatorios son iid y son independientes para los diferentes grupos con S-Plus (o R ya que es muy similar).

Answer 1

Parece que el uso de la nlme paquete. Tal vez valdría la pena intentar R y el lme4 en lugar de ello, aunque no es totalmente comparable wrt. la sintaxis o la llamada a la función.

En su caso, sugiero que para especificar el level cuando usted llama a ranef(), ver ?ranef.lme:

 nivel: opcional, un vector de enteros positivos dando a los niveles de
 agrupación para ser utilizado en la extracción de los efectos aleatorios a partir de una
 objeto de múltiples anidada los niveles de agrupación. Los valores predeterminados para todos
 los niveles de agrupación.

También está presente en el boletín oficial de la documentación para NLME 3.0 (por ejemplo, p. 17).

Echa un vistazo Douglas Bates de ordenadas folletos en GLMM. Él es también la escritura de un libro de texto titulado lme4: efectos Mixtos de modelado con R. Todos ellos están disponibles en R-forge.