Tengo duda respecto a que R-paquete de utilizar para crear una clase latente/modelo de mezcla con tanto categóricas y continuas variables indicadoras. Todavía no he encontrado un buen ejemplo de este uso de R, a pesar de que hay una gran cantidad de la mezcla y latente de la clase de análisis de paquetes en R. El paquete de software Mplus tienen esta capacidad y en la UCLA sitio web no es un ejemplo de cómo crear un modelo con Mplus (http://www.ats.ucla.edu/stat/mplus/seminars/introMplus_part2/lca.htm/ el segundo ejemplo).
La única manera que he sido capaz de replicar este ejemplo en R es el uso de la depmixS4 paquete. Usar el siguiente código:
library(depmixS4)
lca<-read.table("http://www.ats.ucla.edu/stat/mplus/seminars/introMplus_part2/lca.dat", sep=",")
names(lca)<-c( "hm", "he", "voc", "nocol" ,"ach9", "ach10", "ach11","ach12")
mod <- mix(list(hm~1,he~1,voc~1,nocol~1,ach9~1,ach10~1 ,ach11~1 ,ach12~1 ), data=lca, nstates=2,
family=list(multinomial(),multinomial(),multinomial(),multinomial(), gaussian(),gaussian(),gaussian(),gaussian()),
respstart=runif(32))
fmod<-fit(mod)
summary(fmod)
Este código parece hacer el trabajo de replicar el ejemplo con algunas diferencias menores entre la salida en el depmixS4 variante y la Mplus.
Primero, este es el código correcto para la reproducción de la Mplus ejemplo?
Segundo, ¿alguien tiene alguna otra paquetes que recomiendo a la hora de crear este tipo de modelo en R?