psych::principal - explanation for the order and naming of rotated (principal) components

Question

Let x be a sample dataframe.

set.seed(0)
x <- replicate(4, rnorm(10))

A PCA using the principal function from

Answer 1

Here's a partial answer. I looked through the code behind the principal function, and can see clearly where the reordering happens:

 if (nfactors > 1) {
    ev.rotated <- diag(t(loadings) %*% loadings)
    ev.order <- order(ev.rotated, decreasing = TRUE)
    loadings <- loadings[, ev.order]
  }

So the code above is the reason the order changes, but it's less clear what purpose that serves. I don't have enough experience in rotations to be able to discern the package author's intent.

Answer 2

Mark, The logic is to recognize what rotation does. This is more for pedagogical reasons than anything else. I am trying to show the relationship of the original components to the rotated components. To take your example, look at the loadings, not just the variances accounted for.

unrotated:

    PC1   PC2   PC3   PC4 h2       u2
1 -0.77 -0.40  0.39  0.32  1 -6.7e-16
2  0.71 -0.28  0.63 -0.17  1  6.7e-16
3 -0.10  0.93  0.35  0.09  1  6.7e-16
4  0.90 -0.02 -0.13  0.42  1  2.2e-16

Rotated:
    RC4   RC3   RC2   RC1 h2       u2
1  0.95 -0.10 -0.08 -0.29  1 -6.7e-16
2 -0.10  0.97 -0.06  0.22  1  6.7e-16
3 -0.07 -0.06  0.99 -0.05  1  6.7e-16
4 -0.34  0.27 -0.07  0.90  1  2.2e-16

In particular, look at variables 3 and 4. In the unrotated solution, they define PC2 and PC1 respectively. Now look at the rotated solution. These two still mark PC2 and PC1 (and are labeled RC2 and RC1 to reflect that they are rotated), but the variances accounted for have changed as PC4 when rotated to RC4 now soaks up more variance. (This is also true for PC3 and PC4 but not as clear.)

What I am trying to do is represent what happens as you rotate. PC1 is rotated to a simpler structure, and becomes RC1.

Then, because many people like to have their components in order of variance accounted for, I sort by the eigen value (sum squares accounted for).

I believe what other programs do is to rotate and relabel so that the components are always called C1 ... Cn. I just like to see where the components came from.

If you think it is useful, I can (eventually) add this discussion to the documentation for principal as well as fa.

Bill