Hostname: page-component-7bb8b95d7b-qxsvm Total loading time: 0 Render date: 2024-10-04T21:18:48.675Z Has data issue: false hasContentIssue false

Managing the shortfall risk of target date funds by overfunding

Published online by Cambridge University Press:  10 January 2024

Giovanni Barone Adesi
Affiliation:
Swiss Finance Institute at Università della Svizzera Italiana (USI), Institute of Finance, Lugano, Switzerland
Eckhard Platen
Affiliation:
School of Mathematical and Physical Sciences and Finance Discipline Group, University of Technology Sydney, Sydney, Australia
Carlo Sala*
Affiliation:
Department of Financial Management and Control, Universitat Ramon Llull, ESADE, Barcelona, Spain
*
Corresponding author: Carlo Sala; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Is it possible to achieve almost riskless, nonfluctuating investment payoffs in the long run, at a fraction of the traditional funding requirement, using equity investments? What is their shortfall risk? These questions are motivated by the need to increase yields, while limiting the variability of investment results. We show how to use contingent claims, denominated in units of a stock index, to achieve an almost riskless investment outcome. To control the risk of the proposed hedge portfolios, we introduce an overfunded scheme and show its reliability using bootstrapping. Results show that a modest amount of overfunding is an effective risk-management approach that brings the probability of not achieving the target to less than 1 percent. Our results are based on the use of the minimal market model and a change of numeraire. Robustness tests support their validity under different market specifications.

Type
Article
Copyright
Copyright © The Author(s), 2024. Published by Cambridge University Press

1. Introduction

Many investors and pensioners who save for retirement may benefit from the fact that, in the long-term, the equity market tends to outperform the fixed income market (see Dimson et al. (Reference Dimson, March and Staunton2003)).Footnote 1 At retirement, the investors often prefer to target almost riskless payouts, where payouts in units of the savings account are not affected by a severe market drawdown close to retirement. A natural question is, therefore, whether it is possible to use equities to achieve almost riskless payouts of units of a savings account in the long-term. The growing interest in this question is motivated by the presence of low real interest rates, which reduces the real yields of investment policies followed traditionally by insurance and pension funds. Target date fund models show that it is possible to invest in equities less expensively than under risk-neutral assumptions to achieve almost riskless payouts of units of a savings account in the long-term, as is typical for defined-benefit pension lump-sum payouts. Arguably, pension funds should aim to achieve highly predictable payout by means of dynamic investment policies. Therefore, they should reduce the high variability typical of equity investments. The equity component of pension and insurance funds currently leads to highly variable outcomes. As an example, the returns of pension funds varied significantly in European countries in 2019, ranging from 14 percent in the United Kingdom to 4 percent in Switzerland. Furthermore, the performance of such funds has varied as well, as reported in the second and third columns of Table 1. Such performance reflects the trade-off between risk (volatility) and expected return, which is roughly linear for optimal investment policies.

Table 1. Average performance values for the UK and Swiss pension systems for the period 2009–2019

Source: Moneyfacts UK Personal Pension Trends Treasury Report/Lipper Reports for the UK performance, and Swiss Pension Fund Studies 2019, Swiss Canto Ltd, for the Swiss performance.

In this paper, we describe a dynamic approach to manage the shortfall risk of a portfolio that uses equities and risk-free securities to produce almost riskless long-term payouts. The proposed approach has the primary goal of achieving a target outcome of one unit of the risk-free security over a long time horizon by controlling for its variability at maturity. This is achieved by overfunding an innovative basic investment plan introduced by Fergusson and Platen (Reference Fergusson and Platen2014) that invests in the stock market heavily at inception, when investors have time to enjoy its long-term growth, and also whenever the stock market falls. We call our proposed approach ‘the overfunded plan’, and compare it with the ‘basic investment plan’, which is not overfunded. Overfunding and minimum expected shortfall (MES)Footnote 2 investment strategies are two distinct approaches to managing the risk of shortfalls in target date funds (TDFs). While both strategies aim to mitigate the risk of falling short of a TDF's target, they differ in their approach. In fact, while overfunding partially compensates for the shortfall arising in the implementation of the policy of Fergusson and Platen (Reference Fergusson and Platen2014) over a finite time horizon, the MES strategy is an optimization approach that optimizes the expected shortfall over a decision variable. As such, the two approaches differ both theoretically and in their implementations. The proposed overfunded plan is a mathematically founded formalization of the popular concept of investing in risky assets at a young age, and then moving to the risk-free security over time. From a theoretical viewpoint, the approaches we discuss here may be linked to the fact that there are rational asset price bubbles in the market, as theoretically predicted in, for example, Loewenstein and Willard (Reference Loewenstein and Willard2000), Platen (Reference Platen2002), Fernholz et al. (Reference Fernholz, Karatzas and Kardaras2005) and Hugonnier (Reference Hugonnier2012), and empirically confirmed in, for example, Baldeaux et al. (Reference Baldeaux, Grasselli and Platen2015), Baldeaux et al. (Reference Baldeaux, Ignatieva and Platen2018) and Platen and Rendek (Reference Platen and Rendek2020). The benchmark approach, as presented by Platen (Reference Platen, Kohlmann and Tang2001), Platen (Reference Platen2002), and Platen and Heath (Reference Platen and Heath2010), allows a more general modeling world than the risk-neutral modeling world of Black and Scholes (Reference Black and Scholes1973) permits. This does not mean that there is any economically meaningful arbitrage under the benchmark approach. No portfolio can reach an infinite value in finite time. The pricing and hedging methodology we use has been only partially implemented in Fergusson and Platen (Reference Fergusson and Platen2014) and Fergusson and Platen (Reference Fergusson and Platen2023) in order to describe dynamic investment policies that target locally riskless payouts. In contrast to Fergusson and Platen (Reference Fergusson and Platen2014) and Fergusson and Platen (Reference Fergusson and Platen2023), we overfund the respective basic investment plan to achieve with a high probability the desired payoff at a given target date.

The desired final payoff (i.e., one unit of the risk-free security) becomes random when denominated in units of a stock index, which we call the benchmark. This fixed payoff denominated in the benchmark can then be viewed as a contingent claim, to be hedged dynamically in order to achieve the number of units of the stock index necessary to reach the targeted value at maturity. Intuitively, the investment portfolio initially allocates a relatively high weight to equities. By complying with the underlying model, the initially high weight is reduced progressively over time in line with the theoretical hedge ratio in favor of the risk-free security in order to achieve the desired target. The optimal policy is a dynamic hedging strategy that follows a similar logic to the Black and Scholes (Reference Black and Scholes1973) option hedging strategy, which has been used successfully to manage the risk for the payoff of short-term (usually up to three years) options.

Main drivers for the success of the proposed overfunded plan are the mean-reversion and the leverage effect, which are captured by the underlying model, and are missing under the Black–Scholes model. Supported by Summers (Reference Summers1986), Poterba and Summers (Reference Poterba and Summers1988), Fama and French (Reference Fama and French1988), and Spierdijk et al. (Reference Spierdijk, Bikker and van den Hoek2012), mean-reversion is the key feature in the Fergusson and Platen (Reference Fergusson and Platen2014) to model the long-term dynamics of the logarithm of the stock index. The presence of mean-reversion is debated in the literature. It is useful to remark that, while mean-reversion is necessary to asymptotically achieve the desired rate of return with certainty, this paper focuses on finite time horizons. Different robustness tests (Section 4.2) show that, while important, the mean-reversion is not the only driver of the good performance of our proposed overfunded model. Another important feature is the leverage effect, first documented by Black (Reference Black1976). There is an empirical negative correlation between returns and volatility, which needs to be accommodated in a suitable model, and implies that there is an asymmetric reaction of volatility to bad and good news and, hence, to negative or positive return innovations, respectively. Economically, a negative return implies a drop in the firm's equity value. This increases its leverage, which in turn, leads to higher volatility of equity returns.

Aside from the importance of properly modeling these key features, which indeed could potentially be achieved theoretically using various stock index models, it is mostly the strict supermartingale property of the risk-free security denominated in units of the stock index, captured under the employed model, that ensures most of the cost savings of the proposed strategies. This property is of crucial importance for the effective management of long-term assets. More precisely, the supermartingale property states that the expected payoff of the risk-free bond, measured in current units of the stock index, tends to decrease – on average – for increasing time horizons, which is modeled realistically under the model.

