Tutorials › Statistical Learning (Graduate) › The Leave-One-Out Generalization Identity

Localized Complexity and Stability · Tutorial 186 of 989

The Leave-One-Out Generalization Identity

Learn why an expected generalization gap can be written as a difference of out-of-sample and training losses, and how replacement-based and delete-one identities differ.

Advanced 10 min read

What You'll Learn

  • Define expected generalization gap and the assumptions needed for its expectation to exist
  • Derive the fixed-sample-size identity using an independent replacement observation
  • Prove the delete-one identity for an algorithm trained on one fewer observation
  • Distinguish replacement-based evaluation from conventional leave-one-out cross-validation
  • Check both identities with explicit Bernoulli and squared-error calculations
  • Avoid an invalid conditional exchange involving a fixed replacement observation

Why an Expected Generalization Gap Has a Leave-One-Out Representation

A model is evaluated on the same observations used to fit it when we calculate empirical training loss. Its population risk instead averages loss on a new observation, independent of the training data. The difference between these quantities is the generalization gap. Independence allows its expectation to be represented using held-out observations, but the precise identity depends on whether a held-out point is deleted or replaced.

This distinction matters. Replacing one training observation keeps the training sample size fixed; deleting one observation reduces it by one. The corresponding identities are both exact, but they concern algorithms trained on different sample sizes. Neither identity, by itself, says that a model's training loss estimates its risk accurately for a particular realized data set.

Setup and Definitions

Let \(Z_1,\ldots,Z_n\) be independent observations with common distribution \(P\), and write \(S=(Z_1,\ldots,Z_n)\). A deterministic learning algorithm \(A_n\) maps a sample of size \(n\) to a predictor or model. Let \(\ell(a,z)\) be the loss of model \(a\) on observation \(z\). For a fixed model \(a\), define its population risk as \(R(a)=E_{Z\sim P}[\ell(a,Z)]\), when this expectation exists. Its empirical training risk on \(S\) is \(\hat R_S(a)=n^{-1}\sum_{i=1}^n\ell(a,Z_i)\).

Definition: The expected generalization gap of \(A_n\) is \(G_n=E[R(A_n(S))-\hat R_S(A_n(S))]\). The sign convention here is population risk minus training risk. A positive gap means that, on average, the fitted model has larger loss on a fresh observation than on its training observations.

For the identities below, assume the data space and model space are measurable, the algorithm and loss are measurable, and the displayed losses have finite absolute expectations. In particular, these conditions permit writing population risk as an expectation over an independent observation and taking the expectations and finite sums below. We use a deterministic algorithm to keep the notation focused; an algorithm with independent randomization can be handled by including that randomization in the expectations.

The Fixed-Size Replacement Identity

