Exploring Recursive Utility

8. Exploring Recursive Utility#

Authors: Jaroslav Borovicka (NYU), Lars Peter Hansen (University of Chicago) and Thomas J. Sargent (NYU)

8.1. Introduction#

We explore the recursive utility preference specification of [Kreps and Porteus, 1978] and [Epstein and Zin, 1989] including versions of them proposed by [Hansen and Sargent, 2001] and [Anderson et al., 2003] that capture concerns about model misspecification as in. We deploy two distinct approximation approaches. One approach builds on a characterization by [Duffie and Epstein, 1992] and uses a continuous-time limiting approximation to a discrete-time specification in which underlying shocks are normally distributed, though our detailed derivations differ from those of [Duffie and Epstein, 1992]. We represent the limiting approximation with a Brownian motion information structure. Our second approximation m allows macroeconomic uncertainty to have first-order implications. We show how to explore model implications for (nonstandard) first- and second-order approximations to equilibria of dynamic stochastic models. We modify first- and second-order approximations routinely used in the macroeconomics literature in ways designed to focus on macroeconomic uncertainty. Our approximations apply to production-based macro-finance models with opportunities to invest in different kinds of capital. Both of our approximation approaches capture adjustments for uncertainty with a change in probability measure, one that differs from one that leads the so-called risk neutral distribution widely used to price derivative claims. Such a change of measure plays a big role in specifications of preferences like those in [Hansen and Sargent, 2001] and [Anderson et al., 2003] that build on a robust control literature initiated by [Jacobson, 1973] and [Whittle, 1981].

Our approximations allow macroeconomic uncertainty to have first-order implications. We assume risk averse economic decision makers with recursive utility preferences or closely related preferences that also express concerns about model misspecification. We show how to explore model implications for (nonstandard) first and second-order approximations to equilibria of dynamic stochastic models. We modify first and second order approximations routinely used in the macroeconomics literature in ways that focus on macroeconomic uncertainty. Our uncertainty adjustments are captured by a change in probability measure. Furthermore, our approximations apply to production-based macro-finance models with opportunities to invest in diverse forms of capital. We use this framework to advance our understanding of alternative preference specifications and their implications for production and asset pricing.

We extend work by [Schmitt-Grohé and Uribe, 2004] and [Lombardo and Uhlig, 2018] in ways that highlight consequences of uncertainty. We design approximations to make implied stochastic discount factors reside within an exponential linear quadratic class, a class that gives rise tractable formulas for asset valuation over alternative investment horizons. See, for instance, [Ang and Piazzesi, 2003] and [Borovička and Hansen, 2014]. The class is also useful for studying production-based macro-finance models with opportunities to invest in different forms of capital.

Note

We of course recognize that nonlinearities in some models are more accurately captured by global solution methods.

8.2. Recursive utility valuation process#

We construct continuation value and stochastic discount processes, important constituents of many dynamic stochastic models in macroeconomics and finance.

8.2.1. Basic recursion#

A homogeneous of degree one representation of recursive utility is

(8.1)#

V_{t} = {[(1 - β) {(C_{t})}^{1 - ρ} + β {(R_{t})}^{1 - ρ}]}^{\frac{1}{1 - ρ}}

where

(8.2)#

R_{t} = {(E [{(V_{t + 1})}^{1 - γ} ∣ A_{t}])}^{\frac{1}{1 - γ}} .

Value $V_{t}$ defined in equation (8.1) is a homogeneous of degree one function of $C_{t}$ and $R_{t}$ ; equation (8.2) defines $R_{t}$ as a homogeneous of degree one function of another function of $V_{t + 1}$ .

In equation (8.1), $0 < β < 1$ is a subjective discount factor and $\frac{1}{ρ}$ is the elasticity of intertemporal substitution, while $γ$ in equation (8.2) describes attitudes towards risk.

Continuation values are determined only up to an increasing transformation. For computational and conceptual reasons, it is useful to work with the transformation ${\hat{V}}_{t} = \log V_{t}$ . Recursions for ${\hat{V}}_{t}$ expressed in terms of the logarithm of consumption ${\hat{C}}_{t}$ are

(8.3)#

{\hat{V}}_{t} = \frac{1}{1 - ρ} \log [(1 - β) \exp [(1 - ρ) {\hat{C}}_{t}] + β \exp [(1 - ρ) {\hat{R}}_{t}]]

where

(8.4)#

{\hat{R}}_{t} = \frac{1}{1 - γ} \log E [\exp ((1 - γ) {\hat{V}}_{t + 1}) ∣ A_{t}] .

The right side of recursion (8.3) is the logarithm of a constant elasticity of substitution (CES) function of $\exp ({\hat{C}}_{t})$ and $\exp ({\hat{R}}_{t})$ .

Remark 8.1

The limit of ${\hat{R}}_{t}$ as $γ$ approaches $1$ is expected logarithmic utility:

lim_{γ ↓ 1} {\hat{R}}_{t} = lim_{γ ↓ 1} \frac{\log E (\exp [(1 - γ) {\hat{V}}_{t + 1}] | A_{t})}{1 - γ} = E ({\hat{V}}_{t + 1} | A_{t}) .

We shall construct small noise expansions of ${\hat{V}}_{t}$ and ${\hat{R}}_{t}$ separately, then assemble them appropriately. Before doing so, we offer a reinterpretation of our recursions (8.3)-(8.4).

8.2.2. Preference for robustness#

When $γ > 1$ , (8.4) emerges as an indirect utility function for a robust control problem in which $\frac{1}{γ - 1}$ serves as a penalty parameter on entropy of a baseline model relative to an alternative model. The penalty parameter constrains a set of alternative probability models that a decision maker considers when evaluating consumption processes. This interpretation of (8.4) as an indirect utilty function from a minimization problem originated in [Hansen and Sargent, 1995], which built [Jacobson, 1973] and [Whittle, 1981].

Let the random variable $N_{t + 1} \geq 0$ satisfy $E (N_{t + 1} ∣ A_{t}) = 1$ so that it is a likelihood ratio. Think of replacing the expected continuation value $E ({\hat{V}}_{t + 1} ∣ A_{t})$ by

(8.5)#

min_{N_{t + 1} \geq 0, E (N_{t + 1} | A_{t}) = 1} E (N_{t + 1} {\hat{V}}_{t + 1} ∣ A_{t}) + ξ E (N_{t + 1} \log N_{t + 1} ∣ A_{t}),

where $ξ$ is a parameter that penalizes departures of $N_{t + 1}$ from unity as measured by relative entropy. Conditional entropy relative to an alternative conditional probability induced by applying change of measure $N_{t + 1}$ is

E (N_{t + 1} \log N_{t + 1} ∣ A_{t}) \geq 0 .

The inequality follows from Jensen’s inequality because the function $n \log n$ is convex. Relative entropy is evidently zero when $N_{t + 1} = 1$ .

Relative entropy measures a discrepancy between two probability distributions. Think of $N_{t + 1}$ as a relative conditional likelihood ratio of an alternative model vis-a-vis a baseline model. Then $E (N_{t + 1} \log N_{t + 1} ∣ A_{t})$ is an expected (conditional) log-likelihood ratio of the alternative model when the expectation is taken using the alternative probability model. A small expected log likelihood indicates a small discrepancy between two models, i.e., two probability distributions.

Remark 8.2

To solve minimization problem (8.5), we attach a Lagrange multiplier $ℓ$ to the conditional expectation constraint and form a Lagrangian that we want to minimize with respect to the random variable $N$ and maximize with respect to $ℓ$ . This extremization problem separates across states, so to minimze with respect to $N$ we can solve a set of minimum problems

min_{n} n v + ξ n \log n + ℓ (n - 1),

where $n$ is a potential realization of $N_{t + 1}$ and $v$ is a realization of $V_{t + 1} .$ First-order necessary conditions are:

v + ξ + ξ \log n + ℓ = 0,

which implies the minimizer

n^{*} = \exp [- \frac{1}{ξ} (v + ℓ + ξ)],

and the minimized objective

- ξ \exp [- \frac{1}{ξ} (v + ℓ + ξ)] - ℓ .

To determine $ℓ,$ we pose

max_{ℓ} - ξ E (\exp [- \frac{1}{ξ} ({\hat{V}}_{t + 1} + ℓ + ξ)] ∣ A_{t}) - ℓ

whose first-order necessary condition is

E (\exp [- (\frac{1}{ξ}) {\hat{V}}_{t + 1}] ∣ A_{t}) \exp [- (\frac{ℓ + ξ}{ξ})] - 1 = 0.

The maximizing $ℓ$ is

ℓ^{*} = ξ \log E (\exp [- (\frac{1}{ξ}) {\hat{V}}_{t + 1}] ∣ A_{t}) - 1

and the minimized objective function is

- ξ \log E (\exp [- (\frac{1}{ξ}) {\hat{V}}_{t + 1}] ∣ A_{t}) .

The minimizing $N_{t + 1}$ is

(8.6)#

N_{t + 1}^{*} = \frac{\exp (- \frac{1}{ξ} {\hat{V}}_{t + 1})}{E [\exp (- \frac{1}{ξ} {\hat{V}}_{t + 1}) ∣ A_{t}]} .

The minimizer (8.6) of problem (8.5) evidently shifts probabilities toward low continuation values. [Bucklew, 2004] called this a stochastic version of Murphy’s law. Notice that the minimized objective satisfies

- ξ \log E [\exp (- \frac{1}{ξ} {\hat{V}}_{t + 1}) ∣ A_{t}] = {\hat{R}}_{t}

where earlier we described ${\hat{R}}_{t}$ by equation (8.4) after we set $ξ = \frac{1}{γ - 1}$ .
It follows from (8.6) that

N_{t + 1}^{*} = \exp [- \frac{1}{ξ} ({\hat{V}}_{t + 1} - {\hat{R}}_{t})] .

The random variable $N_{t + 1}^{*}$ will appear often below.

8.2.3. Stochastic discount factor process#

A stochastic discount factor (SDF) process $S = {S_{t} : t \geq 0}$ describes a consumer’s attitudes about small changes in uncertainty. SDF processes have several uses. First, they contain shadow prices that tell how a consumer’s attitudes about uncertainty shape marginal valuations of risky assets. Second, they shape first-order necessary conditions for optimally choosing financial and physical investments. Third, they underlie tractable formulas for equilibrium asset prices. Fourth, they can help construct Pigouvian taxes that ameliorate externalities under uncertainty. Fifth, they can help evaluate effects of small (local) changes in government policies.

To deduce an SDF process, we positing that a date zero value of a risky date $t$ consumption payout $χ_{t}$ is

(8.7)#

π_{0}^{t} (χ_{t}) = E [(\frac{S_{t}}{S_{0}}) χ_{t} | A_{0}] .

To compute the ratio $\frac{S_{t}}{S_{0}}$ that appears in formula (8.7) we evaluate the slope of an indifference curve that runs through both a baseline consumption process ${C_{t}}_{t = 0}^{\infty}$ and a perturbed consumption process

(C_{0} - P_{0} (q), C_{1}, C_{2}, \dots, C_{t} + q χ_{t}, C_{t + 1}, . . .) .

We think of $q$ as parameterizing an indifference curve, so $P_{0} (q)$ expresses how much current period consumption must be reduced to keep a consumer on the same indifference curve after we replace $C_{t}$ by $C_{t} + q χ_{t}$ . We set $π_{0}^{t} (χ_{t})$ defined in equation (8.7) equal to the slope of that indifference curve:

π_{0}^{t} (χ_{t}) = {\frac{d}{d q} P_{0} (q) |}_{q = 0} .

For recursive utility, a one-period increment in the stochastic discount factor process is

(8.8)#

\begin{aligned} \frac{S_{t + 1}}{S_{t}} = & β {(\frac{C_{t + 1}}{C_{t}})}^{- ρ} \exp [(1 - γ) ({\hat{V}}_{t + 1} - {\hat{R}}_{t})] \exp [(ρ - 1) ({\hat{V}}_{t + 1} - {\hat{R}}_{t})] \\ = & β N_{t + 1}^{*} \exp ({\hat{S}}_{t + 1} - {\hat{S}}_{t}), \end{aligned}

where

(8.9)#

{\hat{S}}_{t + 1} - {\hat{S}}_{t} \overset{def}{=} - ρ ({\hat{C}}_{t + 1} - {\hat{C}}_{t}) + (ρ - 1) ({\hat{V}}_{t + 1} - {\hat{R}}_{t})

and $N_{t + 1}^{*}$ induces the change of probability measure that equation (8.6) presents as the outcome of a robust valuation problem. We will use the second line in (8.8) in what follows. We will think of $β$ as a subjective discount factor adjustment, $N_{t + 1}^{*}$ as a change-of-measure adjustment for uncertainty, and $\exp ({\hat{S}}_{t + 1} - {\hat{S}}_{t})$ as an adjustment for the elasticity of intertemporal substitution. We interpret the twisted transition probability measure induced by $N_{t + 1}^{*}$ as an adjustment for uncertainty in valuation.
The recursive structure of preferences makes the time $t$ period stochastic discount factor $\frac{S_{t}}{S_{0}}$ be the product of the respective one-period stochastic discount factor increments. Similarly, we can compound one-period transition uncertainty measures into multiple time-horizon measures of uncertainty.

Remark 8.3

To verify formula (8.8), we compute a one-period intertemporal marginal rate of substitution. Given the valuation recursions (8.3) and (8.4), we construct two marginal utilities familiar from CES and exponential utility:

m c = (1 - β) {(c)}^{- ρ} \exp [(ρ - 1) \hat{v}]

m \hat{r} = β \exp [(1 - ρ) (\hat{r} - \hat{v})]

From the certainty equivalent formula, we construct the marginal utility of the next-period logarithm of the continuation value:

m {\hat{v}}^{+} = \exp [(1 - γ) ({\hat{v}}^{+} - \hat{r})]

where the $+$ superscript is used to denote the next-period counterpart. In addition, the next-period marginal utility of consumption is

m c^{+} = (1 - β) {(c^{+})}^{- ρ} \exp [(ρ - 1) {\hat{v}}^{+}]

Putting these four formulas together using the chain rule for differentiation gives a marginal rate of substitution:

\frac{(m r) (m v^{+}) (m c^{+})}{m c} = β {(\frac{c^{+}}{c})}^{- ρ} \exp [(1 - γ) ({\hat{v}}^{+} - \hat{r})] \exp [(ρ - 1) ({\hat{v}}^{+} - \hat{r})] .

Now let ${\hat{v}}^{+} = {\hat{V}}_{t + 1}$ , $c^{+} = C_{t + 1}$ , $C_{t} = c$ and $\hat{r} = {\hat{R}}_{t}$ to obtain the formula for the one-period stochastic discount factor (8.8).

8.3. Continuous-time limit#

We now explore a continuous-time limit that approximates a discrete-time specification. Since we continue to work with normal shocks, the continuous-time counterpart to these shocks are Brownian increments. The continuation value in continuous time will evolve as:

d V_{t} = V_{t} μ_{t}^{V} d t + V_{t} σ_{t}^{V} \cdot d W_{t}

for some drift (local mean), $V_{t} μ_{t}^{V}$ and some local shock exposure vector $V_{t} σ_{t}^{V}$ , where ${W_{t} : t \geq 0}$ is a multivariate Brownian motion. The scaling for the local evolution coefficients by $V_{t}$ is done for convenience, where the continuation value process is presumed to be positive. As in discrete time, it is convenient to work with the logarithm of the continuation value process (a strictly increasing transformation). The implied evolution is

d {\hat{V}}_{t} = {\hat{μ}}_{t}^{V} d t + σ_{t}^{V} \cdot d W_{t}

where ${\hat{μ}}_{t} = μ_{t}^{V} - \frac{1}{2} ∣ σ_{t}^{V} ∣^{2} .$ This adjustment follows from the well known Ito’s formula.

8.3.1. Discrete-time approximation#

To study the utility recursion, start with a discrete-time specification:

\begin{aligned} \frac{1}{1 - ρ} \log [(1 - β_{ϵ}) \exp [(1 - ρ) ({\hat{C}}_{t} - {\hat{V}}_{t})] + β_{ϵ} \exp [(1 - ρ) ({\hat{R}}_{t} - {\hat{V}}_{t})]] = 0 \\ {\hat{R}}_{t} - {\hat{V}}_{t} = \frac{1}{1 - γ} \log E (\exp [(1 - γ) ({\hat{V}}_{t + ϵ} - {\hat{V}}_{t})] ∣ A_{t}) \end{aligned}

where $β_{ϵ} = \exp (- δ ϵ)$ and $δ > 0$ is instantaneous subjective rate of discount. Consider the time derivative of the second recursion:

\begin{aligned} \frac{d}{d ϵ} {({\hat{R}}_{t} - {\hat{V}}_{t}) |}_{ϵ = 0} & = \frac{d}{d ϵ} {\frac{1}{1 - γ} \log E (\exp [(1 - γ) ({\hat{V}}_{t + ϵ} - {\hat{V}}_{t})] ∣ A_{t}) |}_{ϵ = 0} \\ = {\hat{μ}}_{t}^{V} + \frac{(1 - γ)}{2} | σ_{t}^{V} |^{2} \\ = μ_{t}^{V} - \frac{γ}{2} | σ_{t}^{V} |^{2} . \end{aligned}

by local log normality. We turn now to the first recursion and compute time derivatives in three steps. First, we evaluate the term inside the logarithm as $ϵ$ tends to zero :

lim_{ϵ ↓ 0} \exp [(1 - ρ) ({\hat{R}}_{t} - {\hat{V}}_{t})] = 1.

This term is in the denominator as implied by the derivative with respect to a logarithm. Second, we differentiate the term inside the logarithm with respect to $ϵ$ as contributed by $β_{ϵ}$ :

δ \exp [(1 - ρ) ({\hat{C}}_{t} - {\hat{V}}_{t})] - δ .

Third, we differentiate the term inside the logarithm with respect to $ϵ$ as contributed by

\frac{d}{d ϵ} \exp [(1 - ρ) ({\hat{R}}_{t} - {\hat{V}}_{t})] = {\hat{μ}}_{t}^{V} + \frac{(1 - γ)}{2} | σ_{t}^{V} |^{2} .

Putting together the derivative components together gives:

(8.10)#

\begin{aligned} 0 = & \frac{δ [{(\frac{C_{t}}{V_{t}})}^{1 - ρ} - 1]}{1 - ρ} + {\hat{μ}}_{t}^{V} + \frac{(1 - γ)}{2} | σ_{t}^{V} |^{2} \\ = & \frac{δ [{(\frac{C_{t}}{V_{t}})}^{1 - ρ} - 1]}{1 - ρ} + μ_{t}^{V} - \frac{γ}{2} | σ_{t}^{V} |^{2} . \end{aligned}

Notice that this relation imposes a restriction across the local mean, $μ_{t}^{V}$ and local variance, $| σ_{t}^{V} |^{2}$ of the continuation value. [Duffie and Epstein, 1992] refer to $γ$ as a variance multiplier where larger values of $γ$ imply a more substantial adjustment for local volatility. As in discrete-time, the $ρ = 1$ is interesting special case represented as;

0 = - δ ({\hat{V}}_{t} - {\hat{C}}_{t}) + {\hat{μ}}_{t}^{V} + \frac{(1 - γ)}{2} | σ_{t}^{V} |^{2} .

8.3.2. Robustness to misspecification#

To investigate an aversion to model misspecification in continuous time, we now treat the distribution of ${W_{t} : t \geq 0}$ as uncertain. We allow for probability measures that entertain possible Brownian motions with local means or drifts that are history dependent.

We start by considering positive martingales ${M_{t}^{H} : t \geq 0}$ parameterized by alternative ${H_{t} : t \geq 0}$ processes with the same dimension as the underlying Brownian motion.
The martingales have local evolutions:

d M_{t}^{H} = M_{t}^{H} H_{t} \cdot d W_{t},

and we initialize them at $M_{0}^{H} = 1$ . Observe that by applying Ito’s formula, $\log M^{H}$ evolves as:

d \log M_{t}^{H} = - \frac{1}{2} {| H_{t} |}^{2} d t + H_{t} \cdot d W_{t} .

We use these martingales as relative densities or likelihood ratios.

Write the discrete-time counterpart as

\log M_{t + ϵ}^{H} - \log M_{t}^{H} = - \frac{ϵ}{2} {| H_{t} |}^{2} + H_{t} \cdot (W_{t + ϵ} - W_{t})

Let $w$ be a realized $W_{t + ϵ} - W_{t}$ and $h$ be a realization $H_{t}$ . Then $\log M_{t + ϵ} - \log M_{t}$ contributes $- \frac{ϵ}{2} h^{'} h + h \cdot w$ to the log-likelihood. The standard normal density for $\sqrt{ϵ} (W_{t + ϵ} - W_{t})$ contributes $- \frac{1}{ϵ} w^{'} w - \log (2 π ϵ)$ . Put there two components together, we have a log-likelihood:

- \frac{ϵ}{2} h^{'} h + h \cdot w - \frac{1}{2 ϵ} w^{'} w - \log (2 π ϵ) = - \frac{1}{2 ϵ} (w - ϵ h)^{'} (w - ϵ h) - \log (2 π ϵ) .

The altered conditional density has mean $ϵ h$ , which is the realized value of $ϵ H_{t}$ with the same conditional covariance matrix as before. Moreover, the conditional expectation of

E [(\frac{M_{t + ϵ}^{H}}{M_{t}^{H}}) (\log M_{t + ϵ}^{H} - \log M_{t}^{H}) ∣ A_{t}] = \frac{ϵ}{2} {| H_{t} |}^{2},

which measures the statistical divergence or relative entropy between original and altered conditional probabilities.

For the continuous-time limit, under the $H$ change of probability measure:

d W_{t} = H_{t} d t + d {\tilde{W}}_{t}^{H}

where ${\tilde{W}}^{H}$ is a standard Brownian motion. Thus the potential changes of probability measures induce a local means or drift $H$ processes to the Brownian motion. The continuous-time counterpart to conditional relative entropy at time $t$ is $\frac{1}{2} {| H_{t} |}^{2} .$

We may justify focusing on drift distortions for Brownian increments because of our imposition of absolute continuity of the alternative probabilities with respect to the baseline specification of a multivariate standard Brownian motion. This is an implication of the Girsanov Theorem.

We are now in a position to deduce a robustness adjustment in continuous time. Consider formula (8.10) when $γ = 1$ modified for a potential change in the probability measure

0 = \frac{δ [{(\frac{C_{t}}{V_{t}})}^{1 - ρ} - 1]}{1 - ρ} + {\hat{μ}}_{t}^{V} + σ_{T}^{V} \cdot H_{t}

Modify this equation to include minimization over $H_{t}$ subject to a relative entropy penalty $\frac{ξ}{2} | H_{t} |^{2} :$

0 = min_{H_{t}} \frac{δ [{(\frac{C_{t}}{V_{t}})}^{1 - ρ} - 1]}{1 - ρ} + {\hat{μ}}_{t}^{V} + σ_{t}^{V} \cdot H_{t} + \frac{ξ}{2} H_{t} \cdot H_{t} .

The minimizing solution is

H_{t}^{*} = - \frac{1}{ξ} σ_{t}^{V}

with a minimized objective

0 = \frac{δ [{(\frac{C_{t}}{V_{t}})}^{1 - ρ} - 1]}{1 - ρ} + {\hat{μ}}_{t}^{V} - \frac{1}{2 ξ} {| σ_{t}^{V} |}^{2}

Notice that this agrees with formula (8.10) for $γ - 1 = 1 / ξ$ . The explicit link is entirely consistent with our discrete-time equivalence result. By taking the continuous-time limit, we are able to focus our misspecification analysis on changing local means of the underlying Brownian increments.

8.3.3. Uncertainty pricing#

For the purposes of valuation, we compound the equilibrium version of ${H_{t}^{*} : t \geq 0}$ to give a exponential martingale:

M_{t}^{*} = \exp (\int_{0}^{t} H_{τ}^{*} d W_{τ} - \frac{1}{2} \int_{0}^{t} {| H_{τ}^{*} |}^{2} d τ)

provided that the constructed process is a martingale.[1] With this construction we interpret $- H_{t}^{*}$ as the vector of local uncertainty prices that give compensations for exposure to Brownian increment uncertainty. These compensations are expressed as changes in conditional means under the baseline distribution as is typical in continuous-time asset pricing.

8.4. Small noise expansion of dynamic stochastic equilibria#

We next consider a different type of characterization that sometimes gives good approximations for dynamic stochastic equilibrium models. While the approximations build from derivations in [Schmitt-Grohé and Uribe, 2004] and [Lombardo and Uhlig, 2018], and we extend them in a way that features the uncertainty contributions more prominently and are reflected even in first-order contributions. By design, the implied approximations of stochastic discount factors used to represent market or shadow values reside within the exponential linear quadratic class. This class is known to give tractable formulas for asset valuation over alternative investment horizons. See, for instance, [Ang and Piazzesi, 2003] and [Borovička and Hansen, 2014]. Moreover, they are applicable to production-based macro-finance models with investment opportunities in alternative forms of capital.

While these approximations are tractable for the reasons described, researchers may be concerned about some more fundamental aspects of uncertainty that are disguised by these methods. Indeed, for some models nonlinearities are more accurately captured by global solution methods. Nevertheless, the approximations still provide further understanding of the preferences and their implications for asset pricing in endowment and production economies.

8.4.1. Approximate state dynamics#

We follow [Lombardo and Uhlig, 2018] by considering the following class of stochastic processes indexed by a scalar perturbation parameter $q$ :[2]

(8.11)#

X_{t + 1} (q) = ψ [X_{t} (q), q W_{t + 1}, q] .

Here $X$ is an $n$ -dimensional stochastic process and ${W_{t + 1}}$ is an i.i.d.~normally distributed random vector with conditional mean vector $0$ and conditional covariance matrix $I$ . We parameterize this family so that $q = 1$ gives the model of interest.

We denote a zero-order expansion $q = 0$ limit as:

(8.12)#

X_{t + 1}^{0} = ψ (X_{t}^{0}, 0, 0),

and assume that there exists a second-order expansion of $X_{t}$ around $q = 0$ :

(8.13)#

X_{t} \approx X_{t}^{0} + q X_{t}^{1} + \frac{q^{2}}{2} X_{t}^{2}

where $X_{t}^{1}$ is a first-order contribution and $X_{t}^{2}$ is a second-order contribution. In other words, the stochastic processes $X^{j}$ , $j = 0, 1, 2$ are appropriate derivatives of $X$ with respect to the perturbation parameter $q$ evaluated at $q = 0.$

In the remainder of this chapter, we shall construct instances of the second-order expansion (8.13) in which the generic random variable $X_{t}$ is replaced, for example, by the logarithm of consumption, a value function, and so on. I

Processes $X_{t}^{j}, j = 0, 1, 2$ have a recursive structure: first compute the stochastic process $X_{t}^{0},$ then the process $X_{t}^{1}$ next (it depends on $X_{t}^{0}$ ), and finally the process $X_{t}^{2}$ (it depends on both $X_{t}^{0}$ and $X_{t}^{1}$ ).

In this chapter, we use a prime $^{'}$ to denote a transpose of a matrix or vector. When we include $x^{'}$ in a partial derivative of a scalar function it means that the partial derivative is a row vector. Consistent with this convention, let $ψ_{x^{'}}^{i}$ , the $i^{t h}$ entry of $ψ_{x^{'}}$ , denote the row vector of first derivatives with respect to the vector $x$ , and similarly for $ψ_{w^{'}}^{i}$ . Since $q$ is scalar, $ψ_{q}^{i}$ is the scalar derivative with respect to $q$ . Derivatives are evaluated at $X_{t}^{0}$ , which in many examples is invariant over time, unless otherwise stated. This invariance follows when we impose a steady state on the deterministic system.

The first-derivative process obeys a recursion

(8.14)#

\begin{array}{r} X_{t + 1}^{1} = [\begin{array}{c} ψ_{x^{'}}^{1} \\ ψ_{x^{'}}^{2} \\ ⋮ \\ ψ_{x^{'}}^{n} \end{array}] X_{t}^{1} + [\begin{array}{c} ψ_{w^{'}}^{1} \\ ψ_{w^{'}}^{2} \\ ⋮ \\ ψ_{w^{'}}^{n} \end{array}] W_{t + 1} + [\begin{array}{c} ψ_{q}^{1} \\ ψ_{q}^{2} \\ ⋮ \\ ψ_{q}^{n} \end{array}] \end{array}

that we can write compactly as the following a first-order vector autoregression:

X_{t + 1}^{1} = ψ_{x^{'}} X_{t}^{1} + ψ_{w^{'}} W_{t + 1} + ψ_{q}

We assume that the matrix $ψ_{x}^{'}$ is stable in the sense that all of its eigenvalues are strictly less than one in modulus.

It is natural for us to denote second derivative processes with double subscripts. For instance, for the double script used in conjunction with the second derivative matrix of $ψ^{i}$ , the first subscript without a prime ( $^{'}$ ) reports the row location; second subscript with a prime ( $^{'}$ ) reports the column location. Differentiating recursion (8.14) gives:

(8.15)#

\begin{array}{r} \begin{aligned} X_{t + 1}^{2} = ψ_{x^{'}} X_{t}^{2} & + [\begin{array}{c} X_{t}^{1^{'}} ψ_{x x^{'}}^{1} X_{t}^{1} \\ X_{t}^{1^{'}} ψ_{x x^{'}}^{2} X_{t}^{1} \\ ⋮ \\ X_{t}^{1^{'}} ψ_{x x^{'}}^{n} X_{t}^{1} \end{array}] + 2 [\begin{array}{c} X_{t}^{1^{'}} ψ_{x w^{'}}^{1} W_{t + 1} \\ X_{t}^{1^{'}} ψ_{x w^{'}}^{2} W_{t + 1} \\ ⋮ \\ X_{t}^{1^{'}} ψ_{x w^{'}}^{n} W_{t + 1} \end{array}] + [\begin{array}{c} {W_{t + 1}}^{'} ψ_{w w^{'}}^{1} W_{t + 1} \\ {W_{t + 1}}^{'} ψ_{w w^{'}}^{2} W_{t + 1} \\ ⋮ \\ {W_{t + 1}}^{'} ψ_{w w^{'}}^{n} W_{t + 1} \end{array}] \\ + 2 [\begin{array}{c} ψ_{q x^{'}}^{1} X_{t}^{1} \\ ψ_{q x^{'}}^{2} X_{t}^{1} \\ ⋮ \\ ψ_{q x^{'}}^{n} X_{t}^{1} \end{array}] + 2 [\begin{array}{c} ψ_{q w^{'}}^{1} W_{t + 1} \\ ψ_{q w^{'}}^{2} W_{t + 1} \\ ⋮ \\ ψ_{q w^{'}}^{n} W_{t + 1} \end{array}] + [\begin{array}{c} ψ_{q q}^{1} \\ ψ_{q q}^{2} \\ ⋮ \\ ψ_{q q}^{n} \end{array}] . \end{aligned} \end{array}

Recursions (8.14) and (8.15) have a linear structure with some notable properties. The law of motion for $X^{0}$ is deterministic and is time invariant if (8.11) comes from a stationary ${X_{t}}$ process. The dynamics for $X^{2}$ are nonlinear only in $X^{1}$ and $W_{t + 1}$ . Thus, the stable dynamics for $X^{1}$ that prevail when $ψ_{x}$ is a stable matrix imply stable dynamics for $X^{2}$ .

Remark 8.4

Perturbation methods have been applied to many rational expectations models in which partial derivatives of $ψ$ with respect to $q$ are often zero.[3] However, derivatives of $ψ$ with respect to $q$ are not zero in production-based equlibrium models with the robust or recursive utility specifications that we shall study here.

Let $C$ denote consumption and $\hat{C}$ the logarithm of consumption. Suppose that the logarithm of consumption evolves as:

{\hat{C}}_{t + 1} - {\hat{C}}_{t} = κ (X_{t}, q W_{t + 1}, q) .

Approximate this process by:

(8.16)#

{\hat{C}}_{t + 1} - {\hat{C}}_{t} \approx {\hat{C}}_{t + 1}^{0} - {\hat{C}}_{t}^{0} + q ({\hat{C}}_{t + 1}^{1} - {\hat{C}}_{t}^{1}) + \frac{q^{2}}{2} ({\hat{C}}_{t + 1}^{2} - {\hat{C}}_{t}^{2})

where

{\hat{C}}_{t + 1}^{0} - {\hat{C}}_{t}^{0} = κ (X_{t}^{0}, 0, 0) := η_{0}^{c}

{\hat{C}}_{t + 1}^{1} - {\hat{C}}_{t}^{1} = κ_{x^{'}} X_{t}^{1} + κ_{w^{'}} W_{t + 1} + κ_{q}

\begin{aligned} {\hat{C}}_{t + 1}^{2} - {\hat{C}}_{t}^{2} = & κ_{x^{'}} X_{t}^{2} + {X_{t}^{1}}^{'} κ_{x, x^{'}} X_{t}^{1} + 2 κ_{q, x^{'}} X_{t}^{1} + κ_{q q} \\ + 2 {X_{t}^{1}}^{'} κ_{x w^{'}} W_{t + 1} + {W_{t + 1}}^{'} κ_{w w^{'}} W_{t + 1} + 2 κ_{q w^{'}} W_{t + 1} . \end{aligned}

In models with endogenous investment and savings, the consumption dynamics and some of the state dynamics will emerge as the solution to a dynamic stochastic equilibrium model. We use the approximating processes (8.13) and (8.16) as inputs into the construction of an approximating continuation value process and its risk-adjusted counterpart for recursive utility preferences.

8.5. Incorporating preferences with enhanced uncertainty concerns#

To approximate the recursive utility process, we deviate from common practice in macroeconomics by letting the risk aversion or robust parameter in preferences depend on $q :$

ξ = q ξ_{o} γ - 1 = \frac{γ_{o} - 1}{q}

The aversion to model misspecification or the aversion to risk moves inversely with the parameter $q$ when we embed the model of interest within a parameterized family of models. In effect, the variable $q$ is doing double duty. Reducing $q > 0$ limits the overall exposure of the economy to the underlying shocks. This is offset by letting the preferences include a greater aversion to uncertainty. This choice of any expansion protocol has significant and enlightening consequences for continuation value processes and for the minimizing $N$ process used to alter expectations. It has antecedents in the control theory literature, and it has the virtue that implied uncertainty adjustments occur more prominently at lower-order terms in the approximation.

8.5.1. Order-zero#

Write the order-zero expansion of (8.3) as

{\hat{V}}_{t}^{0} = \frac{1}{1 - ρ} \log [(1 - β) \exp [(1 - ρ) {\hat{C}}_{t}^{0}] + β \exp [(1 - ρ) {\hat{R}}_{t}^{0}]]

{\hat{R}}_{t}^{0} = {\hat{V}}_{t + 1}^{0},

where the second equation follows from noting that randomness vanishes in the limit as $q$ approaches $0$ .

For order zero, write the consumption growth-rate process as

{\hat{C}}_{t + 1}^{0} - {\hat{C}}_{t}^{0} = η_{c}^{0} .

The order-zero approximation of (8.3) is:

{\hat{V}}_{t}^{0} - {\hat{C}}_{t}^{0} = \frac{1}{1 - ρ} \log [(1 - β) + β \exp [(1 - ρ) ({\hat{V}}_{t + 1}^{0} - {\hat{C}}_{t + 1}^{0} + η_{c}^{0})]]

We guess that ${\hat{V}}_{t}^{0} - {\hat{C}}_{t}^{0} = η_{v - c}^{0}$ and will have verified the guess once we solve:

\exp [(1 - ρ) (η_{v - c}^{0})] = (1 - β) + β \exp [(1 - ρ) (η_{v - c}^{0})] \exp [(1 - ρ) η_{c}^{0}] .

This equation implies

(8.17)#

\exp [(1 - ρ) (η_{v - c}^{0})] = \frac{1 - β}{1 - β \exp [(1 - ρ) η_{c}^{0}]} .

Equation (8.17) determines $η_{v - c}^{0}$ as a function of $η_{c}^{0}$ and the preference parameters $ρ, β$ , but not the risk aversion parameter $γ$ or its robust counterpart $ξ$ . Specifically,

(8.18)#

η_{v - c}^{0} = \frac{\log (1 - β) - \log (1 - β \exp [(1 - ρ) η_{c}^{0}])}{1 - ρ}

In the limiting $ρ = 1$ case,

η_{v - c}^{0} = (\frac{β}{1 - β}) η_{c}^{0} .

8.5.2. Order-one#

We temporarily take ${\hat{R}}_{t}^{1} - {\hat{C}}_{t}^{1}$ as given. We construct a first-order approximation to the nonlinear utility recursion (8.3)

(8.19)#

{\hat{V}}_{t}^{1} - {\hat{C}}_{t}^{1} = λ ({\hat{R}}_{t}^{1} - {\hat{C}}_{t}^{1})

where

\begin{aligned} λ = & [\frac{β \exp [(1 - ρ) (η_{v - c}^{0} + η_{c}^{0})]}{(1 - β) + β \exp [(1 - ρ) (η_{v - c}^{0} + η_{c}^{0})]}] \\ = & [\frac{β \exp [(1 - ρ) η_{c}^{0}]}{(1 - β) \exp [- (1 - ρ) η_{v - c}^{0}] + β \exp [(1 - ρ) η_{c}^{0}]}] \\ = & [\frac{β \exp [(1 - ρ) η_{c}^{0}]}{1 - β \exp [(1 - ρ) η_{c}^{0}] + β \exp [(1 - ρ) η_{c}^{0}]}] \\ = & β \exp [(1 - ρ) η_{c}^{0}] . \end{aligned}

Notice how the parameter $ρ$ influences the weight $λ$ when $η_{c} \neq 0$ , in which case the log consumption process displays growth or decay. The condition $λ < 1$ restricts the subjective discount rate, $- \log β,$ relative to the consumption growth rate $η_{c}$ since

(1 - ρ) η_{c}^{0} < - \log β .

When $ρ < 1$ , the subjective discount rate has a positive lower bound in contrast to the case in which $ρ \geq 1$ .

To facilitate computing some useful limits we construct:

(8.20)#

{\tilde{V}}_{t} \overset{def}{=} \frac{{\hat{V}}_{t} - {\hat{V}}_{t}^{0}}{q}

(8.21)#

{\tilde{R}}_{t} \overset{def}{=} \frac{{\hat{R}}_{t} - {\hat{V}}_{t + 1}^{0}}{q},

which we assume remain well defined as $q$ declines to zero, with limits denoted by ${\tilde{V}}_{t}^{0}$ and ${\tilde{R}}_{t}^{0}$ . Importantly,

{\tilde{R}}_{t} = (\frac{1}{1 - γ_{o}}) \log E (\exp [(1 - γ_{o}) {\tilde{V}}_{t + 1}] ∣ A_{t}),

Taking limits as $q$ declines to zero:

{\hat{R}}_{t}^{1} = (\frac{1}{1 - γ_{o}}) \log E (\exp [(1 - γ_{o}) {\hat{V}}_{t + 1}^{1}] ∣ A_{t})

Subtracting ${\hat{C}}_{t}^{1}$ from both sides gives:

(8.22)#

\begin{aligned} {\hat{R}}_{t}^{1} - {\hat{C}}_{t}^{1} = \\ (\frac{1}{1 - γ_{o}}) \log E (\exp [(1 - γ_{o}) [({\hat{V}}_{t + 1}^{1} - {\hat{C}}_{t + 1}^{1}) + ({\hat{C}}_{t + 1}^{1} - {\hat{C}}_{t}^{1})]] ∣ A_{t}) . \end{aligned}

Substituting formula (8.22) into the right side of (8.19) gives the recursion for the first-order continuation value:

(8.23)#

\begin{aligned} {\hat{V}}_{t}^{1} - {\hat{C}}_{t}^{1} = \\ (\frac{λ}{1 - γ_{o}}) \log E (\exp [(1 - γ_{o}) [({\hat{V}}_{t + 1}^{1} - {\hat{C}}_{t + 1}^{1}) + ({\hat{C}}_{t + 1}^{1} - {\hat{C}}_{t}^{1})]] ∣ A_{t}) . \end{aligned}

Remark 8.5

We produce a solution by “guess and verify.” Suppose that

(8.24)#

{\hat{V}}_{t}^{1} - {\hat{C}}_{t}^{1} = {υ_{1}}^{'} X_{t}^{1} + υ_{0} .

It follows from (8.23) that

(8.25)#

υ_{1} = λ ({(ψ_{x^{'}})}^{'} υ_{1} + κ_{x}), υ_{0} = λ (υ_{0} + {υ_{1}}^{'} ψ_{q} + κ_{q} + \frac{(1 - γ_{0})}{2} {| {υ_{1}}^{'} ψ_{w^{'}} + κ_{w^{'}} |}^{2}) .

Deduce the second equation by observing that

\exp [(1 - γ_{0}) ({\hat{V}}_{t + 1}^{1} - {\hat{C}}_{t + 1}^{1}) + (1 - γ_{0}) ({\hat{C}}_{t + 1}^{1} - {\hat{C}}_{t}^{1})]

is distributed as a log normal. The solutions to equations (8.25) are:

υ_{1} = λ {(I - λ ψ_{x^{'}})}^{- 1} κ_{x^{'}}, υ_{0} = \frac{λ}{(1 - λ)} ({υ_{1}}^{'} ψ_{q} + κ_{q}) + \frac{λ (1 - γ_{0})}{2 (1 - λ)} {| {υ_{1}}^{'} ψ_{w^{'}} + κ_{w^{'}} |}^{2} .

The continuation value has two components. The first is:

{υ_{1}}^{'} X_{t}^{1} + \frac{λ}{(1 - λ)} ({υ_{1}}^{'} ψ_{q} + κ_{q}) = E [\sum_{j = 1}^{\infty} λ^{j} ({\hat{C}}_{t + j}^{1} - {\hat{C}}_{t + j - 1}^{1}) ∣ A_{t}]

and the second component is a constant long-run risk adjustment given by:

\frac{λ (1 - γ_{0})}{2 (1 - λ)} {| {υ_{1}}^{'} ψ_{w^{'}} + κ_{w^{'}} |}^{2} .

This second term is the variance of

(8.26)#

E [\sum_{j = 1}^{\infty} λ^{j} ({\hat{C}}_{t + j}^{1} - {\hat{C}}_{t + j - 1}^{1}) ∣ A_{t + 1}] = (1 - λ) E [\sum_{j = 1}^{\infty} λ^{j} ({\hat{C}}_{t + j}^{1} - {\hat{C}}_{t}^{1}) ∣ A_{t + 1}]

conditioned on $A_{t}$ scaled by $\frac{λ (1 - γ_{0})}{2 (1 - λ)}$ .

Remark 8.6

The formula for $υ_{1}$ depends on the parameter $ρ$ . Moreover, $υ_{1}$ has a well-defined limit as $λ$ tends to unity as does the variance of (8.26). This limiting variance:

lim_{λ \to 1} {| {υ_{1}}^{'} ψ_{w^{'}} + κ_{w^{'}} |}^{2} .

converges to the variance of the martingale increment of ${\hat{C}}^{1}$ .

Remark 8.7

Consider the logarithm of the uncertainty-adjusted continuation value approximated to the first order. Note that from (8.24),

{\hat{V}}_{t + 1}^{1} - {\hat{C}}_{t}^{1} = {υ_{1}}^{'} X_{t + 1}^{1} + υ_{0} + κ_{x^{'}} X_{t}^{1} + κ_{w^{'}} W_{t + 1} .

Substitute this expression into formula (8.22) and use the formula for the mean of random variable distributed as a log normal to show that

{\hat{V}}_{t + 1}^{1} - {\hat{R}}_{t}^{1} = ({υ_{1}}^{'} ψ_{w}^{'} + κ_{w^{'}}) W_{t + 1} - (\frac{1 - γ_{o}}{2}) | {υ_{1}}^{'} ψ_{w}^{'} + κ_{w^{'}} |^{2} .

Associated with the first-order approximation, we construct:

N_{t + 1}^{0} \overset{def}{=} \frac{\exp [(1 - γ_{o}) {\tilde{V}}_{t + 1}^{0}]}{E (\exp [(1 - γ_{o}) {\tilde{V}}_{t + 1}^{0}] ∣ A_{t})} = \frac{\exp [(1 - γ_{o}) {\hat{V}}_{t + 1}^{1}]}{E (\exp [(1 - γ_{o}) {\hat{V}}_{t + 1}^{1}] ∣ A_{t})} .

Equation (8.22) is a standard risk-sensitive recursion applied to log-linear dynamics. For instance, see [Tallarini, 2000]’s paper on risk-sensitive business cycles and [Hansen et al., 2008]’s paper on measurement and inference challenges created by the presence of long-term risk.[4] Both of those papers assumed a logarithmic one-period utility function, so that for them $ρ = 1.$ Here we have instead obtained the recursion as a first-order approximation without necessarily assuming log utility. Allowing for $ρ$ to be different than one shows up in both the order zero and order one approximations, as reflected in (8.18) and (8.23), respectively. In accordance to (8.23), for the first-order approximation the parameter $λ = β$ when $ρ = 1$ . But otherwise, it is different. Equation (8.22) also is very similar to a first-order approximation proposed in [Restoy and Weil, 2011]. Like formula (8.22), [Restoy and Weil, 2011] allow for $ρ \neq 1.$ In contrast, our equation has an explicit constant term coming from the uncertainty adjustment, and we have explicit formula for $λ$ that depends on preference parameters and the consumption growth rate.

Remark 8.8

The calculation reported in Remark 8.7 implies that

\begin{aligned} \log N_{t + 1}^{0} = (1 - γ_{o}) ({\hat{V}}_{t + 1}^{1} - {\hat{R}}_{t}^{1}) = & (1 - γ_{o}) ({υ_{1}}^{'} ψ_{w^{'}} + κ_{w^{'}}) W_{t + 1} \\ - \frac{(1 - γ_{o})^{2}}{2} | {υ_{1}}^{'} ψ_{w^{'}} + κ_{w^{'}} |^{2} \end{aligned}

As a consequence, under the change in probability measure induced by $N_{t + 1}^{0},$ $W_{t + 1}$ has a mean given by

\begin{array}{r} μ^{0} \overset{def}{=} (1 - γ_{o}) {({υ_{1}}^{'} ψ_{w^{'}} + κ_{w^{'}})}^{'} \\ = (\frac{1}{ξ_{o}}) {({υ_{1}}^{'} ψ_{w^{'}} + κ_{w^{'}})}^{'} \end{array}

and with the same covariance matrix given by the identity. This is an approximation to robustness adjustment expressed as an altered distribution of the underlying shocks. It depends on $γ_{o} - 1 = \frac{1}{ξ_{o}}$ as well as the state dynamics as reflected by $υ_{1}$ and by the shock exposure vectors $ψ_{w^{'}}$ and $κ_{w^{'}}$ . As we will see, this change of measure plays a role in the higher-order approximation, but it also gives a low-order representation of the implied shadow or market one-period compensation for exposure to uncertainty. It captures the following insight from “long-run risk” models: investor concerns of long-term uncertainty impacts short-term asset valuation. In contrast to the “long-run risk” literature, our analysis opens the door to a different interpretation. Instead of aversion to risk it reflects an aversion to the misspecification of to models or simplified perspectives on macroeconomic dynamics.

As we noted in Remark Remark 8.6, $({υ_{1}}^{'} ψ_{w^{'}} + κ_{w^{'}}) W_{t + 1}$ is approximately the martingale component of the logarithm of consumption when $λ$ is close to one. In VAR example we showed that the variance of this component is challenging to estimate, a point originally made by [Hansen et al., 2008]. This finding is part of the reason that we find it important to step back from rational expectations and limit investors confidence in the models they use for decision making.

8.5.3. Order two#

Differentiating equation (8.3) a second time gives:

{\hat{V}}_{t}^{2} = (1 - λ) {\hat{C}}_{t}^{2} + λ {\hat{R}}_{t}^{2} + (1 - ρ) (1 - λ) λ {({\hat{R}}_{t}^{1} - {\hat{C}}_{t}^{1})}^{2} .

Equivalently,

{\hat{V}}_{t}^{2} - {\hat{C}}_{t}^{2} = λ ({\hat{R}}_{t}^{2} - {\hat{C}}_{t}^{2}) + (1 - ρ) (1 - λ) λ {({\hat{R}}_{t}^{1} - {\hat{C}}_{t}^{1})}^{2} .

Rewrite transformations (8.20) and (8.21) as

q {\tilde{V}}_{t} = {\hat{V}}_{t} - {\hat{V}}_{t}^{0}

q {\tilde{R}}_{t} = {\hat{R}}_{t} - {\hat{V}}_{t + 1}^{0} .

Differentiating twice with respect to $q$ and evaluated at $q = 0$ gives:

{2 \frac{d}{d q} {\tilde{V}}_{t} + q \frac{d^{2}}{d q^{2}} {\tilde{V}}_{t} |}_{q = 0} = 2 {\tilde{V}}_{t}^{1} = {\hat{V}}_{t}^{2} .

{2 \frac{d}{d q} {\tilde{R}}_{t} + q \frac{d^{2}}{d q^{2}} {\tilde{R}}_{t} |}_{q = 0} = 2 {\tilde{R}}_{t}^{1} = {\hat{R}}_{t}^{2} .

Differentiating (8.21) with respect to $q$ gives:

\frac{d {\tilde{R}}_{t}}{d q} = \frac{E (\exp [(1 - γ_{o}) {\tilde{V}}_{t + 1}] \frac{d {\tilde{V}}_{t + 1}}{d q} ∣ A_{t})}{E (\exp [(1 - γ_{o}) {\tilde{V}}_{t + 1}] ∣ A_{t})} .

and thus

(8.27)#

{\hat{R}}_{t}^{2} = 2 {\tilde{R}}_{t}^{1} = 2 E (N_{t + 1}^{0} {\tilde{V}}_{t + 1}^{1} ∣ A_{t}) = E (N_{t + 1}^{0} {\hat{V}}_{t + 1}^{2} ∣ A_{t}),

where subtracting ${\hat{C}}_{t}^{2}$ from ${\hat{R}}_{t}^{2}$ gives:

(8.28)#

{\hat{R}}_{t}^{2} - {\hat{C}}_{t}^{2} = E [N_{t + 1}^{0} ({\hat{V}}_{t + 1}^{2} - {\hat{C}}_{t + 1}^{2}) + ({\hat{C}}_{t + 1}^{2} - {\hat{C}}_{t}^{2}) ∣ A_{t}] .

Substituting this formula into (8.27) gives:

(8.29)#

\begin{aligned} {\hat{V}}_{t}^{2} - {\hat{C}}_{t}^{2} = & λ E (N_{t + 1}^{0} [({\hat{V}}_{t + 1}^{2} - {\hat{C}}_{t + 1}^{2}) + ({\hat{C}}_{t + 1}^{2} - {\hat{C}}_{t}^{2})] ∣ A_{t}) \\ + (1 - ρ) (1 - λ) λ {({\hat{R}}_{t}^{1} - {\hat{C}}_{t}^{1})}^{2} . \end{aligned}

Even if the second-order contribution to the consumption process is zero, there will be nontrivial adjustment to the approximation of $\hat{V} - \hat{C}$ because ${({\hat{R}}^{1} - {\hat{C}}^{1})}^{2}$ is different from zero. This term vanishes when $ρ = 1,$ and its sign will be different depending on whether $ρ$ is bigger or smaller than one.

8.6. Stochastic discount factor approximation#

We approximate $[{\hat{S}}_{t + 1} - {\hat{S}}_{t}]$ in formula (8.8) as

{\hat{S}}_{t + 1} - {\hat{S}}_{t} \approx [{\hat{S}}_{t + 1}^{0} - {\hat{S}}_{t}^{0}] + [{\hat{S}}_{t + 1}^{1} - {\hat{S}}_{t}^{1}] + \frac{1}{2} [{\hat{S}}_{t + 1}^{2} - {\hat{S}}_{t}^{2}]

where

{\hat{S}}_{t + 1}^{0} - {\hat{S}}_{t}^{0} \overset{def}{=} \log β - ρ η_{c}^{0}

{\hat{S}}_{t + 1}^{1} - {\hat{S}}_{t}^{1} \overset{def}{=} - {\hat{C}}_{t + 1}^{1} + {\hat{C}}_{t}^{1} + (ρ - 1) ({\hat{V}}_{t + 1}^{1} - {\hat{C}}_{t + 1}^{1}) - (ρ - 1) ({\hat{R}}_{t}^{1} - {\hat{C}}_{t}^{1})

{\hat{S}}_{t + 1}^{2} - {\hat{S}}_{t}^{2} \overset{def}{=} - {\hat{C}}_{t + 1}^{2} + {\hat{C}}_{t}^{2} + (ρ - 1) ({\hat{V}}_{t + 1}^{2} - {\hat{C}}_{t + 1}^{2}) - (ρ - 1) ({\hat{R}}_{t}^{2} - {\hat{C}}_{t}^{2})

We now consider two different approaches to approximating $N_{t + 1}^{*}$ .

8.6.1. Approach 1#

Write

N_{t + 1}^{*} = \frac{\exp [(1 - γ_{o}) {\tilde{V}}_{t + 1}]}{E (\exp [(1 - γ_{o}) {\tilde{V}}_{t + 1}] ∣ A_{t}]}

= \frac{\exp [(1 - γ_{o}) {\tilde{V}}_{t + 1}]}{\exp [(1 - γ_{o}) {\tilde{R}}_{t}]}

Form the ``first-order’’ approximation:

\log N_{t + 1}^{*} \approx (1 - γ_{o}) [({\tilde{V}}_{t + 1}^{0} - {\tilde{R}}_{t}^{0}) + q ({\tilde{V}}_{t + 1}^{1} - {\tilde{R}}_{t}^{1})]

= (1 - γ_{o}) [({\hat{V}}_{t + 1}^{1} - {\hat{R}}_{t}^{1}) + \frac{q}{2} ({\hat{V}}_{t + 1}^{2} - {\hat{R}}_{t}^{2})]

We combine a first-order approximation of $\log N_{t + 1}^{*}$ with a second-order approximation of ${\hat{S}}_{t + 1} - {\hat{S}}_{t}$ :

\log S_{t + 1} - \log S_{t} \approx (1 - γ_{o}) [({\hat{V}}_{t + 1}^{1} - {\hat{R}}_{t}^{1}) + \frac{1}{2} ({\hat{V}}_{t + 1}^{2} - {\hat{R}}_{t}^{2})]

+ ({\hat{S}}_{t + 1}^{0} - {\hat{S}}_{t}^{0}) + ({\hat{S}}_{t + 1}^{1} - {\hat{S}}_{t}^{1}) + \frac{1}{2} ({\hat{S}}_{t + 1}^{2} - {\hat{S}}_{t}^{2})

which preserves the quadratic approximation of $\log S_{t + 1} - \log S_{t} .$ . Note that If we were to use a second-order approximation of $N_{t + 1}^{*}$ , it would push us outside the class of exponentially quadratic stochastic discount factors.

8.6.2. Approach 2#

Next consider an alternative modification of Approach 1 whereby:

(8.30)#

\begin{aligned} N_{t + 1}^{*} \approx & \frac{\exp [(1 - γ_{o}) ({\tilde{V}}_{t + 1}^{0} + {\tilde{V}}_{t + 1}^{1})]}{E (\exp [(1 - γ_{o}) ({\tilde{V}}_{t + 1}^{0} + {\tilde{V}}_{t + 1}^{1})] ∣ A_{t})} \\ = & \frac{\exp [(1 - γ_{o}) [({\hat{V}}_{t + 1}^{1} - {\hat{R}}_{t}^{1}) + \frac{1}{2} ({\hat{V}}_{t + 1}^{2} - {\hat{R}}_{t}^{2})]]}{E (\exp [(1 - γ_{o}) [({\hat{V}}_{t + 1}^{1} - {\hat{R}}_{t}^{1}) + \frac{1}{2} ({\hat{V}}_{t + 1}^{2} - {\hat{R}}_{t}^{2})]] ∣ A_{t})} \\ \overset{def}{=} & {\tilde{N}}_{t + 1} \end{aligned}

and $\log {\tilde{N}}_{t}$ is used in conjunction with

({\hat{S}}_{t + 1}^{0} - {\hat{S}}_{t}^{0}) + ({\hat{S}}_{t + 1}^{1} - {\hat{S}}_{t}^{1}) + \frac{1}{2} ({\hat{S}}_{t + 1}^{2} - {\hat{S}}_{t}^{2}) .

By design, this approximation of $N_{t + 1}^{*}$ will have conditional expectation equal to one in contrast to the approximation used with Approach 1. With a little bit of algebraic manipulation, it may be shown that this approximation induces a distributional change for $W_{t + 1}$ with a conditional mean that is affine in $X_{t + 1}$ and an altered conditional variance matrix that is constant over time.

To understand better this choice of approximation, consider the family of random variables (indexed by $q$ )

(1 - γ_{o}) ({\tilde{V}}_{t + 1}^{0} + q {\tilde{V}}_{t + 1}^{1}) - \log E (\exp [(1 - γ_{o}) ({\tilde{V}}_{t + 1}^{0} + q {\tilde{V}}_{t + 1}^{1})] ∣ A_{t}) .

The corresponding family of exponentials has conditional expectation one and the $q = 1$ member is the proposed approximation for $N_{t + 1}^{*} .$ Differentiate the family with respect to $q$ :

{\tilde{V}}_{t + 1}^{1} - \frac{E (\exp [(1 - γ_{o}) {\tilde{V}}_{t + 1}^{0}] {\tilde{V}}_{t + 1}^{1} ∣ A_{t})}{E (\exp [(1 - γ_{o}) {\tilde{V}}_{t + 1}^{0}] ∣ A_{t})} = {\tilde{V}}_{t + 1}^{1} - {\tilde{R}}_{t}^{1} .

Thus this family of random variables has the same first-order approximation in $q$ as the one we derived previously for $\log N_{t + 1}^{*}$ .

As a change of probability measure, this approximation will induce state dependence in the conditional mean and will alter the covariance matrix of the shock vector. We find this approach interesting because it links back directly to the outcome of the robustness formulation we described in Section 3.1.

8.7. Solving a planner’s problem with recursive utility#

The [Bansal and Yaron, 2004] example along with many others building connections between the macro economy and asset value take aggregate consumption as pre-specified. As we open the door to a richer collection of macroeconomic models, it becomes important to entertain more endogeneity, including investment and other variables familiar to macroeconomics.

Write a triangular system with stochastic growth as:

(8.31)#

\begin{aligned} X_{t + 1} (q) = & ψ^{x} [D_{t} (q), X_{t} (q), q W_{t + 1}, q] \\ \log G_{t + 1} (q) - \log G_{t} (q) = & ψ^{g} [D_{t} (q), X_{t} (q,), q W_{t + 1}, q], \end{aligned}

where $D_{t}$ is a date $t$ decision vector for the planner. Define ${\hat{G}}_{t} = \log G_{t} .$ In addition, we impose

(8.32)#

\begin{aligned} 0 & = ϕ [D_{t} (q), X_{t} (q)] \\ {\hat{C}}_{t} (q) & = κ [D_{t} (q), X_{t} (q)] + {\hat{G}}_{t} (q) . \end{aligned}

where the first equation is a vector of static constraints and the second constructs the measure of consumption that enters preferences.
We extend the approximations by using a co-state formulation. There are two essentially equivalent interpretations of these co-states. One is they function as a set of Lagrange multipliers on the state evolution equations. The other is that they are partial derivatives of value functions. The co-state equations are forward-looking, linking next period’s co-state vector to this period’s co-state vector. Given the recursive utility structure, we must include the implied value functions in the computations as they enter the relations of interest.

The first-order conditions for $D$ are:

(8.33)#

\begin{aligned} (1 - β) \exp [(1 - ρ) ({\hat{C}}_{t} - {\hat{V}}_{t})] κ_{d} (D_{t}, X_{t}) + ϕ_{d^{'}} (D_{t}, X_{t})^{'} M S_{t} \\ + β \exp [(1 - ρ) ({\hat{R}}_{t} - {\hat{V}}_{t})] E (N_{t + 1}^{*} ψ_{d^{'}}^{x} (D_{t}, X_{t}, q W_{t + 1}, q)^{'} M X_{t + 1} ∣ A_{t}) \\ + β \exp [(1 - ρ) ({\hat{R}}_{t} - {\hat{V}}_{t})] E (N_{t + 1}^{*} ψ_{d^{'}}^{g} (D_{t}, X_{t}, q W_{t + 1}, q)^{'} M G_{t + 1} ∣ A_{t}) \\ = 0. \end{aligned}

where $M X_{t + 1}$ and $M G_{t + 1}$ are the co-states, or the implicit multipliers one for each the state evolutions and $M S_{t}$ is a multiplier on the first static constraint in (8.32). Recall that $N_{t + 1}^{*} = \exp [(1 - γ) ({\hat{V}}_{t + 1} - {\hat{R}}_{t})]$ used for making an uncertainty adjustment in valuation.

In addition, we solve a forward-looking co-state equation given by

(8.34)#

\begin{aligned} [\begin{array}{c} M X_{t} \\ M G_{t} \end{array}] = (1 - β) \exp [(1 - ρ) ({\hat{C}}_{t} - {\hat{V}}_{t})] [\begin{array}{c} κ_{x} (D_{t}, X_{t}) \\ 1 \end{array}] + [\begin{array}{c} ϕ_{x} (D_{t}, X_{t})^{'} M S_{t} \\ 0 \end{array}] \\ + β \exp [(1 - ρ) ({\hat{R}}_{t} - {\hat{V}}_{t})] \times \\ E (N_{t + 1}^{*} [\begin{array}{c} ψ_{x^{'}}^{x} (D_{t}, X_{t}, q W_{t + 1}, q)^{'} & ψ_{x^{'}}^{g} (D_{t}, X_{t}, q W_{t + 1}, q)^{'} \\ 0 & 1 \end{array}] [\begin{array}{c} M X_{t + 1} \\ M G_{t + 1} \end{array}] ∣ A_{t}) . \end{aligned}

The approximation formulas, (8.22), (8.23), (8.28), (8.29), that we deduced previously for ${\hat{V}}_{t} - {\hat{C}}_{t}$ and ${\hat{R}}_{t} - {\hat{C}}_{t}$ have and immediate counterparts for ${\hat{V}}_{t} - {\hat{G}}_{t}$ and ${\hat{R}}_{t} - {\hat{G}}_{t}$ , from which can build an approximation of ${\hat{R}}_{t} - {\hat{V}}_{t} .$

From recursive utility updating equation:

\begin{array}{r} \exp [(1 - ρ) {\hat{V}}_{t}] = (1 - β) \exp [(1 - ρ) {\hat{C}}_{t}] + β \exp [(1 - ρ) {\hat{R}}_{t}] . \end{array}

Dividing by both sides of the equation by $\exp [(1 - ρ) {\hat{V}}_{t}]$ gives

\begin{array}{r} 1 = (1 - β) \exp [(1 - ρ) ({\hat{C}}_{t} - {\hat{V}}_{t})] + β \exp [(1 - ρ) ({\hat{R}}_{t} - {\hat{V}}_{t})] . \end{array}

From the relation, we see that $M G_{t} = 1$ for $t \geq 0$ satisfies the second block of co-state equation (8.34). In fact, this is the solution of interest.

8.7.1. An example economy with long-run uncertainty#

Consider an AK model with recursive utility and adjustment costs.

The exogenous state dynamic capture both long-run uncertainty in the mean growth rate and the overall volatility in the economy.

(8.35)#

\begin{aligned} Z_{1, t + ϵ} - Z_{1, t} = & - ν_{1} Z_{1, t} + \exp (\frac{1}{2} Z_{2, t}) σ_{1} W_{t + 1} \\ Z_{2, t + 1} - Z_{2, t} = & - ν_{2} [1 - μ_{2} \exp (- Z_{2, t})] \\ - \frac{1}{2} | σ_{2} |^{2} \exp (- Z_{2, t}) + \exp (- \frac{1}{2} Z_{2, t}) σ_{2} W_{t + 1} . \end{aligned}

The state variable, $Z_{2, t}$ is included to capture stochastic volatility. The discrete-time dynamics for ${\exp (Z_{2, t})}$ approximate a continuous-time version of what is called a square root process due to Feller. Let

X_{t} = [\begin{matrix} Z_{1, t} \\ Z_{2, t} \end{matrix}] .

With these exogenous dynamics, we obtain the following zero and first-order approximations:

Z_{1, t}^{0} = 0 \exp (Z_{2, t}^{0}) = μ_{2},

and

\begin{aligned} Z_{1, t + 1}^{1} - Z_{1, t}^{1} = & - ν_{1} Z_{1, t}^{1} + \sqrt{{\frac{1}{μ}}_{2}} σ_{1} W_{t + 1} \\ Z_{2, t + 1}^{1} - Z_{2, t}^{1} = & - ν_{2} Z_{2, t}^{1} + \sqrt{\frac{1}{μ_{2}}} σ_{2} W_{t + 1} . \end{aligned}

We impose the resource constraint:

C_{t} + I_{t} = α K_{t} .

The endogenous state dynamics are given by:

\begin{aligned} {\hat{K}}_{t + 1} - {\hat{K}}_{t} = & [\frac{1}{ζ} \log (1 + ζ \frac{I_{t}}{K_{t}}) + ν_{k} Z_{1, t} - ι_{k}] \\ - \frac{1}{2} | σ_{k} |^{2} \exp (Z_{2, t}) + \exp (\frac{1}{2} Z_{2, t}) σ_{k} W_{t + ϵ} \end{aligned}

where ${\hat{K}}_{t} = \log K_{t} = {\hat{G}}_{t}$ . The planner choice variable $D_{t} = (\frac{C_{t}}{K_{t}}, \frac{I_{t}}{K_{t}})$ . Rewrite the current-period resource constraints as:

(8.36)#

\begin{aligned} 0 & = α - D_{1, t} - D_{2, t} \\ {\hat{C}}_{t} - {\hat{G}}_{t} & = \log D_{1, t} . \end{aligned}

Express the first-order conditions for the consumption-capital and investment-capital ratios as:

(8.37)#

\begin{aligned} 0 & = - M S_{t} + \frac{1 - β}{D_{1, t}} \exp [(1 - ρ) ({\hat{C}}_{t} - {\hat{G}}_{t})] \\ 0 & = - M S_{t} + \exp [(1 - ρ) ({\hat{R}}_{t} - {\hat{G}}_{t})] (\frac{β}{1 + ζ D_{2, t}}) \end{aligned}

It is convenient to rewrite the first-order conditions for the consumption-capital ratio as:

\log M S_{t} + \log D_{1, t} = \log (1 - β) + (1 - ρ) ({\hat{C}}_{t} - {\hat{G}}_{t}),

which in turn implies that

\log M S_{t} = \log (1 - β) - ρ ({\hat{C}}_{t} - {\hat{G}}_{t}) .

More generally, we will seek to approximate $\log M S_{t}$ , as we expect the multiplier $M S_{t}$ to be positive.
Notice that these first-order conditions do not depend on the co-state process ${M X_{t} : t \geq 0}$ . We may solve this planner’s problem using (8.37) and up-dating the continuation value processes and its uncertainty-adjusted counterpart until convergence. When $ρ = 1$ , ${\hat{R}}_{t} - {\hat{G}}_{t}$ drops out of the first-order conditions and both components of $D_{t}$ are constant since

M S_{t} = \frac{1}{D_{1, t}} = \frac{1}{α - D_{2, t}} .

When $ρ \neq 1$ , $D_{t}$ depends on the exogenous state $X_{t}$ .

8.7.2. First-order approximation when $ρ = 1$ #

We obtain the following zero and first-order approximations for the exogenous dynamics:

Z_{1, t}^{0} = 0 \exp (Z_{2, t}^{0}) = μ_{2},

and

\begin{aligned} Z_{1, t + 1}^{1} - Z_{1, t}^{1} = & - ν_{1} Z_{1, t}^{1} + \sqrt{μ_{2}} σ_{1} W_{t + 1} \\ Z_{2, t + 1}^{1} - Z_{2, t}^{1} = & - ν_{2} Z_{2, t}^{1} + \sqrt{\frac{1}{μ_{2}}} σ_{2} W_{t + 1} . \end{aligned}

The first-order conditions for $D_{2, t}$ imply that:

\frac{1 - β}{α - D} + \frac{β}{1 + ζ D} = 0.

Solving for $D$ gives:

D_{2, t}^{*} = \frac{(β - 1) + β α}{β + (1 - β) ζ},

which is independent of the state, as should be expected since $ρ = 1.$ From the capital evolution it follows from the order zero approximation is

{\hat{K}}_{t + 1}^{0} - {\hat{K}}_{t}^{0} = [\frac{1}{ζ} \log (1 + ζ D_{2, t}^{*}) - ι_{k}] .

The order one approximation is then:

{\hat{K}}_{t + 1}^{1} - {\hat{K}}_{t}^{1} = ν_{k} Z_{1, t}^{1} + \sqrt{μ_{2}} σ_{k} W_{t + 1} .

Stochastic volatility, as in the [Bansal and Yaron, 2004] model of consumption dynamics, will be present in the second-order approximation.

8.7.3. Second-order approximation when $ρ = 1$ #

We next consider the second-order approximations. The second-oder approximation for ${Z_{2, t}}$ does not contribute to the planner’s solution or to the implied shadow prices and thus we drop it from the analysis. For the remaining two state variables, we find that

\begin{aligned} Z_{1, t + 1}^{2} - Z_{1, t}^{2} = & - ν_{1} Z_{1, t}^{2} - μ_{2} | σ_{1} |^{2} + \sqrt{μ_{2}} Z_{2, t}^{1} σ_{1} W_{t + 1} \\ {\hat{K}}_{t + 1}^{2} - {\hat{K}}_{t}^{2} = & - μ_{2} | σ_{k} |^{2} + \sqrt{μ_{2}} Z_{2, t}^{1} σ_{k} W_{t + 1} \end{aligned}

where we previously noted that ${Z_{2, t}^{1}}$ evolves as first-order autoregression.

The combined approximation for $q = 1$ uses:

\begin{aligned} Z_{1, t + 1} - Z_{1, t} \approx & - ν_{1} Z_{1, t} \\ - \frac{μ_{2} | σ_{1} |^{2}}{2} + \sqrt{μ_{2}} σ_{1} W_{t + 1} + \frac{μ_{2} Z_{2, t}}{2} σ_{1} W_{t + 1} \\ Z_{2, t + ϵ} - Z_{2, t} \approx & - ν_{2} Z_{2, t} + \sqrt{\frac{1}{μ_{2}}} σ_{2} W_{t + 1} \\ {\hat{K}}_{t + 1} - {\hat{K}}_{t} \approx & [\frac{1}{ζ} \log (1 + ζ D^{*}) - ι_{k}] + ν_{k} Z_{1, t} \\ - \frac{μ_{2} | σ_{k} |^{2}}{2} + \sqrt{μ_{2}} σ_{k} W_{t + 1} + \frac{μ_{2} Z_{2, t}}{2} σ_{k} W_{t + 1} . \end{aligned}

The approximate dynamics for the exogenous states remains the same for $ρ \neq 1,$ but the solution for $D^{*}$ becomes state dependent and the approximate dynamic evolution for capital is altered.

8.7.4. Shock elasticities#

We use the shock elasticities to explore pricing implications of this recursive utility specification. In what follows, we use exponential/linear/quadratic implementation by [Borovička and Hansen, 2014] and by [Borovička and Hansen, 2016] with the parameter configuration given in Table 1 of [Hansen et al., 2024]. This latter reference combines inputs from other sources including [Schorfheide et al., 2018] and [Hansen and Sargent, 2021].

Fig. 8.1 gives the shock exposure elasticities for consumption to each of the three shocks. They can be interpreted as nonlinear local impulse responses for consumption (in levels not logarithms). The elasticities for the growth rate shock and the stochastic volatility shock start small and increase over the time horizon as dictated by the persistence of the two exogenous state variable processes. The elasticities for the direct shock to capital are flat over the horizon as to be expected since the shock directly impacts log consumption in a manner that is permanent. Notice that while elasticities for the volatility shock are different from zero, their contribution is much smaller than the other shocks. Stochastic volatility does induce state dependence for the other elasticities as reflected by the quantiles.[5]

../_images/expo_010.png — Fig. 8.1 Exposure elasticities for three shocks. $ρ = 1, γ = 8,$ $β = .99 .$ #

Fig. 8.2 gives the shock price elasticities for growth-rate shock when $ρ = 1$ and $γ = 8$ . Stochastic volatility induces state dependence in these plots as reflected by the quantiles. The recursive utility preferences are forward-looking as reflected by the continuation-value contribution to the one-period increment to the stochastic discount factor process (see (8.8)). This forward-looking contribution is reflected in shock price elasticities that are relatively flat for the growth rate shock as is evident from the second column. The shock price elasticity that is induced by uncertainty alone is flat since because the associated multiplicative functional is a martingale.[6]

../_images/price_010.png — Fig. 8.2 Price elasticities for three shocks. $ρ = 1, γ = 8,$ $β = .99 .$ #

Fig. 8.3 and Fig. 8.4 provide the analogous plots for $ρ = 2 / 3, 3 / 2.$ Fig. 8.5 sets $ρ = γ = 8$ which corresponds to preferences that are time separable. The forward-looking component to the stochastic discount factor is shut down as is evident from formula (8.8). Now the shock price elasticities and shock exposure elasticities show a very similar trajectory except that the shock price elasticities are about eight times larger.

../_images/price_006.png — Fig. 8.3 Price elasticities for three shocks. $ρ = 2 / 3, γ = 8,$ $β = .99 .$ #

../_images/price_015.png — Fig. 8.4 Price elasticities for three shocks. $ρ = 3 / 2, γ = 8,$ $β = .99 .$ #

../_images/price_080.png — Fig. 8.5 Price elasticities for three shocks. $ρ = 8, γ = 8,$ $β = .99 .$ #

8.7.5. An alternative model of intertemporal substitution/complementarity.#

We now extend the preference specification of the consumers in the AK model of Section 8.7.1 to explore implications of time nonseparability in preferences. We are motivated to do so by rather substantial previous literetaure. An important earlier contributor is [Ryder and Heal, 1973], who solve a social planner’s problem with stochastic growth. [Sundaresan, 1989], [Constantinides, 1990], [Heaton, 1995], and [Hansen et al., 1999] consider asset pricing implications with internal habit persistence. These papers essentially explore decentralizations of the planner’s problem. The latter paper considers simultaneously a recursive robustness specification as we will do here, although for a different model specification.

To accomplish this, we introduce an additional state variable, which we will call the habit stock, $H_{t}$ . This stock evolves as:

H_{t + 1} = \exp (- ν_{h}) H_{t} + [1 - \exp (- ν_{h})] I_{h, t} .

where $ν_{h} > 0$ is a depreciation rate. We think of this new investment = $I_{h, t}$ as measured consumption.
Notice that the sum of the coefficients on the right side of the evolution equation sum to one. This allows us to interpret $H_{t + 1}$ as a geometric average of current and past values of measured consumption.

We write evolution equation in logarithms as

\log H_{t + 1} = \log [\exp (- ν_{h} + \log H_{t}) + [1 - \exp (- ν_{h})]] I_{h, t}] .

For numerical purposes, we transform this equation to be:

\begin{aligned} (\log H_{t + 1} - \log K_{t + 1}) \\ = \log K_{t} - \log K_{t + 1} + \log [\exp (- ν_{h} + \log H_{t} - \log K_{t}) + [1 - \exp (- ν_{h})] (\frac{I_{h, t}}{K_{t}})] \\ = - [\frac{1}{ζ} \log (1 + ζ \frac{I_{t}}{K_{t}}) - ν_{k} Z_{1, t} + ι_{k}] + \frac{1}{2} | σ_{k} |^{2} \exp (Z_{2, t}) - \exp (\frac{1}{2} Z_{2, t}) σ_{k} W_{t + ϵ} \\ + \log [\exp (- ν_{h} + \log H_{t} - \log K_{t}) + [1 - \exp (- ν_{h})] (\frac{I_{h, t}}{K_{t}})] \end{aligned}

where we treat $\log H_{t} - \log K_{t}$ as one of the components of $X_{t}$ and we substitute the evolution of $\log K_{t + 1}$ to depict the evolution of this component. As previously, ${\hat{G}}_{t} = \log K_{t}$ .

We impose the output constraint as:

I_{h, t} + I_{k, t} = α K_{t}

where $I_{k, t}$ is the investment in the capital stock. Analogous to previous formulation, we transform this equation to be:

\frac{I_{h, t}}{K_{t}} + \frac{I_{k, t}}{K_{t}} = α

where $\frac{I_{h, t}}{K_{t}}$ and $\frac{I_{k, t}}{K_{t}}$ are the two components of $D_{t}$ .

What enters the utility function each date inside the recursive utility preference specification is the CES aggregate:

(8.38)#

C_{t} = {[(1 - λ) {(I_{h, t})}^{1 - τ} + λ {(H_{t})}^{1 - τ}]}^{\frac{1}{1 - τ}}

for $τ \geq 0$ . This specification captures a form of intertemporal complementarity and intertemporal substitution in preferences. In static demand theory, $τ = 0$ implies perfect substitutes, and when $τ > 1$ , the preferences display a form of complementarity as reflected in the cross price effects on demand. In the limiting case when $τ$ approaches $+ \infty$ , preferences treat the two inputs as perfect complements and consumed in fixed proportions. There is a different notion of intertemporal complementarity as originally defined by [Ryder and Heal, 1973]. Under this notion, increasing $H_{t}$ , keeping $I_{h, t}$ fixed should decrease $C_{t}$ . (They impose other restrictions.) This form of intertermporal complementarity can instead be captured by letting $λ < 0$ . Effectively, $H_{t}$ is a “bad’’ not a “good,’’ in the CES aggregator. Setting $0 < λ \leq 1$ captures a form of consumption durability.

For computational purposes, we divide both sides of (8.38)`by $K_{t}$ and take logarithms:

\log C_{t} - \log K_{t} = \frac{1}{1 - τ} \log [(1 - λ) {(\frac{I_{h, t}}{K_{t}})}^{1 - τ} + λ \exp [(1 - τ) (\log H_{t} - \log K_{t})]]

Consistent with the habit-persistence literature, we allow the parameter $λ$ to be negative. When $τ$ is one, this becomes a Cobb-Douglas specification:

\log C_{t} - \log K_{t} = (1 - λ) (\log I_{h, t} - \log K_{t}) + λ (\log H_{t} - \log K_{t})

Remark 8.9

The asset pricing literature often features the computationally more challenging case in which $τ = 0.$ and $λ < 0$ . The local approximation methods we describe here could give particularly poor approximations for this parameter configuration.

For simplicity, consider the case in which $ρ = τ = 1$ . Recall that $ρ$ also impacts intertemporal substitution in preferences so by setting it to one we feature the novel contribution coming from this dynamic specification of preferences. The first-order conditions for $D_{t}$ are:

\begin{aligned} (1 - β) [\begin{array}{c} \frac{1 - λ}{D_{1, t}} \\ 0 \end{array}] - [\begin{array}{c} 1 \\ 1 \end{array}] M S_{t} + E [N_{t + 1}^{*} [\begin{array}{c} \frac{β K_{t} [1 - \exp (- ν_{h})]}{\exp (X_{t + 1}) K_{t + 1}} & 0 & 0 \\ \frac{1}{1 + ζ D_{2, 1}} & 0 & 0 \end{array}] M X_{t + 1} ∣ A_{t}] \\ + β [\begin{array}{c} 0 \\ \frac{1}{1 + ζ D_{2, 1}} \end{array}] E (N_{t + 1}^{*} M G_{t + 1} ∣ A_{t}) = 0. \end{aligned}

Consider now the special case in which $λ = 0.$ In the case the co-state $M X_{t + 1}$ should be zero. From the first-order conditions for $D_{1, t}$ , we find that

M S_{t} = (1 - β) (\frac{1}{D_{1, t}})

More generally,

M S_{t} = (1 - β) (\frac{1 - λ}{D_{1, t}}) + E (N_{t + 1}^{*} [\begin{matrix} \frac{β K_{t} [1 - \exp (- ν_{h})]}{\exp (X_{t + 1}) K_{t + 1}} & 0 & 0 \end{matrix}] M X_{t + 1} ∣ A_{t}) .

This depicts the marginal value of $M S_{t}$ in terms of both a current marginal utility contribution and a forward-looking piece coming from the planner internalizing the intertemporal contribution to preferences. We choose to use $I_{h, t}$ as the date $t$ numeraire instead of $C_{t}$ , as the former is measured consumption. In light of this choice, the logarithm of the one-period stochastic discount factor is

{\hat{S}}_{t + 1} - {\hat{S}}_{t} = \log β + \log M S_{t + 1} - \log M S_{t} + {\hat{K}}_{t + 1} - {\hat{K}}_{t} .

In our calculations, we approximate $\log M S_{t}$ since it enters this formula for the logarithmic stochastic discount factor construction.

Figure Fig. 8.6 explores the median exposure elasticities for the investment-capital ratio as a function of the parameter $λ$ . Recall that these elasticities can be viewed as local impulse responses. The responses are different for the two shocks. For the growth-rate shock, with habit persistence ( $λ = 1$ ), there is an immediate decline followed by an eventual increase. In contrast, for durability ( $λ = .5)$ just the opposite happens. Given the restriction on output, the consumption/capital ratios have an offsetting response. For the direct shock to the technology, with habit persistence, the investment-capital ratio increases and monotonically declines to zero as a function of the investment horizon. Again, under consumption durability, the response is essentially the mirror image of this. Figure Fig. 8.7 provides the $.1$ and $.9$ quantiles as induced by stochastic volatility for the two shocks. Overall, we see that investment/capital and consumption/capital responses are noticeably sensitive to the choice of $λ$ .

../_images/lambda_varying.png — Fig. 8.6 Investment/capital exposure elasticities for the growth and capital shocks with different values of $λ .$ #

../_images/quantiles.png — Fig. 8.7 Investment/capital exposure elasticities for the growth shock with different values of $λ$ and elasticity quantiles.#

We next consider the shock price elasticities for the growth shock. Fig. 8.8 reports the elasticities for three values of $γ$ and $λ = 0$ . Recall that $γ = 1$ is equivalent to investors having no concerns about robustness. Since $λ = 0,$ this figure repeats earlier results and is included here for comparison. As we noted in previous discussion, these elasticities increase in the value of $γ$ ( increase in the concern about robustness) and become very flat for larger values of $γ$ as the martingale component to the stochastic discount factor process comes to dominate. Fig. 8.9 considers the growth shock price elasticities when $λ = .5$ . These elasticities are similar to the ones reported for $λ = 0$ in Fig. 8.8, although there is small drop in the magnitudes relative the $λ = 0$ specifcation. Fig. 8.10 explores the same elasticities when $λ = - 1$ . We now see a modest increase in the shock price elasticities relative to the two previous figures.

../_images/lambda0.png — Fig. 8.8 Price elasticities for the growth shock for $λ = 0$ and different values of $γ$ .#

../_images/lambda09.png — Fig. 8.9 Price elasticities for the growth shock for $λ = 0.5$ and different values of $γ$ .#

../_images/-lambda10.png — Fig. 8.10 Price elasticities for the growth shock for $λ = - 10$ and different values of $γ$ .#

Overall, we see that investment/capital and consumption/capital impulse responses are noticeably sensitive to the choice of $λ$ , which should be expected given the impact of this parameter on intertemporal substitution in preferences. In contrast, $λ$ has a small impact on the shock price elasticities. These elasticities are, however, notably sensitive to the choice of $γ,$ which is consistent with recursive utility construction of uncertainty aversion.

[Pollak, 1970] introduced a version of external habit persistence in consumer demand function. The stock $H_{t}$ enters preferences as a societal or external input and not recognizable as an individual input in the current time period. Thus there is an externality induced by this specification of preferences. [Abel, 1990] and [Campbell and Cochrane, 1999] explore implications along with many subsequent papers explore asset pricing with external habit persistence. Much of the latter literature abstracts from production, but there are exceptions. See for instance, [Lettau and Uhlig, 2000]. The equilibrium solution with external habit persistence deviates from the planner problem in that the co-state for $H_{t}$ is now zero and the co-state equation is omitted from the equations to be solved and approximated.

8.7.6. External habits model#

../_images/external_expo_1.png — Fig. 8.11 Investment/capital exposure elasticities for the growth and capital shocks with different values of $λ .$ #

Please note that the axes in the left panel on the graph below are scaled 0.01 times relative to the right panel.

../_images/external_expo_2.png — Fig. 8.12 Investment/capital exposure elasticities for the growth shock with different values of $λ$ and elasticity quantiles.#

8.8. Solving models#

In this section, we briefly describe one way to extend the approach that builds directly on previous second-order approaches of [Kim et al., 2008], [Schmitt-Grohé and Uribe, 2004], and [Lombardo and Uhlig, 2018]. While such methods should not be viewed as being generically applicable to nonlinear stochastic equilibrium models, we find them useful pedagogically and often as at least initial steps to understanding models that are arguably “smooth.” See [Pohl et al., 2018] for a careful study of nonlinearity in asset pricing models with recursive utility.[7]

We implement these methods for second-order approximation using the following steps.

Solve for $q = 0$ deterministic model.
Take as given first and second-order approximate solutions for ${\hat{C}}_{t} - {\hat{G}}_{t}$ and ${\hat{G}}_{t + ϵ} - {\hat{G}}_{t} .$ . Solve for the approximate solutions for ${\hat{V}}_{t} - {\hat{G}}_{t},$ ${\hat{V}}_{t + ϵ} - {\hat{R}}_{t}$ and $N_{t + ϵ} .$
Compute the first-order expansion and solve the resulting equations following the previous literature for $D_{t},$ ${\hat{C}}_{t} - {\hat{G}}_{t},$ and ${\hat{G}}_{t + ϵ} - {\hat{G}}_{t} .$ When constructing these equations, use expectations computed using the probabilities induced by $N_{t + ϵ}^{0}$ . Substitute the first-order approximation for ${\hat{R}}_{t} - {\hat{G}}_{t} .$
Compute the second-order expansion and solve the resulting equations following the previous literature. Again use the expectations induced by $N_{t + ϵ}^{0}$ . In addition, make another recursive utility adjustment expressed in terms the approximations of ${\hat{R}}_{t} - {\hat{G}}_{t} .$
Return to step 2, and repeat until convergence.

Initialize this algorithm by solving the $γ_{o} = 1$ and $ρ = 1,$ which can be solved without iteration.

See the Appendices that follow for more details and formulas to use in the solution method.

As a second approach we iterate over an $N_{t + ϵ}^{*}$ given by formula (8.30) approximation restricted to induce an alternative probability distribution Call the approximation ${\tilde{N}}_{t + ϵ}$ with an induced distribution for $W_{t + ϵ}$ that is normal with conditional mean ${\tilde{μ}}_{t}$ and covariance matrix $\tilde{Σ}$ .

While we discussed the approximation for resource allocation problems with recursive utility, there is a direct extension of this approach to solve a general class stochastic equilibrium models by stacking a system of expectational-type equations expressed in part using the recursive utility stochastic discount factor that we derived. For resource allocation problems, we expressed the first-order conditions for the planner in utility units, which simplified some formulas. Equilibrium models not derived from a planner’s problems typically use stochastic discount factors expressed in consumption units when representing investment choices. The approximation methods described in this chapter have a direct extension to such models.

8.9. Appendix A: Solving the planner’s problem#

Write the system of interest, including the state equations (8.31), the consumption equation and static constraint (8.32), the first-order conditions (8.33), and the co-state evolution (8.34) as:

\begin{aligned} X_{t + 1} (q) = & ψ^{x} [D_{t} (q), X_{t} (q), q W_{t + 1}, q], \\ \log G_{t + 1} (q) - \log G_{t} (q) = & ψ^{g} [D_{t} (q), X_{t} (q,), q W_{t + 1}, q], \\ {\hat{C}}_{t} (q) = & κ [D_{t} (q), X_{t} (q)] + {\hat{G}}_{t} (q), \\ 0 = & ϕ [D_{t} (q), X_{t} (q)], \\ Q_{t} E (N_{t + 1} H_{t + 1} ∣ A_{t}) + P_{t} L_{t} - M_{t} = & 0 \end{aligned}

where

\begin{aligned} Q_{t} \overset{def}{=} & β \exp [(1 - ρ) ({\hat{R}}_{t} - {\hat{V}}_{t})] \\ P_{t} \overset{def}{=} & (1 - β) \exp [(ρ - 1) ({\hat{V}}_{t} - {\hat{G}}_{t})] \\ H_{t + 1} \overset{def}{=} & [\begin{array}{c} ψ_{d^{'}}^{x} (D_{t}, X_{t}, W_{t + 1})^{'} & ψ_{d^{'}}^{g} (D_{t}, X_{t}, W_{t + 1}) \\ ψ_{x^{'}}^{x} (D_{t}, X_{t}, W_{t + 1})^{'} & ψ_{x^{'}}^{g} (D_{t}, X_{t}, W_{t + 1})^{'} \\ ψ_{g^{'}}^{x} (D_{t}, X_{t}, W_{t + 1})^{'} & 1 \end{array}] [\begin{array}{c} M X_{t + 1} \\ M G_{t + 1} \end{array}] \\ L_{t} \overset{def}{=} & [\begin{array}{c} \exp [(1 - ρ) ({\hat{C}}_{t} - {\hat{G}}_{t})] [\begin{matrix} κ_{d} (D_{t}, X_{t}) \\ κ_{x} (D_{t}, X_{t}) \\ 1 \end{matrix}] \end{array}] \\ M_{t} \overset{def}{=} & [\begin{array}{c} 0 \\ M X_{t} \\ M G_{t} \end{array}] + [\begin{array}{c} ϕ_{d^{'}} (D_{t}, X_{t})^{'} \\ ϕ_{x^{'}} (D_{t}, X_{t})^{'} \\ 0 \end{array}] M S_{t} . \end{aligned}

Here, for computational purposes, we use that ${\hat{C}}_{t} - {\hat{V}}_{t} = ({\hat{C}}_{t} - {\hat{G}}_{t}) - ({\hat{V}}_{t} - {\hat{G}}_{t})$ . We solve for ${\hat{C}}_{t} - {\hat{G}}_{t}, D_{t}, M X_{t}$ as a function of $X_{t}$ . The objects: $\hat{C} - \hat{G}, D$ and $M X$ are sometimes referred to as jump variables since we not impose initial conditions for these variables as part of a solution.

Our solution will entail an iteration. We will impose a specification for $Q^{1}, Q^{2},$ and $N$ and find an approximate solution for the dynamical system. Then given this solution, we will compute a new implied solution for $Q^{1}, Q^{2},$ and $N$ . We then iterate this until we achieve numerical convergence. We use second-order approximations for both steps.

8.9.1. Some steady state calculations#

Observe from the recursive utility updating that:

\begin{aligned} {\hat{V}}_{t}^{0} - {\hat{G}}_{t}^{0} & = \frac{1}{1 - ρ} \log [(1 - β) \exp [(1 - ρ) ({\hat{C}}_{t}^{0} - {\hat{G}}_{t}^{0})] + β \exp [(1 - ρ) ({\hat{R}}_{t}^{0} - {\hat{G}}_{t}^{0})]] \\ {\hat{R}}_{t}^{0} - {\hat{G}}_{t}^{0} & = {\hat{V}}_{t + 1}^{0} - {\hat{G}}_{t + 1}^{0} + {\hat{G}}_{t + 1}^{0} - {\hat{G}}_{t}^{0} . \end{aligned}

In the steady state we view this as two equations in three variables: ${\hat{V}}_{t}^{0} - {\hat{G}}_{t}^{0}, {\hat{R}}_{t}^{0} - {\hat{G}}_{t}^{0}$ and ${\hat{G}}_{t + 1}^{0} - {\hat{G}}_{t}^{0}$ , each of which we assume is time invariant.
We construct

\begin{aligned} Q_{t}^{0} = & β \exp [(1 - ρ) ({\hat{R}}_{t}^{0} - {\hat{V}}_{t}^{0})] \\ = & β \exp [(1 - ρ) ({\hat{G}}_{t + 1}^{0} - {\hat{G}}_{t}^{0})], \end{aligned}

and

P_{t}^{0} = (1 - β) \exp [(ρ - 1) ({\hat{V}}_{t}^{0} - {\hat{G}}_{t}^{0})]

in the steady state, and use them to construct the remaining steady-state equations:

Q_{t}^{0} H_{t + 1}^{0} + P_{t}^{0} - M_{t}^{0} = 0.

In this equation, $H_{t + 1}^{0}, L_{t}^{0},$ and $M_{t}^{0}$ are constructed from the formulas for $H_{t + 1}, L_{t}$ M_t$, defined previously, by setting the shock vector to zero and treating the relevant variables as time invariant.

8.9.2. $Q$ and $P$ derivatives#

For the order one, write

Q_{t}^{1} \overset{def}{=} (1 - ρ) Q_{t}^{0} ({\hat{R}}_{t}^{1} - {\hat{V}}_{t}^{1})

To compute this contribution, we equation (8.19) to write

{\hat{R}}_{t}^{1} - {\hat{C}}_{t}^{1} = \frac{1}{λ} ({\hat{V}}_{t}^{1} - {\hat{C}}_{t}^{1}) .

We then construct

\begin{aligned} {\hat{R}}_{t}^{1} - {\hat{V}}_{t}^{1} = & ({\hat{R}}_{t}^{1} - {\hat{C}}_{t}) + ({\hat{C}}_{t}^{1} - {\hat{V}}_{t}) \\ = & (\frac{1 - λ}{λ}) [({\hat{V}}_{t}^{1} - {\hat{G}}_{t}) - ({\hat{C}}_{t} - {\hat{G}}_{t})] . \end{aligned}

To compute ${\hat{V}}_{t} - {\hat{G}}_{t}$ , we rewrite (8.23) as:

\begin{aligned} {\hat{V}}_{t}^{1} - {\hat{G}}_{t}^{1} = (1 - λ) ({\hat{C}}_{t}^{1} - {\hat{G}}_{t}^{1}) \\ + (\frac{λ}{1 - γ_{o}}) \log E (\exp [(1 - γ_{o}) [({\hat{V}}_{t + 1}^{1} - {\hat{G}}_{t + 1}^{1}) + ({\hat{G}}_{t + 1}^{1} - {\hat{G}}_{t}^{1})]] ∣ A_{t}) \end{aligned}

and solve this equation forward by first computing the $γ_{o} = 1$ answer and then adjusting this answer for $γ_{o} > 1$ analogous to the approach described in Remark 8.5.

For the order two approximation,

Q_{t}^{2} \overset{def}{=} (1 - ρ)^{2} Q_{t}^{0} {({\hat{R}}_{t}^{1} - {\hat{V}}_{t}^{1})}^{2} + (1 - ρ) Q_{t}^{0} ({\hat{R}}_{t}^{2} - {\hat{V}}_{t}^{2}) .

Express:

{\hat{R}}_{t}^{2} - {\hat{V}}_{t}^{2} = ({\hat{R}}_{t}^{2} - {\hat{G}}_{t}^{2}) - ({\hat{V}}_{t}^{2} - {\hat{G}}_{t}^{2}) .

It follows from (8.29) that

\begin{aligned} {\hat{V}}_{t}^{2} - {\hat{G}}_{t}^{2} = & λ E (N_{t + 1}^{0} [({\hat{V}}_{t + 1}^{2} - {\hat{G}}_{t + 1}^{2}) + ({\hat{G}}_{t + 1}^{2} - {\hat{G}}_{t}^{2})] ∣ A_{t}) \\ + (1 - λ) ({\hat{C}}_{t}^{2} - {\hat{G}}_{t}^{2}) + (1 - ρ) (1 - λ) λ {({\hat{R}}_{t}^{1} - {\hat{G}}_{t}^{1} + {\hat{G}}_{t}^{1} - {\hat{C}}_{t})}^{2}, \end{aligned}

which we solve this equation forward under the $N^{0}$ implied change in probability measure. Form:

\begin{aligned} {\hat{R}}_{t}^{2} - {\hat{G}}_{t}^{2} = & (\frac{1}{λ}) ({\hat{V}}_{t}^{2} - {\hat{G}}_{t}^{2}) - (\frac{1 - λ}{λ}) ({\hat{C}}_{t}^{2} - {\hat{G}}_{t}^{2}) \\ - (1 - ρ) (1 - λ) {({\hat{R}}_{t}^{1} - {\hat{G}}_{t}^{1} + {\hat{G}}_{t}^{1} - {\hat{C}}_{t})}^{2} . \end{aligned}

The $P$ approximations use some of these same computations where

\begin{aligned} P_{t}^{0} \overset{def}{=} & (1 - β) \exp [(ρ - 1) ({\hat{V}}_{t}^{0} - {\hat{G}}_{t}^{0})] \\ P_{t}^{1} \overset{def}{=} & (ρ - 1) P_{t}^{0} ({\hat{V}}_{t}^{1} - {\hat{G}}_{t}^{1}) \\ P_{t}^{2} \overset{def}{=} & (ρ - 1)^{2} P_{t}^{0} {({\hat{V}}_{t}^{1} - {\hat{G}}_{t}^{1})}^{2} + (1 - ρ) P_{t}^{0} ({\hat{V}}_{t}^{2} - {\hat{G}}_{t}^{2}) . \end{aligned}

8.9.3. $N$ derivatives#

N_{t + 1}^{0} \overset{def}{=} \exp [(1 - γ_{o}) ({\tilde{V}}_{t + 1}^{0} - {\tilde{R}}_{t}^{0})] = \exp [(1 - γ_{o}) ({\hat{V}}_{t + 1}^{1} - {\hat{R}}_{t}^{1})]

\begin{aligned} N_{t + 1}^{1} \overset{def}{=} & {\frac{d}{d q} \exp [(1 - γ_{o}) ({\tilde{V}}_{t + 1} - {\tilde{R}}_{t})] |}_{q = 0} \\ = & (1 - γ_{o}) N_{t + 1}^{0} ({\tilde{V}}_{t + 1}^{1} - {\tilde{R}}_{t}^{1}) \\ = & (\frac{1 - γ_{o}}{2}) N_{t + 1}^{0} ({\hat{V}}_{t + 1}^{2} - {\hat{R}}_{t}^{2}) . \end{aligned}

Form:

\begin{aligned} {\hat{V}}_{t + 1}^{1} - {\hat{R}}_{t}^{1} = & ({\hat{V}}_{t + 1}^{1} - {\hat{G}}_{t + 1}^{1}) + ({\hat{G}}_{t + 1}^{1} - {\hat{G}}_{t}^{1}) - ({\hat{R}}_{t}^{1} - {\hat{G}}_{t}^{1}) \\ {\hat{V}}_{t + 1}^{2} - {\hat{R}}_{t}^{2} = & ({\hat{V}}_{t + 1}^{2} - {\hat{G}}_{t + 1}^{2}) + ({\hat{G}}_{t + 1}^{2} - {\hat{G}}_{t}^{2}) - ({\hat{R}}_{t}^{2} - {\hat{G}}_{t}^{2}) . \end{aligned}

It may be directly verified that $N_{t + 1}^{1}$ and $N_{t + 1}^{2}$ have conditional expectations equal to zero. Express

(8.39)#

\begin{aligned} {\hat{V}}_{t + 1}^{1} - {\hat{R}}_{t}^{1} = & (\frac{1}{1 - γ_{o}}) [μ^{0} \cdot (W_{t + 1} - μ^{0}) + \frac{1}{2} μ^{0} \cdot μ^{0}] \\ {\hat{V}}_{t + 1}^{2} - {\hat{R}}_{t}^{2} = & \frac{1}{2} {(W_{t + 1} - μ^{0})}^{'} Υ_{2}^{2} (W_{t + 1} - μ^{0}) - \frac{1}{2} t r (Υ_{2}^{2}) \\ - {(W_{t + 1} - μ^{0})}^{'} (Υ_{1}^{2} X_{t}^{1} + Υ_{0}^{2}) . \end{aligned}

In producing these representations, we use that have conditional ${\hat{V}}_{t + 1}^{2} - {\hat{R}}_{t}^{2}$ and ${\hat{V}}_{t + 1}^{2} - {\hat{R}}_{t}^{2}$ have mean zero under the conditional probability distribution induced by $N_{t + 1}^{0} .$

8.10. Appendix B: Approximation formulas (approach one)#

Consider the equation:

Q_{t} E (N_{t + 1} H_{t + 1} ∣ A_{t}) + P_{t} L_{t} - M_{t} = 0 .

not including the state evolution equations.

8.10.1. Order zero#

The order zero approximation of the product: $N_{t + 1} Q_{t} H_{t + 1} + P_{t} L_{t} - M_{0}$ is:

Q_{t}^{0} N_{t + 1}^{0} H_{t + 1}^{0} + P_{t}^{0} L_{t}^{0} - M_{0} = 0

Thus the order zero approximate equation is:

Q_{t}^{0} E [N_{t + 1}^{0} (H_{t + 1}^{0}) ∣ A_{t}] + L_{t + 1}^{0} = Q_{t}^{0} H_{t + 1}^{0} + P_{t}^{0} L_{t}^{0} - M_{0} = 0

since $N_{t + 1}^{0}$ has conditional expectation equal to one. We add to this subsystem the $q = 0$ state dynamic equation inclusive of jump variables, and we compute a stable steady state solution.

8.10.2. Order one#

The order one approximation of the product: $Q_{t} N_{t + 1} H_{t + 1} + P_{t} L_{t} - M_{t}$ is:

Q_{t}^{1} N_{t + 1}^{0} H_{t + 1}^{0} + Q_{t}^{0} N_{t + 1}^{1} H_{t + 1}^{0} + Q_{t}^{0} N_{t + 1}^{0} H_{t + 1}^{1} + P_{t}^{0} L_{t}^{1} + P_{t}^{1} L_{t}^{0} - M_{t}^{1} .

Thus the order one approximate equation is:

\begin{array}{r} \begin{aligned} Q_{t}^{0} E (N_{t + 1}^{1} H_{t + 1}^{0} + N_{t + 1}^{0} H_{t + 1}^{1} ∣ A_{t}) + Q_{t}^{1} H_{t + 1}^{0} + P_{t}^{0} L_{t}^{1} + P_{t}^{1} L_{t}^{0} \\ = Q_{t}^{0} E (N_{t + 1}^{0} H_{t + 1}^{1} ∣ A_{t}) + Q_{t}^{1} H_{t + 1}^{0} + P_{t}^{0} L_{t}^{1} + P_{t}^{1} L_{t}^{0} - M_{t}^{1} \\ = 0 \end{aligned} \end{array}

where we used the implication that $E (N_{t + 1}^{1} ∣ A_{t}) = 0.$

8.10.3. Order two#

The order two approximation of the product: $Q_{t} N_{t + 1} H_{t + 1} + P_{t} L_{t} - M_{t}$ is:

\begin{array}{r} \begin{aligned} Q_{t}^{0} N_{t + 1}^{0} H_{t + 1}^{2} + Q_{t}^{0} N_{t + 1}^{2} H_{t + 1}^{0} + 2 Q_{t}^{0} N_{t + 1}^{1} H_{t + 1}^{1} \\ + P_{t}^{0} L_{t}^{2} + 2 P_{t}^{1} L_{t}^{1} + P_{t}^{2} L_{t}^{0} - M_{t}^{2} \\ + 2 N_{t + 1}^{1} Q_{t}^{1} H_{t + 1}^{0} + 2 N_{t + 1}^{0} Q_{t}^{1} H_{t + 1}^{1} + N_{t + 1}^{0} Q_{t}^{2} H_{t + 1}^{0} \end{aligned} \end{array}

The terms $Q_{t}^{0} N_{t + 1}^{2} H_{t + 1}^{0}$ and $2 N_{t + 1}^{1} Q_{t}^{1} H_{t + 1}^{0}$ have conditional expectation equal to zero. Thus the approximating equation is:

\begin{aligned} Q_{t}^{0} E (N_{t + 1}^{0} H_{t + 1}^{2} ∣ A_{t}) + P_{t}^{0} L_{t}^{2} + 2 P_{t}^{1} L_{t}^{1} + P_{t}^{2} L_{t}^{0} \\ + 2 Q_{t}^{0} E (N_{t + 1}^{1} H_{t + 1}^{1} ∣ A_{t}) + 2 Q_{t}^{1} E (N_{t + 1}^{0} H_{t + 1}^{1} ∣ A_{t}) + H_{t + 1}^{0} Q_{t}^{2} E (N_{t + ϵ}^{0} ∣ A_{t}) \\ = 0. \end{aligned}

To elaborate on the contributions in the second line, express $H_{t + 1}^{1}$ as

(8.40)#

H_{t + 1}^{1} = Θ_{0}^{1} + Θ_{1}^{1} X_{t}^{1} + Θ_{2}^{1} (W_{t + 1} - μ^{0}) .

Then

\begin{array}{r} \begin{aligned} 2 Q_{t}^{0} E (N_{t + 1}^{1} H_{t + 1}^{1} ∣ A_{t}) & = (1 - γ) Q_{t}^{0} E [N_{t + 1}^{0} ({\hat{V}}_{t + 1}^{2} - {\hat{R}}_{t}^{2}) H_{t + 1}^{1} ∣ A_{t}], \\ = (γ_{o} - 1) Q_{t}^{o} Θ_{2}^{1} [Υ_{1}^{2} X_{t}^{1} + Υ_{0}^{2}] \\ 2 Q_{t}^{1} E (N_{t + 1}^{0} H_{t + 1}^{1} ∣ A_{t}) & = 2 Q_{t}^{0} (1 - ρ) ({\hat{R}}_{t}^{1} - {\hat{V}}_{t}^{1}) (Θ_{0}^{1} + Θ_{1}^{1} X_{t}^{1}), \\ H_{t + 1}^{0} Q_{t}^{2} E (N_{t + 1}^{0} ∣ A_{t}) & = H_{t + 1}^{0} Q_{t}^{0} [(1 - ρ)^{2} {({\hat{R}}_{t}^{1} - {\hat{V}}_{t}^{1})}^{2} + (1 - ρ) ({\hat{R}}_{t}^{2} - {\hat{V}}_{t}^{2})] . \end{aligned} \end{array}

The formula for the first of these terms follows from (8.39) and (8.40), along with fact the third central moments of normals are zero.

We add to this second-order subsystem, the second-order approximation of the state dynamics inclusive of the jump variables. We substitute in the solution for the first-order approximation for the jump variables into both the first and second-order approximate state dynamics. In solving the second-order jump variable adjustment we use expectations induced by $N_{t + 1}^{0}$ zero throughout under which $W_{t + 1}$ is conditionally normally distributed with mean $μ^{0}$ and covariance $I$ .

8.11. Appendix C: Approximation formulas (approach two)#

In this approach we use the same order zero approximation. For the order one approximation, we use the formula (8.30) for ${\tilde{N}}_{t + 1},$ which approximates $N_{t + 1}^{*},$ in conjunction with:

(8.41)#

\begin{array}{r} Q_{t}^{0} E ({\tilde{N}}_{t + 1} H_{t + 1}^{1} ∣ A_{t}) + Q_{t}^{1} H_{t + 1}^{0} + P_{t}^{0} L_{t}^{1} + P_{t}^{1} L_{t}^{0} - M_{t}^{1} = 0. \end{array}

From formula (8.39), it follows that under the ${\tilde{N}}_{t + 1},$ induced change in probability, $W_{t + 1}$ is normally distributed with conditional mean

μ^{0} + (1 - γ_{o}) {(Υ_{2}^{2})}^{- 1} (Υ_{1}^{2} X_{t}^{1} + Υ_{0}^{2}),

and conditional precision:

(γ_{o} - 1) Υ_{2}^{2} + I .

For the order two approximation, we use:

Q_{t}^{0} E ({\tilde{N}}_{t + 1} H_{t + 1}^{2} ∣ A_{t}) + P_{t}^{0} L_{t}^{2} + 2 P_{t}^{1} L_{t}^{1} + P_{t}^{2} L_{t}^{0} - M_{t}^{2} + 2 Q_{t}^{1} E ({\tilde{N}}_{t + 1} H_{t + 1}^{1} ∣ A_{t}) + H_{t + 1}^{0} Q_{t}^{2} = 0 .

8.12. Appendix D: Parameter values#

To facilitate a comparison to a global solution method, we write down a discrete-time approximation to a continuous time version of such an economy. (See Section 4.4 of [Hansen et al., 2024] for a continuous-time benchmark model that our discrete-time system approximates. Note that we convert the annual parameters in that paper to quarterly time.) The parameter settings are:

$δ$	$ι_{k}$	$ζ$	$ν_{k}$	$ν_{1}$	$ν_{2}$	$μ_{2}$
0.025	0.01	32	0.01	0.014	0.0485	$6.3 \times 10^{- 6}$

\begin{array}{r} [\begin{array}{c} σ_{k} \\ σ_{1} \\ σ_{2} \end{array}] = \sqrt{3} [\begin{array}{c} .92 & .40 & 0 \\ 0 & 5.7 & 0 \\ 0 & 0 & .00031 \end{array}] \end{array}

The numbers for $η_{k}, ϕ, β_{1},$ $σ_{k}$ and $σ_{1}$ are such that, when multiplied by stochastic volatility, they match the parameters from [Hansen and Sargent, 2021]. In particular, the constant $Z^{2}$ which scales our $σ_{k}$ to match is $7.6 \times 10^{- 6},$ which is the 67th percentile of our $Z^{2}$ distribution. While [Hansen and Sargent, 2021] use a lower triangular representation for the two-by-two right block of

[\begin{matrix} σ_{k} \\ σ_{1} \end{matrix}],

we use an observationally equivalent upper triangular representation for most of the results. Finally, the numbers for $β_{2}$ and $σ_{2}$ come from [Schorfheide et al., 2018], but they are adjusted for approximation purposes as described in Appendix A [Hansen et al., 2024]. In both cases, we use the medians of their econometric evidence as input into our analysis.

For the extension to the habit persistence model in Section 8.7.5, we use a habit persistence of $ν_{h} = 0.025$ .

Exploring Recursive Utility

Contents

8. Exploring Recursive Utility#

8.1. Introduction#

8.2. Recursive utility valuation process#

8.2.1. Basic recursion#

8.2.2. Preference for robustness#

8.2.3. Stochastic discount factor process#

8.3. Continuous-time limit#

8.3.1. Discrete-time approximation#

8.3.2. Robustness to misspecification#

8.3.3. Uncertainty pricing#

8.4. Small noise expansion of dynamic stochastic equilibria#

8.4.1. Approximate state dynamics#

8.5. Incorporating preferences with enhanced uncertainty concerns#

8.5.1. Order-zero#

8.5.2. Order-one#

8.5.3. Order two#

8.6. Stochastic discount factor approximation#

8.6.1. Approach 1#

8.6.2. Approach 2#

8.7. Solving a planner’s problem with recursive utility#

8.7.1. An example economy with long-run uncertainty#

8.7.2. First-order approximation when ρ=1#

8.7.3. Second-order approximation when ρ=1#

8.7.4. Shock elasticities#

8.7.5. An alternative model of intertemporal substitution/complementarity.#

8.7.6. External habits model#

8.8. Solving models#

8.9. Appendix A: Solving the planner’s problem#

8.9.1. Some steady state calculations#

8.9.2. Q and P derivatives#

8.9.3. N derivatives#

8.10. Appendix B: Approximation formulas (approach one)#

8.10.1. Order zero#

8.10.2. Order one#

8.10.3. Order two#

8.11. Appendix C: Approximation formulas (approach two)#

8.12. Appendix D: Parameter values#

8.7.2. First-order approximation when $ρ = 1$ #

8.7.3. Second-order approximation when $ρ = 1$ #

8.9.2. $Q$ and $P$ derivatives#

8.9.3. $N$ derivatives#