Necesito analizar un conjunto de datos de datos de rehabilitación clínica. Estoy interesado en las relaciones basadas en hipótesis entre la "entrada" cuantificada (cantidad de terapia) y los cambios en el estado de salud. Aunque el conjunto de datos es relativamente pequeño (n ~ 70), hemos repetido datos que reflejan cambios temporales en ambos. Estoy familiarizado con el modelado de efectos mixtos no lineales en R, sin embargo, estoy interesado en posibles relaciones "causales" entre la entrada y la salida aquí y, por lo tanto, estoy considerando aplicaciones de medidas repetidas de SEM
Agradecería consejos sobre si alguno de los paquetes SEM para R (sam, lavaan, openmx?) Son los más adecuados para datos de medidas repetidas, y particularmente recomendaciones para libros de texto (¿hay un "Pinheiro y Bates" del campo?) .
Respuestas:
Creo que quieres un modelo de curva de crecimiento latente. Si bien solo lo he usado
LISREL
para esto,lavaan package documentation
indica que se puede usar para adaptarse a este tipo de modelo.No conozco ningún libro que se especialice en este tema, el libro del que estoy trabajando para SEM cubre una variedad de métodos. Quizás alguien más pueda responder ese aspecto de su pregunta.
fuente
lavaan
/ Mx con Mplus, publicaré una respuesta yo mismo.No, no hay "Pinheiro y Bates". Puede encontrar varios libros titulados como "SEM con AMOS / LISREL / Mplus", pero no conozco ninguno con R. El mejor libro, matemáticamente hablando, sobre SEM sigue siendo Bollen (1989) . Está escrito por un sociólogo en lugar de un bioestadístico (¡aunque es muy bueno!), Por lo que está dirigido a científicos sociales y contiene pocas referencias al software (y de todos modos no desea el software de hace un cuarto de siglo) . Bollen también ha sido coautor de un buen artículo recientemente sobre causalidad con Judea Pearl, consulte http://ftp.cs.ucla.edu/pub/stat_ser/r393.pdf . Por lo que puedo decir, Mulaik (2009) también debería ser bueno, pero está escrito por un psicólogo para psicólogos.
No creo que sem package sea lo suficientemente flexible como para ejecutar este tipo de cosas. OpenMx puede manejar datos ordinales (y, por lo tanto, resultados binarios), pero no creo que lavaan pueda hacer esto.
El software que conceptualmente encontrará más fácil de manejar podría ser GLLAMM , un paquete escrito para Stata . Visto de una manera, esto es esencialmente una encarnación de Stata
nlme
. Con un ajuste adicional (que permite que los coeficientes de los efectos aleatorios varíen según los valores de otras variables), se convierte en un paquete de modelado de variables latentes. Todo esto se describe en Skrondal y Rabe-Hesketh (2004) ... que es un gran libro per se que te gustaría tener incluso si lo hacesnlme
.fuente
gllamm
, visto de otra manera, desde la perspectiva de un psicométrico acostumbrado a los modelos IRT: es terriblemente lento :-)polychoric
, por ejemplo, cuando lo necesitaba.Como parece cómodo con los modelos mixtos lineales generalizados, y no parece implicar que le interesan las variables latentes, quizás desee adoptar un enfoque por partes
lmer
que luego pueda evaluar mediante una prueba D-Sep. Ver Shipley, B. (2009). Análisis de ruta confirmatoria en un contexto multinivel generalizado. Ecology, Ecology, 90, 363–368. http://dx.doi.org/10.1890/08-1034.1 para un ejemplo. También proporciona el código R en el apéndice sobre cómo calcular la prueba de D-Separation.Si realmente desea entrar en el modelado de variables latentes y SEM con la máxima probabilidad, consulte http://lavaan.org : hay un gran tutorial que cubre sus capacidades, así como una sección sobre modelos de curva de crecimiento latente que bien puede ser lo que que estas buscando
fuente