g <- rep(c("g1", "g2"), length.out = nrow(GRiPS_raw))
mg <- efa_group(GRiPS_raw, groups = g, n_factors = 1)
# Per-factor congruence against the Lorenzo-Seva & ten Berge bands
plot(mg)
# Per-item cross-group loading-difference heatmap
plot(mg, type = "differences")
Run the code above in your browser using DataLab