Next Article in Journal
The Dynamic Relationship between Investor Attention and Stock Market Volatility: International Evidence
Next Article in Special Issue
Mean–Variance Portfolio Efficiency under Leverage Aversion and Trading Impact
Previous Article in Journal
Customers’ Risk Tolerance and Suppliers’ Investment Inefficiency
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Allocation of Retirement Portfolios

1
Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, NY 11794, USA
2
Department of Finance, University of Illinois at Champaign-Urbana, Champaign, IL 61820, USA
*
Author to whom correspondence should be addressed.
J. Risk Financial Manag. 2022, 15(2), 65; https://doi.org/10.3390/jrfm15020065
Submission received: 14 December 2021 / Revised: 24 January 2022 / Accepted: 26 January 2022 / Published: 1 February 2022
(This article belongs to the Special Issue Dynamic Portfolio Investment with Changing Economic States)

Abstract

:
A retiree with a savings account balance, but without a pension, is confronted with an important investment decision that has to satisfy two conflicting objectives. Without a pension, the function of the savings is to provide post-employment income to the retiree. At the same time, most retirees want to leave an estate to their heirs. Guaranteed income can be acquired by investing in an annuity. However, that decision takes funds away from investment alternatives that might grow the estate. The decision is made even more complicated because one does not know how long one will live. A long life expectancy may require more annuities, and a short life expectancy could promote more risky investments. However there are very mixed opinions about both strategies. A framework has been developed to assess consequences and the trade-offs of alternative investment strategies. We propose a stochastic programming model to frame this complicated problem. The objective is to maximize expected estate value, subject to cash outflow constraints. The model is motivated by the Markowitz mean-variance approach, but with risk measured by CVaR and additional sophisticated constraints. The cash outflow shortages are penalized in the objective function of the problem. We use the kernel method to build position adjustment functions that control how much is invested in each asset. These adjustments nonlinearly depend upon asset returns in previous years. A case study was conducted using two variations of the model. The parameters used in this case study correspond to a typical retirement situation. The case study shows that if the market forecasts are pessimistic, it is optimal to invest in an annuity. The case study results, codes, and data are posted on our website.

1. Introduction

The problem of selecting optimal portfolios for retirement has unique features that are not addressed by more commonly used portfolio selection models used in trading. One distinct feature of a retirement portfolio is that it should incorporate the life span of an investor. The planning horizon depends on the age of investor, or more specifically, on a conditional life expectancy. Another important feature is to guarantee, in some sense, that the individual will be able to withdraw some amount of money every year from a portfolio by selling some predefined amount of assets without injecting external funds. Finally, one of the questions that the model tries to answer is, in what situation is it beneficial to invest in an annuity instead of more risky assets?
Most of the portfolio optimization literature considers portfolios focusing on risk minimization with some budget and expected profit constraints. The famous mean-variance (or Markowitz) portfolio Markowitz (1952) minimizes portfolio variance with constraints on the expected return. There are many directions that extend the original mean-variance portfolio and deal with its shortcomings. One direction is to substitute variance with some other risk measures. Variance does not distinguish positive and negative portfolio returns; however, investors are mostly concerned only with negative returns. Rockafellar and Uryasev (2000, 2002) and Krokhmal et al. (2002) used conditional value-at risk (CVaR) instead of the variance. CVaR is a convex function of its random variable and therefore problems involving CVaR can be solved efficiently in many cases. Another important risk measure, which is frequently used in trading, is drawdown. Drawdown can be optimized with convex and linear programming; see Zabarankin et al. (2014). Another extension of the portfolio theory focuses on dynamic models. In dynamic models, the decision to invest is made over time. The dynamic models can be of two types, continuous-time and discrete-time (multistage). In continuous-time ones, the decision to invest is made continuously, and in discrete-time ones, the investment decisions take place at specific times. For continuous-time portfolio selection, see Merton (1969, 1971). For the discrete-time stochastic control model, see Samuelson (1969). A comprehensive literature review on dynamic models is given in Rizal and Wiryono (2015). Multistage models can be formulated as stochastic optimization problems. Mulvey and Shetty (2004) and Mulvey and Vladimirou (1992) developed a general multistage approach for modeling financial planning problems. Shang et al. (2016) and Bogentoft et al. (2001) used stochastic programming to solve dynamic cash flow matching and asset/liability management problems, respectively. In general, it is very hard to solve multistage stochastic optimization problems formulated with scenario trees, due to the size of the problem (number of variables) growing beyond tractable bounds. It should be mentioned that calibration of such trees is a difficult non-convex optimization problem.
In order to avoid the dimensionality problems, Calafiore (2008) modeled the investment decisions as linear functions that remained the same across all scenarios and produced the investment decision based on previous performance of the asset.
Takano and Gotoh (2014) modeled the investment decisions with the kernel method, resulting in the nonlinear control functions depending upon returns of instruments.
We followed the ideas of Takano and Gotoh (2014) and modeled the multistage portfolio decision process using the kernel method. The investment horizon was 35 years, starting from the retirement of the investor at the age of 65. The objective was to maximize the discounted expected terminal wealth subject to constraints on cash outflows from the portfolio. In every scenario, the discounted weighted portfolio value was calculated; the probabilities of death were used as weights. The probability of death was calculated from the U.S. mortality tables. In our scenario, the investor wants to have predetermined cash outflows obtained by selling a portion of the portfolio. Risk of shortage of these cash outflows was managed by penalizing the cash outflow shortage in the objective function. Furthermore, monotonicity constraints were imposed on the cash outflows from the portfolio. Without the monotonicity constraint, the model could not provide the necessary cash outflow on certain periods, and instead reinvest that amount to increase the expected estate value.
We conducted a case study corresponding to a typical investment decision upon retirement, in order to reveal the conditions leading to investments in annuities. Two types of asset return scenarios were considered. The first type assumes that the asset returns will be similar to the historically observed rates of the asset. The second type of scenario assumes the future asset returns will be significantly lower. These scenarios were created by subtracting 12% from the historical returns of all assets. The case study showed that for scenarios of the first type, where rates are similar to the ones observed in the past, investment in annuities is not optimal. However, when the asset growth rates are significantly lower, our model invests only in annuities.
There have been many attempts to solve the problem of investing for retirement, including how to allocate available funds between risky assets and annuities. See Bayraktar and Young (2009, 2016); Gao and Ulm (2012); Milevsky (1998). While these attempts approach the problem in a variety of ways, including minimizing the chance of lifetime ruin as in Bayraktar and Young (2016) or mandating a target estate value as in Bayraktar and Young (2009), they are not directly relevant to our approach due to the more sophisticated setting involving multi-stage stochastic programming and dynamic risk management.
Summarizing, this study had the following goals:
  • Develop a dynamic mathematical programming model proving an optimal investment strategy for an individual investor at retirement age.
  • Consider two sets of scenarios of indices returns: (1) an optimistic base case set obtained by bootstrapping existing history of returns of indices (this set is called “optimistic”); (2) a “pessimistic” set of scenarios with dramatically reduced returns of indices.
  • Conduct a case study. In particular, identify the percentage of capital which the investor should allocate to annuities.

