Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Portfolio Choice with CRRA Utility (Merton-Samuelson)

Authors
Affiliations
Johns Hopkins University
Econ-ARK
Johns Hopkins University
Econ-ARK

Merton (1969) and Samuelson (1969) study the optimal portfolio choice of a consumer with constant relative risk aversion ρ\CRRA.[1] This consumer has assets at the end of period tt equal to ata_{t} and is deciding how much to invest in a risky asset[2] with a lognormally distributed return factor Rt+1\Risky_{t+1}, whose log is

rt+1=r+θt+1,θt+1∼N(−σr2/2,σr2).\risky_{t+1} = \risky+\ShkMeanOneLog_{t+1}, \qquad \ShkMeanOneLog_{t+1} \sim \mathcal{N}(-\sigma^{2}_{\risky}/2,\sigma^{2}_{\risky}).

The shock is mean-one in levels, Et[eθt+1]=1\Ex_{t}[e^{\ShkMeanOneLog_{t+1}}]=1 (see ELogNormMeanOne), so r\risky is the log of the expected return factor, Et[Rt+1]=er=R\Ex_{t}[\Risky_{t+1}]=e^{\risky}=\Risky (upper-case variables like R\Risky without a subscript are time-invariant means), while the mean of the log return is Et[rt+1]=r−σr2/2\Ex_{t}[\risky_{t+1}]=\risky-\sigma^{2}_{\risky}/2.

The alternative to the risky asset is a riskfree asset that earns return factor R=er\Rfree=e^{\rfree}. Importantly, the consumer is assumed to have no labor income and to face no risk except from the investment in the risky asset.[3][4]

Both papers consider a multiperiod optimization problem, but here we examine a consumer for whom period tt is the second-to-last period of life (the insights, and even the formulas, carry over to the multiperiod case).[5]

If the period-tt consumer invests proportion ς\riskyshare in the risky asset, spending all available resources in the last period of life t+1t+1 will yield:

ct+1=(R(1−ς)+Rt+1ς)at=(R+(Rt+1−R)ς)⏟≡R~t+1at\begin{gathered}\begin{aligned} c_{t+1} & = \left(\Rfree(1-\riskyshare)+\Risky_{t+1}\riskyshare\right)a_{t} \\ & = \underbrace{\left(\Rfree+(\Risky_{t+1}-\Rfree)\riskyshare\right)}_{\equiv \Rport_{t+1}}a_{t} \end{aligned}\end{gathered}

where R~t+1\Rport_{t+1} is the realized arithmetic return factor for the portfolio.

The optimal portfolio share will be the one that maximizes expected utility:

ς=arg max⁡ςEt[u(ct+1)]\begin{gathered}\begin{aligned} \varsigma & = \argmax_{\varsigma} \Ex_{t}[\uFunc(c_{t+1})] \end{aligned}\end{gathered}

and can be calculated numerically for any distribution of returns.

Campbell & Viceira (2002) point out that if we define the realized log excess return

φt+1=rt+1−r\begin{gathered}\begin{aligned} \EpremLog_{t+1} & = \risky_{t+1}-\rfree \end{aligned}\end{gathered}

whose mean is φ−σr2/2\EpremLog-\sigma^{2}_{\risky}/2, where φ≡log⁡φ=r−r\EpremLog \equiv \log \EPrem = \risky-\rfree is the log of the expected excess return factor φ≡R/R\EPrem \equiv \Risky/\Rfree, then for many distributions a good approximation to the log of the portfolio’s return factor is obtained by[6]

r~t+1=r+ςφt+1+ςσr2/2−ς2σr2/2\begin{gathered}\begin{aligned} \rport_{t+1} & = \rfree+ \riskyshare \EpremLog_{t+1}+\riskyshare \sigma^{2}_{\risky}/2 - \riskyshare^{2}\sigma^{2}_{\risky}/2 \end{aligned}\end{gathered}

Using this approximation, the expectation as of date tt of utility at date t+1t+1 is:

Et[u(ct+1)]≈(1−ρ)−1Et[(atereςφt+1+ς(1−ς)σr2/2)1−ρ]=(1−ρ)−1Et[(atR)1−ρ(eςφt+1+ς(1−ς)σr2/2)1−ρ]=(1−ρ)−1(atR)1−ρEt[e(ςφt+1+ς(1−ς)σr2/2)(1−ρ)]=(1−ρ)−1(atR)1−ρ⏟constant <0e(1−ρ)ς(1−ς)σr2/2Et[eςφt+1(1−ρ)]⏟excess return utility factor\begin{split} \Ex_{t}[\uFunc(c_{t+1})] & \approx (1-\CRRA)^{-1}\Ex_{t}\left[\left(a_{t}e^{\rfree}e^{\riskyshare \EpremLog_{t+1}+\riskyshare(1-\riskyshare)\sigma^{2}_{\risky}/2 }\right)^{1-\CRRA}\right] \\ & = (1-\CRRA)^{-1}\Ex_{t}\left[(a_{t}\Rfree)^{1-\CRRA}\left( e^{\riskyshare \EpremLog_{t+1}+\riskyshare(1-\riskyshare)\sigma^{2}_{\risky}/2 }\right)^{1-\CRRA}\right] \\ & = (1-\CRRA)^{-1}(a_{t}\Rfree)^{1-\CRRA}\Ex_{t}\left[e^{(\riskyshare \EpremLog_{t+1}+\riskyshare(1-\riskyshare)\sigma^{2}_{\risky}/2) (1-\CRRA)}\right] \\ & = \underbrace{(1-\CRRA)^{-1}(a_{t}\Rfree)^{1-\CRRA}}_{\text{constant $< 0$}}\underbrace{e^{ (1-\CRRA)\riskyshare(1-\riskyshare)\sigma^{2}_{\risky}/2}\Ex_{t}\left[e^{\riskyshare \EpremLog_{t+1} (1-\CRRA)}\right]}_{\text{excess return utility factor}} \end{split}

where the first term is a negative constant under the usual assumption that relative risk aversion ρ>1\CRRA>1.

Because the return is lognormal, we can make substantial further progress, by obtaining an analytical approximation to the numerical optimum. In this case ς(1−ρ)φt+1∼N(ς(1−ρ)(φ−σr2/2),(ς(1−ρ))2σr2)\riskyshare (1-\CRRA) \EpremLog_{t+1} \sim \mathcal{N}(\riskyshare (1-\CRRA)(\EpremLog - \Evarr/2),(\riskyshare(1-\CRRA))^{2}\Evarr) (using NormTimes). With a couple of extra lines of derivation we can show that the log of the expectation in (6) is

log⁡Et[eςφt+1(1−ρ)]=(1−ρ)ςφ−(1−ρ)ςσr2/2+((1−ρ)ς)2σr2/2=(1−ρ)ςφ−(1−ρ)ς(1−ς(1−ρ))σr2/2=(1−ρ)ςφ−(1−ρ)ς(1−ς)σr2/2−ρ(1−ρ)ς2σr2/2\begin{split} \log \Ex_{t}\left[e^{\riskyshare \EpremLog_{t+1} (1-\CRRA)}\right] & = {(1-\CRRA)\riskyshare \EpremLog-(1-\CRRA)\riskyshare\Evarr/2+ ((1-\CRRA)\riskyshare)^{2}\Evarr/2} \\ & = {(1-\CRRA)\riskyshare \EpremLog-(1-\CRRA)\riskyshare(1-\riskyshare(1-\CRRA))\Evarr/2} \\ & = {(1-\CRRA)\riskyshare \EpremLog-(1-\CRRA)\riskyshare(1-\riskyshare)\Evarr/2-\CRRA (1-\CRRA)\riskyshare^{2}\Evarr/2} \end{split}

Substitute from (7) for the log of the expectation in (6) and note that the resulting expression simplifies because it contains (1−ρ)ς(1−ς)σr2/2−(1−ρ)ς(1−ς)σr2/2=0{(1-\CRRA)\riskyshare(1-\riskyshare)\Evarr/2-(1-\CRRA)\riskyshare(1-\riskyshare)\Evarr/2}=0; thus the log of the “excess return utility factor” in (6) is

−(ρ−1)ςφ−(ρ−1)(−ρς2σr2/2)-(\CRRA-1)\riskyshare \EpremLog - (\CRRA-1)(- \CRRA \riskyshare^{2}\Evarr/2)

and the ς\riskyshare that minimizes the log will also minimize the level; minimizing this when ρ>1\CRRA>1 is equivalent to maximizing the terms multiplied by −(ρ−1)-(\CRRA-1), so our problem reduces to

max⁡ς  ςφ−ρς2σr2/2\begin{gathered}\begin{aligned} \max_{\riskyshare}~~ \riskyshare \EpremLog -\CRRA\riskyshare^{2}\Evarr/2 \end{aligned}\end{gathered}

with FOC

