Tutorials › Econometrics (Graduate) › Inference After Estimating Common Factors

Panel Dependence and Inference · Tutorial 334 of 952

Inference After Estimating Common Factors

See how factor estimation enters panel uncertainty, and compare factor-adjusted covariance estimates with inference based on common-correlated-effects regressions.

Advanced 11 min read

What You'll Learn

  • Separate common-factor variation from idiosyncratic error in a panel model.
  • Derive how a shared shock contributes to the variance of a coefficient estimator.
  • Identify when treating estimated factors and loadings as known understates uncertainty.
  • Work through a numerical adjustment for factor-estimation error.
  • Compare factor-adjusted covariance inference with common-correlated-effects inference.
  • Recognize the assumptions and limitations behind each covariance estimate.

Why Factor Estimation Changes Inference

A panel regression can have a precisely estimated coefficient conditional on the observed regressors and still have an understated standard error. One reason is that outcomes across units may share unobserved shocks. Another is that the common factors and their loadings are estimated rather than known. Inference must address both sources of uncertainty: the dependence left by common shocks and the extra variation introduced by estimating the common component.

This tutorial focuses on inference after common factors have been estimated. It distinguishes a factor-adjusted covariance calculation from common-correlated-effects (CCE) inference. The methods are related in purpose, but they do not make identical assumptions or use identical regressions. Neither method automatically solves every dependence problem: the remaining errors, the panel's asymptotic regime, and the factor-estimation procedure all matter.

Model and Sources of Uncertainty

Consider a panel with units \(i=1,\ldots,N\) and periods \(t=1,\ldots,T\). A linear model with \(r\) common factors is

$$ y_{it}=x_{it}'\beta+\lambda_i'f_t+\varepsilon_{it}. $$

Here \(y_{it}\) is the outcome, \(x_{it}\) is a vector of observed regressors, \(\beta\) is the coefficient vector of interest, \(f_t\) is an \(r\)-vector of common factors, \(\lambda_i\) is unit \(i\)'s loading vector, and \(\varepsilon_{it}\) is the idiosyncratic error. The factor component \(\lambda_i'f_t\) induces dependence across units at a given time whenever units load on the same factor. The factor and loading scales are not separately identified without a normalization; inference concerns the regression coefficient and fitted common component, not an arbitrary factor rotation.

Working conditions for the calculations below: The conditional mean model is correctly specified; regressors have sufficient variation to identify \(\beta\); the common component is adequately represented by the chosen factor dimension; and the idiosyncratic errors satisfy the dependence and moment conditions required by the chosen covariance estimator. Asymptotic claims additionally require a stated growth regime for \(N\) and \(T\), and suitable factor strength and regularity. These conditions are not guaranteed merely by fitting a factor model.

There are two distinct uncertainty questions. First, even if \(f_t\) and \(\lambda_i\) were known, common factors would affect the variance of a coefficient estimator unless their contribution to its score were removed or accounted for. Second, in practice the factors and often the loadings are estimated. Their estimation error can affect the coefficient estimate and its covariance. A covariance calculation that conditions on the estimated factor values as if they were fixed addresses the first question only insofar as its residual covariance model is adequate; it does not, by itself, account for the second.

How a Common Shock Enters the Variance

A short derivation shows why cross-sectional independence is consequential. For illustration, suppose \(x_{it}\) is scalar and treated as fixed, and consider the no-intercept least-squares estimator that ignores the factors:

$$ \hat{\beta}=\frac{\sum_{t=1}^{T}\sum_{i=1}^{N}x_{it}y_{it}} {\sum_{t=1}^{T}\sum_{i=1}^{N}x_{it}^{2}}. $$

Write \(D=\sum_{t,i}x_{it}^{2}\). Substituting the model gives

$$ \hat{\beta}-\beta = \frac{\sum_t f_t\left(\sum_i x_{it}\lambda_i\right) +\sum_{t,i}x_{it}\varepsilon_{it}}{D}. $$

If the scalar factors are independent across time with variance \(\sigma_f^2\), and independent of the idiosyncratic errors, then the common-factor part contributes

$$ \operatorname{Var}_{f}(\hat{\beta}\mid X,\lambda) = \frac{\sigma_f^2\sum_t\left(\sum_i x_{it}\lambda_i\right)^2}{D^2}. $$

This follows because the coefficient on \(f_t\) in the numerator is \(\sum_i x_{it}\lambda_i\), and independent factors have zero cross-time covariances. If the idiosyncratic errors are also independent with variance \(\sigma_\varepsilon^2\), their contribution is \(\sigma_\varepsilon^2/D\). A covariance estimate that treats all \(NT\) observations as independent captures neither the common-factor term nor, generally, other serial or cross-sectional dependence.

Worked Example: A Shared Factor Inflates the Variance

