Why Orthogonal Criteria Differ
An orthogonal rotation changes how a fitted loading matrix allocates loadings across factors while preserving the model's common covariance representation. The rotation does not change a variable's communality, but it can make the loading pattern easier—or harder—to describe. Quartimax, varimax, and equamax are different rules for preferring some rotated patterns over others.
All three belong to the orthomax family. Their differences are easiest to understand by separating two questions: does a criterion reward a variable for loading strongly on just one factor, or does it reward each factor for having a few strong loadings and many weak ones? Quartimax emphasizes the first question, varimax the second, and equamax combines the two weightings.
This states the convention explicitly because criterion formulas can differ in their normalizing constants or in whether an algorithm minimizes the negative of the displayed objective. Here, a rotation with a larger value of the stated objective is preferred by that criterion. The loading matrix is not row-normalized before calculating these formulas.
What the Three Criteria Reward
For variable \(i\), its communality is \(h_i^2=\sum_j\lambda_{ij}^2\). Under an orthogonal rotation, its row of loadings is multiplied by an orthogonal matrix, so its squared length is unchanged:
Thus, a criterion cannot increase a variable's total communality by rotating; it can only change how that fixed amount is distributed across its factor loadings. Quartimax directly rewards concentration within each row through fourth powers. In contrast, varimax evaluates squared loadings column by column, rewarding a factor with a mixture of relatively large and small loadings across variables.
Proof. Put \(x_j=\lambda_{ij}^2\), so \(x_j\geq 0\) and \(\sum_jx_j=h_i^2\). The quartimax contribution from row \(i\) is \(\sum_jx_j^2\). Since \[ \left(\sum_jx_j\right)^2=\sum_jx_j^2+2\sum_{j<k}x_jx_k\geq \sum_jx_j^2, \] it is at most \(h_i^4\). Equality holds when every product \(x_jx_k\) for \(j\ne k\) is zero—that is, when at most one loading in the row is nonzero.
For factor \(j\), let \(x_i=\lambda_{ij}^2\) and \(c_j=\sum_i x_i\), which is fixed for this comparison. Its varimax contribution is \(\sum_i x_i^2-c_j^2/p\). The subtracted term is constant when \(c_j\) is fixed, so maximizing the contribution is equivalent to maximizing \(\sum_i x_i^2\). For nonnegative \(x_i\) with fixed sum, concentrating the sum increases this quantity: transferring an amount \(\delta\) from \(x_b\) to \(x_a\), where \(x_a\geq x_b\), changes the sum of squares by \[ (x_a+\delta)^2+(x_b-\delta)^2-x_a^2-x_b^2 =2\delta(x_a-x_b)+2\delta^2\geq 0 \] whenever \(0\leq\delta\leq x_b\). This proves the stated preference. \(\square\)
The proof describes the criteria's tendencies, not a promise that an entire matrix can attain all of these separate maxima at once. A rotation couples the rows and columns: the same orthogonal transformation must act on every variable. Varimax's columnwise reward commonly produces factors with a few strong loadings and many weak ones. Quartimax can instead favor assigning each variable mainly to one factor, but this can yield a broad general factor on which many variables load. Equamax uses an intermediate orthomax weighting; it does not guarantee a perfect balance or a unique, universally best pattern.
How Equamax Fits the Orthomax Family
The parameter \(\gamma\) controls the penalty for a factor column with a large sum of squared loadings. When \(\gamma=0\), there is no such penalty, giving quartimax. When \(\gamma=1\), the penalty centers the squared loadings within each column, giving varimax. Equamax sets \(\gamma=p/(2m)\), so its weighting depends on the number of variables relative to the number of factors.
To derive the displayed form, substitute \(\gamma=p/(2m)\) into the orthomax definition. Its column penalty coefficient becomes \(\gamma/p=(p/(2m))/p=1/(2m)\). The first term sums the fourth powers over all entries, so the expression follows directly. In particular, equamax does not always define a numerically distinct criterion: when \(p=2m\), its parameter equals the varimax parameter.
Worked Example: Calculating All Three Criteria
Worked Example: Calculating All Three Criteria
Consider the \(6\times2\) loading matrix $$ \Lambda= \begin{pmatrix} 0.8&0.1\\ 0.7&0.2\\ 0.1&0.8\\ 0.2&0.7\\ 0.6&0.3\\ 0.3&0.6 \end{pmatrix}. $$ There are \(p=6\) variables and \(m=2\) factors. Calculate the quartimax, varimax, and equamax objectives.
Quartimax. Sum the fourth powers of every entry. The first four rows contribute \(0.4097, 0.2417, 0.4097,\) and \(0.2417\). The last two each contribute \(0.6^4+0.3^4=0.1377\). Therefore, $$ Q_0=0.4097+0.2417+0.4097+0.2417+0.1377+0.1377=1.5782. $$
Varimax. Each column has squared-loading sum \(1.63\). For example, the first column gives \(0.8^2+0.7^2+0.1^2+0.2^2+0.6^2+0.3^2=1.63\). The sum of its fourth powers is \(0.7891\); the second column has the same values. Each column's contribution is \(0.7891-(1/6)(1.63)^2=0.7891-0.4428167=0.3462833\), so $$ Q_1=2(0.3462833)\approx0.69257. $$
Equamax. Here \(\gamma=p/(2m)=6/4=1.5\). The coefficient of each squared column sum is \(\gamma/p=1.5/6=0.25\). Thus, $$ Q_{1.5}=2\{0.7891-0.25(1.63)^2\} =2(0.124875)=0.24975. $$
The three values differ because the criteria apply different column penalties, not because they describe different fitted covariance matrices. These values are useful for evaluating orientations under their respective objectives; their magnitudes should not be compared as if the criteria shared a common scale.
Worked Example: Quartimax and a Variable's Loading Pattern
Worked Example: Quartimax and a Variable's Loading Pattern
Suppose a variable has communality \(0.64\) and two possible loading patterns are \((0.8,0)\) and \((\sqrt{0.32},\sqrt{0.32})\). Both have squared loadings summing to \(0.64\). Their quartimax row contributions are $$ 0.8^4+0^4=0.4096 $$ and $$ (\sqrt{0.32})^4+(\sqrt{0.32})^4=0.32^2+0.32^2=0.2048, $$ respectively. Quartimax favors the first pattern because the same communality is concentrated in one loading. This rowwise comparison illustrates the criterion's preference; by itself, it does not establish that a single rotation can give every variable its preferred pattern simultaneously.
Worked Example: Varimax and Factor Differentiation
Worked Example: Varimax and Factor Differentiation
Compare two possible sets of squared loadings in a factor column with total squared loading \(0.60\), using \(p=3\). In the first pattern, the squared loadings are \((0.54,0.06,0)\). The sum of the loadings' fourth powers is \(0.54^2+0.06^2=0.2952\), so the varimax column contribution is $$ 0.2952-\frac{1}{3}(0.60)^2=0.2952-0.12=0.1752. $$ In the second pattern, the squared loadings are \((0.20,0.20,0.20)\). Their fourth-power sum is \(3(0.20)^2=0.12\), giving $$ 0.12-\frac{1}{3}(0.60)^2=0. $$ For this fixed column sum of squares, varimax favors the more unequal pattern. In terms of loadings, the first pattern has two relatively strong entries and one zero, whereas the second spreads the column's squared loading evenly across all three variables.
Interpreting and Comparing the Criteria
All three objectives are unchanged by permuting factor columns or changing the sign of any factor: fourth powers and column sums of squares are unaffected. General orthogonal rotations, however, can change the objectives. Although the row communalities are preserved, the entrywise fourth-power sum and the distribution of squared loadings within columns need not be.
The practical choice is about the pattern of simplicity desired, not about which criterion has the largest numerical value. Quartimax puts direct emphasis on simple rows—variables with relatively few substantial loadings—and can encourage a general factor. Varimax focuses on simple columns—factors with relatively few substantial loadings. Equamax retains the quartimax fourth-power term while applying a column penalty whose strength depends on \(p\) and \(m\).
A common mistake is to say that quartimax, varimax, and equamax all simply “maximize sparsity,” then treat their objective values as interchangeable. That skips the essential distinction: quartimax has no column penalty, while varimax and equamax do. The correction is to calculate the specified objective, check its parameter and normalization, and describe the pattern it rewards. A larger quartimax value is not evidence that a solution is better by the varimax criterion.
Check Your Understanding
Use the definitions and examples above to check both the calculation and the interpretation.
- For a \(p\times m\) loading matrix, what value of \(\gamma\) defines equamax in the orthomax convention used here?
- Why does orthogonal rotation preserve each variable's communality even though its individual loadings can change?
- A variable has fixed communality \(h_i^2\). What loading pattern maximizes its quartimax row contribution, and what is that maximum?
- For a fixed factor-column sum of squared loadings, why does the varimax column contribution increase as the squared loadings become more unequal?
- When do equamax and varimax have the same parameter in this convention?
- Why should the numerical value of quartimax not be compared directly with the numerical value of varimax to decide which rotated solution is preferable?