After having performed Factor Analysis (FA), I am now trying to implement the following rotation strategy (Swanson 2017 - details of rotation strategy in appenix A) to interpret the three factors structurally and economically. The strategy is in 3 steps but I will report just the first one as the other steps will follow suit:
1) Quadratic minimization problem subject to a linear constraint:
is the 80x3 matrix of unrotated factors, (3x1) is the first column of λ which is the 7x3 matrix of factor loadings and is the third column of U which is a 3x3 orthogonal matrix. and are known (find data attached: fa_pre = (80x3) and lambda_1 = (1x3)). Obviously, I would like to solve it for to fully have .
Can you help me code the first minimization problem?