φ−ςρσr2=0ς=(φρσr2)\begin{gathered}\begin{aligned} \EpremLog-\riskyshare\CRRA\Evarr & = 0 \\ \riskyshare & = \left(\frac{\EpremLog}{\CRRA \Evarr}\right) \end{aligned}\end{gathered}

Equation (10) says[7] that the consumer allocates a higher proportion of their net worth to the high-risk, high-return asset when

  1. the amount φ\EpremLog by which the risky asset’s return exceeds the riskless return is greater

  2. the consumer is less risk averse (ρ\CRRA is lower)

  3. riskiness σr2\sigma^{2}_{\risky} is less

If there is no excess return, nothing will be put in the risky asset. Similarly, if risk aversion or the variance of the risk is infinity, again nothing will be put in the risky asset.[8]

This formula hints at the existence of an “equity premium puzzle” (Mehra & Prescott (1985)). Interpreting the risky asset as the aggregate stock market, the annual standard deviation of the log of U.S. stock returns has historically been about σr=0.2\sigma_{\risky}=0.2 yielding σr2=0.04\Evarr = 0.04. The equity premium over historical periods has been something like φ=0.08\EpremLog = 0.08 (eight percent). With risk aversion of ρ=2\CRRA=2 this formula implies that the share of risky assets in your portfolio should be 0.08/0.080.08/0.08 or 100 percent! The fact that most people hold far less than 100 percent of their wealth in stocks, and that many hold none, is the “stockholding puzzle,” the microeconomic manifestation of the equity premium puzzle (Haliassos & Bertaut (1995)).

To avoid the problems caused by a prediction of a risky portfolio share greater than one, we can calibrate the model with more modest expectations for the equity premium. Taking evidence from other countries and longer periods into account, Dimson et al. (2008) estimate that investors can expect a premium of about 3 to 3½ percent on a geometric basis, or 4½ to 5 percent on the arithmetic basis that matches φ\EpremLog. Figure 1 shows the relationship between the portfolio share and relative risk aversion for a deliberately modest premium of 3 percent and a large standard deviation of σ=0.2\sigma=0.2. Even when risks are this high and the premium is this low, if relative risk aversion is close to logarithmic (ρ=1\CRRA = 1) the investor wants to put well over half of the portfolio in the risky asset. Only for values of risk aversion greater than 2 does the predicted portfolio share reach plausible small values.

But remember that these calculations are all assuming that the consumer’s entire consumption spending is financed by asset income. If the consumer has other income (for example, labor or transfer income) whose covariance with the risky return is small, they should be willing to take more risk. Since, for most consumers, most of their future consumption will be financed from labor or transfer income, it is not surprising to learn that models calibrated to actual data on capital and noncapital income dynamics imply that people should be investing most of their non-human wealth in the risky asset, with reasonable values of ρ\CRRA Cocco et al., 2005.

Another question is what the expected rate of return on the consumer’s portfolio will be once the portfolio share in risky assets has been chosen optimally. Note first that (16) implies that

log⁡Et[er~t+1−r]=ςφ\begin{gathered}\begin{aligned} \log \Ex_{t}[e^{\rport_{t+1}-\rfree}] & = \riskyshare \EpremLog \end{aligned}\end{gathered}

while the variance of the log of the excess return factor for the portfolio is σr~2=ς2σr2\sigma^{2}_{\rport} = \riskyshare^{2} \sigma^{2}_{\risky}. Substituting the solution (10) for ς\riskyshare into (11), we have

ςφ=(φ2ρσr2)=(φ/σr)2/ρ\begin{gathered}\begin{aligned} \riskyshare \EpremLog & = \left(\frac{\EpremLog^{2}}{\CRRA \Evarr}\right) \\ & = (\EpremLog/\sigma_{\risky})^{2}/\CRRA \end{aligned}\end{gathered}

which is an interesting formula for the excess return of the optimally chosen portfolio because the object φ/σr\EpremLog/\sigma_{\risky} (the excess return divided by the standard deviation) is a well-known tool in finance for evaluating the tradeoff between risk and return (the “Sharpe ratio”). Equation (12) says that the consumer will choose a portfolio that earns an excess return proportional to the square of the Sharpe ratio and inversely proportional to the risk aversion coefficient. Higher reward (per unit of risk) convinces the consumer to take the risk necessary to earn higher returns; but higher risk aversion convinces the investor to sacrifice (risky) return for safety.

Finally, we can ask what effect an exogenous increase in the risk of the risky asset has on the endogenous riskiness of the portfolio once the consumer has chosen optimally. The answer is surprising: The variance of the optimally-chosen portfolio is

