Tutorials › Statistical Theory III - Theories of Inference, Multivariate Methods, and Contingency Tables › Leverage-Based Lindeberg Conditions

Random-Design Regression Asymptotics · Tutorial 579 of 972

Leverage-Based Lindeberg Conditions

Learn to express score contributions in leverage-scaled form and use leverage and conditional error-tail bounds to verify—or diagnose failure of—the conditional Lindeberg condition.

Advanced 10 min read

What You'll Learn

  • Define the conditional Lindeberg condition for regression score contributions.
  • Construct variance-weighted leverage by whitening the score with its conditional covariance.
  • Relate variance-weighted leverage to ordinary hat-matrix leverage under homoskedasticity.
  • Verify the condition using bounded errors or uniformly controlled standardized error tails.
  • Diagnose failure when one observation retains a nonvanishing share of score variance.
  • Distinguish a leverage bound from the additional error-tail control it requires.

Why Leverage Enters a Score Limit

In random-design regression, the score that drives the least-squares coefficient error is a sum of terms \(X_i\varepsilon_i\). A central limit approximation for that sum requires more than a stable average covariance: no small set of unusually influential terms should retain a substantial share of the total variation. The conditional Lindeberg condition makes this requirement precise after conditioning on the observed regressors.

The key quantity is a variance-weighted version of leverage. It measures how much each observation contributes to the covariance of the score, rather than only how unusual its regressor is. A useful sufficient condition combines vanishing maximum weighted leverage with uniform control of the tails of standardized errors.

Conditional Score Contributions and Leverage

For each sample size \(n\), let \(X_{ni}\in\mathbb{R}^p\) be the regressor vector and let \(\varepsilon_{ni}\) be the regression error, for \(i=1,\ldots,n\), with fixed \(p\). Condition on the complete design \(\mathcal X_n=(X_{n1},\ldots,X_{nn})\). Assume that, given \(\mathcal X_n\), the errors are independent, have conditional means zero, and have finite, strictly positive conditional variances \(\sigma_{ni}^2=E[\varepsilon_{ni}^2\mid\mathcal X_n]\). Define the conditional score covariance per observation by

$$ Q_n=\frac{1}{n}\sum_{i=1}^n \sigma_{ni}^2 X_{ni}X_{ni}' , $$

and assume \(Q_n\) is positive definite. The normalized score and its covariance-standardized individual contributions are

$$ S_n=\frac{1}{\sqrt n}\sum_{i=1}^n X_{ni}\varepsilon_{ni}, \qquad V_{ni}=\frac{1}{\sqrt n}Q_n^{-1/2}X_{ni}\varepsilon_{ni}. $$

Here \(Q_n^{-1/2}\) is the symmetric inverse square root. Conditional on the design, the \(V_{ni}\) are independent, have mean zero, and their covariance matrices sum to the identity. Their squared lengths in expectation define the variance-weighted leverages.

Definition: The variance-weighted leverage of observation \(i\) is \(\lambda_{ni}=\sigma_{ni}^2 X_{ni}'Q_n^{-1}X_{ni}/n\). It is the conditional expectation of \(\|V_{ni}\|^2\). Since \(Q_n\) is positive definite, \(\sum_{i=1}^n\lambda_{ni}=p\).

In the homoskedastic case, \(\sigma_{ni}^2=\sigma^2\), this definition reduces to the ordinary diagonal leverage \(h_{ni}=X_{ni}'(\sum_jX_{nj}X_{nj}')^{-1}X_{ni}\). Thus weighted leverage extends familiar design leverage to account for differences in conditional error variance.

The Conditional Lindeberg Condition

Definition: The conditional vector Lindeberg condition holds if, for every \(\delta>0\), \[ \sum_{i=1}^n E\!\left[\|V_{ni}\|^2 \mathbf 1\{\|V_{ni}\|>\delta\}\mid\mathcal X_n\right] \xrightarrow{p}0 . \]

The condition says that the total conditional variance contributed by score terms larger than any fixed threshold tends to zero. Convergence in probability is appropriate when the design, and therefore the conditional quantities, are random. In this setup, write \(\eta_{ni}=\varepsilon_{ni}/\sigma_{ni}\), so \(E[\eta_{ni}^2\mid\mathcal X_n]=1\). Directly from the definitions,