2. Notations

We start with introduction of notation.
  • N := number of assets available for investments,
  • S := number of scenarios,
  • T := portfolio investment horizon,
  • r i , t s := rate of return of asset i = 1 , , N during period t = 1 , , T in scenario s = 1 , , S ; we will call rate of return by just return and denote the vector of returns by r t s = ( r 1 , t s , , r N , t s ) ,
  • v m , t s = { r m s , , r t 1 s } := set of returns observed from period m, until the end of period t 1 (not including the returns r t s ) in scenario s,
  • d t s := discount factor at time t in scenario s; discounting is done using inflation rate ρ t s , d t s = 1 / ( 1 + ρ t s ) t ,
  • p t := probability that a person will die at the age 65 + t (conditional that he is alive at the age of 65),
  • y i := vector of control variables for investment adjustment function,
  • f ( v t s , y i ) := investment adjustment function defining how much investment is made in each scenario s in asset i at the end of period t,
  • G ( y i ) := regularization function of control parameters,
  • K ( v m , t s , v m , t k ) := kernel function, k = 1 , , S ,
  • x i , t s := investment amount to i-th asset at the end of time period t in scenario s,
  • x i := investments to i-th asset at time t = 0 ,
  • u i , t s := adjustment (change in position) of asset i at the beginning of period t in scenario s,
  • R i , t s := cash outflow resulting from selling an asset i at the end of time t in scenario s,
  • V 0 := portfolio value at time t = 0 (initial investment),
  • V t s := portfolio value at time t in scenario s,
  • z := investment in annuity at time t = 0 (in dollars),
  • A t s := yield of annuity at the end of time period t in scenario s,
  • L := amount of money that the investor is planing to withdraw as each time t,
  • λ := regularization coefficient, λ > 0 ,
  • κ t := penalty for the cash flow shortage at time t,
  • α := upper bound on sum of absolute adjustments each year, expressed as a fraction of the portfolio.

3. Model Formulation and Research Hypotheses