Let \(Z_1',\ldots,Z_n'\) be independent observations from \(P\), also independent of \(S\). For each \(i\), define the replacement sample \(S^{(i)}=(Z_1,\ldots,Z_{i-1},Z_i',Z_{i+1},\ldots,Z_n)\). The point \(Z_i\) is held out for evaluation, while \(Z_i'\) takes its place in the training sample. Thus \(A_n(S^{(i)})\) is still trained on \(n\) observations.

Theorem: Under the stated measurability and integrability conditions,
$$ G_n = E\left[\frac{1}{n}\sum_{i=1}^n \left\{ \ell(A_n(S^{(i)}),Z_i) - \ell(A_n(S),Z_i) \right\}\right]. $$
The first loss in each summand is the held-out loss after replacement; the second is the training loss on the original sample.

Proof. Take an independent \(Z'\sim P\). By the definition of risk and independence of \(Z'\) and \(S\), \(E[R(A_n(S))]=E[\ell(A_n(S),Z')]\). For a fixed index \(i\), hold all coordinates of \(S\) other than \(Z_i\) fixed, but do not condition on the replacement observation \(Z_i'\). The pair \((Z_i,Z')\) has distribution \(P\otimes P\), as does \((Z_i',Z_i)\). The first coordinate of the first pair supplies the \(i\)th training entry and the second supplies the evaluation point. In the second pair, those roles are supplied by \(Z_i'\) and \(Z_i\). Consequently,

$$ E[\ell(A_n(S),Z')] = E[\ell(A_n(S^{(i)}),Z_i)]. $$

This equality is unconditional over the two independent draws and the other sample coordinates. It holds for every \(i\). Averaging these equalities over \(i\), and subtracting \(E[\hat R_S(A_n(S))]=n^{-1}\sum_i E[\ell(A_n(S),Z_i)]\), gives the identity. \(\square\)

The proof uses exchangeability of the relevant independent pairs, not a pointwise equality of losses. In general, the two losses inside a summand are dependent, and neither one must have the same value as the other on any particular data set.

Worked Example: A Sample Mean Under Squared Error

Let \(Z_i\sim\operatorname{Bernoulli}(1/2)\), let the algorithm return the sample mean \(A_n(S)=\bar Z\), and use squared-error loss \(\ell(a,z)=(a-z)^2\). Take \(n=2\). We calculate both sides of the fixed-size identity.

Population risk. For a fixed \(a\), \(R(a)=\tfrac12a^2+\tfrac12(1-a)^2=(a-\tfrac12)^2+\tfrac14\). Since \(\operatorname{Var}(\bar Z)=1/8\) and \(E[\bar Z]=1/2\), \(E[R(\bar Z)]=1/8+1/4=3/8\).

Training risk. If \(Z_1=Z_2\), the sample mean equals both observations and the training risk is zero. If they differ, the sample mean is \(1/2\), and each squared residual is \(1/4\), so the average training risk is \(1/4\). The observations differ with probability \(1/2\); therefore \(E[\hat R_S(\bar Z)]=(1/2)(1/4)=1/8\). The expected gap is \(3/8-1/8=1/4\).

Replacement side. For either index \(i\), the replacement-trained mean is based on two independent Bernoulli observations, \(Z_i'\) and the other original observation. The held-out \(Z_i\) is independent of that mean. Its expected squared loss is therefore its population risk averaged over a sample mean, namely \(3/8\). The original training loss at index \(i\) has expectation \(1/8\), since the two original observations are exchangeable and their average loss has expectation \(1/8\). Each summand has expected difference \(3/8-1/8=1/4\), so their average is \(1/4\), agreeing with the expected gap.

The Conventional Delete-One Identity

In conventional leave-one-out cross-validation, observation \(Z_i\) is removed, rather than replaced. Write \(S_{-i}\) for the ordered sample of the other \(n-1\) observations, and let \(A_{n-1}\) be an algorithm trained on \(n-1\) observations. Define the leave-one-out prediction loss and the corresponding training loss by

$$ \hat R_{\mathrm{LOO}} = \frac{1}{n}\sum_{i=1}^n\ell(A_{n-1}(S_{-i}),Z_i), \qquad \hat R_{\mathrm{train,LOO}} = \frac{1}{n}\sum_{i=1}^n \frac{1}{n-1}\sum_{j\ne i} \ell(A_{n-1}(S_{-i}),Z_j). $$

The second quantity averages the training risk of each of the \(n\) models, each fitted after deleting one observation. It is not the empirical risk of one model trained on all \(n\) observations.

Theorem: For \(n\ge 2\), under the same measurability and integrability conditions,
$$ E[\hat R_{\mathrm{LOO}}-\hat R_{\mathrm{train,LOO}}] = E\left[ R(A_{n-1}(T))-\hat R_T(A_{n-1}(T)) \right], $$
where \(T\) is a sample of \(n-1\) independent observations from \(P\).

Proof. For each \(i\), \(Z_i\) is independent of \(S_{-i}\) and has distribution \(P\). Conditioning on \(S_{-i}\), then integrating over \(Z_i\), gives \(E[\ell(A_{n-1}(S_{-i}),Z_i)]=E[R(A_{n-1}(S_{-i}))]\). The sample \(S_{-i}\) has the distribution of \(T\), so each term has expectation \(E[R(A_{n-1}(T))]\).

For the training part, \(S_{-i}\) also has the distribution of \(T\), and its inner average is precisely the empirical risk of \(A_{n-1}(S_{-i})\) on that sample. Hence \(E[(n-1)^{-1}\sum_{j\ne i}\ell(A_{n-1}(S_{-i}),Z_j)]=E[\hat R_T(A_{n-1}(T))]\). Averaging the two equalities over \(i\) and subtracting proves the result. \(\square\)

This identity concerns the expected generalization gap at training size \(n-1\), not \(n\). The fixed-size replacement identity instead concerns \(A_n\), with a replacement observation still included in its training set. Comparing the two expected gaps requires additional information about how the algorithm behaves as its training size changes.

Worked Example: Delete-One Cross-Validation for a Bernoulli Mean

Let \(Z_i\sim\operatorname{Bernoulli}(1/4)\), use \(A_m\) equal to the sample mean of \(m\) observations, and use squared-error loss. Take \(n=3\), so each leave-one-out model is trained on \(m=2\) observations. Write \(p=1/4\) and \(q=3/4\), so \(pq=3/16\).

Risk at training size two. For a sample mean \(\bar Z\) based on \(m\) independent Bernoulli observations, \(E[\bar Z]=p\) and \(\operatorname{Var}(\bar Z)=pq/m\). A fresh \(Z\) is independent of \(\bar Z\), so \(E[(\bar Z-Z)^2]=\operatorname{Var}(\bar Z)+\operatorname{Var}(Z)=pq/m+pq\). At \(m=2\), this equals \((3/16)(3/2)=9/32\).

Training risk at size two. For any two values, their average squared residual from their sample mean is \((Z_1-Z_2)^2/4\). The squared difference is one exactly when the observations differ, an event of probability \(2pq=3/8\). Thus the expected training risk is \((3/8)(1/4)=3/32\). The expected gap is \(9/32-3/32=6/32=3/16\).

For each deleted index, the held-out observation is independent of the two-point training sample, so the expected leave-one-out prediction loss is \(9/32\). The expected training loss of that two-point sample is \(3/32\). Their difference is \(3/16\), exactly the expected gap at training size two.

A Common Pitfall: Conditioning on the Replacement

A tempting but invalid proof conditions on the replacement value \(Z_i'\), treats it as fixed, and then tries to exchange \(Z_i\) with an independent test point. Once \(Z_i'\) is fixed, the replacement sample still contains that fixed value; exchanging \(Z_i\) with the test point does not turn the replacement sample into the original sample in the required way. The valid proof compares the joint distributions of \((Z_i,Z')\) and \((Z_i',Z_i)\) without fixing \(Z_i'\).

Worked Example: Why the Exchange Must Be Unconditional

Take \(n=1\), \(Z,Z'\) independent \(\operatorname{Bernoulli}(1/2)\), algorithm \(A_1(x)=x\), and bounded loss \(\ell(a,z)=a\). The risk of the fitted model is \(R(A_1(Z))=Z\), and its training loss is also \(Z\), so the generalization gap is exactly zero.

The replacement loss is \(\ell(A_1(Z'),Z)=Z'\), with expectation \(1/2\); the original training loss has expectation \(1/2\). Thus the replacement identity gives \(1/2-1/2=0\), as required. But if we condition on \(Z'=0\), the replacement loss becomes zero while the independent risk term \(E[\ell(A_1(Z),Z')]\) remains \(E[Z]=1/2\). The conditional expectations are not equal. Unconditionally, \(Z'\) is again a Bernoulli draw, so both relevant pairs have the same independent product distribution and the exchange is valid.

How to Use the Identities

These identities express an expected gap as a difference between two evaluation schemes. They do not require the algorithm to be an empirical risk minimizer, nor do they assume stability. Their role is to make the comparison exact; further assumptions are needed to bound the size of the difference. In particular, if changing one observation has little effect on the fitted model's loss, the replacement expression can help control the expected gap.

Keep the sample size and evaluation design explicit when applying the result:

  • Fixed-size replacement: train on \(n\) observations, replacing \(Z_i\) by an independent draw, and evaluate on the omitted \(Z_i\). This gives the expected gap for \(A_n\).
  • Delete-one cross-validation: train on the \(n-1\) observations excluding \(Z_i\), then evaluate on \(Z_i\). Comparing its prediction loss with the training loss of those deleted samples gives the expected gap for \(A_{n-1}\).

A frequent mistake is to call these two procedures interchangeable because each involves holding out one observation. They use different training sets and, generally, different fitted models. Another is to interpret an identity in expectation as a statement about each realized sample. The equalities proved here concern expectations under independent sampling; they do not assert equality of the corresponding random quantities for every data set.

Check Your Understanding

Use the definitions and keep track of the training sample size in each answer.

  1. In the fixed-size replacement identity, how many observations are used to train the replacement model, and which observation is used to evaluate it?
  2. Why is \(Z_i\) independent of \(S_{-i}\) in the delete-one identity, and what role does that independence play?
  3. For the delete-one identity with \(n\) original observations, which algorithm's expected generalization gap appears on the right-hand side?
  4. Why does conditioning on a fixed \(Z_i'\) invalidate the pair-exchange step in the replacement proof?
  5. Does either identity imply that the realized training loss and test loss are close for every sample? Explain.
Clinical Biostats