Clustered standard errors

Clustered standard errors (or Liang-Zeger standard errors) are measurements that estimate the standard error of a regression parameter in settings where observations may be subdivided into smaller-sized groups ("clusters") and where the sampling and/or treatment assignment is correlated within each group. Clustered standard errors are widely used in a variety of applied econometric settings, including difference-in-differences or experiments.

Analogous to how Huber-White standard errors are consistent in the presence of heteroscedasticity and Newey–West standard errors are consistent in the presence of accurately-modeled autocorrelation, clustered standard errors are consistent in the presence of cluster-based sampling or treatment assignment. Clustered standard errors are often justified by possible correlation in modeling residuals within each cluster; while recent work suggests that this is not the precise justification behind clustering, it may be pedagogically useful.

Intuitive motivation
Clustered standard errors are often useful when treatment is assigned at the level of a cluster instead of at the individual level. For example, suppose that an educational researcher wants to discover whether a new teaching technique improves student test scores. She therefore assigns teachers in "treated" classrooms to try this new technique, while leaving "control" classrooms unaffected. When analyzing her results, she may want to keep the data at the student level (for example, to control for student-level observable characteristics). However, when estimating the standard error or confidence interval of her statistical model, she realizes that classical or even heteroscedasticity-robust standard errors are inappropriate because student test scores within each class are not independently distributed. Instead, students in classes with better teachers have especially high test scores (regardless of whether they receive the experimental treatment) while students in classes with worse teachers have especially low test scores. The researcher can cluster her standard errors at the level of a classroom to account for this aspect of her experiment.

While this example is very specific, similar issues arise in a wide variety of settings. For example, in many panel data settings (such as difference-in-differences) clustering often offers a simple and effective way to account for non-independence between periods within each unit (sometimes referred to as "autocorrelation in residuals"). Another common and logically distinct justification for clustering arises when a full population cannot be randomly sampled, and so instead clusters are sampled and then units are randomized within cluster. In this case, clustered standard errors account for the uncertainty driven by the fact that the researcher does not observe large parts of the population of interest.

Mathematical motivation
A useful mathematical illustration comes from the case of one-way clustering in an ordinary least squares (OLS) model. Consider a simple model with N observations that are subdivided in C clusters. Let $$Y$$ be an $$n \times 1$$ vector of outcomes, $$X$$ a $$n \times m$$ matrix of covariates, $$\beta$$ an $$m \times 1$$ vector of unknown parameters, and $$e$$ an $$n \times 1$$ vector of unexplained residuals:


 * $$Y = X\beta + e$$

As is standard with OLS models, we minimize the sum of squared residuals $$e$$ to get an estimate $$\hat{\beta}$$:


 * $$ \min_\beta (Y-X\beta)^2$$


 * $$\Rightarrow X'(Y-X\hat{\beta}) = 0$$


 * $$\Rightarrow \hat{\beta} = (X'X)^{-1} X'Y$$

From there, we can derive the classic "sandwich" estimator:


 * $$V(\hat{\beta}) = V((X'X)^{-1}X'Y) = V(\beta + (X'X)^{-1}X'e) = V((X'X)^{-1}X'e) = (X'X)^{-1}X'ee'X(X'X)^{-1} $$

Denoting $$\Omega \equiv ee'$$ yields a potentially more familiar form


 * $$V(\hat{\beta}) = (X'X)^{-1}X'\Omega X (X'X)^{-1}$$

While one can develop a plug-in estimator by defining $$\hat{e} \equiv Y - X\hat{\beta}$$ and letting $$\hat{\Omega} \equiv \hat{e} \hat{e}'$$, this completely flexible estimator will not converge to $$V({\hat{\beta}})$$ as $$N \rightarrow \infty$$. Given the assumptions that a practitioner deems as reasonable, different types of standard errors solve this problem in different ways. For example, classic homoskedastic standard errors assume that $$\Omega$$ is diagonal with identical elements $$\sigma^2$$, which simplifies the expression for $$V(\hat{\beta}) = \sigma^2 (X'X)^{-1}$$. Huber-White standard errors assume $$\Omega$$ is diagonal but that the diagonal value varies, while other types of standard errors (e.g. Newey–West, Moulton SEs, Conley spatial SEs) make other restrictions on the form of this matrix to reduce the number of parameters that the practitioner needs to estimate.

Clustered standard errors assume that $$\Omega$$ is block-diagonal according to the clusters in the sample, with unrestricted values in each block but zeros elsewhere. In this case, one can define $$X_c$$ and $$\Omega_c$$ as the within-block analogues of $$X$$ and $$\Omega$$ and derive the following mathematical fact:


 * $$X'\Omega X = \sum_c X'_c \Omega_c X_c$$

By constructing plug-in matrices $$\hat{\Omega}_c$$, one can form an estimator for $$V(\hat{\beta})$$ that is consistent as the number of clusters $$c$$ becomes large. While no specific number of clusters is statistically proven to be sufficient, practitioners often cite a number in the range of 30-50 and are comfortable using clustered standard errors when the number of clusters exceeds that threshold. Alternatively, finite-sample modifications are also typically used, to reduce downwards bias $$V(\hat{\beta})$$ due to finite C. Thus, the resulting bias corrected estimator follows the form:


 * $$\hat{V}(\hat{\beta}) = (X'X)^{-1} \sum_c X'_c \hat{\Omega}_c X_c (X'X)^{-1} \frac{C}{C-1} \frac{n - 1}{n - k}.$$