Hello,
recently I was looking at the implementation of the "varimax" rotation
procedure from the "stats" package and to me it looks quite different
from the algorithm originally suggested by Kaiser in 1958.
The R procedure iteratively uses singular value decompositions of some
matrices whereas Kaiser proposed to iteratively compute rotation
matrices between all pairs of factors which does not seem to happen ( at
least not explicitely ) in the R version.
My question now is whether R uses a completely different approach than
Kaiser (if so, then could you please point me to a publication or
explanation of the algorithm used since I wasn't able to find any) or if
it is the Kaiser method just well hidden under quite a bit of clever
linear algebra ( explanations on why the methods are equal is also
appreciated).
Thanks for any hints or clarifications!
With best regards
Sebastian Starke
______________________________________________
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