Suppose two units are observed for two periods. In each period \(x_{1t}=1\), \(x_{2t}=2\), and both units have loading \(\lambda_i=1\). Let the factors be independent across periods with variance 1, and let the idiosyncratic errors be independent with variance 1, also independent of the factors. Use the no-intercept estimator above.

The denominator is \(D=2(1^2+2^2)=10\). In each period the factor's score coefficient is \(\sum_i x_{it}\lambda_i=1+2=3\). Therefore the factor contribution to the variance is

$$ \frac{1^2(3^2+3^2)}{10^2} =\frac{18}{100} =0.18. $$

The idiosyncratic contribution is \(1/D=0.10\) if computed using the general variance formula for this estimator with repeated design values? More directly, its numerator variance is \(2(1^2+2^2)=10\), so its contribution is \(10/10^2=0.10\). Thus the total variance is \(0.28\), and the standard error is \(\sqrt{0.28}=0.5292\). By contrast, an independence calculation that ignores the factor reports variance \(0.10\), or standard error \(\sqrt{0.10}=0.3162\). The difference comes from the shared shock, not from a change in the point-estimate formula.

This is a variance illustration, not a recommendation to use the factor-ignorant estimator. If the common component is modeled and removed, inference must instead address whether the estimated factors and loadings have been treated appropriately.

Factor-Adjusted Covariance Inference

A factor-adjusted analysis estimates the common component and uses the adjusted regression to estimate \(\beta\). One common route is to estimate latent factors from the panel, for example by principal components, and then include the estimated factors in the regression. Another uses an estimated common component to adjust the outcome or residuals. In either case, the covariance calculation must match the estimating procedure: a sandwich covariance based on adjusted scores can account for remaining heteroskedasticity or dependence only to the extent its score aggregation and covariance estimator permit.

For a regression with estimated factors \(\hat f_t\), the plug-in calculation treats those factors as fixed regressors. It can be a useful first-order approximation in settings where factor estimates are sufficiently accurate and their effect on \(\hat\beta\) is asymptotically negligible. That conclusion requires conditions; it is not automatic. If factor estimation error affects the coefficient's first-order distribution, the covariance needs an additional correction. The same issue applies to estimated loadings and to uncertainty about the estimated common component.

Key distinction: A factor-adjusted covariance estimate describes uncertainty after accounting for a modeled common component. A plug-in version conditions on estimated factors and loadings. It represents unconditional coefficient uncertainty only when factor-estimation effects are negligible under the relevant asymptotic conditions, or when those effects are explicitly included.

An analytic correction can incorporate the factor-estimation influence into the coefficient's variance. A practical alternative, when justified by the design and asymptotics, is a resampling procedure that repeats the factor estimation, loading estimation, adjustment, and coefficient estimation in every resample. Re-estimating the factors matters: resampling only the final regression residuals while keeping the original factors fixed does not reproduce uncertainty arising from factor estimation. The resampling scheme must also preserve the relevant dependence structure; an ordinary independent-observation bootstrap is not automatically valid for a dependent panel.

Worked Example: Adding a Factor-Estimation Contribution

Return to two units and two periods with \(x_{1t}=1\), \(x_{2t}=2\), and loadings equal to 1. Suppose the true factor has been removed, leaving independent idiosyncratic errors of variance 1. The factor-adjusted coefficient has denominator \(D=10\), so its idiosyncratic variance is \(10/100=0.10\). Now isolate factor-estimation error by assuming the estimated factor is \(\hat f_t=f_t+u_t\), where \(u_t\) is independent across periods, independent of the idiosyncratic errors, and has variance \(0.04\). Treat the loadings as known for this illustration.

The error in the fitted common component contributes a score \(3u_t\) in each period. Its variance contribution is

$$ \frac{0.04(3^2+3^2)}{10^2} =\frac{0.72}{100} =0.0072. $$

The total variance under these simplified assumptions is \(0.10+0.0072=0.1072\), giving standard error \(\sqrt{0.1072}=0.3274\). Treating the factors as known would give standard error \(\sqrt{0.10}=0.3162\), which is smaller. This calculation isolates uncertainty from factor estimation; it is not a universal correction formula. If loadings are estimated, if factor-estimation errors are serially dependent, or if those errors covary with the coefficient score, cross-terms and additional contributions must also be considered.

Comparing Factor Adjustment with CCE Inference

CCE estimation handles unobserved common influences differently. Rather than first estimating latent factors, a CCE regression augments the unit-level regression with cross-sectional averages of the outcome and regressors, such as \(\bar y_t=N^{-1}\sum_i y_{it}\) and \(\bar x_t=N^{-1}\sum_i x_{it}\). Under its identifying conditions, these averages span or proxy the relevant common-factor space as \(N\) grows. Unit-specific CCE regressions can allow heterogeneous factor loadings and heterogeneous regression coefficients; pooled variants impose additional restrictions appropriate to the target parameter.