The strategy we implement here mimics payouts of the risk-free security at maturity. To provide a broad and general analysis, we study different investment plans, with maturities ranging from 10 to 100 years. As a first step, we follow Platen (Reference Platen, Kohlmann and Tang2001), Platen (Reference Platen2002), and Fergusson and Platen (Reference Fergusson and Platen2023), and present novel parametric investment plans. For all of the proposed investment plans, we describe how to dynamically set up the hedge portfolio by determining how much to invest in the stock index, depending on its value and time to maturity. Then, as a second step, we control the risk by overfunding the above strategies. The proposed overfunded plans are again based on the original basic investment plans. Thus, we invest the same fractions of the risky stock index as before, but have higher initial capital (the overfunded initial capital) invested in the hedge portfolio. This change carries relevant results, both theoretically and empirically, that could be important for whoever manages and invests in long-term assets (e.g., pension fund managers, life insurers, and other long-term investors).

In fact, overfunding the basic investment plan allows us to still benefit from the analytical tractability of the underlying model. Furthermore, we can overcome the well-known fact that it might not always be possible to achieve the targeted payoff, owing to model and hedge errors. Our empirical results refer to the US market and show that, using up to a century of data, and irrespective of the time horizon (from 10 to 100 years) and the initial time (from 1920 to 2020), an initial overfunding of 6 percent is enough to keep the probability of not achieving the target below 1 percent.

In theory, managing TDFs this way provides approximately the same benefits of a traditional plan, but at a fraction of the funding requirement. The benefit relative to traditional low-risk investment schemes is related to the equity risk-premium, and increases with the time horizon. It is modest at short maturities, but may become large at long maturities. Specifically, the benefit stems from the faster expected growth of the portion of funds invested in stocks. This is mitigated by the progressive, dynamic shift into the risk-free security, and determined theoretically by assuming appropriate dynamics for the stock index.

To bring the aforementioned theoretical results to fruition in the investment industry, it is necessary to evaluate their robustness and their margin of safety. As such, we test the basic investment plan and different overfunded schemes for 10-, 30-, 50-, and 100-year periods using bootstrapping based on US data for the period from 1920 to 2020. The long-term dynamic hedging strategy derived under the benchmark approach works similarly well, even for long time horizons close to 100 years. It appears that the proposed overfunded scheme can achieve the objectives of typical pension saving with the chosen margin of safety and less expensively than currently practiced.

The results presented in the paper are likely to raise questions for many readers which are only be partly answered in this paper and may require further research.Footnote 3 However, due to the potential of triggering a significant change in the practice of long-term asset management and providing a potential solution to many problems faced by social security, pension plans, and life insurance companies, the authors hope that the main message becomes disseminated and interest in the underlying theory raised. This should generate more interest and research on refined modeling and also on the related institutional aspects.

The remainder of the paper is organized as follows. In Section 2, we review the underpinning theory that allows us to manage and hedge the proposed investment strategies as a contingent claim. In Section 3, we show how to synthetically replicate a risk-free payoff using the stock index and the risk-free security, modeled as a roll-over short-term bond account. In Section 4, we present the overfunded plan and propose various alternatives. Then, we analyze which of the alternatives are best for short (10 year), medium (30 year), and very long (100 year) time horizons and consider different model robustness. Lastly, Section 5 concludes the paper.

2. Proposed investment strategy

To present our model as simply as possible, we denominate, at first, the securities in units of the risk-free security, a roll-over short-term bond account B t, and consider only how to approximately hedge a target payoff V T = 1 of one unit of the risk-free security using a total return stock index S t > 0 and the risk-free security B t = 1 for [0, T], following Fergusson and Platen (Reference Fergusson and Platen2023). This avoids the need to involve a stochastic interest rate process. However, similar dynamic trading strategies can be derived for targeted payoffs of a currency unit or units of a consumer price index. Perhaps more importantly, the proposed strategies can also be extended and generalized to other long-term contingent claims including options.

Under the minimal market model (MMM) introduced in Platen (Reference Platen, Kohlmann and Tang2001), the value S t of the stock index follows a squared Bessel process of dimension four; see (12). By assuming the MMM with linear deterministic market time, Fergusson and Platen (Reference Fergusson and Platen2023) synthetically replicate a risk-free security payoff. As a main drawback, the authors implement a parametric model with an assumed linear market time. As is common for many parametric models, this provides analytical tractability, but might not fit well with real data (e.g., see Campbell et al. (Reference Campbell, Lo and MacKinlay1997), Chapter 2). As a possible resolution to this issue, we propose an overfunded scheme that adds a buffer to the dynamic hedging strategy, thus relaxing the issues linked to the parametric nature of the model, and making the overfunded approach of practical use in the real world. In fact, the overfunded approach allows us to still work with non-numerically intensive closed-form solutions while, at the same time, reducing the probability of not achieving the targeted final value. Future capital, for example, one unit of the risk-free security, becomes a random quantity $\widehat{{V_T}} = 1/S_T$ when measured in units of a well-diversified stock index S T, the numeraire or benchmark. We call any security benchmarked when it is denominated in units of the benchmark. The benchmarked target $\widehat{V}_T$ can then be regarded as the payoff of a contingent claim, to be hedged over time in order to achieve its face value at maturity, where the hedge is investing dynamically in the benchmarked stock index $\widehat{{S_t}} = 1$ and the benchmarked risk-free security $\widehat{{B_t}}$. It is common practice for financial planners to recommend starting investing in risky securities at a young age, and then shifting to risk-free assets over time, where the respective glide path remains rather subjective. The proposed overfunded strategies make this glide path theoretically rigorous and practically feasible by performing the hedging with benchmarked securities, whereas risk-neutral hedging is based on securities denominated in the risk-free security.

The proposed policies invest in the stock index S t and a roll-over short-term bond account B t, which we call the risk-free security. Dividends are reinvested in the stock index. To keep our presentation simple, we discount the values by the risk-free security, and plan to achieve one unit of the risk-free security (i.e., the discounted value V T = 1) in, for example, 30 years (T = 30). For comparison, under the Law of One Price, a traditional pension policy that aims to produce one unit of the risk-free security would prescribe investing the present value of the desired payoff at the risk-free rate, which means it would simply buy and hold one unit of the risk-free security. However, as we mentioned at the beginning, the risk-free rate is, for economic reasons (on average), lower than the long-term average growth rate of a stock index. Thus, under classical assumptions, a pension that targets one unit of the risk-free security would have to buy and hold until maturity one unit of the risk-free security. However, this is rarely followed in practice because pension funds also invest in stocks. The practice of investing in stocks, even when one is targeting a risk-free payout, can be justified theoretically, as we describe in the following.

Under the benchmark approach, which generalizes the classical risk-neutral approach, one has the following fundamental result (see Theorem 10.3.1 in Platen and Heath (Reference Platen and Heath2010)): the expected value of the targeted payout V T = 1, when denominated in the units of the benchmark (the stock index S T) is, at all times t < T < ∞, equal to the investment portfolio value V t denominated in units of the benchmark S t. This is referred to as the real-world pricing formula, because the expectation is taken under the real-world probability measure. Therefore, we invest at time t < T the amount V t, where

(1)$$\widehat{{V_t}} = \displaystyle{{V_t} \over {S_t}} \ge {\rm {\opf E}}_t\left({\displaystyle{{V_T} \over {S_T}}} \right)$$

for all $0\> \le \>t\> \le \>T < \infty$. Here, ${\rm {\opf E}}_t\left( \cdot \right)$ denotes the conditional expectation under the physical or real-word probability measure, given the information available at time t. The quantity $\widehat{{V_t}}$ defined in equation (1) is the investment measured in units of the stock index, also called the benchmarked investment. The benchmarked risk-free security can be interpreted as a stochastic discount factor in the sense of Cochrane (Reference Cochrane2001), if normalized to 1.0 at the initial time. However, under the benchmark approach, the benchmarked risk-free security is not required to form a true martingale, which would be required for risk-neutral prices to coincide with the prices under the real-world pricing formula. This allows for the less expensive production of the targeted payout than possible under the classical risk-neutral assumption, where the benchmarked risk-free security is assumed to follow a martingale. The price S t of the benchmark, the stock index, fluctuates over time. Thus, we need to adjust dynamically, in a self-financing manner, the number of its units in the portfolio to keep equation (1) satisfied. If we can do that continuously over time, similarly to Black–Scholes option hedging, then we can replicate the benchmarked targeted value, exploiting the fact that:

(2)$$\mathop {\lim }\limits_{t\to T} {\rm {\opf E}}_t\left({\displaystyle{{V_T} \over {S_T}}} \right) = \displaystyle{{V_T} \over {S_T}}.$$

Therefore, by equation (1), the sequence of benchmarked holdings $\widehat{V}_t = V_t/S_t$ over time converges to ${\widehat{V}_T} = 1/S_T$. Thus, at maturity, we get $V_T = \widehat{V}_T\cdot S_T = 1$ unit of the roll-over short-term bond account, the risk-free security.

