Tutorials › Statistical Theory III - Theories of Inference, Multivariate Methods, and Contingency Tables › Conditional and Unconditional Regression Limits

Random-Design Regression Asymptotics · Tutorial 577 of 972

Conditional and Unconditional Regression Limits

Learn when conditioning on the observed regressors and averaging over random regressors lead to the same asymptotic regression covariance.

Advanced 11 min read

What You'll Learn

  • Distinguish conditional and unconditional asymptotic distributions for random-design OLS.
  • State conditions for a conditional central limit theorem given the observed regressors.
  • Derive the heteroskedastic sandwich covariance from conditional score variances.
  • Explain why the conditional and unconditional limiting covariance matrices agree under stable iid designs.
  • Identify when finite-sample covariance averages need additional integrability conditions.
  • Diagnose a random-environment example in which the conditional limit remains random.

Why distinguish conditional and unconditional limits?

In random-design regression, the regressors are random before data collection. After observing them, however, an analyst may want to describe the sampling uncertainty in the coefficient estimator while treating that particular design as fixed. This is a conditional question. A different question averages over both the errors and the random regressors, giving an unconditional limit.

These two viewpoints need not agree automatically. The conditional variance can depend on the observed design, and an unconditional distribution can be a mixture of different conditional distributions. Under a stable iid design and suitable moment and central limit conditions, the design-based conditional covariance converges to a fixed matrix. It then agrees with the covariance in the unconditional normal limit. The distinction is about the conditions that justify this agreement, not about two different formulas for OLS.

Setup and precise meanings