$$ \|V_{ni}\|^2=\lambda_{ni}\eta_{ni}^2, \qquad \|V_{ni}\|>\delta \ \Longleftrightarrow\ \sqrt{\lambda_{ni}}|\eta_{ni}|>\delta . $$

Consequently, the Lindeberg sum is exactly

$$ \sum_{i=1}^n \lambda_{ni} E\!\left[ \eta_{ni}^2 \mathbf 1\{\sqrt{\lambda_{ni}}|\eta_{ni}|>\delta\} \mid\mathcal X_n \right]. $$

This expression shows why small leverage alone is not the whole argument: an observation with small weighted leverage can still contribute to the Lindeberg sum if its standardized error has sufficiently heavy tails.

A Leverage and Tail Bound

Theorem: Suppose \(\sum_i\lambda_{ni}=p\), \(\max_i\lambda_{ni}\xrightarrow{p}0\), and there is a deterministic function \(\rho(t)\) with \(\rho(t)\to0\) as \(t\to\infty\) such that, almost surely, \[ \sup_{n,i} E[\eta_{ni}^2\mathbf 1\{|\eta_{ni}|>t\}\mid\mathcal X_n] \leq \rho(t) \quad\text{for every }t>0. \] Then the conditional vector Lindeberg condition holds.

Proof. Fix \(\delta>0\), and let \(m_n=\max_i\lambda_{ni}\). On the event \(m_n\leq a\), where \(a>0\), the indicator \(\{\sqrt{\lambda_{ni}}|\eta_{ni}|>\delta\}\) implies \(|\eta_{ni}|>\delta/\sqrt a\). Therefore,

$$ \begin{aligned} &\sum_i\lambda_{ni} E\!\left[\eta_{ni}^2 \mathbf 1\{\sqrt{\lambda_{ni}}|\eta_{ni}|>\delta\} \mid\mathcal X_n\right]\\ &\qquad\leq \sum_i\lambda_{ni}\rho(\delta/\sqrt a) =p\,\rho(\delta/\sqrt a). \end{aligned} $$

Given any positive tolerance, choose \(a\) sufficiently small that the last bound is below that tolerance; this is possible because \(\rho(t)\to0\). Since \(m_n\xrightarrow{p}0\), the probability that \(m_n>a\) tends to zero. Thus the Lindeberg sum tends to zero in probability. This proves the result. \(\square\)

A standard conditional Lindeberg–Feller theorem now gives a conditional normal limit for the whitened score: under the stated conditional independence and mean-zero assumptions, the conditional distribution of \(\sum_i V_{ni}=Q_n^{-1/2}S_n\) approaches \(N(0,I_p)\) in probability. If, in addition, \(Q_n\xrightarrow{p}Q\) for a deterministic positive-definite matrix \(Q\), the corresponding unwhitened score has limiting distribution \(N(0,Q)\). The theorem above supplies a practical way to check the Lindeberg hypothesis; it does not replace the covariance convergence needed for the latter limit.

Worked Examples

Worked Example: Bounded Errors and Equal Leverage

Take a scalar regression with \(X_{ni}=1\) for every observation and independent errors \(\varepsilon_{ni}\) taking values \(1\) and \(-1\) with equal probability. Then \(\sigma_{ni}^2=1\), \(Q_n=1\), and \(\lambda_{ni}=1/n\). Also \(\eta_{ni}=\varepsilon_{ni}\), so \(|\eta_{ni}|=1\).

For any fixed \(\delta>0\), the indicator in the Lindeberg sum is \(\mathbf 1\{1/\sqrt n>\delta\}\), which equals zero as soon as \(n\geq\delta^{-2}\). The entire sum is therefore eventually zero. Here the maximum leverage tends to zero, and bounded standardized errors make the conditional Lindeberg condition immediate.

Worked Example: Gaussian Errors with Small Maximum Leverage

Suppose the standardized errors are conditionally \(N(0,1)\). Let \(Z\sim N(0,1)\). For \(t>0\), integration by parts gives

$$ E[Z^2\mathbf 1\{|Z|>t\}] =2\int_t^\infty z^2\phi(z)\,dz =2\{t\phi(t)+1-\Phi(t)\}, $$

where \(\phi\) and \(\Phi\) are the standard normal density and distribution function. This quantity tends to zero as \(t\to\infty\). For each \(\delta>0\), the Lindeberg sum is bounded above by