The policy sketched above exploits the fact that our desired final benchmarked payoff corresponds to the price of the risk-free security measured in terms of the stock index value. This allows for the modeling of the investment payout as a contingent claim, which may be hedged analogously, as in Black and Scholes (Reference Black and Scholes1973), to an option payoff.

This approach entails several important differences to the risk-neutral approach. First, under the proposed approach, every security is denominated in units of the stock index, and it aims in expectation at maturity at one unit of the risk-free security, denominated in units of the stock index. Thus, this approach uses the risky asset, rather than the risk-free security, as its numeraire. Risk management evolves then around the stock index, which on average and in the long-term, grows faster than the risk-free security, the traditional numeraire. Second, the chosen model for the dynamics of the stock index, the MMM, captures rather realistically its real-world mean-reversion and the leverage effect. Finally, the physical or real-world probability measure is the chosen pricing measure, whereas risk-neutral option pricing uses the putative risk-neutral probability measure. As a result, targeted payoffs are hedged less expensively than when using the traditional method, as we show later on in the empirical section of the paper. Note that we use the stock index as a proxy for the growth-optimal portfolio, which in the long run, is theoretically the pathwise best performing portfolio for the given stock investment universe, see Kelly (Reference Kelly1956) and Merton (Reference Merton1973), and equivalent to the numeraire portfolio, which was introduced in Long (Reference Long1990) and is central to the benchmark approach; see Platen and Heath (Reference Platen and Heath2010).

Under the real-world pricing formula (1), with our chosen numeraire (the stock index), the investment portfolio value necessary to achieve the desired targeted payoff follows a martingale (when denominated in units of the numeraire, the benchmark) under the physical or real-world probability measure. This means, its current benchmarked value is equal to its expected future benchmarked values. The observed typical dynamics of the risk-free security denominated in units of the stock index can be modeled realistically using the MMM by a strict supermartingale, where the current benchmarked value is strictly greater than the expected future benchmarked values. Assuming here a martingale, would be the risk-neutral assumption. The parametrization of the stock index dynamics through the MMM is preferable for capturing strict supermartingale dynamics and enables the achievement of the benchmarked final payoff at a lower cost. Fergusson and Platen (Reference Fergusson and Platen2023) show how the physical or real-world probability measure replaces the risk-neutral measure. In addition, they show how the strict supermartingale property permits the less expensive production of a payoff approximating one unit of the risk-free security than achieved by buying and holding one unit of the risk-free security. We show in Figure 1 the benchmarked trajectories of the risk-free security and the respective hedge portfolio that targets at maturity one unit of the risk-free security, both denominated in units of the stock index.Footnote 4

Figure 1. Risk-free security (in blue) and investment portfolio (thicker in red); both denominated in stock index.

In this figure, the risk-free security, denominated in units of the stock index, clearly shows on average a declining behavior. In contrast, the hedge portfolio value, denominated in units of the index, looks more like a martingale. Note that if we would use for comparison any risk-neutral model for the parametrization of the index, the risk-free security denominated in units of the stock index would have to be theoretically interpreted as a true martingale, which would make V t more expensive than under the MMM. We emphasize that this is true for any model that uses classical risk-neutral pricing. Under risk neutrality, by the Law of One Price, the theoretical value is equal to one unit of the risk-free security. From a hedging perspective, even in the presence of a working risk-neutral hedging scheme (in this case, a trivial buy and hold strategy), the risk-neutral price is equal to one unit of the risk-free security. This price is higher than that provided by the real-world pricing formula (1) under the MMM, owing to the strict supermartingale property of the benchmarked risk-free security.

3. Parametric synthetic bond strategy

In this section, we propose a parametric approach under the MMM, which we use to synthetically replicate a long-term bond with the discounted payoff V T = 1 by investing dynamically, in a self-financing manner, in a discounted stock index S t and the discounted risk-free security. Throughout this paper, we identify the outcomes of the parametric bond strategy as those of the basic investment plan. This is to differentiate the basic investment plan from the overfunded investment plans presented in the next section.

The desired final risk-free value V T = 1 is approached at time T because, under the MMM we describe below, the investment portfolio tends to glide gradually into the risk-free security when approaching the desired target. It is progressively reducing the fraction of the holdings in the stock index. This fraction is converging to zero, so that at maturity, one holds only units of the risk-free security. More details on this will follow. Fergusson and Platen (Reference Fergusson and Platen2023) provide an explicit formula for the proportion of the portfolio value to be invested in the index at any time. To do so, they assume the following stochastic differential equation (SDE) as the risk-free security discounted stock index process to follow under the MMM:

(3)$$dS_t = \alpha _tdt + \sqrt {S_t\alpha _t} dW_t, \;$$

where $\alpha _t = S_t\theta _t^{ 2} = \alpha \exp ( {\eta \cdot t} )$, with θ t denoting the volatility. Here, α > 0 is a normalization parameter, η > 0 is the (net) growth rate of S t, and W = {W t ≥ 0, t} is the driving canonical Brownian motion. More details about the MMM are described in Appendix A.

In the long-term, when using only constant volatility, as is the case in the Black–Scholes model with log-normal index increments, we would not capture the stock index dynamics sufficiently well. This is because the variance of the logarithm of its value would theoretically continue to increase linearly with time, even after long periods. This is not the case in reality, as can be seen in Figures 2 and 3. Here, we display the monthly log-closing prices and the monthly closing prices, respectively, of the discounted S&P500 total return index for the period from January 1920 until August 2020, denominated in units of the risk-free security. The risk-free security is here taken as the three-month US T-Bill roll-over account. All monthly data are obtained from Global Financial Data.

Figure 2. Discounted S&P500 total return log-index: time series of monthly closing log-prices for the period January 1920 to August 2020.

Figure 3. Discounted S&P500 total return index: time series of discounted monthly closing prices for the period January 1920 to August 2020.

Note that the quadratic variation [X.]t of a process X is approximated by the sum of the squares of its increments. The MMM can be conveniently fitted to historical stock index data, by exploiting the fact that the quadratic variation of the square root of the index, denoted by $\left[{\sqrt {S_.} } \right]_t$, satisfies the equation $\left[{\sqrt {S_.} } \right]_t = 0.25\alpha ( {e^{\eta t}-1} ) /\eta$ owing to the Ito formula. Thus, when α and η are estimated appropriately, theoretically under the MMM, the value $\tau _t = \ln \left({4\eta {\left[{\sqrt {S_.} } \right]}_t/\alpha + 1} \right)/\eta$ of the so-called market time should form an approximately straight line, with slope one. Figure 4 shows the resulting graph and the respective trendline, which obtained a maximum R 2 of 0.9933 when fitting the MMM for the period from January 1920 until August 2020, with α = 0.024 and η = 0.048, with standard errors equal to 0.000613 and 0.002968, respectively. Appendix B describes the estimation approach used to obtain these values. The rather high R 2 value indicates that the MMM parametrization makes good sense. Defined as the expected growth rate minus the short rate, the estimated net growth rate of η = 0.048 is aligned with Table 23 of Damodaran (Reference Damodaran2020).Footnote 5 We use these parameter estimates for the MMM in the remainder of the paper.

Figure 4. Resulting intrinsic time τ t and trendline under the MMM with R 2 of 0.9933 for the period January 1920 to August 2020.

As presented in Fergusson and Platen (Reference Fergusson and Platen2023), it follows that at any time t ∈ [0, T) under the MMM, we have the following discounted investment portfolio value:

(4)$$V_t = 1-{\rm exp}\left\{{-\displaystyle{{2\eta S_t} \over {\alpha ( {{\rm exp}\{ {\eta T} \} -{\rm exp}\{ {\eta t} \} } ) }}} \right\}.$$

Note that $\widehat{V}_T = ( {1/S_T} )$ is the benchmarked payoff. Here, the benchmarked value $\widehat{{V_t}} = V_t/S_t = {\rm {\opf E}}_t( 1/S_T)$ of the investment portfolio V t is its value denominated in units of the benchmark S t, the numeraire and stock index, and forms a martingale. In Figure 5, we show the theoretical discounted value of V t that pays in August 2020 one unit of the risk-free security (the continuous red line). Note that it fluctuates initially, similarly to the index, and approaches the risk-free security when close to maturity. Most importantly, its value is significantly less than one for many years, which makes it less expensive than the risk-free security. Only close to maturity is its value not much different to that of the risk-free security.

Figure 5. Theoretical discounted value of V t that, in August 2020, pays one unit of the risk-free security (red line) and the self-financing discounted hedge portfolio (blue line).