This section develops a model for optimization of a retirement portfolio. We consider a portfolio, including stock indices, bond indices and an annuity. The annuity pays amount A t s z at the end of each period t and does not contribute funds to the expected estate value. The annuity is bought at time t = 0 and can not be bought or sold after that moment. It is also assumed that the tax rate is zero (tax free environment).
Given initial investments in assets x i , the dynamics of investments in stocks and bonds are as follows:
x i , 1 s = ( 1 + r i , 1 s ) x i , x i , t s = ( 1 + r i , t s ) ( x i , t 1 s + u i , t 1 s R i , t 1 s ) t = 2 , , T .
Variables u i , t s and R i , t s control how much is invested at the end of each period in each asset. u i , t s is a position adjustment for asset i at the end of time t in scenario s. R i , t s is cash outflow from the portfolio, generated from selling asset i at time t in scenario s. The variable u i , t s is defined as
u i , t s = f ( v t s , y i ) ,
where v t s is a set of returns for all assets, by time t, in scenario s; and y i are some parameters defining the function f. Therefore, u i , t s , are some nonlinear functions of previous returns of assets. The explicit form of function f is not specified in this section. The only requirement on function f is that it should be linear in y i ; i.e.,
f ( v t s , γ 1 y i 1 + γ 2 y i 2 ) = γ 1 f ( v t s , y i 1 ) + γ 2 f ( v t s , y i 2 ) ,
where γ 1 , γ 2 R . Furthermore, it should be noted that f does not change with t. The linearity of f with respect to y i is introduced to formulate the portfolio optimization problem as a convex programming problem.
The total asset adjustments must sum to 0, this is expressed as a constraint:
i = 1 N u i , t s = 0 .
In addition to (3), the sum of absolute adjustments (over each asset i) in each period t and scenario s is constrained to be less than or equal to some fraction α of the portfolio value in the previous year of the same scenario:
i = 1 N | u i , t s | α V t 1 s .
Constraint (4) serves as additional regularization on the adjustments. Without constraint (4) the values of u i , t s can potentially be very large in absolute value but cancel out due to opposite signs and still satisfy (3).
The value of the portfolio at the end of time period t in scenario s equals
V t s = i = 1 N x i , t s .
The objective is to maximize expected estate value of the portfolio. The expected estate value is the weighted average of the discounted expected portfolio values in each scenario, where the probabilities of death p t are used as weights. For every scenario s the portfolio value V t s , at the end of time period t, is discounted to time 0 using discounting coefficients d t s , defined by inflation; therefore,
discounted estate value in scenario s = t = 1 T p t d t s V t s .
By averaging over scenarios, we obtain the expected estate value:
1 S s = 1 S t = 1 T p t d t s V t s .
In order to avoid over-fitting the data, we included the regularization term G ( y i ) , defined for every instrument i. The total regularization term is
i = 1 N G ( y i ) .
The total cash outflow from selling the assets in the portfolio equals
cash flow from portfolio = i = 1 N R i , t s .
The amount of money that the investor receives from the portfolio and annuity at the end of time period t in scenario s equals A t s z + i = 1 N R i , t s . If A t s z + i = 1 N R i , t s < L then there is a shortage of cash outflow and the resulting amount is penalized in the objective. Let { κ t } t = 1 T be some decreasing sequence of positive numbers, the following function is a penalty term of cash outflow shortages in the objective
t = 1 T κ t L A t s z i = 1 N R i , t s + ,
where [ ] + = max { , 0 } . To illustrate why it is important that { κ t } t = 1 T is a decreasing sequence, consider the case where all κ t are equal. Furthermore, let us assume that there is a shortage of cash outflow, equal to the amount w, at some year t > 0 . As, κ t are all equal in (9), it does not make a difference for that penalty term if there is a shortage equal to w / t during every year until t, or just a single shortage of w at time t. However, if the amount of w / t is reinvested before time t in the portfolio, it will (probably) increase in value by the time t and therefore, it will increase the expected estate value of the portfolio. Thus, if { κ t } t = 1 T is not a decreasing sequence, the model will try to incur penalty as soon as possible, even if there are enough funds in the portfolio at that earlier date, and reinvest that shortage amount in the portfolio. Therefore the penalty from parameters κ t should outweigh any possible benefits from reinvesting at earlier dates. A simple formula for κ t is κ t = κ ( 1 + r ¯ ) T t , where κ > 1 is some constant and r ¯ is some percentage that is significantly greater than the average growth rate of any asset considered in the portfolio. The parameter κ t at time t was chosen to equal 2 times the gain from constant 20% compounding over the period T t , i.e, κ t = 2 × 1 . 2 ( T t ) . Since there is no investment strategy resulting in a larger expected gain than 2 times the constant 20% compounding over T t years for the considered set of scenarios, it is optimal to provide the required cash outflows to the retirees.
The model includes a constraint on monotonicity of the cash outflows from the portfolio
i = 1 N R i , t 1 s i = 1 N R i , t s .
Without the monotonicity constraint, the model might not provide necessary cash outflows at the end of certain years and instead, reinvest that amount to increase the expected estate value of the portfolio. The monotonicty constraint ensures that the cash outflow shortage occurs only in years where the portfolio value drops below the cash outflow amount at the end of the previous year.
We minimize the objective function, containing expected costs with minus signs, a regularization term with penalty coefficient λ > 0 and cash outflow shortage with penalty κ t :
1 S s = 1 S t = 1 T p t d t s V t s + λ i = 1 N G ( y i ) + t = 1 T κ t L A t s z i = 1 N R i , t s + .
The explicit form of function G is not defined in this section, however, it is assumed that the function G ( y ) is a convex function in y . This is important to formulate the problem as a convex optimization. The resulting objective function (11) is a convex function in y i and linear in V t s .
Further we provide the general model formulation.
min u i , t s , R i , t s , V 0 , V t s , y i , x i s , x i , t s , z 1 S s = 1 S t = 1 T p t d t s V t s + λ i = 1 N G ( y i ) + t = 1 T κ t L A t s z i = 1 N R i , t s + s . t . x i , 1 s = ( 1 + r i , 1 s ) x i i = 1 , , N ; s = 1 , , S x i , t s = ( 1 + r i , t s ) ( x i , t 1 s + u i , t 1 s R i , t 1 s ) i = 1 , , N ; t = 2 , , T ; s = 1 , , S i = 1 N x i = V 0 z V t s = i = 1 N x i , t s t = 1 , , T ; s = 1 , , S i = 1 N u i , t s = 0 t = 1 , , T ; s = 1 , , S i = 1 N R i , t s i = 1 N R i , t 1 s t = 2 , , T ; s = 1 , , S u i , t s = f ( v m , t s , y i ) i = 1 , , N ; t = 1 , , T ; s = 1 , , S i = 1 N | u i , t s | α V t 1 s t = 2 , , N ; s = 1 , , S i = 1 N | u i , 1 s | α ( V 0 z ) z 0 R i , t s 0 x i 0 i = 1 , , N x i , t s 0 i = 1 , , N ; t = 1 , , T ; s = 1 , , S

4. Special Case of General Formulation

This section presents a special case of the general problem formulation. We picked functions G ( y i ) and f ( r t s , y i ) similar to the model developed in Takano and Gotoh (2014).
Let m > 0 be some integer and K m ( v t s , v t k ) be the kernel function defined as follows:
K ( v m , t s , v m , t k ) = exp σ m i = 1 N l = t m 1 t 1 ( r i , l k r i , l s ) 2 ,
where σ > 0 is some constant. The parameter m controls how many previous years of information is used by the kernel function to calculate the portfolio adjustments. Given (13), the control function f ( v t s , y i ) is defined as
f ( v t s , y i ) = j = 1 S y i j K ( v m , t s , v m , t j ) , where y i = ( y i 1 , , y i S ) .
Function (14) is linear in y i . By substituting (14) in constraint (2), we get the following adjustment functions:
u i , t s = j = 1 S y i j K ( v m , t s , v m , t j ) i = 1 , , N ; t = 1 , , T ; s = 1 , , S .
We use L2 norm as the regularization function G ( y i ) :
G ( y i ) = y i 2 2 = s = 1 S ( y i s ) 2 .
Substituting (16) in the objective gives
1 S t = 1 T s = 1 S p t d t s V t s + λ i = 1 N y i 2 2 + t = 1 T κ t L A t s z i = 1 N R i , t s + .
This model can be reduced to a convex quadratic problem by linearizing (9). Other formulations are also possible. For example, using the L1 norm instead of the L2 norm in (16) leads to a linear programming formulation after linearization of (9). Another variation of this model could be linear (with respect to rates r i , t s ) adjustment functions instead of the nonlinear kernel adjustment functions. Linear investment adjustments will lead to a lower expected estate value. However, the dimensionality of the problem will be reduced significantly, because the problem size (the number of parameters to be optimized) will increase linearly with the number of scenarios, instead of quadratically, as with kernel functions.

