Hostname: page-component-745bb68f8f-b95js Total loading time: 0 Render date: 2025-02-11T17:28:23.142Z Has data issue: false hasContentIssue false

Lifetime ruin minimization: should retirees hedge inflation or just worry about it?*

Published online by Cambridge University Press:  20 May 2011

HUAXIONG HUANG
Affiliation:
Department of Mathematics and Statistics, York University, 4700 Keele Street, Toronto
MOSHE A. MILEVSKY
Affiliation:
Schulich School of Business, York University, 4700 Keele Street, Toronto (e-mail: milevsky@yorku.ca)
Rights & Permissions [Opens in a new window]

Abstract

Inflation for retirees is different from and mostly higher than the macro-economic (average) inflation rate for the entire population. In the U.S.A, for example, the Consumer Price Index for the Urban population (CPI-U) calculated and reported by the Bureau of Labor Statistics (BLS) has a lesser known cousin called the CPI-E (for the elderly) in which the sub-component weights are based on the consumption patterns of Americans above the age of 62. This suggests that Inflation-Linked Bond Funds (ILBFs) – whose individual component bond adjustments are based on broad population (CPI-U) inflation – might not be the best hedge for individual retirees’ cost of living. But then again, broad shocks to inflation are likely to impact both indices. So, motivated by the question – is it good enough? – the current paper uses lifetime ruin minimization (LRM) techniques to investigate the optimal allocation between an ILBF and a nominal investment fund for a retiree facing an exogenous liability. Our model trades off the benefit of an imperfect hedge against the cost of lower investment growth. However, our numerical results suggest that although ILBFs can be a large part of the optimal retirement portfolio, it should be treated as just another asset class in the broad optimization problem as opposed to a special or unique category.

Type
Articles
Copyright
Copyright © Cambridge University Press 2011

1 Introduction and motivation

The first American baby boomer turned 65 on 1 January 2011, with over 70 million to follow during the next 20 years. As a result, the financial services industry (in the U.S.A. and around the world) is shifting its attention from wealth accumulation to generating a sustainable retirement income. With the ongoing demise of defined benefit (DB) pensions, one of the main concerns in creating an appropriate income stream from a DC-style retirement account is the issue of inflation and price uncertainty over 25–35 potential years of retirement.

Against this backdrop, a number of financial commentators, including writers at the Wall Street Journal and the New York Times, have urged retirees to allocate a substantial portion of their retirement wealth to inflation-linked bonds, such as TIPS and I-BONDS. These investments generate periodic income in real, as opposed to nominal units. In fact, a number of large financial services companies, such as TIAA-CREF, PIMCO, Vanguard and Fidelity actively promote their inflation-linked mutual funds as a hedge against an increasing cost of living, especially for retirees. Holding a bond fund is more cost effective and feasible than holding individual TIPS.

Oddly enough, there seems to be very little discussion as to whether the Consumer Price Index for the Urban population (CPI-U) – which is the underlying adjustment mechanism for most inflation-linked bonds – is an adequate measure of the inflation rate experienced by retirees. In fact, since the early 1980s the U.S. Bureau of Labor Statistics (BLS) has been calculating an experimental index called the CPI-E, for Americans above the age of 62. This index weights its expenditure sub-components based on the spending habits of the elderly, as opposed to the aggregate (urban dwelling or wage earning) population. See Amble and Stewart (Reference Amble and Stewart1994) or Hobijn and Lagakos (Reference Hobijn and Lagakos2003) for more statistical and methodological details about this experimental index. The index was motivated by a group of legislators within the U.S. House of Representatives who wanted to introduce an alternative inflation index for retirees, one that could be used to better adjust Social Security payments, see House of Representatives (2001). Currently, Social Security payments in the U.S.A. are adjusted each year based on the previous year's inflation rate as measured by the CPI-W, which is a sub-component of the CPI-U.Footnote 1

Figure 1 all tables and charts are placed in the appendix exhibits the annualized gap between the CPI-E (elderly) versus the CPI-U (urban, used to adjust inflation-linked bonds) during the last 25 years. On average the gap is positive, by approximately 50 basis points per annum, which implies that the CPI-E increases at a greater rate, compared with the CPI-U. If both of these indices moved in lock-step, the gap would obviously be zero.

Figure 1. Annualized gap: CPI-E inflation minus CPI-U inflation.

The reason for this out performance – or excessive inflation rate – is that retirees have different consumption and spending habits compared with the general population. They spend more of their disposable income on medical care and housing, and of course these prices have increased more than other subcomponents during the last 25 years. Whether this will persist going forward, especially after the 2007/2008 decline in housing prices is debatable. But, once again, our agenda is not to forecast inflation but to alert readers to the implications of different inflation rates for different groups on optimal asset allocation. Table 1 displays the relative importance and weights of the major expenditure categories in the respective indices. Note that medical care is weighted 10.81% in the CPI-E but only 5.27% in the CPI-U, for example. Likewise, food is 12.87% in the CPI-E but is 16.56% higher in the CPI-U. Needless to say, any financial instrument – such as inflation-linked bond funds (ILBFs) – whose returns are linked to the CPI create an imperfect hedge for a typical retiree's personal inflation rate.

Table 1. Components of the CPI, wage-earners (W) versus urban (U) versus elderly (E)

These weights are computed and selected by the BLS, disclosed in April 2008. Note that weights vary over time based on the consumer expenditure survey conducted by the BLS and with inflation adjustments.

Figure 2 illustrates the extent to which the investment returns from holding an ILBF might differ from changes (or shocks) to the inflation rate for the elderly. It displays the ongoing market value of $1 invested in a hypothetical ILBF in December 2003, as well as the total change in the CPI-E starting from December 2003. Two additional curves were traced in Figure 2, showing a constant rate of 3% and 4% applied to the initial $1 investment. This hypothetical fund was constructed (by the authors) by averaging the monthly total returns the four most popular ILBFs available to retail investors in the U.S.A. (they are T. Rowe, Fidelity, TIAA CREF and Vanguard). In other words, at the end of the month the portfolio was rebalanced so that 25% was held in each of the four ILBFs. Note how the value of this investment increases over time, somewhere between 3% and 4% per year, but is more volatile and often falls outside this range. In contrast, the CPI-E-adjusted value of the initial $1 stays within this 3–4% band.

Figure 2. Values of $1 invested ILBF* versus the nominal value of an adjusted $1 liability. *Based on the average returns of four ILBFs.

Figure 3 takes this analysis one step further. It displays the results from regressing the monthly returns from the above-mentioned hypothetical portfolio of ILBFs against the changes in the unique CPI for the elderly. If the fit was perfect, the points would all fall on a 45° line. Clearly, this is not the case. The R 2=0.06, which is an estimate we will use later in our (continuous time) model. We do not have enough data points to conduct a statistically significant analysis for (the lower frequency) annual returns, although the R 2 is higher over longer return windows, as one would expect. On the other hand, we should note that when the analysis underlying Figure 2 and/or Figure 3 is conducted for any of the individual ILBFs, as opposed to the hypothetical average of the four – which diversifies the holder further against risks – the results are even more pronounced. More importantly, each individual retiree will have their own desired weightings for consumer expenditures, which might further reduce the correlations between their cost of living (e.g., CPI-ME) and aggregate inflation. Thus, although any given ILBF (or inflation-adjusted pension) is expected to increase over time as general prices increase, it is important not to confuse an investment that increases over time with an investment that provides perfect insurance against individual-specific shocks.

Figure 3. The fit. Changes in CPI-E versus changes in inflation- linked bond fund total return monthly periods from 12/02 to 11/08.

Once again, we are not the first to point out that personal inflation rates are distinct from macro-economic inflation measures, especially as it pertains to retirees. Just as an example, a group of well-known economists, Ernst et al. (Reference Ernst, Cockburn, Douglas, Epstein and Grilliches1997) used a unique dataset of prescription drug prices to focus on the challenges of measuring inflation for retirees. In fact, some researchers – such as Jennings (Reference Jennings2006), or Aziz et al. (Reference Aziz, Katz and Prisman2001) – have argued that as a result of these imperfect measures, CPI-linked bonds (which are the building blocks of ILBFs) should be disaggregated into individual sub-components so that consumers can reconstruct the bonds using their preferred weights. In this way, a suitably tailored bond fund would generate payments and provide income based on a rate that is linked to the medical care component of the CPI, for example. Retirees and for that matter any consumer with a basket of consumption that differs from the general CPI, would be able to pick and choose their own ILBF that would hedge their own inflation rate.

And, while this idea is great in theory, consumers and current retirees, who are still waiting for access to this engineering innovation, must contend with a generic ILBF that is weakly correlated to their desired consumption's cost of living adjustments. The pragmatic question is: How much wealth should they allocate to an imperfect insurance policy? Are they better off taking their chances on a long-run equity risk premium – or balanced stock/bond fund – that is expected to grow their portfolio over time?