For illustration, recall Figure 1, which depicts the trajectories of the benchmarked risk-free security (in red) and the benchmarked investment portfolio (in blue). Note that the benchmarked risk-free security (red line) seems to fit a strict supermartingale, and the benchmarked investment portfolio (blue line) seems to fit a martingale.

The so-called ‘hedge ratio’, needed for the replication of the desired final target value is, similarly to the Black-Scholes hedge ratio for options, the partial derivative of the investment portfolio value with respect to changes in the stock index value

(5)$$\displaystyle{{\partial V_t} \over {\partial S_t}} = {\rm exp}\left\{{\displaystyle{{-2\eta S_t} \over {\alpha ( {{\rm exp}\{ {\eta T} \} -{\rm exp}\{ {\eta t} \} } ) }}} \right\}\displaystyle{{2\eta } \over {\alpha ( {{\rm exp}\{ {\eta T} \} -{\rm exp}\{ {\eta t} \} } ) }}.$$

Note that the strategy buys stocks when the index value S t falls, and vice versa. This is not what many investment funds do when the stock markets decline severely. However, this is optimal under the given model, and also rational when believing that the stock index will revert back to its average long-term growth behavior.

The classical Law of One Price assigns a discounted price of one to the payout of one unit of the risk-free security. Thus, in theory, managing pension funds the proposed way using hedging, provides with high probability the desired benefits but at a fraction of the traditional (risk-neutral based) funding requirement when targeting one unit of the risk-free security. The benefit relative to choosing a traditional pension scheme under classical risk-neutral assumptions increases with the length of the time horizon. It is modest at short maturities, but becomes substantial at maturities commonly of interest, such as 30 years. For example, Figure 5 shows an approximate value of 0.66 for January 1990, representing a saving of around 33 percent. This is a consequence of the faster expected growth of the portion of funds invested in stocks. Closer to maturity, the progressive shift out of the stock index into the risk-free security (equation (5)), necessary to achieve the risk-free target, reduces the growth of the investment portfolio. When the time horizon is extremely long, say 100 years, the saving increases to about 98 percent, as shown in Figure 5, with V 0 = 0.012 in January 1920.

The fraction π t = S t(dV t/dS t)/V t invested in the stock index at time t, which uses (12), is shown in Figure 6.

Figure 6. Fraction of stock index over time, defined as π t = S t(dV t/dS t)/V t.

Note that, initially, almost all wealth is invested in the index, before gradually shifting to the risk-free security. This appears to be consistent with the traditional popular financial planning advice, where one should invest in equities when young, and then shift to a risk-free security when getting close to retirement. Next, we create a self-financing discounted hedge portfolio, where we reinvest the wealth monthly according to the fraction π t. Figure 5 shows the value of the hedge portfolio (blue continuous line). Note that it evolves similarly to the theoretical value V t (red continuous line), with some minor hedge error at maturity. This dynamic asset allocation is very similar to the well-established risk-neutral hedging of options. A major difference is that it is required to perform successfully over a period of 100 years in this case, whereas most options are hedged over only a few years. In practice, this represents a major mitigation of risk over time, which may be widely exploited in industry to reduce production costs significantly. For illustration, we later examine the hedge error under the MMM for a period of up to 100 years.

Concerning the cost of the strategy, the institution that hedges according to the strategy is providing highly desirable liquidity to the market because it buys more units of the index when the index is falling and vice versa. This can be inferred by the number of units held in the index. Importantly, brokers are usually very supportive to such liquidity provider, and do not charge major transaction costs. As such, we do not add any frictions in our model to keep it mathematically cleaner and easier to follow. Clearly, we are not claiming that the proposed strategy is cost-less, just that its performance and the main message of the paper will not change due to its cost. It is also worth noticing that, from a regulatory and macroeconomic perspective, the strategy also ‘stabilizes’ the market dynamics, because when the index falls and many investors want to sell in a crash, the institution that hedges the bond creates demand.

4. The overfunded plan

The previously presented hedging scheme has been shown to work well under a variety of alternative diffusion models for the index S t by Gnoatto et al. (Reference Gnoatto, Grasselli and Platen2018). However, its implementation in the pension industry requires a careful assessment of its possible shortcomings. For example, an important drawback of Fergusson and Platen (Reference Fergusson and Platen2023) and Gnoatto et al. (Reference Gnoatto, Grasselli and Platen2018) is that their models rely on ad-hoc diffusion processes. Because parametric models are often not realistic, the obtained outcomes might lead to under/overestimated final hedging outcomes. To overcome this problem, we propose an overfunded plan that adds extra funds to the hedging portfolio. As a result, it is possible to control the underestimation risk of the final targeted payout. In this section, we set up the overfunded plan and explain empirically how it works.

The overfunded plan still aims to achieve a long-term payoff V T = 1 using a stock index and a risk-free security, but with a simple overfunding approach. For simplification purposes, we extract Brownian motion increments under the MMM, and bootstrap these to simulate alternative trajectories under the MMM. Let us introduce the discrete observation times t i = i/12, for i = 0, 1, …. We use the Euler scheme (see e.g., Kloeden and Platen (Reference Kloeden and Platen1992)) to discretize the SDE (12) in Appendix B. We do so for two reasons. First, we use it to extract the postulated Brownian motion increments under the MMM, shown in Figure 7, from the observed discounted index increments.

Figure 7. Increments of the MMM extracted Brownian motion: time series of monthly increments for the period January 1920 to August 2020.

Second, we use it to simulate alternative trajectories for the index using the bootstrapped Brownian motion increments. For the SDE (12), the Euler scheme is of the form

(6)$$S_{t_{i + 1}}-S_{t_i} = \alpha e^{\eta t_i}( {t_{i + 1}-t_i} ) + \sqrt {S_{t_i}\alpha e^{\eta t_i}} ( {W_{t_{i + 1}}-W_{t_i}} ) .$$

We rely on a simple, but effective bootstrap approach. However, more complex approaches are indeed applicable (e.g., the filtered historical simulation approach of Barone-Adesi et al. (Reference Barone-Adesi, Giannopoulos and Vosper1999)).

Recall that Figure 3 represents the monthly observed time series of the S&P500 total return index, denominated in units of a monthly rolled-over three-month US T-Bill portfolio, for the period from January 1920 until August 2020. Being a total return index, the time series is computed by reinvesting all dividends paid over time. We first obtain the bootstrapped time series by extracting the historical increments of the Brownian motion under the MMM; (see Figure 7). Then, we generate 10,000 Brownian motion-type trajectories by bootstrapping the extracted increments. Specifically, the historical Brownian motion increments are extracted to provide a pool of about 100⋅12 = 1,200 monthly Brownian motion increments. These are then randomly resampled with replacement to generate the 10,000 simulated (bootstrapped) Brownian motion-type trajectories for the observation period, January 1920 to August 2020. Needless to say, the number of simulations can be increased with almost no limit. These simulated Brownian motion-type trajectories are then used in the Euler scheme (6) to generate the trajectories of an index under the MMM that resemble the path of the discounted S&P500 total return index. The logarithms of 50 of these simulated trajectories are shown in Figure 8.

Figure 8. 50 trajectories of the bootstrapped logarithms of the S&P500 total return index for the period January 1920 to August 2020. Superimposed: the mean and the 5% and 95% confidence intervals.

Note the long-term mean-reversion of these trajectories. Using a straightforward generalization of the Black-Scholes model to bootstrap the observed S&P500 returns would not have generated such mean-reversion.

Following equation (10), we first compute the discounted investment portfolio value, for which only 50 trajectories are depicted for clarity in Figure 9.

Figure 9. Discounted investment portfolio value for the simulated total return index for the period January 1920 to August 2020, computed using equation (10). Superimposed: the mean and the 5% and 95% confidence intervals.

We set the initial value of the hedge portfolio equal to the initial value of the discounted investment portfolio (V 0 = 0.012). Then, following equation (5), and for each of the 10,000 simulated index trajectories, we compute the trajectory of the hedge portfolio in a self-financing manner until maturity, with a monthly reallocation, using the fraction invested in the index involving the hedge ratio (5). From January 1920 until August 2020, the time horizon of all simulations is fixed at 100 years. Indeed, to compute the hedge portfolio, we need to compute both the hedge ratio (dV t/dS t) and the fraction of the investment portfolio invested in the stock index over time (π t); see Figure 10.

Figure 10. Fraction of investment portfolio invested in the stock index over time for the period January 1920 to August 2020, defined as π t = S t(dV t/dS t)/V t. Superimposed: the mean and the 5% and 95% confidence intervals.

Once the hedge ratio and the quantity π t are known, we can compute the investment plan H t, which is defined as follows:

(7)$$H_{t_{i + 1}} = H_{t_i}\cdot \left[{1 + \pi_{t_i}\cdot \left({\displaystyle{{S_{t_{i + 1}}} \over {S_{t_i}}}-1} \right)} \right].$$

We perform the hedge simulation for all bootstrapped index scenarios; 50 trajectories are shown in Figure 11. We interpret this dynamic asset allocation as the basic investment plan. Note that most hedges end up rather close to 1.0 at maturity.

Figure 11. Bootstrapped basic investment plan for the simulated total return index for the period January 1920 to August 2020. Superimposed: the mean and the 5% and 95% confidence intervals.

Taken from another viewpoint, Figure 12 shows the histogram values of the 10,000 bootstrapped investment plans at maturity. As expected, the histogram is centered around one and strongly non-normal.

Figure 12. Histogram of the basic investment plan at maturity. The histogram captures the last value of the investment computed for the period January 1920 to August 2020.

Figures 11 and 12 both show that model errors and discretization errors in the hedge implementation for the proposed basic investment plan mean it is not possible to exactly achieve the target of one unit of the risk-free security at maturity. The main reasons for the hedging errors are that we have no continuous hedging and that the model is not perfect. To overcome this problem, we propose overfunding the basic investment plan. The overfunded plan is defined by following the investment strategy, with the same fractions as before. The only change is that we have a higher initial capital $V_0^{x} = xV_0$, for x > 1, invested in the hedge portfolio, which yields the following overfunded hedge portfolio value:

(8)$$V_t^{x} = xV_t, \;$$

for $0 \le t \le T$, where 100 ⋅ (x − 1) determines the percentage of the overfunding. Note that the overfunded plan simply multiplies the entire basic investment plan by x because the self-financing hedging portfolio evolves such that (8) holds. The returns remain the same, but the initial value is larger. Figure 13 depicts the original basic investment plan with no overfunding x = 1 (top left), and the three overfunded strategies, x 1 = 1.02 (top right), x 2 = 1.04 (bottom left), and x 3 = 1.06 (bottom right). This figure shows how almost all trajectories with x 3 = 1.06 lead to final values greater than the targeted payoff of one.

Figure 13. Bootstrapped overfunded plans: The four panels capture the basic investment plan (top left) with no overfunding (x = 1) and the three overfunded schemes, x 1 = 1.02 (top right), x 2 = 1.04 (bottom left), and x 3 = 1.06 (bottom right). Superimposed in all panels: the mean and the 5% and 95% confidence intervals.

Once more, and to better understand the final outcomes of the overfunded plans, we focus on their final payouts. Figure 14 shows (in clockwise order) histograms of the final values of the basic investment plan and the overfunded plans $V_T^{x} = xV_T$, with x 1, x 2, and x 3 representing an overfunding of 1.02, 1.04, and 1.06, respectively.

Figure 14. Histogram of the final values of the overfunded plans: The four panels capture $V_T^{x} = xV_T$ for the basic investment plan (top left) with no overfunding (x = 1) and the three overfunded plans, with x 1 = 1.02 (top right), x 2 = 1.04 (bottom left), and x 3 = 1.06 (bottom right).

To better understand the magnitude of the still remaining chance of underachieving, Table 2 collects the basic summary statistics of the final values of the overfunded plans, where the first column (overfunding weight x) defines the level of overfunding. Note how the arithmetic means of all schemes confirm equation (8), and the 25 percent and 75 percent quantiles together with the 99.9 percent percentile values confirm the increasing success of the overfunded plans over the basic plan. This is further emphasized by the rate of failure (Ratio), computed as the percentage of the total number of simulations that the plan does not reach the target of one, and the expected shortfall (ES), which accounts for the expected loss in units of the risk-free security when the target has not been achieved. While the basic investment plan misses the target about half the time (50.72%), producing an expected shortfall of 0.0083 units of the risk-free security, overfunding the basic investment plan quickly decreases both the rate of failure and the expected shortfall of the investment. More precisely, with 6 percent overfunding, both the probability of missing the target and the expected shortfall become almost negligible (0.87% and 0.0002, respectively).

Table 2. Summary statistics for the period January 1920 to August 2020 of the basic investment plan with no overfunding (x = 1), and for the three overfunded plans, with x 1 = 1.02, x 2 = 1.04, and x 3 = 1.06

All values are computed for the final value of the investment plan, defined as $V_T^{x} = xV_T$.

To test the robustness of the overfunded plan, we repeat the analysis for three different time horizons of the investment period. Specifically, we analyze three 30-year periods: 1920–1950, 1950–1980, and 1990–2020. For each study, we use the simulated 10,000 stock index trajectories, changing the starting times and maturities for the pricing and hedging accordingly. Note that all hedge portfolios are initialized using their respective theoretical bond values, which are determined fully by the simulated random value of the stock index at that time. This makes the initial value of the hedge portfolio quite random, an important effect that may impact the accuracy of the hedge. However, similarly to Black–Scholes option hedging, if the model is correct, the hedge should theoretically generate the targeted payout independently of the particular scenario for the stock index. Comparing Tables 2 and 3, the basic investment plan performs slightly worse over 30 years than over 100 years in terms of the percentage of failure, but with a lower expected shortfall. Furthermore, the percentage of failure is almost the same for each of the three periods and close to 60 percent. However, the magnitude of the expected shortfall decreases when the stock index value at the initialization of the hedge becomes more random. This result can be explained by the fact that there is more averaging over the stationary density available to diversify over-time the stock index value. Interestingly, overfunding the basic investment plan by 6 percent is again enough to make the likelihood of failure and the expected loss almost negligible.

Table 3. Summary statistics for three 30-year periods of the basic investment plan with no overfunding (x = 1), and of the three overfunded plans with x 1 = 1.02, x 2 = 1.04, and x 3 = 1.06

All values are computed for the final value of the investment plan, defined as $V_T^{x} = xV_T$.

A similar trend is confirmed by keeping the maturity date fixed as August 2020 and considering plans of 10, 20, 30, and 50 years. Table 4 reports the summary statistics of the basic investment plan and the overfunded plans, where the plans end in August 2020 and start in January 2010 (10-year plan), January 2000 (20-year plan), and January 1970 (50-year plan).

Table 4. Summary statistics for 10-, 20-, and 50-year periods of the basic investment plan with no overfunding (x = 1), and of the three overfunded plans with x 1 = 1.02, x 2 = 1.04, and x 3 = 1.06

All values are computed for the final value of the investment plan, defined as $V_T^{x} = xV_T$. The periods end in August 2020 and start in January 2010, January 2000, and January 1970, respectively.

From Table 4, interesting results emerge. First, the 6 percent overfunding rule still applies. Almost irrespective of the time length and period of observations, 6 percent overfunding appears to be high enough to achieve about 99 percent of the targeted value. Second, as the time window widens, the likelihood of failure and the expected shortfall both appear to increase. This is mostly the result of hedging errors that tend to build up over time. However, note that this trend is not monotonic. In fact, as soon as the investment plan covers more than 60 years, a time diversification effect is evident, as demonstrated by the lower failure rate documented for the 100-year analysis presented in Table 2. Overall, the number of positive outcomes increases with the length of the time window. However, this is only true for plans longer than 60 years.Footnote 6

4.1 Model robustness

Given the set-up presented in this paper, an investor is confronted with some basic problems, that is: what happens to the performances of the model in the presence of some model uncertainty? What happens if the model parameters are mis-estimated? For which time horizons is the model more robust? In this subsection we provide answers to these questions and show the robustness of the basic model and the overfunded models. Only results for underestimated parameters are presented because overestimated parameters improve convergence, at the cost of a higher initial investment.

As a first test, we consider the case faced by an investor who estimates the model parameters with an error. In such a case the investor builds the portfolio and hedges it with the net growth rate η = 0.048, which turns out to be different to those we use for the simulation of the testing data. The simulated testing data are generated with a misspecified version of the Euler scheme presented by equation (6):

(9)$$S_{t_{i + 1}}-S_{t_i} = \alpha ^{{\ast} }e^{\eta ^{{\ast} }t_i}( {t_{i + 1}-t_i} ) + \sqrt {S_{t_i}\alpha ^{ {\ast} }e^{ \eta ^{{\ast} }t_i}} ( {W_{t_{i + 1}}-W_{t_i}} ) .$$

where α* and η* are the perturbed versions of the normalization factor and the net growth rate, respectively.

