← Back to Session 3.1  ·  Course contents
Note: This page's design and content were created and enhanced using Claude (Anthropic's AI assistant); the references have been checked. Spot an error? Email jonathan.shock@uct.ac.za.
Session 3.1 • Derivation appendix

The optimizer's curse

A formal derivation of regressional Goodhart, and how the bias grows with optimisation

What this appendix does

Session 3.1 claimed that optimising an unbiased proxy systematically over-states the true value of whatever you select, and that the over-statement grows the harder you optimise. This page proves it. We derive the post-selection bias exactly for the Gaussian case, show how it scales when you take the maximum over many options (the role of "optimisation power"), give the distribution-free version, and map the result back to specification gaming.

Prerequisites: joint Gaussians and conditional expectation, a little order-statistics / extreme-value theory, and comfort with expectations. Nothing beyond a strong second-year probability course. This is optional deep-dive material; the main thread continues in Session 3.1.

① Setup and the claim

The model

We have a population of options indexed by \(i\). Option \(i\) has an unknown true value \(G_i\) (the goal) and we observe only a noisy proxy

\[ P_i = G_i + \varepsilon_i, \qquad \mathbb{E}[\varepsilon_i \mid G_i] = 0, \qquad \varepsilon_i \perp G_i. \]

The proxy is unbiased as a measurement: averaged over the noise, it equals the truth. We optimise by selecting the option with the largest proxy, \(i^* = \operatorname{argmax}_i P_i\).

The claim, precisely

Despite the proxy being unbiased before selection, the selected option's proxy is biased upward as an estimate of its true value: \(\mathbb{E}[P_{i^*} - G_{i^*}] \ge 0\). The distortion is created by the \(\operatorname{argmax}\), not by any error in the measurement, and it grows with how hard you optimise. This post-selection over-statement is the optimizer's curse (Smith & Winkler, 2006); seen through the lens of metric-gaming it is regressional Goodhart (Manheim & Garrabrant, 2018).

② The single-option posterior

First, a sub-question: given that one option scored a proxy value \(p\), what do we expect its true value to be? Take \(G \sim \mathcal{N}(\mu_G, \sigma_G^2)\) and \(\varepsilon \sim \mathcal{N}(0, \sigma_\varepsilon^2)\), independent.

Derivation by precision-weighting (completing the square)

Because \(\varepsilon = P - G\), the conditional density of \(G\) given \(P = p\) is proportional to the product of the prior on \(G\) and the likelihood of the observation:

\[ f(g \mid p) \propto \exp\!\left( -\frac{(g - \mu_G)^2}{2\sigma_G^2} \right) \cdot \exp\!\left( -\frac{(p - g)^2}{2\sigma_\varepsilon^2} \right). \]

The exponent is quadratic in \(g\), so the posterior is Gaussian. Completing the square (equivalently: combining two Gaussian "beliefs" by adding their precisions \(1/\sigma^2\) and taking the precision-weighted mean) gives

\[ \mathbb{E}[G \mid P = p] = \mu_G + \beta (p - \mu_G), \qquad \operatorname{Var}[G \mid P = p] = \sigma_G^2(1 - \beta), \]

where the shrinkage factor is the proxy's reliability, the fraction of its variance that is signal:

\[ \beta = \frac{\sigma_G^2}{\sigma_G^2 + \sigma_\varepsilon^2} \in (0, 1). \]

Two readings of β

This \(\beta\) is exactly the slope of the least-squares regression of \(G\) on \(P\): for jointly Gaussian variables the conditional expectation is the linear regression, with slope \(\operatorname{Cov}(G,P)/\operatorname{Var}(P) = \sigma_G^2/(\sigma_G^2+\sigma_\varepsilon^2)\). Hence "regression to the mean" (Galton, 1886): the estimate of a high scorer is pulled back toward the population mean by the factor \(\beta\). The estimate is also the minimum-mean-squared-error estimator, so no cleverer use of the same data avoids the shrinkage.

③ The over-statement, exactly

Bias of the proxy at a given score

Define the over-statement \(b(p)\) as proxy minus expected truth at proxy value \(p\):

