mgm (no la de las pelis sino la de los modelos gráficos)

Cayeron en mis manos unos datos que no puedo publicar, pero me atreveré a presentar algunos resultados anonimizados. Se trata de una tabla de puntuaciones numéricas (18 en total, cada una en su columna) proporcionadas por unos cuantos centenares de sujetos (filas). Era de interés un estudio cualitativo de las posibles relaciones de dependencia entre las variables.

La manera más rápida de comenzar, un heatmap(cor(dat)), para obtener

Y luego PCA y todas esas cosas.

Pero esos datos fueron la excusa perfecta para ensayar algo que todavía no tengo muy claro lo que hace (porque no he acabado de leer con detalle esto): mgm.

Así que

modelo <- mgm(as.matrix(dat),
    type = rep("g", ncol(dat)),
    level = rep(1, ncol(dat)),
    k = 2,
    lambdaSel = "EBIC", lambdaGam = 0.25)

y luego

qgraph(modelo$pairwise$wadj,
    layout = 'spring', repulsion = 1.3,
    edge.color = modelo$pairwise$edgecolor,
    nodeNames = colnames(dat),
    color = color,
    groups = my_groups,
    legend.mode = "style2", legend.cex=.4,
    vsize = 3.5, esize = 15)

para obtener

que podrá ser muchas cosas (algunas que aún desconozco) pero que nunca dejará de ser la mar de molón.

(Y dejo al lector el ejercicio de compararlo con el heatmap de más arriba).