To better disentangle the role of the misspecification of parameters we modify only one parameter at a time, and we perturb it by a multiple n = 1, 2, …, n < ∞ of their estimated standard error. We start with a misspecified normalization parameter $\alpha ^{{\ast} } = \alpha -( {\varepsilon^{ \alpha }\cdot n} )$, where $\varepsilon ^{\alpha }$ is the standard error of the normalization factor α and n the number of standard errors considered. This first test shows an extreme model robustness for the 100 years 6 percent overfunded plan, as it achieves the target for n = 7. Interestingly, the model becomes more robust as we shorten the time horizons. With a time horizon of 30 years, which is a more realistic time horizon for an investor who saves for her pension, the 6 percent overfunded plan achieves its target even for n = 10 (while the basic plan is 5% off target).

Next, we keep α fixed and misspecify the net growth rate $\eta ^{ {\ast} } = \eta -( {\varepsilon^{\eta }\cdot n} )$, where $\varepsilon ^{ \eta }$ is the standard error of the growth factor η. Empirical results show that, with a one standard error perturbation, the 6 percent overfunded plan still achieves the target at both 100 and 30 years, while the basic plan is 7 percent off target. Interestingly, the 30 years plan is the one that shows more robustness, as with n = 2 the overfunded plan is just 1.5 percent off target.

With an additional robustness test we check what happens if one misspecifies the parameters α and η in both the portfolio construction and hedging. Such a test requires the mis-estimation of

(10)$$V_t = 1-{\rm exp}\left\{{-\displaystyle{{2\eta^{ {\ast} }S_t} \over {\alpha^{ {\ast} }( {{\rm exp}\{ {\eta^{ {\ast} }T} \} -{\rm exp}\{ {\eta^{ {\ast} }t} \} } ) }}} \right\}$$

and the relative passages necessary to build and hedge the portfolio.

Parameters, α* = α − (α ⋅ y) and η* = η − (η ⋅ y) are once more the original parameters, this time perturbed by a percentage y. Once more, the model shows a strong robustness to the normalization parameter α, which still achieves the target with a misspecification of the parameter of y = −5%. For the growth rate, the robustness of the 6 percent overfunded model remains very high up to y = −1%, and it decreases as y increases. Remarkably, and once more, the 30 years plan shows more model robustness than those for most other time horizons.

4.2 The impact of mean-reversion and risk-premium for the overfunding model

Among the tests performed in the previous section we analyzed what happens to the overfunding model in the presence of parameters mis-estimation. Among others, we have tested the robustness of the model in the presence of a contaminated growth factor. In this section, we dig even further into the role of the risk-premium and we shed light on the role of the mean-reversion on the final performance of the overfunding model. We add some robustness simulations where we first eliminate the mean-reversion from our simulations, then we simulate our model with the MMM but only half of the risk-premium, and finally we mix these two effects, thus producing a simulation without mean-reversion and half of the risk-premium.

In particular, we examine the effect of making wrong assumptions on the stochastic process generating our simulated returns. Hedge ratios are computed under the same parameter values we used previously, as if the process generating the data had not changed, but the new simulation data are generated with (i) half of the risk-premium, (ii) from a geometric Brownian motion and (iii) with one half of the risk-premium applied to a Brownian motion in the third panel. In other words, we are hedging incorrectly to reflect our ignorance of the true process. We focus our analysis on finite horizons only, always controlling the risk through overfunding.

For the first test, we simply set η = 0.048/2, which is half of the value used throughout the paper and that we obtained with the numerical approach presented in Appendix A for the period from January 1920 until August 2020. Results of the simulated outcomes are presented in Table 5.

Table 5. Summary statistics for the period January 1990 to August 2020 of the basic investment plan with no overfunding (x = 1), and for the three overfunded plans, with x 1 = 1.02, x 2 = 1.04, and x 3 = 1.06

All values are computed for the final value of the investment plan, defined as $V_T^{ x} = xV_T$ and applying a growth factor equal to η = 0.048/2.

As it emerges from Table 5, the proposed approach remains robust with respect to the effect of a much lower risk-premium, which could be mitigated through a modest overfunding.

To remove the mean-reversion and generate the simulated returns from a geometric Brownian motion, we first need to estimate the parameters entering the Black–Scholes dynamics. We first estimate the volatility of the MMM realized log-returns at 10-year periods. The value of 17 percent for σ is aligned with the literature concerning the annual volatility of the S&P 500. It is theoretically grounded, as the average multi-yearly volatility of returns of the MMM decreases over time because of long-term mean-reversion from 17 percent at 10-year period to 11 percent at 30-year period. This quantity is used as a volatility input (σ) to simulate returns through a geometric Brownian motion:

(11)$$S_{t_{i + 1}} = S_{t_i}exp^{ \left[{\left({\mu -{{\sigma^{ 2}} \over 2}} \right)( {t_{i + 1-t_i}} ) + ( {W_{t_{i + 1}}-W_{t_i}} ) } \right]}.$$

with a drift parameter μ = η + σ 2/2, η = 0.048 and W denoting the stochastic differential of the Brownian motion. The portfolio analysis is then performed as if the simulated returns would come from an MMM dynamics. Results are depicted in Table 6.

Table 6. Summary statistics for the period January 1990 to August 2020 of the basic investment plan with no overfunding (x = 1), and for the three overfunded plans, with x 1 = 1.02, x 2 = 1.04, and x 3 = 1.06

All values are computed for the final value of the investment plan, defined as $V_T^{x} = xV_T$ and generated by a Black–Scholes dynamics.

The effect of removing the mean-reversion from the simulated returns has a minor impact on all variables, and the overall strategy remains robust and an investor can still achieve its target by slightly increasing the overfunding.

Finally, we repeat the same exercise, with a reduced risk-premium, so that η = 0.048/2 and μ = η + σ 2/2. Results are depicted in Table 7.

Table 7. Summary statistics for the period January 1990 to August 2020 of the basic investment plan with no overfunding (x = 1), and for the three overfunded plans, with x 1 = 1.02, x 2 = 1.04, and x 3 = 1.06

All values are computed for the final value of the investment plan, defined as $V_T^{ x} = xV_T$ and generated by a Black–Scholes dynamics and η = 0.048/2.

As from Table 7 the impact of the mean-reversion and the risk-premium are not of major importance, as the overall strategy remains robust.

The intuition behind the above results is that mean-reversion is hard to detect on individual stock paths of finite length. Therefore, it is not surprising that its effects on the proposed strategies are modest at a 30-year time horizons. Overall, we confirm that the role of mean-reversion is potentially very important in the very long run (Pastor and Stambaugh (Reference Pastor and Stambaugh2012)), but much less on the thirty year horizon on which we focus our empirical analysis. It is also worth noticing that our asset allocation tends to reduce stock allocation in later years, reducing the effects of long-term uncertainty in stock parameters. Likewise, also the size risk-premium has just a minor effect on the proposed overfunding scheme, and both the mean-reversion and the size risk-premium are not driving the performances of the overfunding scheme and for more risk averse investors, their effect can be alleviated through a minimal extra overfunding.

5. Conclusion

The strategy we discussed promises to achieve the long-term objectives of traditional pension and target date investment policies at a fraction of their traditional theoretical initial cost. Its practical implementation requires an assessment of its long-term reliability and an evaluation of their safety margin. To fulfill these prudential requirements, the proposed overfunded scheme allows for the choice and the evaluation of the desired safety level. Our results show that an overfunding of 6 percent achieves the desired final value more than 99 percent of the time, with an expected shortfall not greater than 0.2 percent of the targeted value. To corroborate the robustness of our model, similar results are achieved when changing the market scenario by adding market uncertainty, model misspecification, lower risk-premium, and removing the mean-reversion. Therefore, appropriate overfunding appears to allow for the safe use of this strategy in the investment industry. With a wider use of this methodology, a similar gain in the efficiency of long-term risk management for pensions and other long-term investments may become achievable.

Our aim is the introduction of a simple framework for the management of shortfall risk in TDFs. Additional capital or policy changes may be required to improve the robustness of the proposed scheme in order to achieve the desired outcomes with even higher confidence. Where appropriate, remedial measures must be introduced to address potential pitfalls. The effects of randomness in volatility, price jumps, interest rates (when cash is targeted as the payout), and the equity risk-premium on the performance of the proposed policy should also be analyzed. Particular attention should be given to the additional funding requirements under a variety of stress conditions. Finally, the effect of a large-scale implementation on the overall stability of financial markets should be evaluated. As such, although we provide an initial basic overfunded approach to implement a pension scheme, many questions remain for further investigation.

Acknowledgements

Comments by Stephen Brown and the financial support from the grant PID2019-106465GB-I00, funded by MCIN/AEI/10.13039/501100011033 and the AGAUR – SGR 2017-640 grant are gratefully acknowledged.

A. The minimal market model (MMM)

As defined in the main text, we model the stock index dynamics as:

(12)$$dS_t = \alpha _tdt + \sqrt {S_t\alpha _t} dW_t, \;$$