5. Simulation of Return Scenarios and Mortality Probabilities

5.1. Historical Simulations

We simulate return scenarios of considered investment instruments for T years in the future. The simulations are based on end-of-year data of N assets over T ¯ years. Let t ¯ { 1 , , T ¯ } be a year index for a historical dataset and r ¯ i , t ¯ be a historical return of asset i. The returns of the indices are represented as an N × T ¯ matrix:
r ¯ 1 , 1 r ¯ 2 , 1 , , r ¯ N , 1 r ¯ 1 , 2 r ¯ 2 , 2 , , r ¯ N , 2 r ¯ 1 , T ¯ r ¯ 2 , T ¯ , , r ¯ N , T ¯
We generate return sample paths (scenarios) with the historical simulation method, also known as the “Bootstrap” method. The historical simulation method samples a random row from the matrix (18) and uses this row as a possible future realization of returns of instruments. Therefore the future simulation of returns is just a sampling of the matrix (18) with replacement. Each such sample represents a future dynamics of return of the assets. Note that the simulation method samples entire rows from matrix (18); therefore, the correlations among assets were maintained in the random sample.

5.2. Mortality Probabilities p t

Let τ be a random variable that denotes an age of death of the investor. The probability that an investor dies in time interval [ t 1 , t ) since retirement at the age 65 is defined as follows:
p t = P ( t + 64 < τ t + 65 | τ > 65 ) , t = 1 , , T .
It is possible to calculate p t using the mortality table of USA. We use the mortality Table 1, which gives probability p ^ t that t + 64 < τ t + 65 , on the condition that τ > t + 64 :
p ^ t = P ( t + 64 < τ t + 65 | τ > t + 64 ) , t = 1 , , T .
It can be shown that
p t = p ^ t , if t = 1 p ^ t j = 1 t 1 ( 1 p ^ j ) , if t = 2 , , T
Figure 1 shows p t as the function of age t.

6. Case Study

6.1. Case Study Parameters

The case study results, codes and data are posted at the website, see Pertaia and Uryasev (2019).
This case study considers a typical retirement situation in the USA. Two variants of future asset return scenarios are considered. These two variants correspond to an optimistic and pessimistic view regarding the future market dynamics. In the optimistic case, the future returns over 35 years, for all instruments, are sampled from the historical returns over the most recent 30 years. In the pessimistic case, the market is assumed to enter into a stagnation, similar to the Japaneses market, which has approximately zero cumulative return for the most recent 30 years. In the pessimistic case, 12% is subtracted from each asset return, every year, in every scenario.
Here are the parameters of the model, which correspond to typical retirement conditions in the USA.
  • The retiree is 65 years old.
  • Investment horizon is 35 years.
  • Portfolio is re-balanced at the end of each year.
  • Retiree is a male (mortality probabilities for males are used in objective function).
  • $500,000 is available for investment at time t = 0 .
  • Yearly inflation rate is 3% during the entire investment horizon.
  • Yearly rate of return of annuity is 5%.
  • Adjustment rules use kernel functions with parameter σ = 1 .
  • λ = 100 .
  • κ t = 2 × 1 . 2 ( 35 t ) .
  • α = 20 % .
  • m = 5
There are 10 stock and bond indexes available for investment; see Table 2.
For each index, 30 years of yearly returns (from 1985, to 2015) were used to create future scenarios (return sample-paths). Each scenario included 35 yearly returns, sampled from the 30 year historical dataset (see the Historical Simulation method in Section 5). Two-hundred scenarios were generated for both optimistic and pessimistic cases. One-hundred scenarios out of 200, for both optimistic and pessimistic scenario datasets, were used to find optimal investment rules in the model. The remaining 100 scenarios, not included in the optimization, were used for evaluating the out-of-sample performance of the model.

6.2. Optimal Portfolio