$$ \sum_i\lambda_{ni} E[Z^2\mathbf 1\{|Z|>\delta/\sqrt{m_n}\}] = p\,E[Z^2\mathbf 1\{|Z|>\delta/\sqrt{m_n}\}], \qquad m_n=\max_i\lambda_{ni}. $$

If \(m_n\xrightarrow{p}0\), the threshold \(\delta/\sqrt{m_n}\) tends to infinity in probability, so this bound tends to zero in probability. Thus Gaussian errors satisfy the required tail control, and vanishing maximum weighted leverage verifies the condition.

Worked Example: One Observation Retains Half the Score Variance

Consider a scalar design with \(X_{n1}=\sqrt n\) and \(X_{ni}=1\) for \(i=2,\ldots,n\). Let the errors be independent Rademacher variables, so each is \(1\) or \(-1\) with equal probability and has variance one. Then

$$ Q_n=\frac{1}{n}\sum_iX_{ni}^2 =\frac{2n-1}{n}, \qquad \lambda_{n1}=\frac{n}{2n-1}\longrightarrow\frac12, \qquad \lambda_{ni}=\frac{1}{2n-1}\quad(i\geq2). $$

Choose \(\delta=0.6\). For all sufficiently large \(n\), \(\sqrt{\lambda_{n1}}>0.6\), and \(|\eta_{n1}|=1\). The first observation alone contributes \(\lambda_{n1}\) to the Lindeberg sum, which tends to \(1/2\), not zero. Hence the condition fails. This calculation identifies the precise obstruction: the first score contribution retains a nonvanishing share of total variance.

Using the Condition in Regression Limits

In least squares with a correctly specified conditional mean, the coefficient error has a score representation involving \(S_n=n^{-1/2}\sum_iX_{ni}\varepsilon_{ni}\) and the empirical Gram matrix. A conditional normal approximation for the score is one ingredient in an asymptotic normality argument; a stable, nonsingular Gram matrix is another. Leverage-based Lindeberg checks address whether individual score terms are negligible, not whether the design is identified or the covariance matrix converges.

A useful proof strategy is to compute the covariance-weighted leverage, establish that its maximum vanishes, and then verify a uniform tail bound for the standardized errors. Under homoskedasticity, ordinary maximum hat-matrix leverage is the relevant leverage measure. Under heteroskedasticity, the variance weights matter: raw design leverage can be small while one observation still dominates the score covariance.

Common Mistake

A tempting but incomplete argument is: “Every ordinary hat-matrix leverage tends to zero, so the score satisfies Lindeberg.” This overlooks two separate requirements. First, with heteroskedastic errors, the relevant quantity is variance-weighted leverage. Second, even small weighted leverage needs suitable control of the standardized error tails.

For example, take \(X_{ni}=1\) for all \(i\), let \(\sigma_{n1}^2=n\), and let \(\sigma_{ni}^2=1\) for \(i\geq2\). Ordinary leverage is \(1/n\) for every observation, but the weighted leverage of the first observation is

$$ \lambda_{n1}=\frac{n}{n+(n-1)} =\frac{n}{2n-1}\longrightarrow\frac12. $$

If the standardized errors are Rademacher, the same calculation as in the preceding example shows that Lindeberg fails. The correction is to assess leverage using the conditional score covariance and then check the standardized error tails; neither check can be silently omitted.

Key takeaway: The conditional Lindeberg condition rules out large score contributions in total variance. Vanishing maximum variance-weighted leverage, together with uniformly controlled tails of standardized errors, is a transparent sufficient condition.

Check Your Understanding

Use the conditional score and leverage definitions to answer these questions.

  1. Why does \(\sum_i\lambda_{ni}=p\) hold, and what does this identity represent?
  2. In the homoskedastic case, how does variance-weighted leverage relate to the ordinary hat-matrix diagonal?
  3. Why does vanishing maximum weighted leverage not, by itself, establish the conditional Lindeberg condition?
  4. For Gaussian standardized errors, what happens to the second-moment tail \(E[Z^2\mathbf 1\{|Z|>t\}]\) as \(t\) increases?
  5. In the high-leverage example, why does the first observation prevent the Lindeberg sum from tending to zero?
Clinical Biostats