One row per participant, one interval per factor on the bounded
loading scale, with selected flags marked. The dotted rules sit at
the classical descriptive cut-off 1.96 / sqrt(J) (reference).
Usage
plot_loading_posterior(fit, q = 0.05)
# S3 method for class 'bayesqm_fit'
plot(x, ...)