where $\alpha _t = S_t\theta _t^{2} = \alpha \exp ( {\eta \cdot t} )$, with θ t denoting the volatility. Here, α > 0 is a normalization parameter, η > 0 is the (net) growth rate of S t, and W = {W t ≥ 0, t} is the driving canonical Brownian motion.

This is the MMM introduced by Platen (Reference Platen, Kohlmann and Tang2001) and described in Platen and Heath (Reference Platen and Heath2010), where S t models the growth optimal portfolio as a time-transformed squared Bessel process of dimension four. This yields the ergodic Cox–Ingersoll–Ross process Y t = S t/α t when normalized, which is linearly mean-reverting. The index conditioned on previous values follows a noncentral chi-square distribution with four degrees of freedom with respect to market time, which depends on the value of the index and keeps the index going up and down but never reaches zero.

Note that the volatility $\theta _t = \sqrt {\alpha _t/S_t}$ with respect to market time is equal to $1/\sqrt {Y_t}$, which generates the well-known leverage effect. Thus, S t = Y tα t is modeled as the product of the mean-reverting process Y t with stationary density, and an exponential function α t of market time. It is important that we realistically capture with these dynamics the leverage effect and the mean-reversion, which are both typical for the long-term dynamics of a stock index.

The mean-reversion is clearly visible for the logarithm of the stock index shown in Figure 2 and it is difficult to imagine a reason for not assuming that this mean-reversion with respect to market time will continue as long as the market exists. Of course, the mean-reversion is not too important when modeling over shorter time periods. The key property that makes the proposed approach working under the model is the strict supermartingale property of the risk-free security when denominated in units of the stock index.

B. Parameter estimation under the MMM

In this appendix we present the maximum likelihood estimation approach used to estimate the MMM parameters in the paper. Specifically, the goal is to find the parameters α and θ of the SDE presented in equation (12) that maximize the occurrence of the observations under the hypothesis that the MMM holds true. To achieve it we first need to define the transition density function of the discounted stock index (see Revuz and Yor (Reference Revuz and Yor1999)):

(13)$$p_S( {t, \;x_t, \;T, \;x_T} ) = \displaystyle{1 \over {2( {\varphi_T-\varphi_t} ) }}\sqrt {\displaystyle{{x_T} \over {x_t}}} \exp \left({-\displaystyle{{x_t + x_T} \over {2( {\varphi_T-\varphi_t} ) }}} \right)I_t\left({\displaystyle{{\sqrt {x_Tx_t} } \over {\varphi_T-\varphi_t}}} \right)$$

where $I_\upsilon ( z ) = \sum _{m = 0}^{\infty } ( {{( {1/2z} ) }^{ 2m + \upsilon }/m!\Gamma ( {m + 1 + \upsilon } ) } )$ is the modified Bessel function of the first kind with index υ and φt = 1/4η(α(exp (ηt) − 1)) is the quadratic variation of $\sqrt S$.

Given a finite time series t 0 < t 1 < ⋅ ⋅ ⋅ < t n < ∞ of observations of the discounted index $S_{t_0}, \;\;S_{t_1}, \;\;\ldots , \;\;S_{t_n}$ and using the above defined transition density function, we maximize the logarithm of the likelihood function:

(14)$$l( \alpha , \;\eta ) = $$
(15)$$ = \sum\limits_{i = 1}^n {\log } \left({\displaystyle{1 \over {2( \varphi_{t_i}-\varphi_{t_i-1}) }}\sqrt {\displaystyle{{S_{t_i}} \over {S_{t_i-1}}}} \exp \left({-\displaystyle{{S_{t_i-1} + S_{t_i}} \over {2( \varphi_{t_i}-\varphi_{t_i-1}) }}} \right)I_1\left({\displaystyle{{\sqrt {S_{t_i}S_{t_i-1}} } \over {\varphi_{t_i}-\varphi_{t_i-1}}}} \right)} \right)$$
(16)$$ = \sum\limits_{i = 1}^n {\left\{{\log \left({\displaystyle{1 \over {2( \varphi_{t_i}-\varphi_{t_i-1}) }}} \right) + \displaystyle{1 \over 2}\log \left({\displaystyle{{S_{t_i}} \over {S_{t_i-1}}}} \right) + \left({-\displaystyle{{S_{t_i-1} + S_{t_i}} \over {2( \varphi_{t_i}-\varphi_{t_i-1}) }}} \right) + \log \left({I_1\left({\displaystyle{{\sqrt {S_{t_i}S_{t_i-1}} } \over {\varphi_{t_i}-\varphi_{t_i-1}}}} \right)} \right)} \right\}} .$$

To initialize the computation we need to find the initial estimates of both α and θ, which can be found by equating the theoretical quadratic variation of $\sqrt S$:

(17)$$\langle \sqrt S \rangle _{t_j} = \displaystyle{{1}\over {4\eta }} \alpha({\exp ( \eta t_j) -1})$$

with its empirical counterpart

(18)$$\langle \sqrt S \rangle _{t_j} \approx \mathop \sum \limits_{\sigma + j}^{i = 1} \left( {\sqrt {S_{t_i}} -\sqrt {S_{t_{i-1}}} } \right)^{ 2}.$$

For t = t k and t = t 2k, where $k = \lfloor + ; \;\;n/2\rfloor$, the initial estimates of the parameters are defined as:

(19)$$\alpha _0 = \langle \sqrt S \rangle _{t_k}\displaystyle{{4\eta } \over {\exp ( \eta t_k) -1}}$$

and

(20)$$\eta _0 = \log \displaystyle{{{\langle \sqrt S \rangle }_{t_{2k}}/{\langle \sqrt S \rangle }_{t_k}-1} \over {t_k-t_0}}.$$

The first iteration to find the logarithm of the likelihood function is computed using the above defined initial estimates at points α 0 + iδα and η 0 + jδη for i, j = −2,  − 1, 0, 1, 2, δα = α 0/4 and δη = η 0/4 and fitting the quadratic form

(21)$$Q( x ) = x^{T}Ax-2b^{T}x + c$$

where x = (α/η) is a 2-by-1 vector that contains the parameters of interest, A is a 2-by-2 negative definite matrix, b is a 2-by-1 vector and c a scalar. After the algorithm is initialized, the subsequent values (α 1 and η 1) are results of the matrix expression A −1b, which corresponds to the maximum value of the quadratic form. The algorithm stops when it finds the maximum likelihood estimates of both α and η, obtained by applying the Newton-Raphson root-finding method to the first order partial derivatives of the log-likelihood function.

Finally, we obtain the standard errors of the estimated parameters from the Cramer-Rao inequality for the covariance matrix, as this provides the lower bound for the variance of an unbiased estimator of a parameter as the number of observations increases,

(22)$$VAR( \alpha , \;\eta ) \ge -\displaystyle{1 \over {{\rm \bigtriangledown }^2l( \alpha , \;\eta ) }}.$$

Footnotes

1 The outperformance refers to the empirically evidence that the long-term average growth rate of the stock market index is larger than the average interest rate. There seems also to be no rational that would suggest that this empirical fact may reverse.

2 See, e.g.: Acharya et al. (Reference Acharya, Engle and Richardson2012), Brownlees and Engle (Reference Brownlees and Engle2016) for some market data application and Bertsimas et al. (Reference Bertsimas, Lauprete and Samarov2004) for a review of the theoretical properties of expected shortfall.

3 The theory underpinning the basic investment strategy we implement is described in the monograph Platen and Heath (Reference Platen and Heath2010).

4 All details on how to estimate both quantities are presented in Section 3.

5 In the 11th edition of his survey on the risk-premium, using different approaches, Damodaran (Reference Damodaran2020) presents different estimates for the United States. All are between 3.20 percent and 5.20 percent.

6 For ease of space, the graphical outputs of the respective cited simulations have been omitted and are available upon request.

References

