Orthomax Rotation
Produces a simultaneous orthomax rotation of two matrices (using one rotation matrix).
orthmax2(A1, A2, gam1, gam2, conv)
A1 |
First matrix to be rotated with the same number of columns of A2 |
A2 |
Second matrix to be rotated with the same number of columns of A1 |
gam1 |
orthmax parameter for |
gam2 |
orthmax parameter for |
conv |
Optional convergence value (default 1e-6) |
A list including the following components:
B1 |
Rotated version of |
B2 |
Rotated version of |
T |
Rotation matrix |
f |
Orthomax function value |
The function to be maximized is f=sum((A1^2)-1/m1*gam1*sum((sum(A1^2))^2))^2+sum((A2^2)-1/m2*gam2*sum((sum(A2^2))^2))^2.
Maria Antonietta Del Ferraro mariaantonietta.delferraro@yahoo.it
Henk A.L. Kiers h.a.l.kiers@rug.nl
Paolo Giordani paolo.giordani@uniroma1.it
R. Jennrich (1970). Orthogonal rotation algorithms. Psychometrika 35:229–235.
X <- matrix(rnorm(8*3),ncol=3) Y <- matrix(rnorm(6*3),ncol=3) orthXY <- orthmax2(X,Y,1,2) # rotated version of X orthXY$B1 # rotated version of Y orthXY$B2 # rotation matrix orthXY$T
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.