This is the impetus for our paper. In the real-world absence of a perfect hedge, we gauge whether retirees are better-off investing in an asset class that is expected to grow over time – such as a diversified mutual fund – versus investing in an asset class that may (or may not) hedge a retiree's cost of living. We develop a model that illustrates how the correlation or lack thereof, between ILBF and a retiree's personal inflation rate impacts the demand and optimal allocation to this asset.

1.1 Our position within the literature

Most ‘lifecycle choice’ models are built and calibrated assuming the consumer seeks to maximize their end-of-period utility of wealth or their ongoing utility of consumption under a variety of risk-aversion parameters and inter-temporal rates of substitutions. This approach to the dynamic asset allocation problem was initiated by Merton (Reference Merton1971) and has spawned hundreds and possibly thousands of papers using these techniques. See Campbell and Viciera (Reference Campbell and Viciera2002) for an extensive literature review. A recent monograph by Ibbotson et al. (Reference Ibbotson, Milevsky, Chen and Zhu2007) takes the same utility-based approach as well. In related literature, Kothari and Shanken (Reference Kothari and Shanken2004) derive an optimal asset allocation with inflation-linked bonds, although not within the context of a retiree facing exogenous liabilities. We do use a number of their parameter estimates for our inputs, but differ in methodology and objective. The ‘calibrating’ asset allocation literature continues to grow, and due to the lack of space for a comprehensive literature review we refer the interested reader to recent papers by Horneff et al. (Reference Horneff, Maurer and Stamos2008) or Hoevenaars et al. (Reference Hoevenaars, Molenaar, Schotman and Steenkamp2008) as well as Cairns et al. (Reference Cairns, Blake and Dowd2006).

Our primary model, however, is predicated on a retiree who seeks to minimize the probability of lifetime ruin, assuming that their expenditures are stochastic (but exogenous) in nominal terms. The lifetime ruin minimization (LRM) framework has been adopted or discussed in a number of recent papers. Examples within this literature include the original paper by Browne (Reference Browne1995) with a deterministic horizon, Milevsky and Robinson (Reference Milevsky and Robinson2000) in a static framework with a stochastic horizon, Albrecht and Maurer (Reference Albrecht and Maurer2002), Young (Reference Young2004) in a dynamic framework, Moore and Young (Reference Moore and Young2006), Bayraktar and Young (Reference Bayraktar and Young2007) as well as the recent working paper Robinson and Tahani (Reference Robinson and Tahani2007). More broadly, our paper is related to Gupta and Li (Reference Gupta and Li2007) and Haberman and Vigna (Reference Haberman and Vigna2002). In all these articles, normative models are proposed for individuals as they transition into retirement under a variety of objective functions. They do not necessarily maximize discounted utility of lifetime consumption, a laMerton (Reference Merton1971), but instead use a framework that eliminates the need for requiring subjective discount rates and/or coefficients of relative risk aversion. The LRM approach resonates with some fraction of the wealth management community, which is one of the reasons why we have adopted it as a framework for this analysis. More importantly, by stripping away time, risk and substitution preferences we are able to directly measure the impact of not having a perfect hedge on the probability of success and the affordability of a given projected liability. Our paper is also related to recent work by Wang (Reference Wang2009) and Gao (Reference Gao2009) in that we also use the Hamilton–Jacobi–Bellman (HJB) formulation and techniques to analyze and extract the ruin probability from a retirement income perspective. From a methodological point of view the contribution of this paper is to extend the existing literature on ruin-minimizing optimal control to the case where there are three – as opposed to just one – state variables, in an incomplete market.

The additional level of complexity eliminates the ability to derive analytic solutions and expressions for the optimal control, since we are left with a collection of PDEs. Our paper provides numerical results that help shed light on the optimal demand for CPI-linked bonds, by retirees who are concerned with a simple objective function: minimizing the probability of running out of money before running out of life.

The remainder of this paper is organized as follows. Section 2 develops our underlying optimization model. Section 3 presents the derivation of the HJB equation satisfied by the ruin probability and methodologies for solving the equation. Section 4 provides numerical examples and displays relevant results. Section 5 concludes the paper and in the appendix more technical details are provided.

2 The underlying model

Consistent with the LRM references cited above, our model retiree starts with initial (nest egg) wealth denoted by W 0=w, from which he/she spends at a rate of $1 per annum, initially. Practitioners refer to this as an initial spending rate of 1/w. We assume that the retiree's exogenous liabilities or the annual rate of desired consumption – measured in nominal terms – evolve according to the diffusion process:

(1)
{\rm d}L_{t} \equals \pi L_{t} {\rm d}t \plus \xi L_{t} {\rm d}B_{t}^{\pi } \comma \quad L_{\setnum{0}} \equals 1\comma

where the parameter π is the expected cost-of-living adjustment (COLA), which is unique to the retiree, ξ is the volatility rate and B tπ is the Brownian motion driving the uncertainty COLA uncertainty. For example, a retiree might expect his or her spending will increase by π=4% each year, with a volatility of ξ=3%. Think of this as a personal inflation index/rate. Figure 2 – displayed earlier – represents one possible sample path of Lt from the December/03 to December/08 period.

In our simple model, the retiree's net worth of W 0=w can be allocated between two types of investment funds. The first is a lower-risk ILBF that evolves according to the diffusion process:

(2)
{\rm d}I_{t} \equals rI_{t} {\rm d}t \plus \eta I_{t} {\rm d}B_{t}^{r} \comma \quad I_{\setnum{0}} \equals 1\comma

where r is the expected return – in nominal terms – and η is the volatility. In this case Btr denotes the Brownian motion driving the ILBF. It is important to emphasize that It does not represent the price of a particular CPI-linked bond, nor does it represent a coupon rate. Rather, one can best think of It as a unit value of a mutual fundFootnote 2 that invests in a collection of CPI-linked bonds across various maturities, where all coupons and maturity payments are re-invested in the fund. Figure 2 also presents a possible sample path for It , over the same December/03 to December/08 period. It was clearly more volatile than Lt , and had a low correlation, as per Figure 3. Typical parameter values could be an expected r=5‰ nominal investment return, which would consist of 2‰ real-return and perhaps 2.5‰ inflation and a 0.5‰ inflation risk premium, although the exact decomposition is completely unnecessary for our analysis. The typical volatility η would be in the 8% range. These numbers are consistent with Ibbotson Associates estimates of the long-run volatility of a typical CPI-linked mutual fund, again in nominal terms.

The second investment alternative is an equity-based investment mutual fund that evolves according to the diffusion process:

(3)
{\rm d}Y_{t} \equals \mu Y_{t} {\rm d}t \plus \sigma Y_{t} {\rm d}B_{t}^{y} \comma \quad Y_{\setnum{0}} \equals 1\comma

where μ is the expected nominal return and σ>0 is the volatility of the risky fund. Typical parameter values are between 6% and 11% for the expected fund return μ and 10%–30% for the fund volatility σ. Nothing is new here.