The considered optimization problems were reduced to quadratic programming, by linearizing function (9) in the objective with Portfolio Safeguard (PSG) Package, AORDA (2021). Gurobi version 8.1.0 and Pyomo version 5.5.0 were used for solving the resulting quadratic programming problem, which were automatically called by PSG. The case study link (Pertaia and Uryasev (2019)) contains the corresponding code.
The coefficients of the adjustment functions y i were obtained by solving the quadratic optimization problem corresponding to the optimal portfolio problem (12). Next, the adjustment values for the out-of-sample dataset were evaluated, according to the formula (14). The adjustment functions, for end of the time moment t, took previous m rates of returns of all assets in the portfolio, observed in time interval [ t m , t 1 ] , and produced an asset adjustment for that time moment. Note that returns that go into these functions were different on each scenario; therefore, the adjustment values were different in each scenario as well.
In order to calculate the portfolio values on the out-of-sample data, the cash outflows R i , t s are required. The model does not provide the cash outflow R i , t s for the out-of-sample scenarios, as those values are calculated for the in-sample scenarios. Therefore, it is unclear what values of R i , t s should be used in the out-of-sample scenarios. Additionally, despite the constraint on positivity of asset positions in the in-sample optimization problems, a small portion of the assets may be allocated to short positions in out-of-sample runs. Usually, retirement portfolios do not have short positions, since that is considered a risky strategy and is therefore not suitable for a risk averse retiree investor. Next, we show how to circumvent these problems for the out of sample datasets.
Let P + s , t and P s , t be the total dollar investment in long and short positions, in a portfolio at the end of time period t in scenario s:
P + s , t = i = 1 N [ x i , t s ] + , P s , t = i = 1 N [ x i , t s ] + .
The cash outflows are calculated as follows:
R i , t s = L [ x i , t 1 s ] + P + s , t 1 .
Thus, the cash outflows originate only from the long positions and are proportional to P + s , t 1 .
All short positions, at the end of time period t in scenario s, are set to 0. As a result, the amount of money equal to P s , t has to be subtracted from the remaining (long) part of the portfolio. To shrink the portfolio by P s , t , each long asset position is reduced in a proportion to P + s , t . Thus, the new positions x ¯ i , t s are
x ¯ i , t s = 0 , if x i , t s 0 x i , t s x i , t s P + s , t P s , t , otherwise .
Table 3, Table 4, Table 5, Table 6 and Table 7 show the average (over scenarios) investments in assets over time for optimistic out-of-sample scenarios, corresponding to the model (12), with the minimum cash flow requirements L { $ 10,000 ; $ 30,000 ; $ 50,000 ; $ 70,000 ; $ 90,000 } . Table 8, Table 9, Table 10 and Table 11 show the average (over scenarios) investments in assets over time for pessimistic out-of-sample scenarios, corresponding to the model (12), with the minimum cash flow requirements L { $ 10,000 ; $ 25,000 ; $ 30,000 ; $ 50,000 } . Table 8, Table 9 and Table 10, show that, in the pessimistic case, for L = $ 10,000 , the model invests 30% of funds in the annuity and for L = $ 25,000 , 100% of investment goes into the annuities. However, for L = $ 30,000 the model decreases the annuity investment to 56%. As for L = $ 50,000 (and higher) nothing is invested in the annuities and the model selects the stock/bond indexes. Figure 2 shows the average (taken over scenarios) portfolio values through time, constructed using the adjustment functions, corresponding to the model (12) with the minimum cash flow requirements of L { $ 10,000 ; $ 30,000 ; $ 50,000 ; $ 70,000 ; $ 90,000 } . However, in the optimistic scenarios, the model does not invest in annuities at any minimum cash outflow requirement L.
The main conclusion from the numerical experiments is that the optimal policy does not invest in the annuities in the base-case optimistic set of scenarios based on historical bootstrapping. Even for the pessimistic set of scenarios, investment in annuities is optimal only for low cash outflows not exceeding $30,000. Therefore, investment in annuities can be recommended only to conservative investors anticipating very pessimistic market conditions and only in the case that the investor needs small cash outflows.

6.3. Expected Shortage Times for Different Cash Outflows L

When the investor demands higher cash outflows from the portfolio, the expected estate value of the portfolio should decrease. Furthermore, with higher cash outflow demands, there are higher chances that there will not be enough money in the portfolio, at some point, to finance these outflows.
To measure the cash outflow shortage resulting from the different values of L, the following measure, named expected time shortage (or ETS) is defined:
E T S ( L ) = 1 S s = 1 S t = 1 T p t ( T t ) L t = 1 T R i , t s + L
ETS is measured in years and calculates the amount of time the retiree will spend without the necessary cash outflow L.
The parameters of the case study are used to construct the ETS values for the optimistic and pessimistic sets of scenarios. ETS is calculated on the in-sample data, for the cash outflow values of L { $ 10,000 ; $ 15,000 ; $ 20,000 ; ; $ 100,000 } . The resulting ETS values are shown in Figure 3 and Figure 4 for optimistic and pessimistic sets of scenarios, respectively.
Figure 3 shows that, with the optimistic set of scenarios, the retiree can have cash outflows up to $50,000, without having any shortage at any time. For the values of L greater than $50,000, the ETS grows roughly linearly. For L = $ 100,000 the retiree will spend most of their expected life without necessary cash outflow, because the portfolio cannot provide this much cash outflow, given the initial investment of $500,000.
It should be noted that, in the pessimistic case, if L $ 25,000 , the annuities can fully cover the cash flow requirements, and therefore ETS = 0. However, if L > $ 25,000 , the investment in the annuities can no longer cover the cash outflow requirements. Even if the entire initial investment goes into the annuities, it will provide only A · z = 5 % · $ 500,000 = $ 25,000 . Therefore, for L values higher than $25,000, the model starts to invest in stock and bond indexes and the ETS is greater than 0.
For the pessimistic case, if the cash flow requirement is L = $ 100,000 , the ETS is almost equal to the life expectancy of the retiree. This happens because, in most pessimistic scenarios, the portfolio shrinks to 0 in 3 or 4 years for L = $ 100,000 . However, if L = $ 30,000 , in the pessimistic case, the retiree still has relatively small ETS values of around 3 years.
Higher values of expected estate result in lower values of ETS. Figure 5 illustrates the relationship between expected estate and ETS for the optimistic case. Figure 5 is constructed by solving problem (12) for cash outflow values of L { $ 10,000 ; $ 15,000 ; ; $ 100,000 } and plotting the resulting values of ETS and expected estate.
We can conclude that in the optimistic case, the retiree can withdraw up to $50,000 with zero ETS. For the pessimistic case, the retiree can withdraw up to $25,000 with zero ETS.

7. Summary