Let \(X_i\) be a \(p\)-dimensional regressor vector, and suppose \(Y_i=X_i'\beta+u_i\), for \(i=1,\ldots,n\). The pairs \((X_i,u_i)\) are iid. Write \(X_1,\ldots,X_n\) for the observed design, and suppose the conditional errors are independent given these regressors. Define \(\sigma^2(x)=E[u_i^2\mid X_i=x]\).

Definition: The conditional asymptotic distribution of an estimator given the design describes the limiting distribution of its centered, scaled value under the conditional law given \(X_1,\ldots,X_n\). Convergence of these conditional laws in probability means that, for every bounded continuous function \(g\), the conditional expectation of \(g\) converges in probability to its expectation under the stated limiting law. The unconditional asymptotic distribution is the limit under the joint law of the regressors and errors.

Assume \(E[u_i\mid X_i]=0\), \(E\|X_i u_i\|^2<\infty\), and that \(A=E[X_iX_i']\) is positive definite. Put \(B=E[X_iX_i'u_i^2]\), and define \(A_n=n^{-1}\sum_{i=1}^n X_iX_i'\). The conditional Lindeberg condition below is one convenient precise condition for the conditional central limit step:

$$ \frac{1}{n}\sum_{i=1}^n E\!\left[ \|X_i u_i\|^2 \mathbf{1}\{\|X_i u_i\|>\epsilon\sqrt{n}\} \mid X_1,\ldots,X_n \right]\xrightarrow{p}0 \quad\text{for every }\epsilon>0. $$

Also require that the conditional score covariance \(B_n=n^{-1}\sum_{i=1}^n X_iX_i'\sigma^2(X_i)\) converges in probability to \(B\). Under iid sampling this follows from the law of large numbers when \(E\|X_i u_i\|^2<\infty\). The conditional Lindeberg condition follows, for example, if \(E\|X_i u_i\|^{2+\delta}<\infty\) for some \(\delta>0\): the expectation of the nonnegative left side is \(E[\|X_i u_i\|^2\mathbf{1}\{\|X_i u_i\|>\epsilon\sqrt n\}]\), which tends to zero by dominated convergence, so the left side tends to zero in probability by Markov's inequality.

Conditional and unconditional OLS limits

When \(A_n\) is invertible, OLS satisfies the exact identity \(\hat\beta-\beta=A_n^{-1}n^{-1}\sum_i X_i u_i\). The following result makes the two limiting viewpoints precise.

Theorem: Under the stated iid, exogeneity, moment, positive-definiteness, conditional covariance convergence, and conditional Lindeberg conditions, \(\sqrt n(\hat\beta-\beta)\xrightarrow{d}N(0,V)\) unconditionally, and its conditional distribution given \(X_1,\ldots,X_n\) converges in probability to the same \(N(0,V)\), where \(V=A^{-1}BA^{-1}\).

Unconditional derivation. The score vectors \(X_i u_i\) are iid. Their mean is zero because \(E[X_i u_i]=E[X_iE[u_i\mid X_i]]=0\), and their covariance is \(B\). The multivariate central limit theorem therefore gives \(n^{-1/2}\sum_i X_i u_i\xrightarrow{d}N(0,B)\). The law of large numbers gives \(A_n\xrightarrow{p}A\); because \(A\) is positive definite, \(A_n^{-1}\xrightarrow{p}A^{-1}\). Substitution into the exact OLS identity and Slutsky's theorem give the claimed unconditional limit, with covariance \(A^{-1}BA^{-1}\).

Conditional derivation. Given the regressors, the score terms \(X_i u_i\) are independent, have conditional means zero, and have conditional covariance matrices \(X_iX_i'\sigma^2(X_i)\). The conditional Lindeberg condition and \(B_n\xrightarrow{p}B\) give the conditional Lindeberg–Feller limit: the conditional law of \(n^{-1/2}\sum_i X_i u_i\) converges in probability to \(N(0,B)\). Also \(A_n\xrightarrow{p}A\), so \(A_n^{-1}\xrightarrow{p}A^{-1}\). Applying this convergence to the exact OLS identity yields the conditional limit \(N(0,A^{-1}BA^{-1})\).

This also verifies that the unconditional limit follows from the conditional one. For bounded continuous \(g\), conditional convergence in probability implies convergence in \(L^1\) of the conditional expectations of \(g\): their differences are bounded, and a bounded sequence converging to zero in probability has expectations tending to zero. Taking expectations over the design then gives convergence of the unconditional expectations to those under \(N(0,V)\). Thus the two limits agree under the assumptions; this does not assert that the finite-sample conditional law is normal.

Where the covariance formula comes from

The conditional covariance calculation is particularly direct. Conditional on the design, the scaled score has covariance \(B_n\). Thus the conditional covariance of the scaled OLS error is \(A_n^{-1}B_nA_n^{-1}\), whenever \(A_n\) is invertible. Its probability limit is the same \(V\) as in the unconditional central limit argument.

$$ \operatorname{Var}\!\left( \sqrt n(\hat\beta-\beta)\mid X_1,\ldots,X_n \right) =A_n^{-1}B_nA_n^{-1} \xrightarrow{p}A^{-1}BA^{-1}=V. $$

This agreement does not require homoskedasticity. If instead \(E[u_i^2\mid X_i]=\sigma^2\) is constant, then \(B=\sigma^2 A\), so the common limit simplifies to \(V=\sigma^2 A^{-1}\). With heteroskedastic errors, \(B\) retains the variation of the conditional error variance across regressor values.

There is an important qualification about the word “covariance.” The theorem identifies the covariance matrix of the limiting normal distribution and the probability limit of the conditional covariance matrix. It does not, by itself, prove that the expectation of the finite-sample conditional covariance converges to \(V\). That stronger assertion needs sufficient integrability of those random covariance matrices. In particular, convergence in probability alone does not justify exchanging a limit and an expectation.

Worked Example: Heteroskedastic Errors with Two Design Values

Consider scalar regression through the origin, \(Y_i=X_i\beta+u_i\). Let \(P(X_i=1)=P(X_i=2)=1/2\), with \(E[u_i\mid X_i]=0\), conditional error variance \(1\) when \(X_i=1\), and \(4\) when \(X_i=2\). Assume the iid and conditional central limit conditions above.

The population Gram quantity is \(A=E[X_i^2]=(1+4)/2=5/2\). The score variance is \(B=E[X_i^2u_i^2]=(1\cdot1+4\cdot4)/2=17/2\). Consequently the limiting variance of \(\sqrt n(\hat\beta-\beta)\) is \[ V=A^{-2}B=(2/5)^2(17/2)=34/25=1.36. \]

The conditional covariance calculation reaches the same value: \(A_n\) converges to \(5/2\), and \(B_n\) converges to \(17/2\), since each is an average of iid variables with the displayed finite expectations. Hence \(A_n^{-2}B_n\xrightarrow{p}34/25\). The variance is not obtained by treating the errors as homoskedastic; the larger conditional variance at \(X_i=2\) contributes \(4\cdot4\) to \(B\).

Worked Example: Homoskedastic Simplification

Keep the same scalar design, \(P(X_i=1)=P(X_i=2)=1/2\), but suppose \(E[u_i\mid X_i]=0\) and \(E[u_i^2\mid X_i]=4\) for both design values. Then \(A=5/2\), and \(B=E[4X_i^2]=4(5/2)=10\). The limiting variance is \[ V=A^{-2}B=(2/5)^2(10)=8/5=1.6. \]

Equivalently, under homoskedasticity the formula is \(\sigma^2 A^{-1}=4/(5/2)=8/5\). At finite \(n\), the conditional covariance of the scaled error is \(4/A_n\). As the observed design accumulates, \(A_n\xrightarrow{p}5/2\), so this conditional variance converges to \(8/5\), matching the unconditional limit.

When agreement can fail

The stable iid design assumptions matter. If the design is governed by a random quantity shared across the whole sample, conditioning can reveal a persistent feature that does not average away. The conditional limiting covariance may then remain random, and the unconditional limit may be a mixture rather than a single normal distribution.

Worked Example: A Shared Random Design Scale

Consider regression through the origin with \(X_i=Z\) for every \(i\), where \(P(Z=1)=P(Z=2)=1/2\). Let \(u_i\) be iid \(N(0,1)\), independent of \(Z\), and let \(Y_i=X_i\beta+u_i\). Conditional on \(Z\), \[ \sqrt n(\hat\beta-\beta) =\frac{1}{Z\sqrt n}\sum_{i=1}^n u_i \sim N(0,1/Z^2) \] for every \(n\). Thus the conditional limit has variance \(1\) when \(Z=1\) and \(1/4\) when \(Z=2\); it does not converge to one fixed conditional covariance.

Unconditionally, the distribution is the mixture \(\tfrac12N(0,1)+\tfrac12N(0,1/4)\), not a single normal distribution. Its variance is \(\tfrac12(1)+\tfrac12(1/4)=5/8\). The unconditional variance is the average of the two conditional variances, but it does not describe either conditional limit. This example is not an iid random-design sample: all regressors share the same \(Z\). It shows why a stable design limit is needed for the agreement theorem.

Common mistake

  • Wrong reasoning: “Once the design is observed, OLS has covariance \(\sigma^2(X'X)^{-1}\), so that must be the right covariance in every model.” This uses the homoskedastic formula even when the conditional error variance changes with the regressors.
  • Correction: Under conditional heteroskedasticity, use the conditional score covariance. For the scaled estimator it is \(A_n^{-1}B_nA_n^{-1}\), which converges to \(A^{-1}BA^{-1}\) under the stated conditions. Also distinguish this probability limit from the expectation of the finite-sample conditional covariance; the latter requires additional control of integrability.
Key takeaway: In a stable iid random-design regression, conditional and unconditional OLS limits have the same covariance \(A^{-1}BA^{-1}\) when the conditional score variances converge and a conditional central limit condition holds. If the design retains a random feature at the sample limit, the conditional covariance can remain random and the unconditional limit can be a mixture.

Check Your Understanding

Use the definitions and conditions in this tutorial to answer the following questions.

  1. What is conditioned on in the conditional asymptotic distribution of OLS?
  2. For scalar regression through the origin, write \(A\), \(B\), and \(V\) in terms of \(X_i\) and \(u_i\).
  3. Why does \(E[u_i\mid X_i]=0\) imply that the score \(X_i u_i\) has mean zero?
  4. Under homoskedasticity with conditional variance \(\sigma^2\), simplify \(A^{-1}BA^{-1}\).
  5. In the shared-scale example, why is the unconditional limiting law not a single normal distribution?
Clinical Biostats