-
Notifications
You must be signed in to change notification settings - Fork 19
/
Copy path_sectn-method-of-moderation.tex
739 lines (636 loc) · 49.2 KB
/
_sectn-method-of-moderation.tex
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
\input{@resources/tex-add-search-paths}\documentclass[SolvingMicroDSOPs]{subfiles}
\input{subfile-start}
\begin{document}
\hypertarget{the-method-of-moderation}{}
\subsection{The Method of Moderation}\label{sec:method-of-moderation}
\begin{verbatimwrite}{./cctwMoM/EndogGptsProbs.tex}
Unfortunately, this endogenous gridpoints solution is not very
well-behaved outside the original range of gridpoints targeted by
the solution method. (Though other common solution methods are no
better outside their own predefined ranges).
Figure~\ref{fig:ExtrapProblem} demonstrates the point by plotting
the amount of precautionary saving implied by a linear extrapolation
of our approximated consumption rule (the consumption of the perfect
foresight consumer $\cFuncAbove_{\prd-1}$ minus our approximation to
optimal consumption under uncertainty, $\Aprx{\cFunc}_{\prd-1}$).
Although theory proves that precautionary saving is always positive,
the linearly extrapolated numerical approximation eventually
predicts negative precautionary saving (at the point in the figure
where the extrapolated locus crosses the horizontal axis).
\hypertarget{ExtrapProblemPlot}{}
\begin{figure}
\includegraphics[width=6in]{./Figures/ExtrapProblemPlot}
\caption{For Large Enough $m_{\prd-1}$, Predicted Precautionary Saving is Negative (Oops!)}
\label{fig:ExtrapProblem}
\end{figure}
This error cannot be fixed by extending the upper gridpoint; in the presence of serious uncertainty, the consumption rule will need to be evaluated outside of \textit{any} prespecified grid (because starting from the top gridpoint, a large enough realization of the uncertain variable will push next period's realization of assets above that top; a similar argument applies below the bottom gridpoint). While a judicious extrapolation technique can prevent this problem from being fatal (for example by carefully excluding negative precautionary saving), the problem is often dealt with using inelegant methods whose implications for the accuracy of the solution are difficult to gauge.
\end{verbatimwrite}
\input{./cctwMoM/EndogGptsProbs.tex}\unskip
%\renewcommand{\prd}{t} % For the rest of the doc, use generic t vs t+1
\begin{verbatimwrite}{./cctwMoM/MoM-Prelims.tex}
As a preliminary to our solution, define $\hNrm_{\EndStg}$ as end-of-period human wealth (the present discounted value of future labor income) for a perfect foresight version of the problem of a `risk optimist:' a period-$t$ consumer who believes with perfect confidence that the shocks will always take their expected value of \PermShkOn{1, $\tranShkEmp_{t+n} = \Ex[\tranShkEmp]=1~\forall~n>0$ and $\permShk_{t+n} = \Ex[\permShk]=1~\forall~n>0$.} {1, $\tranShkEmp_{t+n} = \Ex[\tranShkEmp]=1~\forall~n>0$.} The solution to a perfect foresight problem of this kind takes the form\footnote{For a derivation, see \cite{BufferStockTheory}; $\MPCmin_{\prd}$ is defined therein as the MPC of the perfect foresight consumer with horizon $\trmT-t$.}
\end{verbatimwrite}
\input{./cctwMoM/MoM-Prelims.tex}\unskip
\begin{verbatimwrite}{./Equations/cFuncAbove.tex}
\begin{equation}\begin{gathered}\begin{aligned}
\cFuncAbove_{\prd}(\mNrm_{\prd}) & = (\mNrm_{\prd} + \hNrm_{\EndStg})\MPCmin_{\prd} \label{eq:cFuncAbove}
\end{aligned}\end{gathered}\end{equation}
for a constant minimal marginal propensity to consume $\MPCmin_{\prd}$ given below.
\end{verbatimwrite}
\input{./Equations/cFuncAbove.tex}\unskip
\begin{verbatimwrite}{./cctwMoM/MoM-Words.tex}
We similarly define $\hEndMin_{\EndStg}$ as `minimal human wealth,' the
present discounted value of labor income if the shocks were to take on
their worst possible value in every future period \PermShkOn
{$\tranShkEmp_{t+n} = \tranShkEmpMin ~\forall~n>0$ and $\permShk_{t+n} =
\permShkMin ~\forall~n>0$} {$\tranShkEmp_{t+n} = \tranShkEmpMin
~\forall~n>0$} (which we define as corresponding to the beliefs of a
`pessimist').
\ctw{}{We will call a `realist' the consumer who correctly perceives the true
probabilities of the future risks and optimizes accordingly.}
A first useful point is that, for the realist, a lower bound for the
level of market resources is $\Min{m}_{\prd} = -\hEndMin_{\EndStg}$, because
if $m_{\prd}$ equalled this value then there would be a positive finite
chance (however small) of receiving \PermShkOn
{$\tranShkEmp_{t+n}=\tranShkEmpMin$ and $\permShk_{t+n}=\permShkMin$}
{$\tranShkEmp_{t+n}=\tranShkEmpMin$}
in
every future period, which would require the consumer to set $c_{\prd}$
to zero in order to guarantee that the intertemporal budget constraint
holds\ctw{.}{~(this is the multiperiod generalization of the discussion in
section \ref{subsec:LiqConstrSelfImposed} explaining the derivation of the `natural borrowing constraint' for period $\trmT-1$,
$\Min{a}_{\prd-1}$).} Since consumption of zero yields negative
infinite utility, the solution to realist consumer's problem is not well
defined for values of $m_{\prd} < \Min{m}_{\prd}$, and the limiting
value of the realist's $c_t$ is zero as $m_{\prd} \downarrow \Min{m}_{\prd}$.
Given this result, it will be convenient to define `excess' market
resources as the amount by which actual resources exceed the lower
bound, and `excess' human wealth as the amount by which mean expected human wealth
exceeds guaranteed minimum human wealth:
\begin{equation*}\begin{gathered}\begin{aligned}
\aboveMin \mNrm_{\prd} & = m_{\prd}+\overbrace{\hEndMin_{\EndStg}}^{=-\Min{m}_{\prd}}
\\ \aboveMin \hNrm_{\EndStg} & = \hNrm_{\EndStg}-\hEndMin_{\EndStg}.
\end{aligned}\end{gathered}\end{equation*}
We can now transparently define the optimal
consumption rules for the two perfect foresight problems, those of the
`optimist' and the `pessimist.' The `pessimist' perceives human
wealth to be equal to its minimum feasible value $\hEndMin_{\EndStg}$ with certainty, so
consumption is given by the perfect foresight solution
\begin{equation*}\begin{gathered}\begin{aligned}
\cFuncBelow_{\prd}(m_{\prd}) & = (m_{\prd}+\hEndMin_{\EndStg})\MPCmin_{\prd}
\\ & = \aboveMin \mNrm_{\prd}\MPCmin_{\prd}
.
\end{aligned}\end{gathered}\end{equation*}
The `optimist,' on the other hand, pretends that there is no uncertainty
about future income, and therefore consumes
\begin{equation*}\begin{gathered}\begin{aligned}
\cFuncAbove_{\prd}(m_{\prd}) & = (m_{\prd} +\hEndMin_{\EndStg} - \hEndMin_{\EndStg} + \hNrm_{\EndStg} )\MPCmin_{\prd}
\\ & = (\aboveMin \mNrm_{\prd} + \aboveMin \hNrm_{\EndStg})\MPCmin_{\prd}
\\ & = \cFuncBelow_{\prd}(m_{\prd})+\aboveMin \hNrm_{\EndStg} \MPCmin_{\prd}
.
\end{aligned}\end{gathered}\end{equation*}
It seems obvious that the spending of the realist will be strictly greater
than that of the pessimist and strictly less than that of the
optimist. Figure~\ref{fig:IntExpFOCInvPesReaOptNeedHiPlot} illustrates the proposition for the consumption rule in period $\trmT-1$.
\end{verbatimwrite}
\input{./cctwMoM/MoM-Words.tex}\unskip
\begin{verbatimwrite}{./Figures/IntExpFOCInvPesReaOptNeedHiPlot.tex}
\hypertarget{IntExpFOCInvPesReaOptNeedHiPlot}{}
\begin{figure}
\includegraphics[width=6in]{./Figures/IntExpFOCInvPesReaOptNeedHiPlot}
\caption{Moderation Illustrated: $\Min{\cFunc}_{\prd-1} < \Aprx{\cFunc}_{\prd-1} < \bar{\cFunc}_{\prd-1}$}
\label{fig:IntExpFOCInvPesReaOptNeedHiPlot}
\end{figure}
\end{verbatimwrite}
\input{./Figures//IntExpFOCInvPesReaOptNeedHiPlot.tex}\unskip
\begin{verbatimwrite}{./cctwMoM/MoM-Words-Rest.tex}
\indent The proof is more difficult than might be imagined, but
the necessary work is done in \cite{BufferStockTheory} so we will take
the proposition as a fact and proceed by manipulating the inequality:
\end{verbatimwrite}
\input{./cctwMoM/MoM-Words-Rest.tex}\unskip
\begin{verbatimwrite}{./Equations/MoM-Inequalities.tex}
\begin{center}
\begin{tabular}{rcl}
$ \aboveMin \mNrm_{\prd} \MPCmin_{\prd} < $ & $ \cFunc_{\prd}(\Min{m}_{\prd}+\aboveMin \mNrm_{\prd}) $ $< (\aboveMin \mNrm_{\prd}+\aboveMin \hNrm_{\EndStg})\MPCmin_{\prd} $
\\ $- \aboveMin \mNrm_{\prd} \MPCmin_{\prd} > $ & $ -\cFunc_{\prd}(\Min{m}_{\prd}+\aboveMin \mNrm_{\prd}) $ & $> -(\aboveMin \mNrm_{\prd}+\aboveMin \hNrm_{\EndStg})\MPCmin_{\prd} $
\\ $ \aboveMin \hNrm_{\EndStg} \MPCmin_{\prd} > $ & $ \bar{\cFunc}_{\prd}(\Min{m}_{\prd}+\aboveMin \mNrm_{\prd})-\cFunc_{\prd}(\Min{m}_{\prd}+\aboveMin \mNrm_{\prd}) $ & $> 0$
\\ $1 > $ & $ \underbrace{\left(\frac{\bar{\cFunc}_{\prd}(\Min{m}_{\prd}+\aboveMin \mNrm_{\prd})-\cFunc_{\prd}(\Min{m}_{\prd}+\aboveMin \mNrm_{\prd})}{\aboveMin \hNrm_{\EndStg} \MPCmin_{\prd}}\right)}_{\equiv \Max{\koppa}_{\prd}} $ & $> 0$
\end{tabular}
\end{center}
\end{verbatimwrite}
\input{./Equations/MoM-Inequalities.tex}\unskip
\begin{verbatimwrite}{./cctwMoM/MoM-Inequalities-Describe.tex}
\noindent where the fraction in the middle of the last inequality is the ratio
of actual precautionary saving (the numerator is the difference
between perfect-foresight consumption and optimal consumption in the
presence of uncertainty) to the maximum conceivable amount of
precautionary saving (the amount that would be undertaken by the
pessimist who consumes nothing out of any future income beyond the perfectly certain component).
\end{verbatimwrite}
\input{./cctwMoM/MoM-Inequalities-Describe.tex}\unskip
\begin{verbatimwrite}{./Equations/MoM-KoppaOfMu.tex}
Defining $\mu_{\prd} =
\log \aboveMin \mNrm_{\prd}$ (which can range from $-\infty$ to $\infty$), the object in the middle of the last inequality is
\begin{equation}\begin{gathered}\begin{aligned}
\Max{\koppa}_{\prd}(\mu_{\prd}) & \equiv \left(\frac{\bar{\cFunc}_{\prd}(\Min{m}_{\prd}+e^{\mu_{\prd}})-\cFunc_{\prd}(\Min{m}_{\prd}+e^{\mu_{\prd}})}{\aboveMin \hNrm_{\EndStg} \MPCmin_{\prd}}\right), \label{eq:koppa}
\end{aligned}\end{gathered}\end{equation}
and we now define
\begin{equation}\begin{gathered}\begin{aligned}
\Max{\chiFunc}_{\prd}(\mu_{\prd}) & = \log \left(\frac{1-\Max{\koppa}_{\prd}(\mu_{\prd})}{\Max{\koppa}_{\prd}(\mu_{\prd})}\right)
\\ & = \log \left(1/\Max{\koppa}_{\prd}(\mu_{\prd})-1\right) \label{eq:chi}
\end{aligned}\end{gathered}\end{equation}
\end{verbatimwrite}
\input{./Equations/MoM-KoppaOfMu.tex}\unskip
\begin{verbatimwrite}{./cctwMoM/MoM-KoppaOfMu-Describe.tex}
which has the virtue that it is linear in the limit as $\mu_{\prd}$ approaches $+\infty$.
Given $\Max{\chiFunc}$, the consumption function can be recovered from
\end{verbatimwrite}
\input{./cctwMoM/MoM-KoppaOfMu-Describe.tex}\unskip
\begin{verbatimwrite}{./Equations/cFuncHi.tex}
\begin{equation}\begin{gathered}\begin{aligned}
\Max{\cFunc}_{\prd} & = \bar{\cFunc}_{\prd}-\overbrace{\left(\frac{1}{1+\exp(\Max{\chiFunc}_{\prd})}\right)}^{=\Max{\koppa}_{\prd}} \aboveMin \hNrm_{\EndStg} \MPCmin_{\prd}. \label{eq:cFuncHi}
\end{aligned}\end{gathered}\end{equation}
\end{verbatimwrite}
\input{./Equations/cFuncHi.tex}\unskip
\begin{verbatimwrite}{./cctwMoM/MoM-End.tex}
Thus, the procedure is to calculate $\Max{\chiFunc}_{\prd}$ at the points
$\vctr{\mu}_{\prd}$ corresponding to the log of the $\aboveMin
\vctr{m}_{\prd}$ points defined above, and then using these to construct an
interpolating approximation $\Aprx{\Max{\chiFunc}}_{\prd}$ from which we indirectly obtain our
approximated consumption rule $\Aprx{\Max{\cFunc}}_{\prd}$ by substituting $\Aprx{\Max{\chiFunc}}_{\prd}$ for $\Max{\chiFunc}$ in equation \eqref{eq:cFuncHi}.
Because this method relies upon the fact that the problem is easy to
solve if the decision maker has unreasonable views (either in the
optimistic or the pessimistic direction), and because the correct
solution is always between these immoderate extremes, we call our
solution procedure the `method of moderation.'
Results are shown in Figure~\ref{fig:ExtrapProblemSolved}; a reader
with very good eyesight might be able to detect the barest hint of a
discrepancy between the Truth and the Approximation at the far
righthand edge of the figure\ctw{.}{ -- a stark contrast with the calamitous
divergence evident in Figure~\ref{fig:ExtrapProblem}.}{}
\hypertarget{ExtrapProblemSolvedPlot}{}
\begin{figure}
\includegraphics[width=6in]{./Figures/ExtrapProblemSolvedPlot}
\caption{Extrapolated $\Aprx{\Max{\cFunc}}_{\prd-1}$ Constructed Using the Method of Moderation}
\label{fig:ExtrapProblemSolved}
\end{figure}
\end{verbatimwrite}
\input{./cctwMoM/MoM-End.tex}\unskip
\hypertarget{approximating-the-slope-too}{}
\subsection{Approximating the Slope Too}
Until now, we have calculated the level of consumption at various different gridpoints and used linear interpolation\ctw{.}{ (either directly for $\cFunc_{\prd-1}$ or indirectly for, say, $\Max{\chiFunc}_{\prd-1}$).} But the resulting piecewise linear approximations have the unattractive feature that they are not differentiable at the `kink points' that correspond to the gridpoints where the slope of the function changes discretely.
\cite{BufferStockTheory} proves that the true consumption function for
this problem
is `smooth:' It
exhibits a well-defined unique marginal propensity to consume at every
positive value of $m$. This suggests that we should calculate, not
just the level of consumption, but also the marginal propensity to
consume (henceforth $\MPC$) at each gridpoint, and then find an
interpolating approximation that smoothly matches both the level and the slope
at those points.
This requires us to differentiate \eqref{eq:koppa} and \eqref{eq:chi}, yielding
\begin{equation}\begin{gathered}\begin{aligned}
\Max{\koppa}_{\prd}^{\mu}(\mu_{\prd}) & = (\aboveMin \hNrm_{\EndStg} \MPCmin_{\prd})^{-1}e^{\mu_{\prd}}\left(\MPCmin_{\prd}-\overbrace{\cFunc^{\mNrm}_{\prd}(\Min{m}_{\prd}+e^{\mu_{\prd}})}^{\equiv \MPCFunc_{\prd}(\mNrm_{\prd})}\right) \label{eq:koppaPrime}
\\ \Max{\chiFunc}_{\prd}^{\mu}(\mu_{\prd}) & = \left(\frac{-\Max{\koppa}_{\prd}^{\mu}(\mu_{\prd})/\Max{\koppa}_{\prd}^{2}}{1/\Max{\koppa}_{\prd}(\mu_{\prd})-1}\right)
\end{aligned}\end{gathered}\end{equation}
and (dropping arguments) with some algebra these can be combined to yield
\begin{equation}\begin{gathered}\begin{aligned}
\Max{\chiFunc}_{\prd}^{\mu} & = \left(\frac{\MPCmin_{\prd} \aboveMin \mNrm_{\prd} \aboveMin \hNrm_{\EndStg} (\MPCmin_{\prd}-\MPC_{\prd})}
{(\cFuncAbove_{\prd}-\cFunc_{\prd})(\cFuncAbove_{\prd}-\cFunc_{\prd} - \MPCmin_{\prd} \aboveMin \hNrm_{\EndStg})}\right).
\end{aligned}\end{gathered}\end{equation}
To compute the vector of values of \eqref{eq:koppaPrime} corresponding
to the points in $\vctr{\mu}_{\prd}$, we need the marginal propensities to
consume (designated $\MPC$) at each of the gridpoints,
$\cFunc^{\mNrm}_{\prd}$ (the vector of such values is
$\vctr{\MPC}_{\prd}$). These can be obtained by differentiating the
Euler equation \eqref{eq:upEqbetaOp} (where we define
$\mFunc_{\EndStg}(a) \equiv \cFunc_{\EndStg}(a)+{a}$, and drop the (a) arguments to reduce clutter):
\begin{equation}\begin{gathered}\begin{aligned}
\uFunc^{c}(\cFunc_{\EndStg}) & = \hat{\vFunc}_{\EndStg}^{\aNrm}(\mFunc_{\EndStg}-\cFunc_{\EndStg}),
\end{aligned}\end{gathered}\end{equation}
yielding a marginal propensity to
\textit{have consumed} $\cFunc_{\EndStg}^{\aNrm}$ at each gridpoint:
\begin{equation}\begin{gathered}\begin{aligned}
\uPP(\cEndStg)\cEndStg^\aNrm & = \hat{\vFunc}_{\EndStg}^{\aNrm}(\mFunc_{\EndStg}-\cFunc_{\EndStg})
\\ \cEndStg^\aNrm & = \hat{\vFunc}^{\aNrm}(\mFunc_{\EndStg}-\cFunc_{\EndStg})/\uPP(\cEndStg)
\end{aligned}\end{gathered}\end{equation}
and the marginal propensity to consume at the beginning of the period is obtained from the marginal propensity to have consumed by differentiating the identity with respect to $\aNrm$:
\begin{equation*}\begin{gathered}\begin{aligned}
\cEndStg & = \mFunc_{\EndStg} - \aNrm
\\ \cEndStg^{\aNrm}+1 & = \mFunc_{\EndStg}^{\aNrm}
\end{aligned}\end{gathered}\end{equation*}
which, together with the chain rule $\cEndStg^\aNrm = \cFunc^{\mNrm}_{\MidStg}\mFunc_{\EndStg}^{\aNrm}$, yields the MPC from
\begin{equation}\begin{gathered}\begin{aligned}
\cFunc^{\mNrm}(\overbrace{\cEndStg^{\aNrm}+1}^{=\mFunc_{\EndStg}^{\aNrm}}) & = \cEndStg^{\aNrm}
\\ \cFunc^\mNrm & = \cEndStg^{\aNrm}/(1+\cEndStg^{\aNrm}) \label{eq:MPCfromMPTHC}.
\end{aligned}\end{gathered}\end{equation}
Designating $\Aprx{\Max{\cFunc}}_{\prd-1}$ as the approximated consumption rule obtained using an interpolating polynomial approximation to $\Max{\chiFunc}$ that matches both the level and the first derivative at the gridpoints, Figure~\ref{fig:IntExpFOCInvPesReaOptGapPlot} plots the difference between this latest approximation and the true consumption rule for period $\trmT-1$ up to the same large value (far beyond the largest gridpoint) used in prior figures. Of course, at the gridpoints the approximation will exactly match the true function; but this figure illustrates that the approximation is quite accurate far beyond the last gridpoint (which is the last point at which the difference touches the horizontal axis). (We plot here the difference between the two functions rather than the level plotted in previous figures, because in levels the difference between the approximate and the exact function would not be detectable even to the most eagle-eyed reader.)
\hypertarget{IntExpFOCInvPesReaOptGapPlot}{}
\begin{figure}
\includegraphics[width=6in]{./Figures/IntExpFOCInvPesReaOptGapPlot}
\caption{Difference Between True $\cFunc_{\prd-1}$ and $\Aprx{\Max{\cFunc}}_{\prd-1}$ Is Minuscule}
\label{fig:IntExpFOCInvPesReaOptGapPlot}
\end{figure}
\hypertarget{value}{}
\subsection{Value}
\begin{verbatimwrite}{./cctwMoM/value-Intro.tex}
Often it is useful to know the value function as well as the consumption rule. Fortunately, many of the tricks used when solving for the consumption rule have a direct analogue in approximation of the value function.
Consider the perfect foresight (or ``optimist's'') problem in period $\trmT-1$. Using the fact that in a perfect foresight model the growth factor for consumption is $(\Rfree \DiscFac)^{1/\CRRA}$, we can use the fact that $\cNrm_{\prd} = (\Rfree \DiscFac)^{1/\CRRA} \cNrm_{\prd-1}$ to calculate the value function in period $\trmT-1$:
\begin{equation*}\begin{gathered}\begin{aligned}
\bar{\vFunc}_{\prd-1}(m_{\prd-1}) & \equiv \uFunc(\cNrm_{\prd-1})+\DiscFac \uFunc(\cNrm_{\prd})
\\ & = \uFunc(\cNrm_{\prd-1})\left(1+\DiscFac ((\DiscFac\Rfree)^{1/\CRRA})^{1-\CRRA}\right)
% \\ & = \uFunc(\cNrm_{\prd-1})\left(1+\DiscFac (\DiscFac\Rfree)^{1/\CRRA-1}\right)
\\ & = \uFunc(\cNrm_{\prd-1})\left(1+(\DiscFac\Rfree)^{1/\CRRA}/\Rfree\right)
\\ & = \uFunc(\cNrm_{\prd-1})\underbrace{\mbox{PDV}_{\prd}^{T}(\cNrm)/\cNrm_{\prd-1}}_{\equiv \PDVCoverc_{\prd-1}^{T}}
\end{aligned}\end{gathered}\end{equation*}
where $\PDVCoverc_{\prd}^{T}=\mbox{PDV}_{\prd}^{T}(\cNrm)$ is the present discounted value of consumption, normalized by current consumption. Using the fact demonstrated in \cite{BufferStockTheory} that $\PDVCoverc_{\prd}=\MPC^{-1}_{\prd}$, a similar function can be constructed recursively for earlier periods, yielding the general expression \hypertarget{vFuncPF}{}
\end{verbatimwrite}
\input{./cctwMoM/value-Intro.tex}\unskip
\begin{verbatimwrite}{./Equations/vFuncPF.tex}
\begin{equation}\begin{gathered}\begin{aligned}
\bar{\vFunc}_{\prd}(m_{\prd}) & = \uFunc(\bar{\cNrm}_{\prd})\PDVCoverc_{\prd}^{T}\label{eq:vFuncPF}
\\ & = \uFunc(\bar{c}_{\prd}) \MPCmin_{\prd}^{-1} % 20190820
\\ & = \uFunc((\aboveMin \mNrm_{\prd}+\aboveMin \hNrm_{\EndStg})\MPCmin_{\prd}) \MPCmin_{\prd}^{-1} % 20190820
\\ & = \uFunc(\aboveMin \mNrm_{\prd}+\aboveMin \hNrm_{\EndStg})\MPCmin_{\prd}^{1-\CRRA} \MPCmin_{\prd}^{-1} % 20190820
\\ & = \uFunc(\aboveMin \mNrm_{\prd}+\aboveMin \hNrm_{\EndStg})\MPCmin_{\prd}^{-\CRRA} % 20190820
\end{aligned}\end{gathered}\end{equation}
This can be transformed as
\begin{equation*}\begin{gathered}\begin{aligned}
\bar{\vInv}_{\prd} & \equiv \left((1-\CRRA)\bar{\vFunc}_{\prd}\right)^{1/(1-\CRRA)}
\\ & = \cNrm_{\prd}(\PDVCoverc_{\prd}^{T})^{1/(1-\CRRA)}
\\ & = (\aboveMin \mNrm_{\prd}+\aboveMin \hNrm_{\EndStg})\MPCmin_{\prd}^{-\CRRA/(1-\CRRA)} % 20190820
\end{aligned}\end{gathered}\end{equation*}
\end{verbatimwrite}
\input{./Equations/vFuncPF.tex}\unskip
\begin{verbatimwrite}{./cctwMoM/value-Rest.tex}
\MPCMatch{with derivative
\begin{equation*}\begin{gathered}\begin{aligned}
\bar{\vInv}_{\prd}^m & = (\mathbb{C}_{\prd}^{T})^{1/(1-\CRRA)}\MPCmin_{\prd},
\\ & = \MPCmin_{\prd}^{-\CRRA/(1-\CRRA)} % 20190820
\end{aligned}\end{gathered}\end{equation*}}{}
and since $\PDVCoverc_{\prd}^{T}$ is a constant while the consumption
function is linear, $\bar{\vInv}_{\prd}$ will also be linear.
We apply the same transformation to the value function for the problem with uncertainty (the ``realist's'' problem)\MPCMatch{ and differentiate}:
\begin{equation*}\begin{gathered}\begin{aligned}
\bar{\vInv}_{\prd} & = \left((1-\CRRA)\bar{\vFunc}_{\prd}(m_{\prd})\right)^{1/(1-\CRRA)}
\MPCMatch{\\ \bar{\vInv}^{m}_{\prd} & = \left((1-\CRRA)\bar{\vFunc}_{\prd}(m_{\prd})\right)^{-1+1/(1-\CRRA)}\bar{\vFunc}_{\prd}^{m}(m_{\prd})}{}
\end{aligned}\end{gathered}\end{equation*}
and an excellent approximation to the value function can be obtained by
calculating the values of $\bar{\vInv}$ at the same gridpoints used by the
consumption function approximation, and interpolating among those points.
However, as with the consumption approximation, we can do even better if we
realize that the $\bar{\vInv}$ function for the optimist's problem is
an upper bound for the ${\vInv}$ function in the presence of uncertainty, and the value function
for the pessimist is a lower bound. Analogously to \eqref{eq:koppa}, define an upper-case
\begin{equation}\begin{gathered}\begin{aligned}
\hat{\Koppa}_{\prd}(\mu_{\prd}) & = \left(\frac{\bar{\vInv}_{\prd}(\Min{m}_{\prd}+e^{\mu_{\prd}})-\vInv_{\prd}(\Min{m}_{\prd}+e^{\mu_{\prd}})}{\aboveMin \hNrm_{\EndStg} \MPCmin_{\prd} (\PDVCoverc_{\prd}^{T})^{1/(1-\CRRA)}}\right) \label{eq:Koppa}
\end{aligned}\end{gathered}\end{equation}
\MPCMatch{with derivative (dropping arguments)
\begin{equation}\begin{gathered}\begin{aligned}
\hat{\Koppa}_{\prd}^{\mu} & = (\aboveMin \hNrm_{\EndStg} \MPCmin_{\prd} (\PDVCoverc_{\prd}^{T})^{1/(1-\CRRA)})^{-1}e^{\mu_{\prd}}\left(\bar{\vInv}^{m}_{\prd}-\vInv^{m}_{\prd}\right) \label{eq:KoppaPrime}
% \\ & = (\aboveMin \hNrm_{\EndStg} \MPCmin_{\prd})^{-1}e^{\mu_{\prd}}\left((\PDVCoverc_{\prd}^{T})^{1/(1-\CRRA)}\MPCmin_{\prd}-\left((1-\CRRA)\vFunc_{\prd}(m_{\prd})\right)^{-1+1/(1-\CRRA)}\vFunc_{\prd}^{m}(m_{\prd})\right) \notag
\end{aligned}\end{gathered}\end{equation}}{}
and an upper-case version of the $\chiFunc$ equation in \eqref{eq:chi}:
\begin{equation}\begin{gathered}\begin{aligned}
\hat{\Chi}_{\prd}(\mu_{\prd}) & = \log \left(\frac{1-\hat{\Koppa}_{\prd}(\mu_{\prd})}{\hat{\Koppa}_{\prd}(\mu_{\prd})}\right)
\\ & = \log \left(1/\hat{\Koppa}_{\prd}(\mu_{\prd})-1\right) \label{eq:Chi}
\end{aligned}\end{gathered}\end{equation}
\MPCMatch{with corresponding derivative
\begin{equation}\begin{gathered}\begin{aligned}
\hat{\Chi}_{\prd}^{\mu} & = \left(\frac{-\hat{\Koppa}_{\prd}^{\mu}/\hat{\Koppa}_{\prd}^{2}}{1/\hat{\Koppa}_{\prd}-1}\right)
\end{aligned}\end{gathered}\end{equation}}{}
and if we approximate these objects then invert them (as above with
the $\Max{\koppa}$ and $\Max{\chiFunc}$ functions) we obtain a very high-quality
approximation to our inverted value function at the same points for
which we have our approximated value function:
\begin{equation}\begin{gathered}\begin{aligned}
\hat{\vInv}_{\prd} & = \bar{\vInv}_{\prd}-\overbrace{\left(\frac{1}{1+\exp(\hat{\Chi}_{\prd})}\right)}^{=\hat{\Koppa}_{\prd}} \aboveMin \hNrm_{\EndStg} \MPCmin_{\prd} (\PDVCoverc_{\prd}^{T})^{1/(1-\CRRA) }
\end{aligned}\end{gathered}\end{equation}
from which we obtain our approximation to the value function\MPCMatch{ and its derivatives~}~as \hypertarget{vHatFunc}{}
\begin{equation}\begin{gathered}\begin{aligned}
\hat{\vFunc}_{\prd} & = \uFunc(\hat{\vInv}_{\prd})
\\ \hat{\vFunc}^{m}_{\prd} & = \uFunc^{c}(\hat{\vInv}_{\prd}) \hat{\vInv}^{m}
\MPCMatch{\\ \hat{\vFunc}^{mm}_{\prd} & = \uFunc^{c{c}}(\hat{\vInv}_{\prd}) (\hat{\vInv}^{m})^{2} + \uFunc^{c}(\hat{\vInv}_{\prd})\hat{\vInv}^{mm}}{}
.
\end{aligned}\end{gathered}\end{equation}
Although a linear interpolation that matches the level of $\vInv$ at the gridpoints is simple, a Hermite interpolation that matches both the level and the derivative of the $\bar{\vInv}_{\prd}$ function at the gridpoints has the considerable virtue that the $\bar{\vFunc}_{\prd}$ derived from it numerically satisfies the envelope theorem at each of the gridpoints for which the problem has been solved.
\MPCMatch{If we use the double-derivative calculated above to produce a higher-order Hermite polynomial, our approximation will also match
marginal propensity to consume at the gridpoints; this would
guarantee that the consumption function generated from the value
function would match both the level of consumption and the
marginal propensity to consume at the gridpoints; the numerical
differences between the newly constructed consumption function and
the highly accurate one constructed earlier would be negligible
within the grid.}{}
\end{verbatimwrite}
\input{./cctwMoM/value-Rest.tex}\unskip
\hypertarget{refinement-a-tighter-upper-bound}{}
\subsection{Refinement: A Tighter Upper Bound}
\begin{verbatimwrite}{./cctwMoM/Tighter.tex}
\cite{BufferStockTheory} derives an upper limit $\MPCmax_{\prd}$ for the MPC as $m_{\prd}$
approaches its lower bound. Using this
fact plus the strict concavity of the consumption function yields the
proposition that
\begin{equation}\begin{gathered}\begin{aligned}
\cFunc_{\prd}(\Min{m}_{\prd}+\aboveMin \mNrm_{\prd}) & < \MPCmax_{\prd} \aboveMin \mNrm_{\prd}.
\end{aligned}\end{gathered}\end{equation}
The solution method described above does not guarantee that
approximated consumption will respect this constraint between gridpoints, and a failure to
respect the constraint can occasionally cause computational problems in solving
or simulating the model. Here, we
describe a method for constructing an approximation that always
satisfies the constraint.
\begin{comment} % Old text needs to be revised or eliminated
That is, the realist's consumption function is bounded from above by both
the \textit{unconstrained} optimist's problem already treated, as well as
by the \textit{constrained} optimist's problem, which is a 45 degree line
originating from $\Min{m}_{\prd}$ on the $m$-axis, as shown in
Figure~\ref{fig:IntExpFOCInvPesReaOptNeed45Plot}. The same is true for
the value function, as illustrated in Figure
\ref{fig:IntExpFOCInvPesReaOptNeed45ValuePlot}.
\hypertarget{IntExpFOCInvPesReaOptNeed45Plot}{}
\begin{figure}
\includegraphics[width=6in]{./Figures/IntExpFOCInvPesReaOptNeed45Plot}
\caption{45 Degree Line as Another Upper Bound}
\label{fig:IntExpFOCInvPesReaOptNeed45Plot}
\end{figure}
\hypertarget{IntExpFOCInvPesReaOptNeed45ValuePlot}{}
\begin{figure}
\includegraphics[width=6in]{./Figures/IntExpFOCInvPesReaOptNeed45ValuePlot}
\caption{A Constrained Optimist's Value Function as Another Upper Bound}
\label{fig:IntExpFOCInvPesReaOptNeed45ValuePlot}
\end{figure}
\end{comment}
\newcommand{\mtCusp}{\ensuremath{\mNrm_{\prd}^{\#}}}
% \newcommand{\aboveMin \mtCusp}{\ensuremath{\aboveMin \mNrm_{\prd}^{\#}}}
\end{verbatimwrite}
\input{./cctwMoM/Tighter.tex}\unskip
\begin{verbatimwrite}{./Equations/mtCusp.tex}
Defining $\mtCusp$ as the `cusp' point where the two upper bounds
intersect:
\begin{equation*}\begin{gathered}\begin{aligned}
\left(\aboveMin \mtCusp+\aboveMin \hNrm_{\EndStg}\right)\MPCmin_{\prd} & = \MPCmax_{\prd} \aboveMin \mtCusp \\
\aboveMin \mtCusp & = \frac{\MPCmin_{\prd}\aboveMin \hNrm_{\EndStg}}{(1-\MPCmin_{\prd})\MPCmax_{\prd}} \\
\mtCusp & = \frac{\MPCmin_{\prd}\hNrm_{\EndStg}-\hEndMin_{\EndStg}}{(1-\MPCmin_{\prd})\MPCmax_{\prd}},
\end{aligned}\end{gathered}\end{equation*}
\end{verbatimwrite}
\input{./Equations/mtCusp.tex}\unskip
\begin{verbatimwrite}{./Equations/TighterUpperBound.tex}
we want to construct a consumption function for $m_{\prd} \in (\Min{m}_{\prd}, \mtCusp]$ that respects the
tighter upper bound:
\begin{center}
\begin{tabular}{rcl}
$ \aboveMin \mNrm_{\prd} \MPCmin_{\prd} < $ & $ \cFunc_{\prd}(\Min{m}_{\prd}+\aboveMin \mNrm_{\prd}) $ $< \MPCmax_{\prd} \aboveMin \mNrm_{\prd} $
% \\ $-\aboveMin \mNrm_{\prd} \MPCmin_{\prd} > $ & $ -\cFunc_{\prd}(\Min{m}_{\prd}+\aboveMin \mNrm_{\prd}) $ & $> -\aboveMin \mNrm_{\prd} $
\\ $ \aboveMin \mNrm_{\prd}(\MPCmax_{\prd}- \MPCmin_{\prd}) > $ & $ \MPCmax_{\prd} \aboveMin \mNrm_{\prd}-\cFunc_{\prd}(\Min{m}_{\prd}+\aboveMin \mNrm_{\prd}) $ & $> 0$
\\ $1 > $ & $ \left(\frac{\MPCmax_{\prd} \aboveMin \mNrm_{\prd}-\cFunc_{\prd}(\Min{m}_{\prd}+\aboveMin \mNrm_{\prd})}{\aboveMin \mNrm_{\prd}(\MPCmax_{\prd}- \MPCmin_{\prd})}\right) $ & $> 0$.
\end{tabular}
\end{center}
\end{verbatimwrite}
\input{./Equations/TighterUpperBound.tex}\unskip
\begin{verbatimwrite}{./Equations/koppaLo.tex}
Again defining $\mu_{\prd} =\log \aboveMin \mNrm_{\prd}$, the object in the middle of the inequality is
\begin{equation*}\begin{gathered}\begin{aligned}
\Min{\koppa}_{\prd}(\mu_{\prd}) & \equiv \frac{\MPCmax_{\prd}-\cFunc_{\prd}(\Min{m}_{\prd}+e^{\mu_{\prd}})e^{-\mu_{\prd}}}{\MPCmax_{\prd}-\MPCmin_{\prd}} \label{eq:koppaL}
\MPCMatch{\\ \Min{\koppa}^{\mu}_{\prd}(\mu_{\prd}) & = \frac{\cFunc_{\prd}(\Min{m}_{\prd}+e^{\mu_{\prd}})e^{-\mu_{\prd}}-\MPCFunc_{\prd}^{m}(\Min{m}_{\prd}+e^{\mu_{\prd}})}{\MPCmax_{\prd}-\MPCmin_{\prd}}}{} .
\end{aligned}\end{gathered}\end{equation*}
\end{verbatimwrite}
\input{./Equations/koppaLo.tex}\unskip
\begin{verbatimwrite}{./cctwMoM/TighterThreeFuncs.tex}
As $m_{\prd}$ approaches
$-\Min{m}_{\prd}$, $\Min{\koppa}_{\prd}(\mu_{\prd})$ converges to zero, while as $m_{\prd}$
approaches $+\infty$, $\Min{\koppa}_{\prd}(\mu_{\prd})$ approaches $1$.
As before, we can derive an approximated consumption function; call it $\Aprx{\Min{\cFunc}}_{\prd}$. This function will clearly do a better job approximating the consumption function for low values of $\mNrm_{\prd}$ while the previous approximation will perform better for high values of $\mNrm_{\prd}$.
For middling values of $\mNrm$ it is not clear which of these functions will perform better. However, an alternative is available which performs well. Define the highest gridpoint below $\mtCusp$ as $\bar{\check{\mNrm}}_{\prd}^{\#}$ and the lowest gridpoint above $\mtCusp$ as $\Min{\hat{\mNrm}}_{\prd}^{\#}$. Then there will be a unique interpolating polynomial that matches the level and slope of the consumption function at these two points. Call this function $\tilde{\cFunc}_{\prd}(\mNrm)$.
Using indicator functions that are zero everywhere except for specified intervals,
\end{verbatimwrite}
\input{./cctwMoM/TighterThreeFuncs.tex}\unskip
\begin{verbatimwrite}{./Equations/TighterThreeEqns.tex}
\begin{equation*}\begin{gathered}\begin{aligned}
\vctr{1}_{\text{Lo}}(\mNrm) & = 1 \text{~if $ \mNrm \leq \bar{\check{\mNrm}}_{\prd}^{\#} \phantom{< \mNrm < \Min{\hat{\mNrm}}_{\prd}^{\#} \leq \mNrm}$}
\\ \vctr{1}_{\text{Mid}}(\mNrm) & = 1 \text{~if $\phantom{ \mNrm \leq}~ \bar{\check{\mNrm}}_{\prd}^{\#} < \mNrm < \Min{\hat{\mNrm}}_{\prd}^{\#} \phantom{\leq \mNrm}$}
\\ \vctr{1}_{\text{Hi}}(\mNrm) & = 1 \text{~if $\phantom{ \mNrm \leq ~\bar{\check{\mNrm}}_{\prd}^{\#} < \mNrm < } \Min{\hat{\mNrm}}_{\prd}^{\#} \leq \mNrm$}
\end{aligned}\end{gathered}\end{equation*}
we can define a well-behaved approximating consumption function
\begin{equation}\begin{gathered}\begin{aligned}
\Aprx{\cFunc}_{\prd} & = \vctr{1}_{\text{Lo}} \Aprx{\Min{\cFunc}}_{\prd} + \vctr{1}_{\text{Mid}} \Aprx{\tilde{\cFunc}}_{\prd}+\vctr{1}_{\text{Hi}} \Aprx{\Max{\cFunc}}_{\prd}.
\end{aligned}\end{gathered}\end{equation}
\end{verbatimwrite}
\input{./Equations/TighterThreeEqns.tex}\unskip
\begin{verbatimwrite}{./cctwMoM/TighterThreeFuncsExplain.tex}
This just says that, for each interval, we use the approximation that
is most appropriate. The function is continuous and
once-differentiable everywhere, and is therefore well behaved for
computational purposes.
\begin{comment}
In practice, in our problem the difference due to this refinement is displayed in Figure \ref{fig:IntExpFOCInvPesReaOpt45GapPlot}.
\hypertarget{IntExpFOCInvPesReaOpt45GapPlot}{}
\begin{figure}
\includegraphics[width=6in]{./Figures/IntExpFOCInvPesReaOpt45GapPlot}
\caption{Difference Between $\Aprx{\Max{\cFunc}}_{L, T-1}$ and $\Aprx{\Max{\cFunc}}_{H,T-1}$ is Small}
\label{fig:IntExpFOCInvPesReaOpt45GapPlot}
\end{figure}
\end{comment}
We now construct an upper-bound value function implied for a consumer whose spending behavior is consistent with the refined upper-bound consumption rule.
For $\mNrm_{\prd} \geq \mNrm_{\prd}^{\#}$, this consumption rule is the same as before,
so the constructed upper-bound value function is also the same. However, for
values $\mNrm_{\prd} < \mNrm_{\prd}^{\#}$ matters are slightly more complicated.
Start with the fact that at the cusp point,
\begin{equation*}\begin{gathered}\begin{aligned}
\bar{\vFunc}_{\prd}(\mtCusp) & = \uFunc(\bar{\cNrm}_{\prd}(\mtCusp))\PDVCoverc_{\prd}^T \\
& = \uFunc(\aboveMin \mtCusp \MPCmax_{\prd})\PDVCoverc_{\prd}^{T}
.
\end{aligned}\end{gathered}\end{equation*}
But for \textit{all} $\mNrm_{\prd}$,
\begin{equation*}\begin{gathered}\begin{aligned}
\bar{\vFunc}_{\prd}(\mNrm) & = \uFunc(\bar{\cNrm}_{\prd}(\mNrm))+ \bar{\EndPrd}(\mNrm-\bar{\cNrm}_{\prd}(\mNrm)),
\end{aligned}\end{gathered}\end{equation*}
and we assume that for the consumer below the cusp point consumption is given by $\MPCmax \aboveMin \mNrm_{\prd}$ so for $\mNrm_{\prd}< \mtCusp$
\begin{equation*}\begin{gathered}\begin{aligned}
\bar{\vFunc}_{\prd}(\mNrm) & = \uFunc( \MPCmax_{\prd} \aboveMin \mNrm)+ \bar{\EndPrd}((1-\MPCmax_{\prd})\aboveMin \mNrm),
\end{aligned}\end{gathered}\end{equation*}
which is easy to compute because $\bar{\EndPrd}(\aNrm_{\prd}) = \DiscFac \bar{\vFunc}_{\prd+1}(\aNrm_{\prd}\RNrmByG+1)$ where $\bar{\vFunc}_{\prd}$ is as defined above because a consumer who ends the current period with assets exceeding the lower bound will not expect to be constrained next period. (Recall again that we are merely constructing an object that is guaranteed to be an \textit{upper bound} for the value that the `realist' consumer will experience.) At the gridpoints defined by the solution of the consumption problem can then construct
\begin{equation*}\begin{gathered}\begin{aligned}
\bar{\vInv}_{\prd}(\mNrm) & = ((1-\CRRA)\bar{\vFunc}_{\prd}(\mNrm))^{1/(1-\CRRA)}
\end{aligned}\end{gathered}\end{equation*}
\MPCMatch{and its derivatives}{} which yields the appropriate vector for constructing $\check{\Chi}$ and $\check{\Koppa}$. The rest of the procedure is analogous to that performed for the consumption rule and is thus omitted for brevity.
\end{verbatimwrite}
\input{./cctwMoM/TighterThreeFuncsExplain.tex}\unskip
\hypertarget{extension-a-stochastic-interest-factor}{}
\subsection{Extension: A Stochastic Interest Factor}
Thus far we have assumed that the interest factor is constant at $\Rfree$. Extending the
previous derivations to allow for a perfectly forecastable time-varying interest factor $\Rfree_{\prd}$
would be trivial. Allowing for a stochastic interest factor is less trivial.
The easiest case is where the interest factor is i.i.d.,
\begin{verbatimwrite}{./Equations/distRisky.tex}
\begin{equation}\begin{gathered}\begin{aligned}
\log \Risky_{t+n} & \sim & \Nrml(\rfree + \eprem - \sigma^{2}_{\risky}/2,\sigma^{2}_{\risky}) ~\forall~n>0 \label{eq:distRisky}
\end{aligned}\end{gathered}\end{equation}
\end{verbatimwrite}
\input{./Equations/distRisky.tex}\unskip
where $\eprem$ is the risk premium and the $\sigma^{2}_{\risky}/2$ adjustment to the mean log return
guarantees that an increase in $\sigma^{2}_{\risky}$ constitutes a mean-preserving spread in the level of the return.
This case is reasonably straightforward because \cite{merton:restat} and \cite{samuelson:portfolio} showed
that for a consumer without labor income (or with perfectly forecastable labor income) the consumption
function is linear, with an infinite-horizon MPC\footnote{See \handoutC{CRRA-RateRisk} for a derivation.}
\begin{equation}\begin{gathered}\begin{aligned}
\MPC & = 1- \left(\DiscFac \Ex_{\BegStg}[\Risky_{\prd+1}^{1-\CRRA}]\right)^{1/\CRRA} \label{eq:MPCExact}
\end{aligned}\end{gathered}\end{equation}
and in this case the previous analysis applies once we substitute this MPC for the one that characterizes
the perfect foresight problem without rate-of-return risk.
The more realistic case where the interest factor has some serial correlation is more complex. We consider
the simplest case that captures the main features of empirical interest rate dynamics: An AR(1) process. Thus
the specification is
\begin{equation}\begin{gathered}\begin{aligned}
\risky_{\prd+1}-\risky & = (\risky_{\prd}-\risky) \gamma + \epsilon_{\prd+1}
\end{aligned}\end{gathered}\end{equation}
where $\risky$ is the long-run mean log interest factor, $0 < \gamma < 1$ is the AR(1) serial correlation
coefficient, and $\epsilon_{\prd+1}$ is the stochastic shock.
The consumer's problem in this case now has two state variables, $\mNrm_{\prd}$ and $\risky_{\prd}$, and
is described by
\begin{equation}\begin{gathered}\begin{aligned}
\vFunc_{\prd}(m_{\prd},\risky_{\prd}) & = \max_{{c}_{\prd}} ~ \uFunc(c_{\prd})+
\Ex_{\BegStg}[{\DiscFac}_{\prd+1}\PermGroFacAdjV{\vFunc}_{\prd+1}(m_{\prd+1},\risky_{\prd+1})] \label{vNormedRisky}
\\ & \text{s.t.} \nonumber \\
a_{\prd} & = m_{\prd}-c_{\prd} \nonumber
\\ \risky_{\prd+1}-\risky & = (\risky_{\prd}-\risky)\gamma + \epsilon_{\prd+1} \notag
\\ \Risky_{\prd+1} & = \exp(\risky_{\prd+1}) \notag
\\ m_{\prd+1} & = \underbrace{\left(\Risky_{\prd+1}/\PermGroFac_{\prd+1}\right)}_{\equiv \Rprod_{\prd+1}}a_{\prd}+\tranShkEmp_{\prd+1} \nonumber.
\end{aligned}\end{gathered}\end{equation}
% Kiichi: I will need you to read the literature and figure out how exactly we want to choose the Markov points and transition probabilities.
% When done, you will fill in the [how] text below.
We approximate the AR(1) process by a Markov transition matrix using standard techniques. The stochastic interest factor is allowed to take
on 11 values centered around the steady-state value $\risky$. Given this Markov transition matrix, \textit{conditional} on the Markov AR(1) state the consumption functions for the `optimist' and the `pessimist' will still be linear,
with identical MPC's that are computed numerically. Given these MPC's, the (conditional) realist's consumption function can be computed for each Markov state, and the converged consumption rules constitute the solution contingent on the dynamics of the stochastic
interest rate process.
In principle, this refinement should be combined with the previous one;
further exposition of this combination is omitted here because no new
insights spring from the combination of the two techniques.
\hypertarget{imposing-artificial-borrowing-constraints}{}
\subsection{Imposing `Artificial' Borrowing Constraints}
Optimization problems often come with additional constraints that must
be satisfied. Particularly common is an `artificial' liquidity constraint that
prevents the consumer's net worth from falling below some value, often
zero.\footnote{The word artificial is chosen only because of its clarity in distinguishing
this from the case of the `natural' borrowing constraint examined above; no derogation is
intended -- constraints of this kind certainly exist in the real world.} The problem then becomes
\begin{equation*}\begin{gathered}\begin{aligned}
\vFunc_{\prd-1}(m_{\prd-1}) & = \max_{\cNrm_{\prd-1}} ~~ \uFunc(c_{\prd-1}) + \Ex_{\prd-1} [\DiscFac \PermGroFacAdjV{\vFunc}_{\cntn(T)}(m_{\prd})] \label{eq:ConstrArt}
\\ & \mbox{s.t.} \nonumber
\\ a_{\prd-1} & = m_{\prd-1} - c_{\prd-1}
\\ m_{\prd} & = \RNrmByG_{\prd} a_{\prd-1} + \tranShkEmp_{\prd}
\\ a_{\prd-1} & \geq 0 .
\end{aligned}\end{gathered}\end{equation*}
By definition, the constraint will bind if the unconstrained consumer
would choose a level of spending that would violate the constraint.
Here, that means that the constraint binds if the $c_{\prd-1}$
that satisfies the unconstrained FOC
\begin{equation}\begin{gathered}\begin{aligned}
c_{\prd-1}^{-\CRRA} & = \vFunc^{a}_{({\prd-1})_\cntn}(m_{\prd-1}-c_{\prd-1}) \label{eq:cUnc}
\end{aligned}\end{gathered}\end{equation}
is greater than $m_{\prd-1}$. Call $\grave{\cFunc}^{\ast}_{\prd-1}$ the approximated function
returning the level of $c_{\prd-1}$ that satisfies \eqref{eq:cUnc}.
Then the approximated constrained optimal consumption function will be
\begin{verbatimwrite}{./Equations/LiqCons.tex}
\begin{equation}\begin{gathered}\begin{aligned}
\grave{\cFunc}_{\prd-1}(m_{\prd-1}) & = \min[{m}_{\prd-1},\grave{\cFunc}^{\ast}_{\prd-1}(m_{\prd-1})] \label{eq:LiqCons}.
\end{aligned}\end{gathered}\end{equation}
\end{verbatimwrite}
\input{./Equations/LiqCons.tex}\unskip
The introduction of the constraint also introduces a sharp
nonlinearity in all of the functions at the point where the constraint
begins to bind. As a result, to get solutions that are anywhere close
to numerically accurate it is useful to augment the grid of values of
the state variable to include the exact value at which the constraint
ceases to bind. Fortunately, this is easy to calculate. We know that
when the constraint is binding the consumer is saving nothing, which
yields marginal value of $\vFunc^{a}_{({\prd-1})_\cntn}(0)$. Further, when the
constraint is binding, $c_{\prd-1} = m_{\prd-1}$. Thus, the largest
value of consumption for which the constraint is binding will be the
point for which the marginal utility of consumption is exactly equal
to the (expected, discounted) marginal value of saving 0. We know
this because the marginal utility of consumption is a downward-sloping
function and so if the consumer were to consume $\tinyAmount$ more,
the marginal utility of that extra consumption would be \textit{below}
the (discounted, expected) marginal utility of saving, and thus the
consumer would engage in positive saving and the constraint would no
longer be binding. Thus the level of $m_{\prd-1}$ at which the
constraint stops binding is:\footnote{The logic here repeats an insight from \cite{deatonLiqConstr}.}
\begin{equation}\begin{gathered}\begin{aligned}
\uFunc^{c}(m_{\prd-1}) & = \vFunc^{a}_{({\prd-1})_\cntn}(0) \nonumber \\
m_{\prd-1} & = (\vFunc^{a}_{({\prd-1})_\cntn}(0))^{(-1/\CRRA)} \nonumber
\\ & = \cFunc_{({\prd-1})_\cntn}(0). \label{eq:LCbindsTm1}
\end{aligned}\end{gathered}\end{equation}
\hypertarget{cVScCon}{}
\begin{figure}
\includegraphics[width=6in]{./Figures/cVScCon}
\caption{Constrained (solid) and Unconstrained (dashed) Consumption}
\label{fig:cVScCon}
\end{figure}
The constrained problem is solved in section ``Artifical Borrowing Constraint''
of the notebook, where the variable
\texttt{constrained} is set to be a boolean type object. If the value of \texttt{constrained}
is true, then the constraint is binding and their consumption behavior is computed to match
\eqref{eq:LiqCons}. The resulting consumption rule is shown in Figure \ref{fig:cVScCon}. For comparison purposes,
the approximate consumption rule from Figure \ref{fig:cVScCon} is
reproduced here as the solid line; this is accomplished by setting the boolean value
of \texttt{constrained} to false.
The presence of the liquidity
constraint requires three changes to the procedures outlined above:
\begin{enumerate}
\item We redefine
$\hEndMin_{\EndStg}$, which now is the PDV of receiving
$\tranShkEmp_{\prd+1}=\tranShkEmpMin$ next period and
$\tranShkEmp_{t+n}=0~\forall~n>1$ -- that is, the pessimist believes he
will receive nothing beyond period $t+1$
\item We augment the end-of-period \code{aVec} with zero and with a point with a small positive value so that the generated
{\mVec} will the binding point $\mNrm^{\#}$ and a point just above it (so that we can better capture the curvature
around that point)
\item We redefine the optimal consumption rule as
in equation (\ref{eq:LiqCons}). This ensures that the
liquidity-constrained `realist' will consume more than the redefined
`pessimist,' so that we will have $\koppa$ still between $0$ and $1$
and the `method of moderation' will proceed smoothly.
\end{enumerate}
As expected, the liquidity constraint only causes a divergence between the two functions at the point where the optimal unconstrained consumption rule runs into the 45 degree line.
\hypertarget{recursion}{}
\section{Recursion}\label{sec:recursion}
\hypertarget{theory}{}
\subsection{Theory}
Before we solve for periods earlier than $\trmT-1$, we assume for
convenience that in each such period a liquidity constraint exists of
the kind discussed above, preventing $c$ from exceeding $m$. This
simplifies things a bit because now we can always consider an
\code{aVec} that starts with zero as its smallest element.
Recall now equations~(\ref{eq:vEndPrimeTm1}) and (\ref{eq:upEqbetaOp}):
\begin{equation*}\begin{gathered}\begin{aligned}
\vPEndStg(a_{\prd}) & = \Ex_{\BegStg}[\DiscFac \Rfree \PermGroFac_{\prd+1}^{-\CRRA}
\uFunc^{c}(\cFunc_{\prd+1}(\RNrmByG_{\prd+1} a_{\prd}+{\tranShkEmp}_{\prd+1}))]
\\\uFunc^{c}(c_{\prd}) & = \vEndStg^{a}(m_{\prd}-c_{\prd}).
\end{aligned}\end{gathered}\end{equation*}
Assuming that the problem has been solved up to period $t+1$ (and thus assuming that we have an approximated $\Aprx{\cFunc}_{\prd+1}(m_{\prd+1})$), our solution method essentially involves using these two equations in succession to work back progressively from period $\trmT-1$ to the beginning of life. Stated generally, the method is as follows. (Here, we use the original, rather than the ``refined,'' method for constructing consumption functions; the generalization of the algorithm below to use the refined method presents no difficulties.)
\begin{enumerate}
\item For the grid of values $a_{t,i}$ in \texttt{aVec\_eee}, numerically calculate the values
of $\cFunc_{\Max{t}}(a_{t,i})$ and $\cFunc_{\Max{t}}^{a}(a_{t,i})$,
% \begin{verbatimwrite}{./Equations/vEndeq.tex}
\begin{equation}\begin{gathered}\begin{aligned}
\cFunc_{\Max{t},i} & = \left(\vEndStg^{a}(a_{t,i})\right)^{-1/\CRRA},
\\ & = \left(\DiscFac \Ex_{\BegStg} \left[\Rfree \PermGroFac_{\prd+1}^{-\CRRA}(\grave{\cFunc}_{\prd+1}(\RNrmByG_{\prd+1} a_{t,i} + {\tranShkEmp}_{\prd+1}))^{-\CRRA}\right]\right)^{-1/\CRRA}, \label{eq:vEndeq}
\MPCMatch{\\ \cFunc^{a}_{\Max{t},i} & = -(1/\CRRA)\left(\vEndStg^{a}(a_{t,i})\right)^{-1-1/\CRRA} \vEndStg^{a{a}}(\aNrm_{t,i}),}{}
\end{aligned}\end{gathered}\end{equation}
% \end{verbatimwrite}
% \input{./Equations/vEndeq.tex} \unskip
generating vectors of values $\vctr{\cFunc}_{\prd}$\MPCMatch{ and $\vctr{\cFunc}^{a}_{\Max{t}}$.}{.}
\item Construct a corresponding vector of values of $\vctr{m}_{\prd}=\vctr{\cNrm}_{\prd}+\vctr{\aNrm}_{\prd}$\MPCMatch{; similarly construct a corresponding list of MPC's $\vctr{\MPC}_{\prd}$ using equation \eqref{eq:MPCfromMPTHC}.}{.}
\item Construct a corresponding vector $\vctr{\mu_{\prd}}$, the levels\MPCMatch{ and first derivatives}{} of $\vctr{\koppa}_{\prd}$, and the levels\MPCMatch{ and first derivatives}{} of $\vctr{\chi}_{\prd}$.
\item Construct an interpolating approximation $\Aprx{\chi}_{\prd}$ that\MPCMatch{ smoothly matches both the level and the slope}{the level} at those points.
\item If we are to approximate the value function, construct a corresponding list of values of $\vctr{v}_{\prd}$, the levels\MPCMatch{ and first derivatives of $\vctr{\Koppa}_{\prd}$,}{,} and the levels\MPCMatch{ and first derivatives}{} of $\hat{\vctr{\Chi}}_{\prd}$; and construct an interpolating approximation function $\hat{\Chi}_{\prd}$ that matches those points.
\end{enumerate}
With $\Aprx{\chi}_{\prd}$ in hand, our approximate consumption function
is computed directly from the appropriate substitutions in \eqref{eq:cFuncHi}
and related equations. With this consumption
rule in hand, we can continue the backwards recursion to period $t-1$
and so on back to the beginning of life.
Note that this loop does not contain an item for constructing $\hat{\vFunc}_{\prd}^{a}(m_{\prd})$. This is because with $\Aprx{\Max{\cFunc}}_{\prd}(m_{\prd})$ in hand, we simply \textit{define} $\hat{\vFunc}^{m}_{\prd}(m_{\prd}) = \uFunc^{c}(\Aprx{\Max{\cFunc}}_{\prd}(m_{\prd}))$ so there is no need to construct interpolating approximations - the function arises `free' (or nearly so) from our constructed $\Aprx{\Max{\cFunc}}_{\prd}(m_{\prd})$ via the usual envelope result (cf.\ \eqref{eq:envelope}).
\end{document}