@@ -625,43 +625,68 @@ We start from a benchmark with $\hat\beta R = 1$, so $\hat\beta = \beta$.
625625---
626626mystnb:
627627 figure:
628- caption: The two observational-equivalence loci
628+ caption: |
629+ Two observational-equivalence experiments. Locus I (below $\beta$) holds the
630+ *non-robust* agent fixed at $\beta R=1$ and reports the *robust* twin's
631+ discount factor $\hat\beta(\sigma)$; locus II (above $\beta$) holds the
632+ *robust* agent fixed at $\beta R=1$ and reports the *non-robust* twin's
633+ discount factor $\tilde\beta(\hat\sigma)$.
629634 name: fig-lqcs-oe-loci
630635---
631-β_bench = β # benchmark with β̂ R = 1
632-α2 = σ1**2 + (1 - β)**2 * σ2**2 # two-factor α² (see eq:bew_alpha2)
636+β_bench = β # benchmark with β R = 1
637+α2 = σ1**2 + (1 - β)**2 * σ2**2 # two-factor α^2 (see eq:bew_alpha2)
633638634639σ_hat_vals = np.linspace(0.0, -0.16, 60)
635640636-# Locus II (eq:obsequivn2): robustness ⟺ an *increase* in β (σ = 0)
641+# Locus I (eq:obseq / eq:bew_locus): non-robust agent fixed at βR=1;
642+# report the robust twin's discount factor β̂(σ) < β
643+β_hat = β_bench + σ_hat_vals * α2 * β_bench / (1 - β_bench)
644+645+# Locus II (eq:obsequivn2): robust agent fixed at βR=1;
646+# report the non-robust twin's discount factor β̃(σ̂) > β
637647disc = 1 - 4 * β_bench * (1 + σ_hat_vals * α2) / (1 + β_bench)**2
638648β_tilde = (β_bench * (1 + β_bench)) / (2 * (1 + σ_hat_vals * α2)) \
639649 * (1 + np.sqrt(disc))
640650641-# Locus I (eq:obseq / eq:bew_locus): robustness ⟺ a *decrease* in β (σ = 0)
642-β_hat = β_bench + σ_hat_vals * α2 * β_bench / (1 - β_bench)
643-644651fig, ax = plt.subplots()
645-ax.plot(-σ_hat_vals, β_tilde, lw=2, color='C0',
646- label=r'locus II: $\tilde\beta(\hat\sigma)$ (upward drift)')
647652ax.plot(-σ_hat_vals, β_hat, lw=2, color='C3',
648- label=r'locus I: $\hat\beta(\sigma)$ (downward drift)')
653+ label=r'locus I: robust twin $\hat\beta(\sigma)<\beta$'
654+ '\n(non-robust agent fixed at $\\beta R=1$)')
655+ax.plot(-σ_hat_vals, β_tilde, lw=2, color='C0',
656+ label=r'locus II: non-robust twin $\tilde\beta(\hat\sigma)>\beta$'
657+ '\n(robust agent fixed at $\\beta R=1$)')
649658ax.axhline(β_bench, color='k', linestyle=':', lw=1,
650659 label=r'benchmark $\beta$ ($\beta R = 1$)')
651-ax.set_xlabel(r'robustness concern $-\hat\sigma$')
652-ax.set_ylabel('observationally equivalent discount factor')
653-ax.legend()
660+ax.set_xlabel(r'robustness concern $|\sigma|$')
661+ax.set_ylabel('discount factor of the equivalent agent')
662+ax.legend(fontsize=8.5)
654663plt.show()
655664656665print(f"at σ̂ = {σ_hat_vals[-1]:.3f}: β̃ = {β_tilde[-1]:.4f} > β = {β_bench}")
657666print(f" β̂ = {β_hat[-1]:.4f} < β = {β_bench}")
658667```
659668660-The two loci pass through the benchmark $\beta$ at $\hat\sigma = 0$ and separate as the robustness concern grows.
669+Both loci pass through the benchmark $\beta$ at $\sigma = 0$ and separate as the robustness concern grows.
670+671+The key to reading the figure is that the two loci hold *different* agents fixed, so the discount factor plotted on the vertical axis refers to a different agent on each curve.
672+673+Locus I, from {prf:ref}`thm-lqcs-oe1`, holds the **non-robust** agent fixed at the benchmark $(\sigma = 0, \beta)$ with $\beta R = 1$ and reports the discount factor $\hat\beta(\sigma) < \beta$ of the **robust** agent that mimics it.
674+675+This is the sense in which HST call a concern for robustness observationally equivalent to a *lower* discount factor: because robustness already makes the agent save more, its discount factor must be lowered to hold the allocation at the benchmark.
676+677+Because the non-robust benchmark has $\beta R = 1$, its optimal consumption is a martingale, $\mathbb{E}_t c_{t+1} = c_t$.
678+679+The robust twin chooses the identical consumption process, so it too satisfies $\mathbb{E}_t c_{t+1} = c_t$.
680+681+The lower $\hat\beta$, which has $\hat\beta R < 1$, would on its own impart a downward drift, but the robust agent's precautionary saving offsets it exactly, leaving expected consumption flat.
682+683+Locus II, from {prf:ref}`thm-lqcs-oe2`, instead holds the **robust** agent fixed at $(\hat\sigma, \beta)$ with $\beta R = 1$ and reports the discount factor $\tilde\beta(\hat\sigma) > \beta$ of the **non-robust** agent that mimics it.
684+685+Here there is no impatience offset, so the common allocation inherits the robust agent's precautionary *upward* drift, which the non-robust twin reproduces through $\tilde\beta R > 1$.
661686662-Locus I, from {prf:ref}`thm-lqcs-oe1`, lies *below* $\beta$: activating robustness looks like an increase in impatience, which imparts a downward drift to expected consumption.
687+The two experiments encode the *same* economics: a concern for robustness adds precautionary saving that acts like extra patience.
663688664-Locus II, from {prf:ref}`thm-lqcs-oe2`, lies *above* $\beta$: the same robustness concern, viewed from a benchmark with $\beta R = 1$, looks like an increase in patience, which imparts an upward drift.
689+They differ only in which agent is anchored at $\beta R = 1$, and hence in whether the common saving motive shows up as an exactly-offsetting impatience adjustment (locus I, expected consumption flat) or as an upward drift in expected consumption (locus II).
665690666691### A robust LQ Bewley model
667692@@ -670,7 +695,7 @@ Locus II, from {prf:ref}`thm-lqcs-oe2`, lies *above* $\beta$: the same robustnes
670695671696We now synthesise the lecture by embedding the Bewley economy of {doc}`lq_bewley_complete_markets` into the HST framework and applying the observational-equivalence theorem.
672697673-In this way, we construct a family of **robust Bewley economies**, parameterised by a robustness level $\sigma \leq 0$, whose equilibrium quantities are identical to those of the plain vanilla Bewley model.
698+We shall construct a family of **robust Bewley economies**, parameterised by a robustness level $\sigma \leq 0$, whose equilibrium quantities are identical to those of the plain vanilla Bewley model.
674699675700We first map the Bewley economy into HST notation, specialising the robust model to
676701$\lambda = \delta_h = 0$ (no habits, no durable goods) and to a