Is it possible to retrieve sums of squared loadings after applying varimax
rotation?

Here's the setup to my problem:

I ran PCA using prcomp(). I then applied the Kaiser criterion to retain
only the components having eigenvalues >= 1. (I know there's debate about
the wisdom of that criterion, but I don't want to get sucked into that.) I
then fed the reduced set of components to varimax() to simplify the
components. (I know there's debate about that as well, but that's also
irrelevant in this case.)

What I'd like to do is obtain the post-rotation sums of squared loadings
(eigenvalues). I know the original eigenvalues can be obtained by squaring
what's reported in the sdev vector, but I'm lost on how to get them after
rotation. I've been searching for hours to no avail. Any ideas?

Thanks :)

Eric

	[[alternative HTML version deleted]]

