orthmax2 {ThreeWay}

Orthomax Rotation

Description

Produces a simultaneous orthomax rotation of two matrices (using one rotation matrix).

Usage

 orthmax2(A1, A2, gam1, gam2, conv)

Arguments

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 A1

gam2

orthmax parameter for A2

conv

Optional convergence value (default 1e-6)

Value

A list including the following components:

B1

Rotated version of A1

B2

Rotated version of A2

T

Rotation matrix

f

Orthomax function value

Note

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.

Author(s)

Maria Antonietta Del Ferraro [e-mail]
Henk A.L. Kiers [e-mail]
Paolo Giordani [e-mail]

References

R. Jennrich (1970). Orthogonal rotation algorithms. Psychometrika 35:229–235.

See Also

varim

Examples

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