Tutorials › Multivariate Statistics (Graduate) › Common Mistakes in Rotation, Scores, and Fit Assessment

Factor Rotation and Assessment · Tutorial 579 of 975

Common Mistakes in Rotation, Scores, and Fit Assessment

Use the factor model’s covariance structure to check rotated loadings, factor comparisons, score interpretations, and claims about fit.

Advanced 10 min read

What You'll Learn

  • Distinguish pattern loadings from structure loadings in an oblique solution.
  • Calculate communalities using the factor correlation matrix.
  • Align factors before interpreting cross-solution comparisons.
  • Describe factor scores as uncertain estimates rather than observed latent values.
  • Separate rotation criteria from covariance-model fit assessment.
  • Check which quantities remain unchanged under a factor transformation.

Why These Errors Matter

A rotated factor solution can be easy to read but easy to overinterpret. A simple-looking loading matrix does not, by itself, establish good model fit; an estimated factor score is not the person’s known latent value; and two sets of loadings can look different merely because their factor axes are ordered or signed differently. These distinctions matter because rotation changes the coordinates used to describe common variation, while fit concerns how well the model reproduces the observed covariance structure.

Consider the common factor model for a \(p\)-dimensional observation \(X\) and \(q\) factors:

$$ X=\mu+\Lambda F+\varepsilon,\qquad \operatorname{Cov}(F)=\Phi,\qquad \operatorname{Cov}(\varepsilon)=\Psi. $$

Here \(\Lambda\) is the \(p\times q\) pattern loading matrix, \(\Phi\) is the factor covariance matrix, and \(\Psi\) is the uniqueness covariance matrix, commonly taken to be diagonal. Assume \(F\) and \(\varepsilon\) are uncorrelated and have finite second moments. Then the model-implied covariance matrix is \(\Sigma=\Lambda\Phi\Lambda^\mathsf{T}+\Psi\). In an orthogonal solution, \(\Phi=I\); in an oblique solution, factors may be correlated.

Definition: The communality of variable \(i\), the variance attributed to the common factors, is \(h_i^2=(\Lambda\Phi\Lambda^\mathsf{T})_{ii}=\lambda_i^\mathsf{T}\Phi\lambda_i\), where \(\lambda_i^\mathsf{T}\) is row \(i\) of \(\Lambda\). For an orthogonal solution this reduces to the sum of squared pattern loadings in that row. For an oblique solution, the factor correlations also contribute.

Oblique Loadings: Do Not Square the Pattern Row and Stop

In an oblique solution, the pattern matrix describes the conditional contribution of each factor to a variable, given the other factors. The structure matrix describes the variable–factor covariances; with standardized variables and factors, its entries are correlations and it is \(S=\Lambda\Phi\). The two matrices are not interchangeable when factors correlate. In particular, summing squared entries of a pattern row omits cross-products involving the factor correlations.

Worked Example: Computing an Oblique Communality

Suppose a standardized variable has pattern loadings \(\lambda_i^\mathsf{T}=(0.65,0.25)\), and the two standardized factors have correlation \(0.50\). Thus

$$ \Phi= \begin{pmatrix} 1&0.50\\ 0.50&1 \end{pmatrix}. $$

The squared-pattern shortcut gives \(0.65^2+0.25^2=0.4225+0.0625=0.485\). The correct communality includes the cross-product terms:

$$ h_i^2 =(0.65,0.25) \begin{pmatrix} 1&0.50\\ 0.50&1 \end{pmatrix} \begin{pmatrix} 0.65\\ 0.25 \end{pmatrix} =0.65^2+0.25^2+2(0.50)(0.65)(0.25) =0.485+0.1625=0.6475. $$

If the variable is standardized and the model has a nonnegative uniqueness variance, the uniqueness is \(1-0.6475=0.3525\), so the quantities sum to one as required. The shortcut would instead report a uniqueness of \(1-0.485=0.515\), which incorrectly assigns too much variance to uniqueness.

Why Communality Survives a Change of Factor Coordinates