The correlation structure between the Brownian motions driving the uncertainty is as follows. First, the correlation between Bt π (driving the retiree's personal inflation rate or COLA) and Btr (driving the ILBF), is relatively high but not 100%. It is denoted by ρπr. Likewise, the correlation between Bt π and Bty (the mutual fund) is obviously lower than ρπr, and will be denoted by ρπy. Finally, the correlation between Bty and Btr is denoted by ρry. Visually, the correlation structure is as follows:

which – to be invertible – imposes some natural conditions on the relationship between the three available correlations, which we will address later. In terms of calibration, Kothari and Shanken (Reference Kothari and Shanken2004) estimated the correlation between equity market returns and (hypothetical) CPI-linked funds returns is close to zero, so that ρ ry ≈0. For the most part, we will use this parameter estimate in our numerical examples. Likewise, the correlation ρπr between the CPI-linked fund and the retiree's COLA is a matter of debate – and the core concern within our paper – so we will use two extreme cases; the first being a low correlation of 40% and the second being a high correlation of 95%. Finally, the correlation between the investment fund and the liabilities is partially constrained by the invertability of the correlation matrix (4), which we discuss later.

Putting all this together, the dynamics for the retiree's investment portfolio will obey:

(5)
{\rm d}W_{t} \equals \alpha _{t} W_{t} {{{\rm d}I_{t} } \over {I_{t} }} \plus \lpar 1 \minus \alpha_{t} \rpar W_{t} {{{\rm d}Y_{t} } \over {Y_{t} }} \minus L_{t} {\rm d}t\comma \quad W_{\setnum{0}} \equals w\comma

where αt is the fraction of the portfolio that is allocated to the ILBF and (1−αt) is the fraction allocated to the investment or equity-based mutual fund. Again, Lt is defined as the consumption rate, i.e., consumption per unit of the time and is assumed exogenous as opposed to optimized.

With the LRM approach, the retiree's problem is to minimize the probability of ruin over a retirement horizon T=∞, but with a random time of death τd, under a given law of mortality. The infinite horizon might sound odd in the context of a human life, but obviously the force of mortality will eventually kill the retiree well before T=∞. So in some sense it is just a technicality. Also, the uncertainty in the length of human life can be hedged using pension annuities and our model can always be framed in terms of the choice between CPI-linked payout annuities and market-linked variable annuities. In other words, irrespective of whether they pose the problem as optimal allocations to non-annuitized investments or to income annuities, the underlying issue is the same.

Formally, the objective function can be stated as

(6)
\phi \lpar w\comma l\comma t\rpar \equals \mathop {\min}\limits_{\setnum{0} \leqslant \alpha _{t} \leqslant \setnum{1}} \Pr \lsqb \tau _{\setnum{0}} \lt \tau _{\rm d} \vert W_{t} \equals w\comma L_{t} \equals l\comma \tau_{\rm d} \gt t\rsqb \comma

where τ0 is the time at which Wt first becomes zero. The additional constraint αt⩾0 is imposed to preclude any short positions in the ILBF and αt⩽1 is imposed to preclude any leveraged position. Our technical objective is to locate the optimal control αt*(w, t) as a function of the drift and diffusion parameters μ, σ, π, ξ, r, η as well as the correlation parameters ρπr, ρπy and ρry. As we have stated earlier, we are particularly interested in how αt*(w, t) varies with respect to changes in the correlation ρπr between the CPI and retiree's COLA or personal inflation rate.

Before we proceed to solve for the optimal control and the corresponding minimal ruin probability, we note the following. The qualitative interpretation of ξ=0 is that the retiree's consumption liabilities will not fluctuate in nominal terms, although they can certainly fluctuate in CPI-adjusted terms. Likewise, the interpretation of η=0 is that the ILBF grows deterministically over time at a rate or r, which means that there is no CPI risk in the system. Finally, when both ξ=0 (i.e., no liability volatility) and η=0 (no CPI risk), our problem boils down to minimizing the probability of ruin in a complete market, which most recently was examined by Moore and Young (Reference Moore and Young2006). Keep in mind, however, that the absence of any risk-free asset in either nominal or personal-inflation adjusted terms, implies that the retiree can never reduce the ruin probability φ(w, t) to zero by immunizing (e.g., annuitizing) the consumption liability stream. The problem boils down to an economic trade-off between: (i) investing retirement wealth in an asset Yt that is expected to earn μ>π, and is greater than the projected increase in the retiree's cost of living, and (ii) investing in the relatively safe (and highly correlated) asset It , which is not expected to earn as much as Yt , since r<μ. Hence, the title of this paper, should retirees insure against their inflation rate by allocating a substantial portion of their wealth to an imperfect hedge or should they just ‘worry’ about it, and instead invest more aggressively.

3 The HJB equation

In this paper, we use Gompertz mortality for the hazard rate and the time variable is equivalent to biological age. The ruin probability satisfies the following HJB equation:

(7)
\lambda _{a} \psi \equals \psi _{a} \plus\! \mathop {\min}\limits_{\setnum{0}\leqslant \alpha \leqslant \setnum{1}} H\comma

where

(8)
\eqalign{H := \tab {1 \over 2}\left\{ {\left[ {\alpha ^{\setnum{2}} \eta ^{\setnum{2}} \plus \lpar 1 \minus \alpha \rpar ^{\setnum{2}} \sigma ^{\setnum{2}} \plus 2\alpha \lpar 1 \minus \alpha \rpar \rho _{ry} \eta \sigma } \right] \plus \xi ^{\setnum{2}} \minus 2\left[ {\alpha \rho _{\pi r} \eta \plus \lpar 1 \minus \alpha \rpar \rho _{\pi y} \sigma } \right]\xi } \right\}x^{\setnum{2}} \psi _{xx} \cr \tab \plus \left\{ {\left[ {\alpha rx \plus \lpar 1 \minus \alpha \rpar \mu x \minus 1 \minus \pi x} \right] \plus \xi ^{\setnum{2}} x \minus \left[ {\alpha \rho _{\pi r} \eta \plus \lpar 1 \minus \alpha \rpar \rho _{\pi y} \sigma } \right]\xi x} \right\}\psi _{x} \comma \cr}\hskip-12pt

with boundary conditions

(9)
\psi \lpar 0\comma a\rpar \equals 1\comma \quad \psi \lpar \infty \comma a\rpar \equals 0\comma \quad \psi \lpar x\comma \infty \rpar \equals 0.

Here we collapse w and l into a single variable x=w/l with ψ(x, a)=φ(w, l, a), and reduce the dimension of the HJB by one. The detailed derivation can be found in the appendix.

On a related note, when the investment portfolio consists of a risk-free asset and consumption rate is a constant, it can be shown that the ruin can be avoided under certain conditions, and the probability of ruin, when it is not zero, is a convex function of the wealth, as shown by Young (Reference Young2004). However, the situation considered in this paper is more complicated. When the liability and the risky assets are completely correlated, it is possible to avoid ruin when the wealth level is sufficiently high. On the other hand, when the risky assets and liability are not perfectly correlated, the probability of ruin is no longer a convex function of wealth. This adds an additional complexity to the solution procedure. We refer the readers to the appendix for more details.

4 Numerical solutions and examples

4.1 Parameter values

Before we present the results we provide the rationale for the choice of relevant parameter values. The first risky asset Yt , is a classical investment mutual fund that is expected to earn a nominal μ between 6% and 11% with a volatility of σ=20%. The second risky asset is the ILBF It , whose expected (nominal) return is r=5% with a volatility of η=8%. In most of our numerical examples we assume that the two risky assets (ILBF and equity/stock fund) are uncorrelated (even though our methodology can readily handle any correlation between the two). In fact, one of the benefits of using two uncorrelated assets is that this deliberately biases our results away from investing in the equity-based fund Yt as a hedge against the retiree's liabilities.

The retiree starts with a consumption liability rate of L 0=$1 per year at retirement but at a varying level of initial wealth W 0, allocated among Yt and It . The retiree withdraws a stochastic Lt , t>0 per year. As discussed within the derivation, we assume that 0⩽αt⩽1, and thus prohibit leverage and short sales. We assume that Lt increases by π=4% per annum, under two separate cases: volatility of ξ=3% (high) and 1% (low). We assume that the CPI-linked fund has a higher expected return r, relative to the liabilities and with greater volatility η=8%. Ultimately, we are interested in the minimal probability of ruin assuming the correlation between liabilities Lt and bond funds returns It is ρπr=95% (high) and 40% (low). Moreover, we are interested in how the decreasing return of the equity-based asset will force a greater allocation to the bond fund. We are also interested in the minimal ruin probability as a function of the correlation between the CPI-based bond fund and the consumption liability ρπr. The value of the remaining correlation coefficient is taken to be ρπy=0.25, which is chosen so that the positivity of the correlation matrix is guaranteed.

All the computations are carried out with an initial retirement age of a=65 and a time horizon of T=60 (a maximum age of 125). We have experimented with a longer time horizon and the results are virtually identical, since obviously the retiree dies well before T=60. We used the so-called Gompertz mortality rate in the numerical results, which implies that the actuarial force of mortality satisfies the following equation:

(10)
\lambda _{a} \equals \lambda _{\setnum{0}} \plus {1 \over b} \exp \left( {{{a \minus m} \over b}} \right).

The hazard/mortality rates become sufficiently large for a large T, so the survival rate is negligible after 125 years of life. For the result we displayed, we used the set of parameter values: m=86.3, b=9.5 and λ0=0.003. The source of these numbers is Milevsky (Reference Milevsky2006). Also, to carry out the computations, we have replaced the infinite domain 0⩽x<∞ by a finite one 0⩽xx max. We found that as long as x max is sufficiently large, the effect on the solution is small. We have chosen x max=200. This is equivalent to assuming that the ruin probability is zero when the wealth to liability ratio is 200. A total of 1,000 grid points are used for the state variable x=w/l and 240 time steps in time. We also experimented with larger number of grid points in x and a, and almost identical results were obtained. Finally, we have also compared the solutions obtained by explicit and implicit methods. Again, the results are virtually the same.

4.2 Cases and examples

In Table 2, we present the numerically computed values of the minimum ruin probability and the asset allocation percentage α* for a high correlation ρπr=95% between (shocks to) the COLA and (innovations to) the CPI-linked fund, together with a high volatility for the consumption liability ξ=3%. We are interested in the interaction of the assumed equity return μ and the initial spending rate 1/w.

Table 2. Optimal allocation to ILBF assuming high (95%) correlation with personal inflation

Note: Assets that are not invested in the ILBF are allocated to Equity Fund with mean return μ, and volatility σ=20%. The expected return from ILBF is r=5%, with volatility η=8%. Correlation between ILBF and Equity Fund return is ρ ry =0%. Correlation between personal inflation and Equity Fund return is ρπy=25%. Note that there is a constraint on the third correlation figure, ρπy needed for the positivity of the covariance matrix. For ρ ry =0 and ρπr=0.95, we find that ρπy<31.225%. In the computations, we have chosen ρπy=25%. All these numbers are consistent with – although not identical to – the Ibbotson Associates estimates for the behavior of TIPS-based funds, which invest in long-maturity CPI-linked bonds. The underlying model assumes the retiree starts with W 0 in liquid assets and spends L 0=1 per year, which increases by an expected π=4% with a volatility of ξ per year. The column headed by α* denotes the optimal allocation and the column headed by ψ denotes the minimum ruin probability.

High (ρπr=95%) correlation with personal inflation rate

High (ξ=3%) volatility of personal inflation rate; expected π=4%

Here is how to interpret the results. A 65-year-old retiree with initial investment wealth (nest egg) of $1,000,000 wants to spend/consume $40,000 per year adjusted for personal inflation. In the language of our model, this is a desired spending rate of 4% and is equivalent to an initial wealth to consumption ratio of W 0=25. The retiree estimates that this liability Lt will increase on average by π=4% per year (using continuous compounding) with a volatility of ξ=3%. The retiree has a choice of investing wealth between a CPI-linked mutual fund that is expected to earn r=5% per year with a volatility of η=8% and a riskier equity (or balanced) fund. The investment fund has a volatility of σ=20% and the expected return is denoted by μ. Tables 2 and 3 display the optimal allocation to the CPI-linked fund depending on the retiree's subjective estimate of the expected return μ. For example, if the equity (balanced) fund is expected to earn a mere μ=6%, which is 100 basis points more than the CPI-linked fund, the initial optimal allocation to the CPI-linked fund is 72.6%. The other 27.4% is allocated to the equity fund. Obviously, since the expected return μ from the equity fund is low the majority of the nest egg should rationally be placed in the CPI-linked fund. In this case, under the optimal strategy, the lifetime probability of ruin is 24.4%. This, recall, is the lowest possible ruin probability – from the set of all possible dynamic investment strategies – given the specified parameters. Now, if the retiree expects the return on the equity fund to be μ=11%, which is 500 basis points higher than the previous example, then the optimal allocation of the CPI-linked fund is a (much) lower than 56%. This means that 64% should be allocated to equity. Naturally, if one is more optimistic about equity (μ), the optimal allocation is greater. In this case, the minimal ruin probability is a mere 6.5% as opposed to 24.4%.

Table 3. Optimal allocation to ILBF assuming low (40%) correlation with personal inflation

Note: Assets that are not invested in the ILBF are allocated to Equity Fund with mean return μ, and volatility σ=20%. The expected return from ILBF is r=5%, with volatility η=8%. Correlation between ILBF and Equity Fund return is ρry=0%. Correlation between personal inflation and Equity Fund return is ρπy=25%. Note that the third correlation figure, ρπy<91.66% is forced by the two other values ρry=0 and ρπr=0.40. For consistency, we have chosen the same value as before, ρπy=25%.

Low (ρπr=40%) correlation with personal inflation rate

High (ξ=3%) volatility of personal inflation rate; expected π=4%

All the numbers in Table 2 assume that the correlation between the (shocks to) retirement liabilities Lt and the CPI-linked fund returns is quite high, at ρπr=95%. To contrast this case, in Table 3 we display the results of the minimum ruin probability ψ and optimal asset allocation percentage α* for a lower correlation value of ρπr=40% between the personal inflation rate shocks and the returns to the CPI-linked fund. We still assume the (high) volatility of personal inflation, ξ=3%. In this case, notice that for the same earlier-mentioned W 0=25, which is a spending rate of 4%, the optimal allocations to the CPI-linked fund are (generally, although not always) lower. The minimal ruin probabilities ψ are uniformly higher in Table 3 compared to Table 2. The intuition is simple. If the CPI-linked fund is a ‘worse’ hedge for the retirement liabilities, the lifetime probability of ruin is higher, even under the most efficient or optimal investment strategy.

In Table 4, we present a different perspective on the results, by comparing the allocations as a function of perceived equity-market volatility. For example, a retiree who expects equities (investment fund) to earn μ=10% and market volatility to be of the order of σ=20%, should allocate 72% of his nest egg to CPI-linked bonds. However, if the volatility can be controlled and is reduced to a mere σ=10% annualized, then the optimal allocation to the CPI-linked fund falls to 48%. And, these results are obtained when the correlation between personal inflation and returns to the CPI-linked fund are high at ρπr=95%. When the correlation is much weaker at ρπr=25%, the optimal allocations to the CPI-linked fund drop to between 32% and 65%. Note that despite the desire to minimize lifetime ruin probability, the allocation to the index-linked bond fund (ILBF) can be less than a third of the optimal portfolio.

Table 4. Optimal allocation to ILBF

Cost of living to increase by π=4% per year, with a volatility of ξ=3%. Expected return from the ILBF is r=5% per annum with a volatility of η=8%. The correlation between the equity fund and the ILBF is ρry=0%. Correlation between the liability and the investment fund is ρπy=25%. This table examines the low (3%) spending rate case.

Finally, in order to provide insight on how a dynamic (versus static) asset allocation strategy helps to lower the ruin probability ψ, we have plotted several key quantities in Figure 4. The underlying parameter values are the same as those in Table 2 with an expected equity return of μ=8%. First, the minimum ruin probability ψ as well as the ruin probabilities assuming a 100% allocation to the ILBF (α=1) or 100% allocation to the equity-based fund (α=0) are presented in Figure 4(a). Note how the optimal strategy results in a ruin probability that is lower than either of the extreme allocations. Then, the variation of the minimum ruin probability with age is plotted in Figure 4(b) and the dynamics of the optimal asset allocation ratio α* at different ages is plotted in Figure 4(c). In order to gain some intuition for Figures 4(b) and 4(c) one should focus on the x-axis point where w/l is 20 units. This represents an individual with assets worth 20 times their current level of consumption. At the age of 80, this is a very fortunate situation because the ruin probability is less than 5% in Figure 4(b). However, at the age of 65, this wealth-to-liability ratio induces a ruin probability of approximately 20%. Obviously, as the wealth-to-liability ratio increases, the entire vector of ruin probabilities declines, and reaches zero at approximately 50-to-1, Tables 5 and 6.

Figure 4. Ruin probability as function of wealth-to-liability ratio. (a) Impact of strategy. (b) Impact of age. (c) Allocation to CPI fund as function of wealth-to-liability ratio.

Table 5. Optimal allocation to ILBF (β=0.5)

Note: Assets that are not invested in the ILBF are allocated to Investment Fund with mean return μ, and volatility σ=20%. The expected return from ILBF is r=5%, with volatility η=8%. Correlation between ILBF and Investment Fund return is ρry=0%. Correlation between liability and Investment Fund return is ρπy=25%.

Table 6. Optimal allocation to ILBF (β=1.5)

Note: Assets that are not invested in the ILBF are allocated to Investment Fund with mean return μ and volatility σ=20%. The expected return from ILBF is r=5%, with volatility η=8%. Correlation between ILBF and Investment Fund return is ρry=0%. Correlation between liability and Investment Fund return is ρπy=25%.

Note in Figure 4(c) that the optimal allocation to the CPI-linked fund, which is represented in the y-axis, depends on age but asymptotes to about 80% when the wealth ratios are very high. At lower levels of wealth, the allocation to the CPI-linked fund is zero, especially at younger ages.

One of the main objectives of this study is to investigate the ‘correlation impact.’ In order to provide a graphical illustration, we have plotted α* (asset allocation) and ψ (ruin probability) against the wealth-to-liability ratio w/l for three different correlation scenarios: ρπr=0, 0.4 and 0.95. The results are displayed in Figure 5. This allows us to directly observe (only) the impact of the correlation between the personal inflation rate and the returns to the CPI-linked fund. One of the more noticeable results is the reversal of the asset allocation strategy – from safety to gambling – at a critical value of w/l. This is typical of all ruin minimization models, see Browne (Reference Browne1999), for example, where at some point the retiree decides to gamble (and hence volatility is welcomed.).

Figure 5. (a) CPI fund allocation as function of W/L ratio: high, medium and zero correlation. (b) Ruin probability as function of W/L ratio: 96%, 40% and 0% liability correlation.

In sum, we do not expect the normative model developed here to drive actual asset allocation recommendations at retirement, especially given the rather restricted nature of the assumed price dynamics (geometric Brownian motion) and investment choices. Furthermore, our selection of the minimum ruin criteria, as opposed to a more economically driven utility maximization criteria is meant to flush out the following point. Even the most risk averse preferences – i.e. ‘my main goal is to never run out of money’ – do not necessarily induce a complete allocation to the relatively safest asset, when the hedging and insurance are not perfect. To some readers this result might be obvious in an incomplete market, but we believe that it is worth emphasizing and illustrating in the context of inflation and the cost of living in retirement. There is no guarantee that an ILBF will hedge personal retirement liabilities. In our opinion, it should be treated as just another asset class in a liability-driven investment portfolio strategy, but with a unique set of risk, return and correlation parameters.

5 Final words

It is not widely known outside the specialized world of inflation connoisseurs and U.S. fed-watchers that inflation as measured for retirees (a.k.a. the elderly) is different from and mostly higher than the average inflation rate for the aggregate population. Indeed, in the U.S.A. the CPI has a lesser known relative, called the CPI-E (for the elderly) in which the sub-component weights are based on the consumption patterns of Americans above the age of 62. The CPI-E has consistently outpaced the broader CPI over the last 25 years that it has been measured, partially due to the greater weight placed on categories such as health-care expenses, long-term care facilities, the utilization of nursing services, etc. Although housing prices are also allocated greater weight in the CPI-E, it remains to be seen whether the recent decline in U.S. real estate will reverse the historical trend of CPI-E exceeding CPI-U (urban) and CPI-W (wage earner). Likewise, we do not address the persistent concern by some researchers and commentators that the CPI itself is biased downwards by hedonic adjustments and imputed rents, which would likely impact the CPI-E as well. One thing is for certain, there is a growing awareness that personal inflation rates for retirees are quite different from year-to-year changes in the CPI.

All of this suggests that investments linked to the CPI – such as retail mutual funds which contain various inflation-linked bonds (ILBF) – might not be a perfect hedge for individual retirees’ unique cost of living. In order to focus on this idea in a more rigorous manner, this paper extends LRM techniques – recently employed by Moore and Young (Reference Moore and Young2006), Young (Reference Young2004), Albrecht and Maurer (Reference Albrecht and Maurer2002) as well as Browne (Reference Browne1995) – to investigate the optimal asset allocation between a generic ILBF and a generic investment fund, for a retiree facing an exogenous liability stream that is imperfectly correlated to the returns on the CPI-linked fund. Our simple model trades off the benefits of an imperfect insurance hedge against the risky potential for investment growth.

Our results suggest that currently available CPI-linked funds – although they certainly are a large part of the optimal portfolio, especially for risk-averse investors as argued by Bodie and Clowes (Reference Bodie and Clowes2005) – should in fact be treated as just another asset class in the broad optimization problem, as opposed to a special or unique category. Practically speaking, retirees should (obviously) include CPI-linked investments in their portfolio universe, but not disproportionately so. In other words, it is very important to remember that inflation-linked assets might not keep up with the liability's inflation rate. So, in response to the question, which is the title of this paper, retirees should both hedge inflation risk and also worry about it.

Appendix – Detailed derivation and solution

A.1 Derivation of the HJB

The derivation is based on the method in Björk (1998, chapter 14). Assume an individual of age a at time zero. At any time t the individual aged a+t continues the optimal asset allocation strategy within the range 0⩽α⩽1 if they survive until t+h. On the other hand, if the individual aged a+t dies before time t+h, then the probability of ruin is zero. Let φ(w, l, t) denote the minimum ruin probability, i.e.

(11)
\phi \lpar w\comma l\comma t\rpar \equals \mathop {min}\limits_{\setnum{0}\leqslant \alpha _{t} \leqslant \setnum{1}} {\rm Pr}\lsqb \tau _{\setnum{0}} \lt \tau _{\rm d} \vert W_{t} \equals w\comma L_{t} \equals l\comma \tau _{\rm d} \gt t\rsqb.

Therefore,

(12)
\eqalign{\phi \lpar w\comma l\comma t\rpar \leqslant E\hskip2pt\lsqb \phi \lpar W_{t \plus h} \comma L_{t \plus h} \comma t \plus h\rpar \vert W_{t} \equals w\comma L_{t} \equals l\hskip1pt \rsqb _{h} \; p_{a \plus t}\plus E\hskip1pt\lsqb 0\vert W_{t} \equals w\comma L_{t} \equals l\hskip2pt \rsqb_{\hskip1pt h} q_{a \plus t} \comma \cr}

where the equality holds if and only if the optimal strategy is used. Here the symbols hpa+t and hqa+t ≔1− hpa+t are standard actuarial notation for the probabilities of surviving and dying between the time period h for an individual aged a+t, given by

(13)
_{h} p_{a} \equals {\rm e}^{ \minus \int_{\setnum{0}}^{h} \lambda _{{a \plus s}} {\rm d}s} \comma

where λt is the instantaneous hazard rate of the force of mortality. Now, let A denote the operator defined as follows:

(14)
\eqalign{ Af\hskip 1\lpar w\comma l\comma t\rpar \equals \tab f_{t} \plus \left[ {\alpha rw \plus \lpar 1 \minus \alpha \rpar \mu w \minus l} \hskip1pt \right]\hskip2pt f_{w} \cr \tab \plus {1 \over 2}\left\{ {\lpar \alpha \eta \rpar ^{\setnum{2}} \plus \lsqb \lpar 1 \minus \alpha \rpar \sigma \rsqb ^{\setnum{2}} \plus 2\alpha \lpar 1 \minus \alpha \rpar \rho _{ry} \eta \sigma } \right\}w^{\setnum{2}} f_{ww} \cr \tab \plus \pi lf_{l} \plus {1 \over 2}\lpar \xi l\rpar ^{\setnum{2}} f_{ll} \plus \lsqb \alpha \rho _{\pi r} \eta \plus \lpar 1 \minus \alpha \rpar \rho _{\pi y} \sigma \rsqb \xi wlf_{wl} \comma \cr}

where ft , fw , fl , etc. denote derivatives with respect to time, wealth and liability, respectively.

Using Ito's lemma on the function φ, we have

(15)
\eqalign{ \phi \lpar W_{t \plus h} \comma L_{t \plus h} \comma t \plus h\rpar \equals \tab \phi \lpar w\comma l\comma t\rpar \plus \int\limits_{t}^{t \plus h} {A\phi \lpar w_{s} \comma l_{s} \comma s\rpar {\rm d}s} \cr \tab \plus \int\limits_{t}^{t \plus h} {\phi _{w} \lpar w_{s} \comma l_{s} \comma s\rpar \lsqb \alpha \eta w_{s} {\rm d}B_{s}^{r} \plus \lpar 1 \minus \alpha \rpar \sigma w_{s} {\rm d}B_{s}^{y} \rsqb } \cr \tab \plus \int\limits_{t}^{t \plus h} {\phi _{l} \lpar w_{s} \comma l_{s} \comma s\rpar \xi l{\rm d}B_{s}^{\pi }.} \cr}

Combining equations (12), (14) and (15), we obtain

(16)
\phi \lpar w\comma l\comma t\rpar \leqslant _{h} p_{a \plus t} E^{w\comma l\comma t} \left[ {\phi \lpar w\comma l\comma t\rpar \plus \int\limits_{t}^{t \plus h} {A\phi \lpar W_{s} \comma L_{s} \comma s\rpar {\rm d}s} } \right]\comma

or

(17)
\phi \lpar w\comma l\comma t\rpar _{h} q_{a \plus t} \leqslant _{h} p_{a \plus t} E^{w\comma l\comma t} \left[ {\int\limits_{t}^{t \plus h} {A\phi \lpar W_{s} \comma L_{s} \comma s\rpar {\rm d}s} } \right].

Taking the limit of h→0 and using limh→0h −1 hqa+t a+t , we obtain

(18)
\phi \lpar w\comma l\comma t\rpar \lambda _{a \plus t} \leqslant A\phi \lpar w\comma l\comma t\rpar

and the equality holds when the optimal strategy is used, which is either given by the first-order condition or the boundary values α=0 or 1. Thus, we obtain the HJB equation

(19)
\phi \lpar w\comma l\comma t\rpar \lambda _{a \plus t} \equals \mathop {\min}\limits_{\setnum{0}\leqslant \alpha \leqslant \setnum{1}} \lcub A\phi \lpar w\comma l\comma t\rpar \rcub.

The boundary conditions are

(20)
\phi \lpar \infty \comma l\comma t\rpar \equals 0\comma \quad \phi \lpar 0\comma l\comma t\rpar \equals 1\comma \quad \phi \lpar w\comma 0\comma t\rpar \equals 0\comma \quad \phi \lpar w\comma \infty \comma t\rpar \equals 1.

Assuming an infinite time horizon, the terminal condition is

(21)
\phi \lpar w\comma l\comma \infty \rpar \equals 0.

When the mortality rate λ a+t is a constant, the infinite time horizon case leads to a time-independent version of (19) by dropping the term φ t . However, this is not the case (for realistic biological mortality), where the variable λ a+t is time dependent and the solution φ is therefore also time dependent. In this case, time is equivalent to age, and we can rewrite the HJB equation (19) using age a as the time variable

(22)
\lambda _{a} \phi \equals \phi _{a} \plus \mathop {\min}\limits_{\setnum{0}\leqslant \alpha \leqslant \setnum{1}} H\comma

where

(23)
\eqalign{ H:= \tab \left[ {\alpha rw \plus \lpar 1 \minus \alpha \rpar \mu w \minus l} \hskip2pt \right]\phi _{w} \plus \pi l\phi _{l} \cr \tab \plus {1 \over 2}\left\{ {\lpar \alpha \eta \rpar ^{\setnum{2}} \plus \lsqb \lpar 1 \minus \alpha \rpar \sigma \rsqb ^{\setnum{2}} \plus 2\alpha \lpar 1 \minus \alpha \rpar \rho _{ry} \eta \sigma } \right\}w^{\setnum{2}} \phi _{ww} \cr \tab \plus {1 \over 2}\lpar \xi l\rpar ^{\setnum{2}} \phi _{ll} \plus \lsqb \alpha \rho _{\pi r} \eta \plus \lpar 1 \minus \alpha \rpar \rho _{\pi y} \sigma \rsqb \xi wl\phi _{wl}. \cr}

The boundary conditions are still:

(24)
\eqalign{\phi \lpar \infty \comma l\comma a\rpar \equals 0\comma \quad \phi \lpar 0\comma l\comma a\rpar \equals 1\comma \quad \phi \lpar w\comma 0\comma a\rpar \equals 0\comma \quad \phi \lpar w\comma \infty \comma a\rpar \equals 1\comma \quad \phi \lpar w\comma l\comma \infty \rpar \equals 0.

A.2 Solution procedure

By introducing a similarity variable x=w/l, which is the ratio of wealth to the liabilities, it can be shown that the now two-dimensional equation ψ(x, a)=φ(w, l, a) satisfies

(25)
\lambda _{a} \psi \equals \psi _{a} \plus \mathop {\min}\limits_{\setnum{0}\leqslant \alpha \leqslant \setnum{1}} H\comma

where

(26)
\eqalign{ H:= \tab {1 \over 2}\left\{ {\left[ {\alpha ^{\setnum{2}} \eta ^{\setnum{2}} \plus \lpar 1 \minus \alpha \rpar ^{\setnum{2}} \sigma ^{\setnum{2}} \plus 2\alpha \lpar 1 \minus \alpha \rpar \rho _{ry} \eta \sigma } \right] \plus \xi ^{\setnum{2}} \minus 2\left[ {\alpha \rho _{\pi r} \eta \plus \lpar 1 \minus \alpha \rpar \rho _{\pi y} \sigma } \right]\xi } \right\}x^{\setnum{2}} \psi _{xx} \cr \tab \plus \left\{ {\left[ {\alpha rx \plus \lpar 1 \minus \alpha \rpar \mu x \minus 1 \minus \pi x} \right] \plus \xi ^{\setnum{2}} x \minus \left[ {\alpha \rho _{\pi r} \eta \plus \lpar 1 \minus \alpha \rpar \rho _{\pi y} \sigma } \right]\xi x} \right\}\psi _{x} \comma \cr}

or

(27)
\lambda _{a} \psi \equals \psi _{a} \plus H{ \ast } \comma

where

(28)
\scale95%\eqalign{\!\! H{ \ast } := \tab {1 \over 2}\left\{ {\left[ {\alpha { \ast }^{\setnum{2}} \eta ^{\setnum{2}} \plus \lpar 1 \minus \alpha { \ast } \rpar ^{\setnum{2}} \sigma ^{\setnum{2}} \plus 2\alpha { \ast } \lpar 1 \minus \alpha { \ast } \rpar \rho _{ry} \eta \sigma } \right] \plus \xi ^{\setnum{2}} \minus 2\left[ {\alpha { \ast } \rho _{\pi r} \eta \plus \lpar 1 \minus \alpha { \ast } \rpar \rho _{\pi y} \sigma } \right]\xi } \right\}x^{\setnum{2}} \psi _{xx} \cr \tab \plus \left\{ {\left[ {\alpha { \ast } rx \plus \lpar 1 \minus \alpha { \ast } \rpar \mu x \minus 1 \minus \pi x} \right] \plus \xi ^{\setnum{2}} x \minus \left[ {\alpha { \ast } \rho _{\pi r} \eta \plus \lpar 1 \minus \alpha { \ast } \rpar \rho _{\pi y} \sigma } \right]\xi x} \right\}\psi _{x}. \cr}

The boundary conditions are

(29)
\psi \lpar 0\comma a\rpar \equals 1\comma \quad \psi \lpar \infty \comma a\rpar \equals 0\comma \quad \psi \lpar x\comma \infty \rpar \equals 0.

To determine α*, we first apply the first-order condition Hα =0, which yields

(30)
\alpha { \ast } \equals \minus {{\lpar \rho _{\pi y} \sigma \minus \rho _{\pi r} \eta \rpar \xi \plus r \minus \mu } \over {\eta ^{\setnum{2}} \plus \sigma ^{\setnum{2}} \minus 2\rho _{ry} \eta \sigma }}{{\psi _{x} } \over {x\psi _{xx} }} \minus {{\sigma \lpar \rho _{ry} \eta \minus \sigma \rpar \plus \lpar \rho _{\pi y} \sigma \minus \rho _{\pi r} \eta \rpar \xi } \over {\eta ^{\setnum{2}} \plus \sigma ^{\setnum{2}} \minus 2\rho _{ry} \eta \sigma }}.

As we will demonstrate in the appendix, under certain conditions, ψ(x, a) will not be a strictly convex or concave function of x. Thus, the value of α* cannot be determined only by the first-order optimality condition. Since

(31)
H_{\alpha \alpha } \equals \lpar \eta ^{\setnum{2}} \plus \sigma ^{\setnum{2}} \minus 2\rho _{ry} \eta \sigma \rpar ^{\setnum{2}} x^{\setnum{2}} \psi _{xx}

and ψ is not convex, H αα can change signs. Therefore, we need both the first- and second-order conditions to determine the value of α* as follows:

  • If H αα>0 and α* given by (30) satisfies 0⩽α*⩽1, then that is the value for α*.

  • If H αα>0 and α* given by (30) satisfies α*>1, then α*=1.

  • If H αα>0 and α* given by (30) satisfies α*<0, then α*=0.

  • If H αα⩽0, we need to compare H at α=0 and α=1 and choose the α* value that gives the smaller H.

In order to solve equation (27) numerically, we need to truncate the domain in x as well as y and replace the boundary conditions in (29) by

(32)
\psi \lpar 0\comma a\rpar \equals 1\comma \quad \psi \lpar x_{{\rm max}} \comma a\rpar \equals 0\comma \quad \psi \lpar x\comma a_{{\rm max}} \rpar \equals 0\comma

where x max and a max are relatively large numbers. The computational results presented in this paper are obtained by choosing x max=200 and a max=125, based on the observation that further increasing the size of the domain produces almost identical solutions.

There are two possible ways of solving (27), one by discretizing the time derivative ψ a explicitly, whereas the other by using an implicitly method. In both cases, we discretize a as an =a 0+nδa for n=0, 1, …, N, where δa=(a maxa 0)/N, a 0 is the current age of the retiree. Furthermore, we set up a grid in x, using xj =jδx for j=0, 1, …, J, where δx=x max/J. Let ψ j (n) be an approximation of ψ(xj , an ), in the explicit method, we use

(33)
\psi _{a} \lpar x_{j} \comma a_{n} \rpar \approx {{\psi _{j}^{\lpar n \plus \setnum{1}\rpar } \minus \psi _{j}^{\lpar n\rpar } } \over {\delta a}}\comma

while

(34)
\psi _{a} \lpar x_{j} \comma a_{n \plus \setnum{1}} \rpar \approx {{\psi _{j}^{\lpar n \plus \setnum{1}\rpar } \minus \psi _{j}^{\lpar n\rpar } } \over {\delta a}}

is used for the implicit method.

A.2.1 Explicit method

In this case, equation (27) is approximated by

(35)
\hskip56pt\psi _{j}^{\lpar n \plus \setnum{1}\rpar } \equals \psi _{j}^{\lpar n\rpar } \plus \delta aF\lpar \psi _{j}^{\lpar n\rpar } \rpar \comma\hfill
(36)
\eqalign{ F\lpar \psi _{j}^{\lpar n\rpar } \rpar \equals \tab \minus \lambda _{a} \psi _{j}^{\lpar n\rpar } \cr \tab \plus {1 \over 2}\left\{ {\left[ {\alpha { \ast \comma j}^{\setnum{2}} \eta ^{\setnum{2}} \plus \lpar 1 \minus \alpha { \ast \comma j} \rpar ^{\setnum{2}} \sigma ^{\setnum{2}} \plus 2\alpha { \ast \comma j} \lpar 1 \minus \alpha { \ast \comma j} \rpar \rho _{ry} \eta \sigma } \right]} \right. \cr \tab \left. { \plus \xi ^{\setnum{2}} \minus 2\left[ {\alpha { \ast \comma j} \rho _{\pi r} \eta \plus \lpar 1 \minus \alpha { \ast \comma j} \rpar \rho _{\pi y} \sigma } \right]\xi } \right\}x^{\setnum{2}} \psi _{xx}^{\lpar n\rpar } \lpar x_{j} \rpar \cr \tab \plus \left\{ {\left[ {\alpha { \ast \comma j} rx \plus \lpar 1 \minus \alpha { \ast \comma j} \rpar \mu x_{j} \minus 1 \minus \pi x_{j} } \right] \plus \xi ^{\setnum{2}} x_{j} } \right. \cr \tab \left. { \minus \left[ {\alpha { \ast \comma j} \rho _{\pi r} \eta \plus \lpar 1 \minus \alpha { \ast \comma j} \rpar \rho _{\pi y} \sigma } \right]\xi x_{j} } \right\}\psi _{x}^{\lpar n\rpar } \lpar x_{j} \rpar \comma \cr}

where

(37)
\psi _{x}^{\lpar n\rpar } \lpar x_{j} \rpar \equals {{\psi _{j \plus \setnum{1}}^{\lpar n\rpar } \minus \psi _{j \minus \setnum{1}}^{\lpar n\rpar } } \over {2\delta x}}\comma \quad \psi _{xx}^{\lpar n\rpar } \lpar x_{j} \rpar \equals {{\psi _{j \plus \setnum{1}}^{\lpar n\rpar } \plus \psi _{j \minus \setnum{1}}^{\lpar n\rpar } \minus 2\psi _{j}^{\lpar n\rpar } } \over {\delta x^{\setnum{2}} }}.

In general, the value of α*,j can be determined using both the first- and second-order optimality conditions, as described earlier. However, owing to the explicit nature of the method, we can use a simpler approach as follows.

Let ψ j (n, *) denote the solution obtained using (36) and the first-order condition (30), ψ j (n, 0) and ψ j (n, 1) the solution of (36) using α=0 or α=1. We choose the final solution as

(38)
\psi _{j}^{\lpar n\rpar } \equals {\rm min}\left\{ {\psi _{j}^{\lpar n\comma} \vskip-1pt\ast ^{\rpar } \comma \psi _{j}^{\lpar n\comma \setnum{0}\rpar } \comma \psi _{j}^{\lpar n\comma \setnum{1}\rpar } } \right\}.

This method is conceptually simple and easy to implement. However, it is not efficient as small δa has to be used because of stability constraint associated with the explicit method.

A.2.2 Implicit method

Alternatively, we can use an implicit method by approximating equation (27) as

(39)
\psi _{j}^{\lpar n \plus \setnum{1}\rpar } \equals \psi _{j}^{\lpar n\rpar } \plus \delta aF\lpar \psi _{j}^{\lpar n \plus \setnum{1}\rpar } \rpar \comma
(40)
\eqalign{ F\lpar \psi _{j}^{\lpar n \plus \setnum{1}\rpar } \rpar \equals \tab\minus \lambda _{a} \psi _{j}^{\lpar n \plus \setnum{1}\rpar }\plus {1 \over 2}\left\{ {\left[ {\alpha { \ast \comma j}^{\setnum{2}} \eta ^{\setnum{2}} \plus \lpar 1 \minus \alpha { \ast \comma j} \rpar ^{\setnum{2}} \sigma ^{\setnum{2}} \plus 2\alpha { \ast \comma j} \lpar 1 \minus \alpha { \ast \comma j} \rpar \rho _{ry} \eta \sigma } \right]} \right. \cr \tab \left. { \plus \xi ^{\setnum{2}} \minus 2\left[ {\alpha { \ast \comma j} \rho _{\pi r} \eta \plus \lpar 1 \minus \alpha { \ast \comma j} \rpar \rho _{\pi y} \sigma } \right]\xi } \right\}x_{j}^{\setnum{2}} \psi _{xx}^{\lpar n \plus \setnum{1}\rpar } \lpar x_{j} \rpar \cr \tab \plus \left\{ {\left[ {\alpha { \ast \comma j} rx_{j} \plus \lpar 1 \minus \alpha { \ast } \rpar \mu x_{j} \minus 1 \minus \pi x_{j} } \right] \plus \xi ^{\setnum{2}} x_{j} } \right. \cr \tab \left. { \minus \left[ {\alpha { \ast \comma j} \rho _{\pi r} \eta \plus \lpar 1 \minus \alpha { \ast \comma j} \rpar \rho _{\pi y} \sigma } \right]\xi x_{j} } \right\}\psi _{x}^{\lpar n \plus \setnum{1}\rpar } \lpar x_{j} \rpar \comma \cr}

where

(41)
\psi _{x}^{\lpar n \plus \setnum{1}\rpar } \lpar x_{j} \rpar \equals {{\psi _{j \plus \setnum{1}}^{\lpar n \plus \setnum{1}\rpar } \minus \psi _{j \minus \setnum{1}}^{\lpar n \plus \setnum{1}\rpar } } \over {2\delta x}}\comma \quad \psi _{xx}^{\lpar n \plus \setnum{1}\rpar } \lpar x_{j} \rpar \equals {{\psi _{j \plus \setnum{1}}^{\lpar n \plus \setnum{1}\rpar } \plus \psi _{j \minus \setnum{1}}^{\lpar n \plus \setnum{1}\rpar } \minus 2\psi _{j}^{\lpar n \plus \setnum{1}\rpar } } \over {\delta x^{\setnum{2}} }}.

The simple approach used in the explicit method is not applicable here. Instead, we have to use both the first and second optimality condition and apply the procedure described earlier to determine α*,j. In a fully implicit method, we need to apply the iterative method to find the solution owing to the non-linear nature of the problem. In this paper, we use a semi-implicit version by computing α*,j using ψ j (n) instead of ψ j (n+1). By comparison, this semi-implicit method is simpler than the fully implicit one and much more efficient than the explicit method since a relatively large δa can be used for the computation.

Finally, as we alluded to earlier, when the investment portfolio consists of a risk-free asset and consumption rate is a constant, it can be shown that the ruin can be avoided under certain conditions, and the probability of ruin, when it is not zero, is a convex function of the wealth, as shown by Young (Reference Young2004). However, when the portfolio consists of only risky assets, it may not be possible to eliminate ruin completely. Furthermore, the probability of ruin may not be a convex function of the wealth; see Huang et al. (Reference Huang, Milevsky and Wang2004) for the case with one risky asset and a constant rate of consumption. The situation considered in this paper is more complicated. When the liability and the risky assets are completely correlated, it is possible to avoid ruin when the wealth level relative to liability is sufficiently high. On the other hand, when the risky assets and liability are not perfectly correlated, the situation is not as simple. In order to see this, we restrict ourselves to the case when the age if relatively small, so that both the mortality rate and the change of ruin probability with time/age is small. In other words, the Hamiltonian H*≈0 in equation (28). When x is small, from equation (28), it can be seen that parameter of the second-order derivative remains positive while that of first derivative is approximately −1. Since ψ x <0, we must have ψ xx <0, i.e., ψ is convex for small x. For large x, on the other hand, the coefficient of ψx is approximately

\left\{ {\alpha { \ast } r \plus \lpar 1 \minus \alpha { \ast } \rpar \mu \plus \xi ^{\setnum{2}} \minus \left[ {\alpha { \ast } \rho _{\pi r} \eta \plus \lpar 1 \minus \alpha { \ast } \rpar \rho _{\pi y} \sigma } \right]\xi \minus \pi } \right\}x\comma

which could be positive or negative. When

\pi \plus \lsqb \alpha { \ast } \rho _{\pi r} \eta \plus \lpar 1 \minus \alpha { \ast } \rpar \rho _{\pi y} \sigma \rsqb \xi \lt \alpha { \ast } r \plus \lpar 1 \minus \alpha { \ast } \rpar \mu \plus \xi ^{\setnum{2}}

for any 0⩽α*⩽1, the coefficient of ψ x is positive. Since the coefficient of ψ xx is also positive, ψ xx must be positive, i.e., ψ is now a concave function of x (scaled wealth).

Footnotes

1 We obviously don't want to get caught-up in the minutia of inflation measurement and forecasting. Technically though, the CPI-W is for a representative household where more than one-half of the household's income must come from clerical or wage occupations and at least one of the household's earners must have been employed for at least 37 weeks during the previous 12 months. The CPI-W's population represents about 32 percent of the total U.S. population and is a subset or part of the CPI-U's (which stands for Urban) population. The CPI-U is used for computing the coupons and inflation-adjustments on Treasury Inflation Protection Securities (TIPS) as well as I-Bonds. From a practical perspective the CPI-W and the CPI-U are much closer to each other in composition, relative to the experimental CPI-E.

2 We obviously do not attempt to offer or develop a complete model that prices all CPI-linked bonds in the most advanced term-structure sense. For something along those lines, we refer interested readers to Jarrow and Yildirim (Reference Jarrow and Yildirim2003) for example. We are focused on a basic portfolio of CPI-linked bonds within a mutual fund structure.

References

Albrecht, P. and Maurer, R. (2002) Self-annuitization, consumption shortfall in retirement and asset allocation: the annuity benchmark. Journal of Pension Economics and Finance, 1(3): 269288.CrossRefGoogle Scholar
Amble, N. and Stewart, K. (1994) Experimental price index for elderly consumers. Bureau of Labor Statistics Monthly Labor Review, 117(5): 1116.Google Scholar
Aziz, A., Katz, E. and Prisman, E. Z. (2001) Managing the risk of relative price changes by splitting index-linked bonds. Journal of Risk, 3(4): 6988.CrossRefGoogle Scholar
Bayraktar, E. and Young, V. R. (2007) Correspondence between lifetime minimum wealth and utility of consumption. Finance and Stochastics, 11(2): 213236.CrossRefGoogle Scholar
Bjork, T. (1998) Arbitrage Theory in Continuous Time. Oxford University Press, UK.CrossRefGoogle Scholar
Bodie, Z. and Clowes, M. J. (2005) Worry-Free Investing: A Safe Approach to Achieving Your Lifetime Financial Goals. 1st edition, Financial Times/Prentice Hall, New Jersey.Google Scholar
Browne, S. (1995) Optimal investment policies for a firm with a random risk process: exponential utility and minimizing the probability of ruin. Mathematics of Operations Research, 20(4): 937958.CrossRefGoogle Scholar
Browne, S. (1999) The risk and reward of minimizing shortfall probability. Journal of Portfolio Management, 25(4): 7685.CrossRefGoogle Scholar
Cairns, A. J. G., Blake, D. and Dowd, K. (2006) Stochastic lifestyling: optimal dynamic asset allocation for defined contribution pension plans. Journal of Economic Dynamics and Control, 30: 843877.CrossRefGoogle Scholar
Campbell, J. and Viciera, L. (2002) Strategic Asset Allocation: Portfolio Choice for Long Term Investors. Oxford University Press, Oxford, UK.CrossRefGoogle Scholar
Ernst, B., Cockburn, I., Douglas, C., Epstein, A. and Grilliches, Z. (1997) Is price inflation different for the elderly? A empirical analysis of prescription drugs. National Bureau of Economic Research (NBER) Working Paper, No. 6182.Google Scholar
Gao, J. (2009) Optimal investment strategy for annuity contracts under the constant elasticity of variance (CEV) model. Insurance: Mathematics and Economics, 45(1): 9–18.Google Scholar
Gupta, A. and Li, Z. (2007) Integrated optimal annuity planning with consumption-investment selections in retirement planning. Insurance: Mathematics and Economics, 41(1): 96–110.Google Scholar
Haberman, S. and Vigna, E. (2002) Optimal investment strategies and risk measures in defined contribution pension schemes. Insurance: Mathematics and Economics, 31(1): 3569.Google Scholar
Hobijn, B. and Lagakos, D. (2003) Social Security and the Consumer Price Index for the Elderly. Current Issues in Economics and Finance, Volume 9, No. 5, Federal Reserve Bank of New York, pp. 1–6. http://www.newyorkfed.org/research/current_issues/ci9-5.pdfGoogle Scholar
Hoevenaars, P. M. M., Molenaar, R. D. J., Schotman, P. C. and Steenkamp, T. B. M. (2008) Strategic asset allocation with liabilities: beyond stocks and bonds. Journal of Economic Dynamics and Control, 32: 29392970.CrossRefGoogle Scholar
Horneff, W. J., Maurer, R. H. and Stamos, M. Z. (2008) Lifecycle asset allocation with annuity markets. Journal of Economic Dynamics and Control, 32: 35903612.CrossRefGoogle Scholar
House of Representatives (2001) Consumer Price Index for the Elderly Consumers Act, H.R. 2035 (107th Congress).Google Scholar
Huang, H., Milevsky, M. A. and Wang, J. (2004) Ruined moments in your life: how good are the approximation. Insurance: Economics and Mathematics, 34: 421447.Google Scholar
Ibbotson, R. G., Milevsky, M. A., Chen, P. and Zhu, K. (2007) Lifetime Financial Advice: Human Capital, Asset Allocation and Insurance. CFA Institute, Charlottesville, VA, USA.Google Scholar
Jarrow, R. and Yildirim, Y. (2003) Pricing treasury inflation protected securities and relative derivatives using HJM model. Journal of Financial and Quantitative Analysis, 38(2): 337358.CrossRefGoogle Scholar
Jennings, W. (2006) Disaggregated TIPS: the case for disaggregating inflation linked bonds into bonds linked to narrower CPI components. Journal of Pension Economics and Finance, 5(3): 325343.CrossRefGoogle Scholar
Kothari, S. P. and Shanken, J. (2004) Asset allocation with inflation-protected bonds. Financial Analysts Journal, 60(1): 5470.CrossRefGoogle Scholar
Merton, R. C. (1971) Optimum consumption and portfolio rules in a continuous time model. Journal of Economic Theory, 3(December): 373413.CrossRefGoogle Scholar
Milevsky, M. A. (2006) The Calculus of Retirement Income: Financial Models for Life Insurance and Pension Annuities. Cambridge: Cambridge University Press.CrossRefGoogle Scholar
Milevsky, M. A. and Robinson, C. (2000) Self-annuitization and ruin in retirement. North American Actuarial Journal, 4(4): 112129.CrossRefGoogle Scholar
Moore, K. S. and Young, V. R. (2006) Optimal and simple: nearly optimal rules for minimizing the probability of financial ruin in retirement. North American Actuarial Journal, 10(4): 145162.CrossRefGoogle Scholar
Robinson, C. and Tahani, N. (2007) Sustainable retirement income for the socialite, the gardener and the uninsured. Atkinson College Working Paper. SSRN. Available at http://ssrn.com/abstract=989165 CrossRefGoogle Scholar
Wang, Y. (2009) Quantile hedging for guaranteed minimum death benefits. Insurance: Mathematics and Economics, 45(3): 449458.Google Scholar
Young, V. R. (2004) Optimal investment strategy to minimize the probability of lifetime ruin. North American Actuarial Journal, 8(4): 106126.CrossRefGoogle Scholar
Figure 0

Figure 1. Annualized gap: CPI-E inflation minus CPI-U inflation.

Figure 1

Table 1. Components of the CPI, wage-earners (W) versus urban (U) versus elderly (E)

Figure 2

Figure 2. Values of $1 invested ILBF* versus the nominal value of an adjusted $1 liability. *Based on the average returns of four ILBFs.

Figure 3

Figure 3. The fit. Changes in CPI-E versus changes in inflation- linked bond fund total return monthly periods from 12/02 to 11/08.

Figure 4

Table 2. Optimal allocation to ILBF assuming high (95%) correlation with personal inflation

Figure 5

Table 3. Optimal allocation to ILBF assuming low (40%) correlation with personal inflation

Figure 6

Table 4. Optimal allocation to ILBF

Figure 7

Figure 4. Ruin probability as function of wealth-to-liability ratio. (a) Impact of strategy. (b) Impact of age. (c) Allocation to CPI fund as function of wealth-to-liability ratio.

Figure 8

Table 5. Optimal allocation to ILBF (β=0.5)

Figure 9

Table 6. Optimal allocation to ILBF (β=1.5)

Figure 10

Figure 5. (a) CPI fund allocation as function of W/L ratio: high, medium and zero correlation. (b) Ruin probability as function of W/L ratio: 96%, 40% and 0% liability correlation.