\[ b(p) = p - \mathbb{E}[G \mid P = p] = (1 - \beta)(p - \mu_G) = \frac{\sigma_\varepsilon^2}{\sigma_G^2 + \sigma_\varepsilon^2}(p - \mu_G). \]

Equivalently, the expected noise of an option conditioned on its score is positive for above-average scores: \(\mathbb{E}[\varepsilon \mid P = p] = b(p) > 0\) for \(p > \mu_G\). Selecting a high proxy preferentially selects options that "got lucky".

Two consequences of the formula

(i) The harder you optimise, the worse it is. \(b(p)\) is proportional to \((p - \mu_G)\), how far into the tail your winner sits. Selecting an extreme proxy value inflates the over-statement without bound.

(ii) The worse the proxy, the worse it is. \(b(p)\) scales with \((1 - \beta) = \sigma_\varepsilon^2/(\sigma_G^2+\sigma_\varepsilon^2)\), the noise fraction. A perfect proxy (\(\sigma_\varepsilon = 0\), \(\beta = 1\)) gives zero bias; a noisy one gives a large one.

Equal-variance example

If signal and noise have equal variance, \(\sigma_G^2 = \sigma_\varepsilon^2\), then \(\beta = \tfrac{1}{2}\). An option whose proxy sits two standard deviations above the mean has an expected true value only one standard deviation above the mean: half of its apparent excellence is luck that will not reproduce. Optimise on a metric that is half noise and you systematically harvest options half as good as they look.

④ Selecting the best of many: optimisation power

"Optimising harder" often means searching more options. Taking the \(\operatorname{argmax}\) over \(N\) independent draws pushes the winning proxy deep into the tail, and §③ then says the over-statement grows with it.

Extreme-value scaling

Let \(P_1,\ldots,P_N\) be i.i.d. \(\mathcal{N}(\mu_P, \sigma_P^2)\) with \(\sigma_P^2 = \sigma_G^2 + \sigma_\varepsilon^2\), and \(M_N = \max_i P_i\). The normal distribution lies in the Gumbel domain of attraction, with centering constant (Embrechts, Klüppelberg & Mikosch, 1997)

\[ b_N = \sqrt{2 \ln N} - \dfrac{\ln \ln N + \ln 4\pi}{2\sqrt{2 \ln N}}, \qquad a_N = \frac{1}{\sqrt{2 \ln N}}. \]

Then \((M_N - b_N)/a_N\) converges in distribution to a standard Gumbel, whose mean is the Euler–Mascheroni constant \(\gamma \approx 0.5772\). Hence the expected maximum is

\[ \mathbb{E}[M_N] = \mu_P + \sigma_P\!\left( b_N + \gamma \cdot a_N \right) + o\!\left(1/\sqrt{\ln N}\right) \approx \mu_P + \sigma_P\sqrt{2 \ln N} \]

(leading order). (The leading term is \(\sqrt{2 \ln N}\); \(b_N\) is the precise centering constant and the \(\gamma \cdot a_N\) term is the same order as the correction, so we keep both rather than mixing them.)

The winner's shortfall

The selected option's expected over-statement is, to leading order, \(b\) evaluated at the winning score:

\[ \mathbb{E}[P_{i^*} - G_{i^*}] \approx (1 - \beta)\!\left( \mathbb{E}[M_N] - \mu_P \right) \approx (1 - \beta) \cdot \sigma_P \cdot \sqrt{2 \ln N}. \]

So the gap between what the winner scored and what it is actually worth grows like \(\sqrt{\ln N}\) in the number of options searched: slowly, but without bound. More search, more capability, more compute → a winner further into the tail → a larger illusory inflation.

A worked number

Search \(N = 1000\) options, equal signal/noise (\(\beta = \tfrac{1}{2}\)), \(\sigma_P = 1\): \(\sqrt{2 \ln 1000} \approx 3.7\), so the winner's proxy sits ≈3.7 σ above the mean, of which the expected over-statement is ≈ \((1-\tfrac{1}{2}) \cdot 3.7 \approx 1.9\) standard deviations of pure illusion. Push to \(N = 10^6\) and \(\sqrt{2 \ln N} \approx 5.3\): the curse worsens, but only logarithmically.