A useful check is that the communality does not depend on the particular admissible factor coordinates used to represent the common covariance. Let \(T\) be any nonsingular \(q\times q\) matrix, and define \(F^*=T^{-1}F\). The corresponding loadings and factor covariance are \(\Lambda^*=\Lambda T\) and \(\Phi^*=T^{-1}\Phi T^{-\mathsf{T}}\).

Theorem: Under this nonsingular change of factor coordinates, the common covariance and each variable’s communality are unchanged.

Proof. Using \(\Lambda^*=\Lambda T\) and \(\Phi^*=T^{-1}\Phi T^{-\mathsf{T}}\), multiply the transformed quantities:

$$ \Lambda^*\Phi^*(\Lambda^*)^\mathsf{T} =(\Lambda T)(T^{-1}\Phi T^{-\mathsf{T}})(\Lambda T)^\mathsf{T} =\Lambda T T^{-1}\Phi T^{-\mathsf{T}}T^\mathsf{T}\Lambda^\mathsf{T} =\Lambda\Phi\Lambda^\mathsf{T}. $$

The diagonal entries of these equal matrices are the communalities, so they are equal as well. The uniqueness matrix \(\Psi\) is unchanged in this coordinate transformation, and therefore the full model covariance is unchanged. \(\square\)

This result also explains why squaring an oblique pattern row is not a coordinate-invariant way to calculate communality: it discards \(\Phi\), which must change along with the loadings when the factor coordinates change.

Compare Factors Only After Addressing Alignment

Factor labels have no inherent order, and the sign of a factor can be reversed together with the signs of its loadings. Consequently, raw column-by-column comparisons may confuse an equivalent relabeling with a substantive difference. One standard summary for comparing two loading columns \(a\) and \(b\) is Tucker’s congruence coefficient:

$$ \phi(a,b)=\frac{a^\mathsf{T}b}{\sqrt{(a^\mathsf{T}a)(b^\mathsf{T}b)}}. $$

It is a normalized inner product, so reversing one column’s sign reverses the coefficient’s sign. A negative coefficient can therefore signal a sign mismatch rather than a different factor. Permutations and sign changes should be considered before interpreting column matches. For oblique solutions, matching may require more than these simple adjustments: distinct oblique coordinate systems can mix factors. Any alignment should respect the model and substantive constraints, not just make columns look similar.

Worked Example: A Low Congruence Caused by Swapped Factors

Suppose one solution has columns \(a_1=(0.8,0.7,0.1)^\mathsf{T}\) and \(a_2=(0.1,0.2,0.75)^\mathsf{T}\). A second solution lists the same factors in the order \(b_1=a_2,\ b_2=a_1\). Comparing the first columns without checking their order gives

$$ \phi(a_1,b_1) =\frac{(0.8)(0.1)+(0.7)(0.2)+(0.1)(0.75)} {\sqrt{(0.8^2+0.7^2+0.1^2)(0.1^2+0.2^2+0.75^2)}} =\frac{0.295}{\sqrt{(1.14)(0.6125)}} \approx 0.353. $$

That modest value might look like poor agreement. But \(b_2=a_1\) exactly, and \(b_1=a_2\) exactly. After matching the columns in the correct order, both congruence coefficients are 1. The apparent mismatch was only a permutation.

If two solutions are not related by a simple permutation or sign change, a low congruence may reflect a different orientation, sampling variation, or genuinely different structure. It is not, by itself, a verdict about fit. In particular, aligning factors is an interpretive comparison; it does not make two fitted covariance models equivalent.

Factor Scores Are Estimates, Not Observed Latent Values

A factor score estimate summarizes information about an unobserved factor given observed variables and a chosen scoring method. It is not a direct measurement of the factor. Score indeterminacy and conditional uncertainty mean that multiple plausible latent values can be compatible with the same observed data. A point estimate should therefore be reported as an estimate, and its uncertainty should not be suppressed when it affects interpretation.

Worked Example: A Score Estimate With Remaining Uncertainty