We developed a multi-period investment model for retirement portfolios. The parameters of the model represent a typical portfolio selection problem solved in the beginning of retirement. The model maximizes the expected estate value of an investor subject to constraints on minimum cash outflows from the portfolio. Investment decisions are based on adjustment rules implemented with kernel functions.
The case study showed the performance of the model in pessimistic and optimistic asset return scenarios. In the pessimistic case, the market was assumed to enter a long term stagnation modeled by subtracting 12% from all rates of returns of the stock/bond indexes considered for investment. In this case, it is optimal to invest a considerable portion of initial capital in annuities. However, investment in annuities is optimal only for low cash outflows not exceeding $30,000. In the optimistic case, the returns of stock/bond indexes are expected to remain similar to past observations. In this case, it is not beneficial to invest in annuities for the given model parameters. Summarizing, investment in annuities can be recommended only to conservative investors anticipating very pessimistic market conditions and only in the case that the investor needs small cash outflows.
We defined a new cash outflow shortage measure called expected time shortage (ETS). The ETS calculates the number of years with a shortage of cash outflow, given the retiree’s minimum cash outflow requirements. The case study shows that even in the pessimistic asset return set of scenarios, a retiree can have zero ETS for some small cash outflows, due to significant investment in the annuities. In particular, in the optimistic and pessimistic cases, the retiree could withdraw up to $50,000 and $25,000, accordingly, with zero ETS.

Author Contributions

Methodology, K.M., M.L., M.M. and S.U.; formal analysis, K.M., M.L., M.M. and S.U.; investigation, K.M., M.L., M.M. and S.U.; data curation, K.M., M.L., M.M. and S.U.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is available at the website: Pertaia and Uryasev (2019).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. AORDA. 2021. American optimal decisions (aorda), Portfolio Safeguard (PSG). Available online: http://www.aorda.com/ (accessed on 13 December 2021).
  2. Bayraktar, Erhan, and Virginia R. Young. 2009. Optimal deferred life annuities to minimize the probability of lifetime ruin. North American Actuarial Journal 13: 141–54. [Google Scholar] [CrossRef]
  3. Bayraktar, Erhan, and Virginia R. Young. 2016. Optimally investing to reach a bequest goal. Insurance: Mathematics and Economics 70: 1–10. [Google Scholar] [CrossRef] [Green Version]
  4. Bogentoft, Erik, H. Edwin Romeijn, and Stan Uryasev. 2001. Asset/liability management for pension funds using cvar constraints. The Journal of Risk Finance 3: 57–71. [Google Scholar] [CrossRef] [Green Version]
  5. Calafiore, Giuseppe Carlo. 2008. Multi-period portfolio optimization with linear control policies. Automatica 44: 2463–73. [Google Scholar] [CrossRef]
  6. Gao, Jin, and Eric R. Ulm. 2012. Optimal consumption and allocation in variable annuities with guaranteed minimum death benefits. Insurance: Mathematics and Economics 51: 586–98. [Google Scholar] [CrossRef]
  7. Krokhmal, Pavlo, Jonas Palmquist, and Stan Uryasev. 2002. Portfolio optimization with conditional value-at-risk objective and constraints. The Journal of Risk 4: 11–27. [Google Scholar] [CrossRef] [Green Version]
  8. Markowitz, Harry. 1952. Portfolio selection. The Journal of Finance 7: 77–91. [Google Scholar]
  9. Merton, Robert C. 1969. Lifetime portfolio selection under uncertainty: The continuous-time case. The Review of Economics and Statistics 51: 247–57. [Google Scholar] [CrossRef] [Green Version]
  10. Merton, Robert C. 1971. Optimum consumption and portfolio rules in a continuous-time model. Journal of Economic Theory 3: 373–413. [Google Scholar] [CrossRef] [Green Version]
  11. Milevsky, Moshe A. 1998. Optimal asset allocation towards the end of the life cycle: To annuitize or not to annuitize? Journal of Risk and Insurance 65: 401–26. [Google Scholar] [CrossRef]
  12. Mulvey, John M., and Bala Shetty. 2004. Financial planning via multi-stage stochastic optimization. Computers and Operations Research 31: 1–20. [Google Scholar] [CrossRef]
  13. Mulvey, John M., and Hercules Vladimirou. 1992. Stochastic network programming for financial planning problems. Management Science 38: 1642–64. [Google Scholar] [CrossRef]
  14. Pertaia, Giorgi, and Stan Uryasev. 2019. Retirement Portfolio Selection Case Study. Available online: http://uryasev.ams.stonybrook.edu/index.php/research/testproblems/financial_engineering/case-study-retirement-portfolio-selection/ (accessed on 13 December 2021).
  15. Rizal, Nora A., and Sudarso K. Wiryono. 2015. A literature review: Modelling dynamic portfolio strategy under defaultable assets with stochastic rate of return, rate of inflation and credit spread rate. GSTF Journal on Business Review (GBR) 4: 80. [Google Scholar] [CrossRef]
  16. Rockafellar, R. Tyrrell, and Stan Uryasev. 2000. Optimization of conditional value-at-risk. The Journal of Risk 2: 21–41. [Google Scholar] [CrossRef] [Green Version]
  17. Rockafellar, R. Tyrrell, and Stan Uryasev. 2002. Conditional value-at-risk for general loss distributions. Journal of Banking and Finance 26: 1443–71. [Google Scholar] [CrossRef]
  18. Samuelson, Paul A. 1969. Lifetime portfolio selection by dynamic stochastic programming. The Review of Economics and Statistics 51: 239–46. [Google Scholar] [CrossRef]
  19. Shang, Danjue, Victor Kuzmenko, and Stan Uryasev. 2016. Cash flow matching with risks controlled by buffered probability of exceedance and conditional value-at-risk. Annals of Operations Research, 1–14. [Google Scholar]
  20. Takano, Yuichi, and Jun-ya Gotoh. 2014. Multi-period portfolio selection using kernel-based control policy with dimensionality reduction. Expert Systems with Applications 41: 3901–14. [Google Scholar] [CrossRef]
  21. Zabarankin, Michael, Konstantin Pavlikov, and Stan Uryasev. 2014. Capital asset pricing model (capm) with drawdown measure. European Journal of Operational Research 234: 508–17. [Google Scholar] [CrossRef]
