Skip to content

Eigenfactor risk adjustment

How eigenfactor risk adjustment calibrates selected covariance directions, why volatility multipliers are squared, and why corrected variances do not recover the true covariance.

Eigenfactor risk adjustment recalibrates the variances assigned to the eigenvector directions of an estimated factor covariance matrix to reduce risk forecast bias caused by sampling error. Its practical target is a portfolio optimiser's tendency to concentrate in directions whose estimated risk is unusually low. It changes the risk matrix used for construction, rather than adding one blanket buffer to the final portfolio volatility.
Menchero, Wang and Orr develop this approach in Improving Risk Forecasts for Optimized Portfolios. Their eigenportfolios reveal strong underprediction in low estimated variance directions and some overprediction at the high end. These are empirical findings under their research design, not universal multipliers.

Why an unbiased matrix can mislead an optimiser

For a covariance estimate F^\widehat F satisfying E[F^]=F\mathbb E[\widehat F]=F, a fixed, nonrandom exposure vector a\mathbf a has
E[a⊤F^a]=a⊤Fa.\mathbb E[\mathbf a^\top\widehat F\mathbf a]=\mathbf a^\top F\mathbf a.

That equality does not establish unbiased risk forecasts for a direction u^\widehat{\mathbf u} selected using the same estimated matrix. The direction and estimation errors then depend on each other. Selecting the smallest estimated variance preferentially finds errors that make risk look cheap. Unbiased variance estimation also does not imply unbiased volatility estimation, because taking a square root is nonlinear.

F^=Udiag⁡(d1,…,dK)U⊤,d1≤⋯≤dK.\widehat F=U\operatorname{diag}(d_1,\ldots,d_K)U^\top, \qquad d_1\le\cdots\le d_K.
The columns uk\mathbf u_k are unit-length factor exposure directions. Their estimated variances are dkd_k, and they are mutually uncorrelated under this estimated matrix. The smallest eigenvector minimises u⊤F^u\mathbf u^\top\widehat F\mathbf u subject to u⊤u=1\mathbf u^\top\mathbf u=1. This is a unit-norm problem; a budget constraint such as weights summing to one defines a different portfolio problem. The authors' May 2011 research paper, “Eigenfactors”, explains this distinction.

What is estimated in the actual method

The historical USE4 methodology, Appendix B, describes a plug-in Monte Carlo procedure:
  1. Use the observed covariance as the generating covariance inside the simulation.
  2. Generate synthetic factor histories; refit the same covariance estimator.
  3. Diagonalise each refit; compare each eigenportfolio's forecast variance with its variance under the generating covariance.
  4. Aggregate comparisons by eigenfactor rank into volatility correction factors vkv_k.

The reconstruction is

Fadj=Udiag⁡(v12d1,…,vK2dK)U⊤.F_{\mathrm{adj}} =U\operatorname{diag}(v_1^2d_1,\ldots,v_K^2d_K)U^\top.
A volatility multiplier must be squared when applied to a variance. The original eigenvectors are retained. Known simulation truth is not knowledge of the real population covariance.

Two corrected directions still leave a missing covariance

Consider a stipulated oracle illustration, with all entries in common squared-return units at one fixed horizon:
F^=(0.005000.050),Ftrue=(0.0250.0150.0150.025).\widehat F=\begin{pmatrix}0.005&0\\0&0.050\end{pmatrix}, \qquad F_{\mathrm{true}}=\begin{pmatrix}0.025&0.015\\0.015&0.025\end{pmatrix}.

The estimated eigenvectors are the coordinate axes, so U=IU=I. Supply the true matrix only to demonstrate the arithmetic; no return sample or Monte Carlo fit produces these numbers.

Both coordinate directions have true variance 0.0250.025. Their variance ratios and volatility multipliers are therefore

r1=0.025/0.005=5,v1=5≈2.2361,r2=0.025/0.050=0.5,v2=0.5≈0.7071.\begin{aligned} r_1&=0.025/0.005=5,&v_1&=\sqrt5\approx2.2361,\\ r_2&=0.025/0.050=0.5,&v_2&=\sqrt{0.5}\approx0.7071. \end{aligned}

The first direction's variance rises fivefold; the second's halves. Applying vk2=rkv_k^2=r_k gives

Fadj=(0.025000.025).F_{\mathrm{adj}}=\begin{pmatrix}0.025&0\\0&0.025\end{pmatrix}.
Stipulated oracle example. Direction 1 variance moves from 0.005 to 0.025 and direction 2 from 0.050 to 0.025. Adjusted cross-covariance remains zero although stipulated truth is 0.015. Along the unit-length sum direction, adjusted variance is 0.025 versus true variance 0.040.
Stipulated oracle example. Direction 1 variance moves from 0.005 to 0.025 and direction 2 from 0.050 to 0.025. Adjusted cross-covariance remains zero although stipulated truth is 0.015. Along the unit-length sum direction, adjusted variance is 0.025 versus true variance 0.040.
Original stipulated oracle example: the two directional variances match, while the cross-covariance remains missing. These are supplied calibration ratios, not estimates from a vendor simulation.
This adjustment does not recover FtrueF_{\mathrm{true}}. Its off-diagonal entry stays zero. For the unit-length combined direction a=(1,1)⊤/2\mathbf a=(1,1)^\top/\sqrt2,
a⊤Fadja=0.025,a⊤Ftruea=0.025+0.025+2(0.015)2=0.040.\mathbf a^\top F_{\mathrm{adj}}\mathbf a=0.025, \qquad \mathbf a^\top F_{\mathrm{true}}\mathbf a =\frac{0.025+0.025+2(0.015)}2=0.040.

Matching the two chosen directional variances does not match every combination. The true matrix has eigenvalues 0.0100.010 and 0.0400.040; they are not the two 0.0250.025 directional variances measured along the estimated axes. The axes themselves differ from the true eigendirections. The trace also changes, from 0.0550.055 to 0.0500.050: trace preservation is not imposed here.

Read the calibration in context

This is distinct from covariance eigenvalue clipping, which can pool a designated noise group at one average, and from fitted nonlinear covariance shrinkage. Equal adjusted eigenvalues in this example are a consequence of the supplied ratios, not a clipping rule.
The USE4 document assumes normal, stationary simulated returns and discusses the limits of that approximation. It also distinguishes an additional empirically scaled correction from the milder simulated correction adopted in that historical model. These August 2011 choices are not a specification or performance guarantee for today's products. Appendix B, equations B7–B10, gives the details.

A correction trained on one estimated covariance and sampling scheme cannot establish that every future portfolio forecast is calibrated. Constraints, changing covariance and factor-model omissions still matter. Evaluate risk forecasts on portfolios formed with information available at the time.

For the portfolio construction step, Mean–Variance Optimisation derives inverse-covariance weights and the fully invested constraint. Its public page is a preview; the full lesson requires membership.

Sources

Last updated
  • factor risk model
  • covariance estimation
  • eigenvectors
  • portfolio optimisation
  • risk calibration
← Back to the glossaryBack to top ↑