Let a single factor satisfy \(F\sim N(0,1)\), and let one standardized observation follow \(X=0.8F+\varepsilon\), where \(\varepsilon\sim N(0,0.36)\) is independent of \(F\). For an observed \(x=1.2\), the conditional distribution is normal. Completing the square in the joint density gives conditional variance \(0.36/(0.8^2+0.36)\) and conditional mean \(0.8x/(0.8^2+0.36)\). Since \(0.8^2+0.36=1\), these are \(0.36\) and \(0.96\), respectively:

$$ F\mid X=1.2\sim N(0.96,0.36). $$

The conditional mean is a score estimate, but the conditional standard deviation is \(\sqrt{0.36}=0.60\). Reporting only “the person’s factor value is 0.96” overstates what the observation establishes: the latent value remains uncertain. This example assumes a specified Gaussian model; it does not imply that every scoring procedure produces a conditional mean or the same uncertainty.

Scores also depend on the factor coordinates. If a model is transformed consistently, the latent coordinates change with it; estimates should be interpreted in the coordinates of the reported solution. Do not compare score values across solutions until factor alignment, scaling, and the scoring method have been checked.

A Rotation Criterion Is Not a Fit Test

A rotation criterion evaluates an orientation of the loading matrix, often rewarding a feature such as simple structure. Fit assessment instead asks how adequately the model reproduces the data, for example through the discrepancy between the observed covariance matrix and \(\Sigma=\Lambda\Phi\Lambda^\mathsf{T}+\Psi\), or through an inferential criterion appropriate to the fitted models. These are different questions. A low rotation-criterion value is not a test statistic for covariance fit and does not establish that the model is adequate.

Theorem: An orthogonal rotation changes loading coordinates but leaves the model-implied covariance unchanged. In particular, an exactly equivalent orthogonal rotation cannot improve or worsen covariance fit.

Proof. Let \(Q\) be orthogonal, so \(Q Q^\mathsf{T}=I\), and set \(\Lambda^*=\Lambda Q\). For an orthogonal factor model, \(\Phi=I\); the rotated common covariance is

$$ \Lambda^*(\Lambda^*)^\mathsf{T} =\Lambda Q Q^\mathsf{T}\Lambda^\mathsf{T} =\Lambda\Lambda^\mathsf{T}. $$

Keeping \(\Psi\) fixed therefore gives the same full covariance matrix before and after rotation. Any fit measure depending only on the fitted covariance and the data has the same value for these equivalent representations. A rotation criterion may still change because it evaluates the loading orientation, not covariance reproduction. \(\square\)

The corresponding oblique coordinate change also preserves the common covariance when \(\Phi\) is transformed consistently, as shown above. By contrast, comparing models with different restrictions, factor counts, or fitted covariance matrices is a model-comparison problem. A rotation criterion cannot substitute for a fit diagnostic or an appropriate comparison procedure.

Common Mistakes and Corrections

  • Wrong: “The oblique communality is the sum of squared pattern loadings.” Correction: Use \(\lambda_i^\mathsf{T}\Phi\lambda_i\), including the factor-correlation cross-products.
  • Wrong: “The first factor in each solution must represent the same construct.” Correction: Check permutations and signs, then assess whether a more general alignment is warranted.
  • Wrong: “A person’s estimated score is their known latent value.” Correction: Call it an estimate and acknowledge the uncertainty and scoring assumptions.
  • Wrong: “The rotation with the best simple-structure criterion fits the data best.” Correction: Assess covariance fit separately; a criterion for orientation is not a fit test.

The practical lesson is to keep four objects distinct: the loading coordinates, the common covariance they represent, the estimated scores in those coordinates, and the diagnostics used to judge fit. Confusing them can turn an equivalent representation into an apparent disagreement, or a convenient interpretation into an unjustified model claim.

Check Your Understanding

Use the model equations and distinctions in this tutorial to answer the following.

  1. For an oblique two-factor solution, what extra terms appear in a communality beyond the sum of squared pattern loadings?
  2. Why can two identical factor structures have low column-by-column congruence before alignment?
  3. In the worked Gaussian score example, what does the conditional variance say that the conditional mean alone does not?
  4. What changes, and what remains unchanged, under an orthogonal rotation of an orthogonal factor model?
  5. Why is a rotation criterion not evidence that a model reproduces the observed covariance well?
Clinical Biostats