Figure 1. Probability that a person dies while he/she is t + 64 years old ( t = 1 , , T ), conditional that he/she is alive at the age of 65.
Figure 1. Probability that a person dies while he/she is t + 64 years old ( t = 1 , , T ), conditional that he/she is alive at the age of 65.
Jrfm 15 00065 g001
Figure 2. The average (over scenarios) portfolio values for the optimistic out-of-sample dataset, constructed using adjustment functions, corresponding to model (12) with minimum cash outflow requirements L { $ 10,000 ; $ 30,000 ; $ 50,000 ; $ 70,000 ; $ 90,000 } .
Figure 2. The average (over scenarios) portfolio values for the optimistic out-of-sample dataset, constructed using adjustment functions, corresponding to model (12) with minimum cash outflow requirements L { $ 10,000 ; $ 30,000 ; $ 50,000 ; $ 70,000 ; $ 90,000 } .
Jrfm 15 00065 g002
Figure 3. ETS values for required cash flows L { $ 20,000 ; $ 30,000 ; $ 40,000 ; ; $ 100,000 } , calculated for the optimistic set of scenarios.
Figure 3. ETS values for required cash flows L { $ 20,000 ; $ 30,000 ; $ 40,000 ; ; $ 100,000 } , calculated for the optimistic set of scenarios.
Jrfm 15 00065 g003
Figure 4. ETS values for required cash flows L { $ 10,000 ; $ 15,000 ; $ 20,000 ; ; $ 100,000 } , calculated for the pessimistic set of scenarios.
Figure 4. ETS values for required cash flows L { $ 10,000 ; $ 15,000 ; $ 20,000 ; ; $ 100,000 } , calculated for the pessimistic set of scenarios.
Jrfm 15 00065 g004
Figure 5. Relationship between expected estate value and the ETS, for the optimistic set of scenarios with cash outflow values of L { $ 10,000 ; $ 15,000 ; ; $ 100,000 } .
Figure 5. Relationship between expected estate value and the ETS, for the optimistic set of scenarios with cash outflow values of L { $ 10,000 ; $ 15,000 ; ; $ 100,000 } .
Jrfm 15 00065 g005
Table 1. USA mortality table for the year 2017 with probabilities of death for male and female USA citizens. This table gives the conditional probability of death at some age, given that person is alive and one year younger than that.
Table 1. USA mortality table for the year 2017 with probabilities of death for male and female USA citizens. This table gives the conditional probability of death at some age, given that person is alive and one year younger than that.
p ^ (Age)
AgeMaleFemale
650.01580.0098
660.01700.0107
1190.88200.8820
Table 2. The list of assets in the retirement portfolio.
Table 2. The list of assets in the retirement portfolio.
Index NameIndex Abbreviation
Barclays MuniFI-MUNI
Barclays AggFI-INVGRD
Russell 2000USEQ-SM
Russell 2000 ValueUSEQ-SMVAL
Russell 2000 GrowthUSEQ-SMGRTH
S&P 500USEQ-LG
S&P 400 Mid CapUSEQ-MID
S&P Citi 500 ValueUSEQ-LGVAL
S&P Citi 500 GrowthUSEQ-LGGRTH
MSCI EAFENUSEQ
Table 3. Average investment (in thousand dollars) in assets over time for the optimistic out-of-sample scenarios with L = $ 10,000 . Results are averages across scenarios.
Table 3. Average investment (in thousand dollars) in assets over time for the optimistic out-of-sample scenarios with L = $ 10,000 . Results are averages across scenarios.
Asset Investmentt = 0t = 5t = 10t = 15t = 20t = 25t = 30t = 35
Annuity00000000
FI-MUNI00000000
FI-INVGRD03467111625
USEQ-SM00000000
USEQ-SMVAL028541041713606351177
USEQ-SMGRTH0112471319
USEQ-LG00000000
USEQ-MID50077914752791499310,59320,18336,797
USEQ-LGVAL00000000
USEQ-LGGRTH048153072139380
NUSEQ050801532684447621186
Table 4. Average investment (in thousand dollars) in assets over time for the optimistic out-of-sample scenarios with L = $ 30,000 . Results are averages across scenarios.
Table 4. Average investment (in thousand dollars) in assets over time for the optimistic out-of-sample scenarios with L = $ 30,000 . Results are averages across scenarios.
Asset Investmentt = 0t = 5t = 10t = 15t = 20t = 25t = 30t = 35
Annuity00000000
FI-MUNI00000000
FI-INVGRD33444516495138194
USEQ-SM00000000
USEQ-SMVAL6928571051923665921121
USEQ-SMGRTH01247142742
USEQ-LG00000000
USEQ-MID402612102518183136664212,57422,657
USEQ-LGVAL00000000
USEQ-LGGRTH2530541021814489202594
NUSEQ04569124209365628996
Table 5. Average investment (in thousand dollars) in assets over time for the optimistic out-of-sample scenarios with L = $ 50,000 . Results are averages across scenarios.
Table 5. Average investment (in thousand dollars) in assets over time for the optimistic out-of-sample scenarios with L = $ 50,000 . Results are averages across scenarios.
Asset Investmentt = 0t = 5t = 10t = 15t = 20t = 25t = 30t = 35
Annuity00000000
FI-MUNI07777101421
FI-INVGRD330244202206246328492680
USEQ-SM00000000
USEQ-SMVAL5713719428142469311631875
USEQ-SMGRTH00000000
USEQ-LG00000000
USEQ-MID36475884108224416820
USEQ-LGVAL00000000
USEQ-LGGRTH776566921573868572515
NUSEQ0333574104154255349
Table 6. Average investment (in thousand dollars) in assets over time for the optimistic out-of-sample scenarios with L = $ 70,000 . Results are averages across scenarios.
Table 6. Average investment (in thousand dollars) in assets over time for the optimistic out-of-sample scenarios with L = $ 70,000 . Results are averages across scenarios.
Asset Investmentt = 0t = 5t = 10t = 15t = 20t = 25t = 30t = 35
Annuity00000000
FI-MUNI00000000
FI-INVGRD195117674032354456
USEQ-SM00000000
USEQ-SMVAL46666748436599132
USEQ-SMGRTH00000000
USEQ-LG00000000
USEQ-MID107118736988170320596
USEQ-LGVAL00000000
USEQ-LGGRTH1369277901423507482300
NUSEQ166748354278162300
Table 7. Average investment (in thousand dollars) in assets over time for the optimistic out-of-sample scenarios with L = $ 90,000 . Results are averages across scenarios.
Table 7. Average investment (in thousand dollars) in assets over time for the optimistic out-of-sample scenarios with L = $ 90,000 . Results are averages across scenarios.
Asset Investmentt = 0t = 5t = 10t = 15t = 20t = 25t = 30t = 35
Annuity00000000
FI-MUNI00000000
FI-INVGRD65541765567
USEQ-SM00000000
USEQ-SMVAL70835130294676115
USEQ-SMGRTH00000000
USEQ-LG00000000
USEQ-MID1648530283367133302
USEQ-LGVAL00000000
USEQ-LGGRTH1401075648762044391522
NUSEQ615826129142342
Table 8. Average investment (in thousand dollar) in assets over time for the pessimistic out-of-sample scenarios L = $ 10,000 . Results are averages across scenarios.
Table 8. Average investment (in thousand dollar) in assets over time for the pessimistic out-of-sample scenarios L = $ 10,000 . Results are averages across scenarios.
Asset Investmentt = 0t = 5t = 10t = 15t = 20t = 25t = 30t = 35
Annuity147147147147147147147147
FI-MUNI00000000
FI-INVGRD13322211
USEQ-SM00000000
USEQ-SMVAL23232110
USEQ-SMGRTH00010000
USEQ-LG00000000
USEQ-MID350350360378384355311303
USEQ-LGVAL00000000
USEQ-LGGRTH04777554
NUSEQ04433322
Table 9. Average investment (in thousand dollar) in assets over time for the pessimistic out-of-sample scenarios L = $ 25,000 . Results are averages across scenarios.
Table 9. Average investment (in thousand dollar) in assets over time for the pessimistic out-of-sample scenarios L = $ 25,000 . Results are averages across scenarios.
Asset Investmentt = 0t = 5t = 10t = 15t = 20t = 25t = 30t = 35
Annuity500500500500500500500500
FI-MUNI00000000
FI-INVGRD00000000
USEQ-SM00000000
USEQ-SMVAL00000000
USEQ-SMGRTH00000000
USEQ-LG00000000
USEQ-MID00000000
USEQ-LGVAL00000000
USEQ-LGGRTH00000000
NUSEQ00000000
Table 10. Average investment (in thousand dollar) in assets over time for the pessimistic out-of-sample scenarios L = $ 30,000 . Results are averages across scenarios.
Table 10. Average investment (in thousand dollar) in assets over time for the pessimistic out-of-sample scenarios L = $ 30,000 . Results are averages across scenarios.
Asset Investmentt = 0t = 5t = 10t = 15t = 20t = 25t = 30t = 35
Annuity282282282282282282282282
FI-MUNI00000000
FI-INVGRD4315100000
USEQ-SM00000000
USEQ-SMVAL3518200000
USEQ-SMGRTH00000000
USEQ-LG00000000
USEQ-MID6133410000
USEQ-LGVAL00000000
USEQ-LGGRTH5427300000
NUSEQ2411100000
Table 11. Average investment (in thousand dollar) in assets over time for the pessimistic out-of-sample scenarios L = $ 50,000 . Results are averages across scenarios.
Table 11. Average investment (in thousand dollar) in assets over time for the pessimistic out-of-sample scenarios L = $ 50,000 . Results are averages across scenarios.
Asset Investmentt = 0t = 5t = 10t = 15t = 20t = 25t = 30t = 35
Annuity00000000
FI-MUNI00000000
FI-INVGRD6732610000
USEQ-SM00000000
USEQ-SMVAL95642463100
USEQ-SMGRTH00000000
USEQ-LG00000000
USEQ-MID14810542135100
USEQ-LGVAL00000000
USEQ-LGGRTH128853083100
NUSEQ62371331000
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Maritato, K.; Lane, M.; Murphy, M.; Uryasev, S. Optimal Allocation of Retirement Portfolios. J. Risk Financial Manag. 2022, 15, 65. https://doi.org/10.3390/jrfm15020065

AMA Style

Maritato K, Lane M, Murphy M, Uryasev S. Optimal Allocation of Retirement Portfolios. Journal of Risk and Financial Management. 2022; 15(2):65. https://doi.org/10.3390/jrfm15020065

Chicago/Turabian Style

Maritato, Kevin, Morton Lane, Matthew Murphy, and Stan Uryasev. 2022. "Optimal Allocation of Retirement Portfolios" Journal of Risk and Financial Management 15, no. 2: 65. https://doi.org/10.3390/jrfm15020065

Article Metrics

Back to TopTop