I would like to add a constraint in my EFA model to ensure that Z is diagonal. Z is a function of other matrices. I used an mxConstraint, where J is a conformable identity matrix. My question is, do you see a way of specifying this using labels instead of an mxConstraint?
algebraZ = mxAlgebra(expression = t(L)%%G%%solve(F), name = "Z") # should be diag
constraint1 = mxConstraint(Z == J*Z, name = "oblique")
Thanks in advance, I can give my full code if needed.