⑤ The distribution-free version

The Gaussian gave us exact numbers, but the bias does not depend on normality. Smith & Winkler's (2006) general statement needs only unbiased estimates and a one-line convexity argument.

The optimizer's curse, in general

Suppose options \(1,\ldots,N\) have fixed true values \(V_i\) and we hold unbiased estimates \(\hat{V}_i\) (so \(\mathbb{E}[\hat{V}_i] = V_i\)). Select \(i^* = \operatorname{argmax}_i \hat{V}_i\). Because the maximum is a convex function, Jensen's inequality gives

\[ \mathbb{E}[\hat{V}_{i^*}] = \mathbb{E}[\max_i \hat{V}_i] \ge \max_i \mathbb{E}[\hat{V}_i] = \max_i V_i \ge V_{i^*}. \]

Taking expectations, \(\mathbb{E}[\hat{V}_{i^*} - V_{i^*}] \ge 0\): the estimate of the selected option is biased high, regardless of the distribution of the estimation errors. Normality only fixes the size of the gap (§②–④); its sign is forced by selection alone.

⑥ Scope of the result

Regressional Goodhart

The shrinkage by \(\beta\) and the noise-driven over-statement are Manheim & Garrabrant's (2018) regressional variant: the unavoidable consequence of selecting on (signal + noise) with a fixed joint distribution. It needs no adversary and no breakdown of the proxy.

The gentle case

The derivation assumes the \(G\)–\(P\) relationship stays linear-Gaussian out into the tail. Real optimisers also reach regions where the proxy stops tracking the goal (extremal Goodhart), where the correlation was never causal (causal), or where another agent games the metric (adversarial). Those are worse, and not captured here.

The closed form as a lower bound

\((1-\beta) \cdot \sigma_P \cdot \sqrt{2 \ln N}\) is the over-statement you pay even in the friendliest possible world. Reality adds the extremal/causal/adversarial terms on top. That is why §⑥ matters: regressional Goodhart proves the problem is real in principle; the other variants are why it is often severe in practice.

⑦ Back to specification gaming

The dictionary

  • Options ↔ the policies / behaviours a model can reach.
  • Proxy \(P\) ↔ the reward or training metric we actually optimise.
  • True value \(G\) ↔ the intended objective we cannot write down.
  • argmax over N ↔ training: searching the policy space for high-reward behaviour. A more capable optimiser explores a larger effective \(N\).

Read in these terms, the result gives a quantitative basis for Session 3.1's empirical claims. The winner's-shortfall formula says the gap between measured reward and true quality grows with optimisation power, the mechanism behind Pan, Bhatia & Steinhardt's (2022) finding that more capable agents game misspecified rewards more. And it sharpens Skalse et al.'s (2022) impossibility result: a non-trivial proxy can be optimised in a value-lowering direction, and the expected over-statement of the best-by-proxy option is strictly positive and increasing in search. "Just write a better reward" reduces \((1-\beta)\); it cannot reach the \(\beta = 1\) that would be needed to make the curse vanish.

References

Primary

Smith, J. E. & Winkler, R. L. (2006), "The Optimizer's Curse: Skepticism and Postdecision Surprise in Decision Analysis", Management Science 52(3), 311–322. DOI 10.1287/mnsc.1050.0451.

Manheim, D. & Garrabrant, S. (2018), "Categorizing Variants of Goodhart's Law", arXiv:1803.04585.

Background

Galton, F. (1886), "Regression towards Mediocrity in Hereditary Stature", Journal of the Anthropological Institute 15, 246–263: the origin of "regression to the mean".

Embrechts, P., Klüppelberg, C. & Mikosch, T. (1997), Modelling Extremal Events for Insurance and Finance, Springer (ISBN 978-3-540-60931-5): §3.3 for the Gaussian maximum and Gumbel norming constants.

• Pan, Bhatia & Steinhardt (2022, arXiv:2201.03544) and Skalse et al. (2022, arXiv:2209.13085): the empirical and impossibility results this derivation underpins (see Session 3.1).