CCE inference is based on the augmented regression's scores and a covariance estimator suited to the panel's remaining dependence and asymptotic regime. It does not estimate and then plug in a particular latent factor series in the same way as principal-components factor adjustment. However, cross-sectional averages are themselves data-derived regressors, and their use does not make inference assumption-free. Validity relies on adequate factor proxies, sufficient cross-sectional information, identification, and a covariance method appropriate to residual dependence. A CCE covariance estimate and a factor-adjusted covariance estimate can therefore differ even when both aim to account for common shocks.

FeatureFactor-adjusted inferenceCCE inference
Common-component representationEstimated latent factors and loadings, or an estimated common componentCross-sectional averages used as proxies for common influences
Primary uncertainty concernResidual dependence after adjustment and uncertainty from estimating factors/loadingsValidity of the proxy augmentation and covariance for the augmented regression
Inference requirementJustify a plug-in approximation or include factor-estimation effectsJustify the proxy and coefficient assumptions, plus the selected covariance estimator

Worked Example: Comparing Two Covariance Calculations

Consider an illustrative panel regression summary with one coefficient, \(\hat\beta=0.40\), and scalar bread \(H=10\). For a simple unit-clustered sandwich calculation, let the unit-level score contributions after factor adjustment be \(q_i=(1,-1,2)\). For the CCE-augmented regression, suppose the corresponding score contributions are \(q_i^{CCE}=(1.2,-0.8,1.6)\), with the same bread for easy comparison. Use the uncorrected sandwich formula \(H^{-2}\sum_i q_i^2\).

For factor adjustment, the score meat is \(1^2+(-1)^2+2^2=6\), so

$$ \widehat{\operatorname{Var}}(\hat\beta)_{\mathrm{FA}} =\frac{6}{10^2}=0.06, \qquad SE_{\mathrm{FA}}=\sqrt{0.06}=0.2449. $$

For CCE, the meat is \(1.2^2+(-0.8)^2+1.6^2=1.44+0.64+2.56=4.64\), giving

$$ \widehat{\operatorname{Var}}(\hat\beta)_{\mathrm{CCE}} =\frac{4.64}{10^2}=0.0464, \qquad SE_{\mathrm{CCE}}=\sqrt{0.0464}=0.2154. $$

The smaller CCE standard error in this arithmetic example is not evidence that CCE is generally more efficient or preferable. The two score sets come from different adjusted regressions, and validity depends on the assumptions behind each method. Moreover, three clusters are far too few for this calculation to support reliable large-sample inference; the numbers illustrate how sandwich meat affects a standard error, not a recommended sample size or a calibrated p-value.

Choosing and Reporting the Inference

The central practical question is not simply whether a factor adjustment was performed. It is whether the reported covariance estimate reflects the complete estimation procedure and the dependence that remains. Before interpreting a confidence interval or test, report the factor or proxy method, the factor dimension or augmentation, the covariance estimator, the relevant clustering or time-dependence treatment, and the asymptotic justification.

A plug-in factor-adjusted standard error is most defensible when the factor-estimation contribution is asymptotically negligible for the coefficient under the stated growth regime and regularity conditions. When that is not justified, use a covariance correction or a suitable procedure that re-estimates the common component. CCE can be attractive when cross-sectional averages proxy the factor space under its conditions, but those conditions and the covariance treatment still require scrutiny. The two approaches are alternatives with distinct assumptions, not interchangeable labels for the same calculation.

  • Mistake: “The estimated factors explain common shocks, so an ordinary regression standard error is valid.”
  • Correction: The adjustment must be paired with a covariance estimator valid for the remaining errors, and factor-estimation uncertainty must be negligible or included.
  • Mistake: “CCE uses cross-sectional averages, so it automatically handles every form of cross-sectional dependence.”
  • Correction: CCE relies on those averages adequately proxying the relevant common factors under specified conditions; unmodeled dependence or unsuitable asymptotics can still invalidate inference.
Key takeaway: Common-factor adjustment changes both the regression's dependence structure and the uncertainty calculation. State what was estimated, what covariance estimator was used, and why factor-estimation effects and residual dependence are covered.

Check Your Understanding

Answer each question by referring to the model, assumptions, or calculations above.

  1. In the common-factor variance derivation, what part of the coefficient's score creates dependence across units?
  2. Why can a plug-in covariance estimate that treats estimated factors as fixed understate uncertainty?
  3. In the second worked example, what extra assumption allowed factor-estimation and idiosyncratic variance contributions to be added directly?
  4. How does CCE represent common influences differently from latent-factor adjustment?
  5. Why does the smaller standard error in the CCE arithmetic example not establish that CCE is generally preferable?
Clinical Biostats