ς2σr2=(φρσr2)2σr2=((φ/ρ)2σr2)\begin{gathered}\begin{aligned} \riskyshare^{2} \sigma^{2}_{\risky} & = \left(\frac{\EpremLog}{\CRRA \Evarr}\right)^{2} \sigma^{2}_{\risky} \\ & = \left(\frac{(\EpremLog/\CRRA)^{2}}{\Evarr}\right) \end{aligned}\end{gathered}

which is actually smaller when σr2\sigma^{2}_{\risky} is larger. The reason: a riskier asset pays less per unit of variance. Cutting the share only enough to restore the portfolio’s variance would leave the same risk at a lower expected return, so the investor cuts it further.

The Approximate Risky Portfolio Share \riskyshare Declines as Relative Risk Aversion \CRRA Increases

Figure 1:The Approximate Risky Portfolio Share ς\riskyshare Declines as Relative Risk Aversion ρ\CRRA Increases

Note: The approximation ς=φ/(ρσr2)\riskyshare = \EpremLog/(\CRRA\Evarr) with a premium of φ=0.03\EpremLog = 0.03 and σr=0.2\sigma_{\risky} = 0.2, drawn by the notebook.

The Approximation Error for the Portfolio Share in Risky Assets \riskyshare Is Small

Figure 2:The Approximation Error for the Portfolio Share in Risky Assets ς\riskyshare Is Small

Note: The exactly optimal share minus the approximation, for the same parameters. The notebook finds the exact share numerically, and its sliders vary the premium and the volatility; run it in your browser with nothing to install, in JupyterLite, or open it in Colab or Binder.

1Appendix: The Campbell & Viceira (2002) Approximation

For mathematical analysis (especially under the assumption of CRRA utility) it would be convenient if we could approximate the realized arithmetic portfolio return factor R~t+1\Rport_{t+1} by the geometric one, R1−ςRt+1ς\Rfree^{1-\riskyshare}\Risky_{t+1}^{\riskyshare}, because the log of the geometric factor is r(1−ς)+rt+1ς=r+(rt+1−r)ς=r+φt+1ς\rfree (1-\riskyshare)+\risky_{t+1}\riskyshare = \rfree + (\risky_{t+1}-\rfree)\riskyshare = \rfree + \EpremLog_{t+1} \riskyshare, so the realized “portfolio excess return” would be simply ςφt+1\riskyshare\EpremLog_{t+1}. Unfortunately, for ς\riskyshare values well away from 0 and 1 (that is, for any interesting values of portfolio shares), the log of the geometric mean is a badly biased approximation to the log of the arithmetic mean when the variance of the risky asset is substantial.

Campbell & Viceira (2002) propose instead

r~t+1≈r+ςφt+1+ςσr2/2−ς2σr2/2\begin{gathered}\begin{aligned} \rport_{t+1} & \approx \rfree+ \riskyshare \EpremLog_{t+1}+\riskyshare \sigma^{2}_{\risky}/2 - \riskyshare^{2}\sigma^{2}_{\risky}/2 \end{aligned}\end{gathered}

To see one virtue of this approximation,[9] note (using NormTimes and SumNormsIsNorm) that since the mean and variance of φt+1ς\EpremLog_{t+1} \riskyshare are respectively ς(r−σr2/2−r)\riskyshare(\risky - \sigma^{2}_{\risky}/2- \rfree) and ς2σr2\riskyshare^{2} \sigma^{2}_{\risky}, fact LogELogNormTimes implies that

log⁡Et[eφt+1ς]=ς(r−r−σr2/2)+ς2σr2/2\begin{gathered}\begin{aligned} \log \Ex_{t}[e^{\EpremLog_{t+1} \riskyshare}] & = \riskyshare (\risky - \rfree - \sigma^{2}_{\risky}/2) + \riskyshare^{2}\sigma^{2}_{\risky}/2 \end{aligned}\end{gathered}

which means that exponentiating then taking the expectation then taking the logarithm of (5) gives

log⁡Et[er~t+1]=log⁡er+log⁡Et[eςφt+1]+log⁡eςσr2/2−ς2σr2/2=r+ς(r−r−σr2/2)+ς2σr2/2+ςσr2/2−ς2σr2/2log⁡Et[er~t+1]−r=ς(r−r)\begin{gathered}\begin{aligned} \log \Ex_{t}[e^{\rport_{t+1}}] & = \log e^{\rfree} + \log \Ex_{t}[e^{\riskyshare \EpremLog_{t+1}}] + \log e^{\riskyshare\sigma^{2}_{\risky}/2 - \riskyshare^{2}\sigma^{2}_{\risky}/2} \\ & = \rfree + \riskyshare (\risky-\rfree-\sigma^{2}_{\risky}/2)+\riskyshare^{2}\sigma^{2}_{\risky}/2 +\riskyshare \sigma^{2}_{\risky}/2 - \riskyshare^{2}\sigma^{2}_{\risky}/2 \\ \log \Ex_{t}[e^{\rport_{t+1}}] - \rfree & = \riskyshare (\risky-\rfree) \end{aligned}\end{gathered}