Acharya, V, Engle, R and Richardson, M (2012) Capital shortfall: a new approach to ranking and regulating systematic risk. American Economic Review 10, 5964.CrossRefGoogle Scholar
Baldeaux, J, Grasselli, M and Platen, E (2015) Pricing currency derivatives under the benchmark approach. Journal of Banking & Finance 53, 3448.CrossRefGoogle Scholar
Baldeaux, J, Ignatieva, K and Platen, E (2018) Detecting money market bubbles. Journal of Banking & Finance 87, 369379.CrossRefGoogle Scholar
Barone-Adesi, G, Giannopoulos, K and Vosper, L (1999) Backtesting derivative portfolios with filtered historical simulation (FHS). Journal of Futures Market 8, 3158.Google Scholar
Bertsimas, D, Lauprete, GJ and Samarov, A (2004) Shortfall as a risk measure: properties, optimization and applications. Journal of Economic Dynamics & Control 28, 13531381.CrossRefGoogle Scholar
Black, F (1976) Studies of stock price volatility changes. Proceedings of the 1976 Meetings of the American Statistical Association, pp. 171–181.Google Scholar
Black, F and Scholes, M (1973) The pricing of options and corporate liabilities. Journal of Political Economy 81, 637654.CrossRefGoogle Scholar
Brownlees, C and Engle, RF (2016) Srisk: a conditional capital shortfall measure of systemic risk. The Review of Financial Studies 30, 4879.CrossRefGoogle Scholar
Campbell, Y, Lo, A and MacKinlay, A (1997) The Econometrics of Financial Markets. Princeton, NJ: Princeton University Press.CrossRefGoogle Scholar
Cochrane, J (2001) Asset Pricing. Princeton: Princeton University Press.Google Scholar
Damodaran, A (2020) Equity Risk Premiums: Determinants, Estimation and Implications – The 2020 Edition. NYU Stern School of Business. http://dx.doi.org/10.2139/ssrn.3550293Google Scholar
Dimson, E, March, P and Staunton, M (2003) Triumph of the optimists: 101 years of global investment returns. Journal of Pension Economics and Finance 2, 9195.Google Scholar
Fama, EF and French, KR (1988) Permanent and temporary components of stock prices. Journal of Political Economy 96, 246273.CrossRefGoogle Scholar
Fergusson, K and Platen, E (2014) Hedging long-dated interest rate derivatives for Australian pension funds and life insurers. Australian Journal of Actuarial Practice 1, 2944.Google Scholar
Fergusson, K and Platen, E (2023) Less-expensive long-term annuities linked to mortality, cash and equity. Annals of Actuarial Science 17, 170207.CrossRefGoogle Scholar
Fernholz, R, Karatzas, I and Kardaras, C (2005) Diversity and relative arbitrage in equity markets. Finance and Stochastics 9, 127.CrossRefGoogle Scholar
Gnoatto, A, Grasselli, M and Platen, E (2018) A penny saved is a penny earned: less expensive zero coupon bonds. ArXiv 1608.04683. https://arxiv.org/abs/1608.04683Google Scholar
Hugonnier, J (2012) Rational asset pricing bubbles and portfolio constraints. Journal of Economic Theory 147, 22602302.CrossRefGoogle Scholar
Kelly, JL (1956) A new interpretation of information rate. Bell System Technical Journal 35, 917926.CrossRefGoogle Scholar
Kloeden, PE and Platen, E (1992) Numerical Solution of Stochastic Differential Equations. Berlin, Heidelberg: Springer.CrossRefGoogle Scholar
Loewenstein, M and Willard, GA (2000) Local martingales, arbitrage, and viability. Economic Theory 16, 135161.CrossRefGoogle Scholar
Long, JB (1990) The numeraire portfolio. Journal of Financial Economics 26, 2969.CrossRefGoogle Scholar
Merton, RC (1973) Theory of rational option pricing. Bell Journal of Economics and Management Science 4, 141183.Google Scholar
Pastor, L and Stambaugh, RF (2012) Are stocks really less volatile in the long run? Journal of Finance 67, 431478.CrossRefGoogle Scholar
Platen, E (2001) A minimal financial market model. In Kohlmann, M and Tang, S (eds), Mathematical Finance. Trends in Mathematics. Basel: Birkhäuser, pp. 293301. https://doi.org/10.1007/978-3-0348-8291-0_27.CrossRefGoogle Scholar
Platen, E (2002) Arbitrage in continuous complete markets. Advances in Applied Probability 34, 540558.CrossRefGoogle Scholar
Platen, E and Heath, D (2010) A Benchmark Approach to Quantitative Finance. Berlin, Heidelberg: Springer.Google Scholar
Platen, E and Rendek, R (2020) Approximating the growth optimal portfolio and stock price bubbles. International Journal of Theoretical and Applied Finance 23, 2050048.CrossRefGoogle Scholar
Poterba, JM and Summers, LH (1988) Mean reversion in stock prices: evidence and implications. Journal of Financial Economics 22, 2759.CrossRefGoogle Scholar
Revuz, D and Yor, M (1999) Continuous Martingales and Brownian Motion, Volume 3rd Edn. Berlin, Heidelberg: Springer.CrossRefGoogle Scholar
Spierdijk, L, Bikker, JA and van den Hoek, P (2012) Mean reversion in international stock markets: an empirical analysis of the 20th century. Journal of International Money and Finance 31, 228249.CrossRefGoogle Scholar
Summers, LH (1986) Does the stock market rationally reflect fundamental values? Journal of Finance 41, 591601.CrossRefGoogle Scholar
Figure 0

Table 1. Average performance values for the UK and Swiss pension systems for the period 2009–2019

Figure 1

Figure 1. Risk-free security (in blue) and investment portfolio (thicker in red); both denominated in stock index.

Figure 2

Figure 2. Discounted S&P500 total return log-index: time series of monthly closing log-prices for the period January 1920 to August 2020.

Figure 3

Figure 3. Discounted S&P500 total return index: time series of discounted monthly closing prices for the period January 1920 to August 2020.

Figure 4

Figure 4. Resulting intrinsic time τt and trendline under the MMM with R2 of 0.9933 for the period January 1920 to August 2020.

Figure 5

Figure 5. Theoretical discounted value of Vt that, in August 2020, pays one unit of the risk-free security (red line) and the self-financing discounted hedge portfolio (blue line).

Figure 6

Figure 6. Fraction of stock index over time, defined as πt = St(dVt/dSt)/Vt.

Figure 7

Figure 7. Increments of the MMM extracted Brownian motion: time series of monthly increments for the period January 1920 to August 2020.

Figure 8

Figure 8. 50 trajectories of the bootstrapped logarithms of the S&P500 total return index for the period January 1920 to August 2020. Superimposed: the mean and the 5% and 95% confidence intervals.

Figure 9

Figure 9. Discounted investment portfolio value for the simulated total return index for the period January 1920 to August 2020, computed using equation (10). Superimposed: the mean and the 5% and 95% confidence intervals.

Figure 10

Figure 10. Fraction of investment portfolio invested in the stock index over time for the period January 1920 to August 2020, defined as πt = St(dVt/dSt)/Vt. Superimposed: the mean and the 5% and 95% confidence intervals.

Figure 11

Figure 11. Bootstrapped basic investment plan for the simulated total return index for the period January 1920 to August 2020. Superimposed: the mean and the 5% and 95% confidence intervals.

Figure 12

Figure 12. Histogram of the basic investment plan at maturity. The histogram captures the last value of the investment computed for the period January 1920 to August 2020.

Figure 13

Figure 13. Bootstrapped overfunded plans: The four panels capture the basic investment plan (top left) with no overfunding (x = 1) and the three overfunded schemes, x1 = 1.02 (top right), x2 = 1.04 (bottom left), and x3 = 1.06 (bottom right). Superimposed in all panels: the mean and the 5% and 95% confidence intervals.

Figure 14

Figure 14. Histogram of the final values of the overfunded plans: The four panels capture $V_T^{x} = xV_T$ for the basic investment plan (top left) with no overfunding (x = 1) and the three overfunded plans, with x1 = 1.02 (top right), x2 = 1.04 (bottom left), and x3 = 1.06 (bottom right).

Figure 15

Table 2. Summary statistics for the period January 1920 to August 2020 of the basic investment plan with no overfunding (x = 1), and for the three overfunded plans, with x1 = 1.02, x2 = 1.04, and x3 = 1.06

Figure 16

Table 3. Summary statistics for three 30-year periods of the basic investment plan with no overfunding (x = 1), and of the three overfunded plans with x1 = 1.02, x2 = 1.04, and x3 = 1.06

Figure 17

Table 4. Summary statistics for 10-, 20-, and 50-year periods of the basic investment plan with no overfunding (x = 1), and of the three overfunded plans with x1 = 1.02, x2 = 1.04, and x3 = 1.06

Figure 18

Table 5. Summary statistics for the period January 1990 to August 2020 of the basic investment plan with no overfunding (x = 1), and for the three overfunded plans, with x1 = 1.02, x2 = 1.04, and x3 = 1.06

Figure 19

Table 6. Summary statistics for the period January 1990 to August 2020 of the basic investment plan with no overfunding (x = 1), and for the three overfunded plans, with x1 = 1.02, x2 = 1.04, and x3 = 1.06

Figure 20

Table 7. Summary statistics for the period January 1990 to August 2020 of the basic investment plan with no overfunding (x = 1), and for the three overfunded plans, with x1 = 1.02, x2 = 1.04, and x3 = 1.06