or, in words: The expected excess portfolio return is equal to the proportion invested in the risky asset times the expected excess return of the risky asset.[10]

Footnotes
  1. u(c)=(1−ρ)−1c1−ρ\uFunc(\cRat) = (1-\CRRA)^{-1}\cRat^{1-\CRRA}.

  2. Both papers present the solution in the case with multiple risky assets; for the two-asset case, see two risky assets.

  3. A common interpretation is that this is the problem of a retired investor who expects to receive no further labor income. Note however that all risks other than the returns from financial investments have been ruled out; for example, health expense risk is not possible in this model, though recent research has argued such risk is important (maybe even dominant) later in life (cf. Ameriks et al. (2011)).

  4. Riskless labor income can trivially be added to the problem, because its risklessness means that (in the absence of liquidity constraints) it is indistinguishable from a lump sum of extra current wealth with a value equal to the present discounted value (using the riskless rate) of the (riskless) future labor income. Of course, in practice, labor income is not riskless, but when labor income is risky the problem no longer has the tidy analytical solution described here and must be solved numerically. See Carroll (2023) for an introduction to numerical solution methods.

  5. This expression differs slightly from that derived by Campbell & Viceira (2002), because we adjust the mean logarithmic return of the risky investment for its variance in order to keep the mean return factor constant for different values of the variance (cf. (4)), which makes comparisons of alternative levels of risk more transparent.

  6. Figure 2 shows the quality of the approximation.

  7. The approximation is motivated by the continuous-time solution, which is obtained using Ito’s lemma.

  8. We use the word “return” always to mean the logarithm of the corresponding “factor”; and when not explicitly specified, we always take the expectation before taking the log; if we wanted to refer to Et[r~t+1]\Ex_{t}[\rport_{t+1}] we would call it the expected log portfolio return (to distinguish it from the expected portfolio return, log⁡Et[er~t+1]\log \Ex_{t}[e^{\rport_{t+1}}]).

References
  1. Merton, R. C. (1969). Lifetime Portfolio Selection under Uncertainty: The Continuous-Time Case. Review of Economics and Statistics, 51(3), 247–257. 10.2307/1926560
  2. Samuelson, P. A. (1969). Lifetime Portfolio Selection by Dynamic Stochastic Programming. Review of Economics and Statistics, 51, 239–246.
  3. Campbell, J. Y., & Viceira, L. M. (2002). Appendix to Strategic Asset Allocation: Portfolio Choice for Long-Term Investors. Oxford University Press. https://scholar.harvard.edu/files/campbell/files/bookapp.pdf
  4. Mehra, R., & Prescott, E. C. (1985). The Equity Premium: A Puzzle. Journal of Monetary Economics, 15(2), 145–161. 10.1016/0304-3932(85)90061-3
  5. Haliassos, M., & Bertaut, C. C. (1995). Why Do So Few Hold Stocks? The Economic Journal, 105(432), 1110–1129. 10.2307/2235407
  6. Dimson, E., Marsh, P., & Staunton, M. (2008). The Worldwide Equity Premium: A Smaller Puzzle. In R. Mehra (Ed.), Handbook of the Equity Risk Premium (pp. 467–514). Elsevier.
  7. Cocco, J. F., Gomes, F. J., & Maenhout, P. J. (2005). Consumption and Portfolio Choice over the Life Cycle. Review of Financial Studies, 18(2), 491–533.
  8. Ameriks, J., Caplin, A., Laufer, S., & Van Nieuwerburgh, S. (2011). The Joy Of Giving Or Assisted Living? Using Strategic Surveys To Separate Public Care Aversion From Bequest Motives. The Journal of Finance, 66(2), 519–561.
  9. Carroll, C. D. (2023). Solving Microeconomic Dynamic Stochastic Optimization Problems. Econ-ARK REMARK. https://llorracc.github.io/SolvingMicroDSOPs
  10. Samuelson, P. A. (1979). Why we should not make mean log of wealth big though years to act are long. Journal of Banking and Finance, 3(4), 305–307. 10.1016/0378-4266(79)90023-2