Table of Contents

Table of Contents

Links to posts in chronological order -- use browser back arrow to return to this post

Plans for this blog

Why Scratch?

Longevity Graphs

Retirement Income Scenarios

Video on Longevity

Savings and Income Scenario Settings

The X% Rule

Investment Returns and Inflation

RMD Accounts

Analyzing Multiple Scenarios

Year/Year Incomes

Present Values

Pricing Kernels 

The RISMAT E-book and Software

Using the Scratch Software

     Go to scratch.mit.edu/projects/20867413


    Notes:
         To put Scratch in turbo mode, hold the shift key then click the green flag
         To see the program full-screen, click the icon in the upper left corner
         To select an item from the main menu, click the corresponding screen button
         To turn context-sensitive help on or off, press the keyboard up arrow key
         To return to the main menu, press the keyboard left arrow key
         The release date of the software will be included in the title (RIS-yyyymmdd)

The RISMAT E-book and Software

The astute reader of this blog will note that this post was written more than three years after its predecessor. This was not due to lethargy on the part of the author. Rather, for the last three years I have devoted my research and writing to a project titled:


            Retirement Income Scenario Matrices (RISMAT)




The project has produced both an e-book and a suite of software written in the MATLAB language. Both are available at:


      www.stanford.edu/~wfsharpe/RISMAT


A compressed (.zip) file with the matlab software is available at:

https://drive.google.com/file/d/1GTrNOkDrKnbE4kdmVJ95GKrhjUKD44FM/view?usp=sharing 

A compressed (.zip) file with the ebook chapters is also available at:

https://drive.google.com/file/d/14sTIRbzfCF0ChyXo8lQJeJTs8DX361cB/view?usp=sharing

The material may be used for any purpose without the payment of any fees, subject to the terms of the Create Commons Attribution 4.0 License, for which a link is provided on the RISMAT table of contents page.



The overall approach is similar to that taken in the previous chapters of this blog, which can be considered a precursor for RISMAT. However, the MATLAB software can handle many more scenarios at speeds that are orders of magnitude faster than those of the Scratch software used for this blog. Moreover, the RISMAT e-book covers many more topics, each in considerable depth.



My hope is that this new material and the accompanying programs will be utilized by financial engineers and financial advisors who offer their services to retirees. I believe the RISMAT e-book and software can provide a useful base for comprehensive analyses of alternative approaches for providing retirement income.



It seems fitting to conclude this blog with the last sentence from the RISMAT e-book:



Following tradition, we conclude with the admonition that more research is needed (and probably more programming) -- tasks that the author, having reviewed actuarial tables, chooses to leave to others.

Pricing Kernels

In June of 2011, I was invited to give the keynote speech at the annual meeting of the European Financial Management Association in Braga, Portugal. The talk was to be followed by an article for the EFMA journal. I of course chose to speak about my favorite topic. Here is the reference to the published version:

"Post-retirement Financial Strategies: Forecasts and Valuation", European Financial Management, Vol. 18, No. 3, 2012, pp. 324-351.

A pre-publication version is available here.

In the paper I tried to succinctly and somewhat formally summarize the approach to asset pricing that I and my coauthors Jason Scott and John Watson took in our earlier publications, that I have taken in my subsequent research and that I have used for the present value calculations in the RIS software. Here I'll use portions of the EFMA paper with added comments to provide a semi-formal description. 

Warning: some of this will be technical. Feel free to skim or ignore it, as needed.

Here goes. 

------------
Building on the approach of [Arrow 1952] and [Debreu 1959], I consider first a one- period setting in which there are n possible states of the world. Assume that securities are priced based on their payoffs in the states of the world, using a set of state prices, where p_s is the price today of $1 if and only if state s occurs at the end of the period. If the probability of state s is π_s, define the pricing kernel value for state s at time t as m_st≡ p_st/π_st,  a value that I called in [Sharpe, 2007] the price per chance (PPC). The set of n values of m_st is the pricing kernel for time t. To say that only the market portfolio is priced in a single period is to assert that all states with the same market return have the same price per chance. Moreover, societal risk-aversion implies that the higher a state’s market return, the lower should be its price per chance. More fundamentally, if markets are to clear, prices must adjust so that income in a state of scarcity (low market return) costs more than income in a state of plenty (high market return).

 -------------

The references are as follows:

--------------
Arrow, Kenneth J., 1952, “Le Role de valeurs boursiers pour la repartition le meillure des risques,” Econometrie, Colloques Internationaux du Centre National de la Recherche Scientifique 11, pp. 41-47.

Debreu, Gerard, 1959, The Theory of Value, Wiley and Sons, New York.

Sharpe, William F., 2007, Investors and Markets: Portfolio Choices, Asset Prices and Investment Advice, Princeton University Press, 2007. 
---------------

There are two major approaches to asset pricing in Financial Economics. One, starting with the Capital Asset Pricing Model (CAPM), builds on the mean/variance portfolio theory of Harry Markowitz. The other, using Pricing Kernels, builds on the state/preference theory developed by Kenneth Arrow and Gerard Debreu. It may surprise some that I favor the latter for applications such the analysis of retirement income strategies; after all, I was awarded a Nobel Prize in Economics for my work on the CAPM. But over the years I have found the state/preference approach better suited to the task of modeling the determination of asset prices in a one-period setting and even moreso in settings involving many periods, as with retirement income analysis.

I labored long and hard to justify this preference in a series of lectures given at Princeton and the subsequent 2007 book referenced above. However, the focus of the book was on a traditional one-period setting in which investment is made at the present time based on estimates of the probabilities of various returns over a single period such as a year. As I showed there, those who have been brought up on the CAPM need not fear obsolescence, since its major qualitative conclusions hold in the state/preference setting. The market portfolio is still an efficient investment for the average investor. Moreover, security and portfolio expected returns are linearly related to a measure of return sensitivity to the returns on the market. There is of course a difference. In the CAPM the appropriate measure (beta) is based on the covariance of a security or portfolio's return with that of the market portfolio. In the Pricing Kernel version, the measure is based on the covariance of a security or portfolio's return with a function of the return on the market portfolio, a result that I called its kernel beta. But the main qualitative message holds in both approaches: securities and portfolios are priced so that higher expected returns are associated with greater probabilities of doing badly in bad times. 

The CAPM can be considered a special case of the more general pricing kernel approach in which the kernel for a single investment period is a linear function of the return on the market portfolio over that period. In the book I argued that it is much  more plausible to assume that the relationship is non-linear, with a positive price for every possible future state of the world.

The case for such a relationship is far stronger when one attempts to model the determination of asset prices in a multi-period world -- a necessity for analysis of retirement income strategies. This requires some assumption about the efficiency of alternative strategies for investors with different horizons. Moreover, the assumption must be consistent with market clearing -- that is, the collective demands for available securities must equal the supplies. To return to the EFMA paper:

------------

Consistent with the focus on strategies that utilize only the market portfolio and a riskless asset, I assume that only market returns are priced, both in any single period and also for any multi-period horizon. This could be consistent with a model of multi-period equilibrium in the capital markets, although I have no aspirations to develop one here (or elsewhere). In any event, as I will show, the assumption greatly restricts the characteristics of the pricing function. 
------------

This may be far too simple a characterization of equilibrium in capital markets. But to value different combinations of income over future years one must make some assumptions and to be credible they should be consistent with market clearing. The approach taken in the EFMA paper seems reasonable. And some such model is better than none. Interestingly, in a world of this sort, it is quite simple to determine the manner in which assets should be priced. Here is the formal analysis from the paper:

------------

Now, consider a two-period case. Letting m_t and r_t represent vectors of pricing kernel values and market total returns (value-relatives) respectively for time t, the pricing kernels for periods 1 and 2 can be written as: 

        m₁ = f₁(r₁ )
        m₂ = f₂(r₂)

The pricing kernel for a horizon that includes both periods 1 and 2 will be the (dot) product of the two kernels. Thus: 

       m₁⋅ m₂ =  f₁(r₁ )⋅ f₂(r₂)
    

In order for (1) the market to be priced in the same manner for each period and (2) for only the market to be priced for any multi-year horizon, it must be the case that: 

       f(r₁ ) ⋅ f(r₂) =  g(r₁⋅ r₂)

A necessary and sufficient condition for this to be the case is that the one-period pricing function be isoelastic: 

      m_t  = Ar_t^-b 

More generally, if M_t represents the pricing kernel for payments t periods hence and V_t the cumulative market return over that horizon: 

     M_t  = A^t V_t ^-b
------------

While this may look formidable, it is quite sensible. Back to the paper:

------------
Taking the logarithms of both sides of the equation:

    log(M_t ) = log(A^t ) − b⋅log(V_t )

Clearly, the b coefficient indicates the elasticity of the pricing kernel with respect to cumulative market return – for every one percent increase in the latter, the pricing kernel decreases by approximately b percent. As is well known, this can be interpreted as indicating that a “representative investor” who holds the market portfolio has a utility function with a constant relative risk-aversion coefficient of b (for a further discussion, see Sharpe, 2007).
------------

The relationship is log-linear. Among other things, this means that every state price is positive, no matter how large or small the return on the market portfolio. This is not the case if the kernel is linear, which can be if one assumes that all investors have quadratic utility functions (and hence care only about the mean and variance of portfolio returns). 

This is not the place to pick a fight between proponents of traditional mean/variance portfolio theory with a strict interpretation of the of the CAPM and those (including me) who prefer the state/preference pricing kernel approach. The key point is that we now have a pricing kernel with desirable properties that could be consistent with a multi-period equilibrium. For any period in the future, we can use it to determine the price per chance (PPC) for $1 to be received at that future time if and only if the state occurs.  PPC is a function of the cumulative compounded return on the market portfolio from the present to that future time and the function is of the form shown above. And clearly, the price for $1 to be received at a future time and in a given state is equal to the PPC times the chance that the state will occur at the time in question.

To return to the paper:
------------
I will use the pricing kernel primarily to estimate the values of the distributions of spending (also called “payments” or “paychecks”) provided by a strategy. This requires a set of state prices – each of which represents the cost today of obtaining $1 at a given future time and state. Since the pricing kernel is simply a set of ratios of state prices to state probabilities: 

     P_t = M_t⋅∏_t

Where P_t  is a vector of state prices for payments at a future time t, and M_t and ∏_t are, respectively, vectors for the pricing kernel and probabilities of the states for that time. In the simulations ... market returns were drawn randomly from the underlying probability distributions for n multi-year scenarios. Considering each scenario as a state, the probabilities all equal 1/n so that:

P_t  = M_t /n
----------

But how to determine the key parameters in the equation (A and b)? In the simulations for the EFMA paper, I ran a million simulations for each year for every strategy analyzed. Moreover, my Matlab program made it possible to run all the results for each year at one time, then find the pricing kernel that best fit the resulting million simulated compound market returns. For each future year, I found values for the two variables (A and b) that would be consistent with valuations of the market portfolio and the riskless asset for that horizon. As I indicated in a footnote in the paper:

------------
In the simulations ... to reduce sampling errors of draws from the underlying return distribution, the values of A and b used for each horizon were found iteratively based on the requirement that the implied present values of both the cumulative market returns and the cumulative risk-free return were both within a very small distance from 1.0.
------------

John Watson is currently doing research on the efficacy of this and other possible approaches to valuation when using a simulated sample to approximate the universe of all possible scenarios. This procedure produces pricing kernels that may fit the simulated scenarios but may violate the underlying assumption about the multi-period equilibrium. I shall have more to say about this in subsequent posts. But the discussion is not particularly germane for those using the RIS software for two reasons. First, it is infeasible to run a million scenarios. And second, since the software generates results a scenario at a time, this fitting procedure is not feasible since it would require storing at least 40 or 50 million results (a matrix with a row for each scenario and a column for each year). 

Fortunately, there is another way, using a formula derived by John Watson in our earlier paper:

Scott, Jason S. with William F. Sharpe and John G. Watson, 2009, "The 4% Rule -- At What Price?", Journal of Investment Management, Vol. 7, No. 3, Third Quarter 2009, pp. 31-48 

 
It is available here.

As derived in the paper, for the special case in which the return on the market portfolio is assumed to be lognormally distributed in each period (year), the coefficients of the asset pricing equation can be calculated directly. Here is the EFMA paper's description:

-----------
The values of the coefficients in the pricing equation can be computed directly from the parameters of the assumed distributions of annual returns as follows:

In the above equations, R_f is the total risk-free return, E_m = E[R_t] is the yearly expected total market return, and S_m = (Var[R_t])^1/2 is the annual market volatility.

------------

The RIS software uses these formulas to compute the parameters for the pricing kernel (computing b first, then A). While producing multiple scenarios, the formula is used to compute the state price for each income and fee payment, obtain the resulting present value and update the appropriate cumulative sum.

Unfortunately, with feasible numbers of scenarios the results will be imperfect since the sample scenarios will not be fully representative of the universe of possible future scenarios. As indicated in the previous post, this is reflected in the design of the RIS  software. First, present values are shown only when 5,000 or more scenarios have been analyzed. Second, the computed values are not shown, only the percentages of the total value associated with the prospects for three recipients (the household, the estate and those who receive fees). Third, the actual percentages are not shown, only a pie chart indicating their magnitudes. And finally, the user is encouraged to generate additional sets of multiple scenarios to determine the variation in the allocation of the present values from case to case. 

Enough caveats. Despite these limitations, the breakdown of present values in the RIS software should provide valuable information in most cases. And some of the analyses that I will describe in subsequent posts will be obtained using Matlab software,  employing large numbers of scenarios and, in some cases using different and hopefully more accurate approximations of present values. 

I will have much to say about the present values of future possible outcomes in the future analyses of alternative retirement income strategies. You should expect nothing less from an economist.

Present Values

A retirement income strategy is designed to provide incomes at different times to different recipients. However, the outcomes will depend on future events -- in our scenarios, longevity, investment returns and inflation. But it is helpful to consider the present value of prospective incomes and how it is distributed among the key recipients.

In this post I'll describe the basic approach used for the Present Value Analysis included in the RIS software; more details will follow in the next post.

To start, consider the problem of valuation of an investment in a share of Apple common stock. This is a claim that provides the owner with any future dividends that Apple might pay as well as the ability to sell the share at any future date. Since shares are traded publicly it is relatively straightforward to determine the market price, which represents a consensus opinion of the present value of Apple's future prospects. If you think the share is worth more, you may well want to buy more shares; if less, you may want to sell any shares you hold. But the market price reflects the overall evaluation of other investors. Economists define this as its present value and consider it a highly relevant measure.

Note that the price of an Apple share is based on its uncertain future prospects. In some possible scenario, future earnings and dividends will be disappointing and the prior price will turn out to have been too high. In another scenario, the future earnings and dividends will be spectacular and the prior price will turn out to have been too low. But the price today reflects investors' assessment of many different possible scenarios.

The RIS software takes the same approach when calculating the present values of prospective incomes. In particular, it generates a number of scenarios, each covering many years. 

Now, imagine a large matrix (table) with 5,000 rows (one for each scenario) and 40 columns (one for each future year). In each cell of the matrix there is an income payment and an indication of the recipient (you and your partner, you alone, your partner alone, or your estate).

Next, imagine that for every cell in this matrix there is a number representing the present value of $1. For example, in row (scenario) 123, column (year) 10 there is a number representing the amount that investors would pay today to receive $1 in year 10 if and only if scenario 123 (and no other) happened. What would determine this value? First, the probability that the scenario would take place. And second, the state of the world in that scenario.

The probability is easily determined. In our simulation there are 5,000 scenarios so the probability that any one will occur is 1/5000.  But there is good reason to believe that the prices will differ across scenarios.

To approach this in a more general context, economists use a concept that I like to term the price per chance (PPC) for payment at a given future time in a given future state of the world. And there is good reason to assume that this value will depend on (1) the number of years before the state occurs and (2) the economic health of the economy in that state. In particular, it makes sense to assume that other things equal, the present (market) value of a future unit of purchasing power is greater for states of scarcity than for states of plenty. In our simulations the best available measure of the state of the economy is the future level of the real value of the market portfolio of world bonds and stocks. So we assume that the present value of $1 in a future state divided by the probability of that outcome depends on the number of years before the state occurs and on the cumulative real value of $1 invested in the market portfolio today and held until that future year. In sum, other things equal, the present value of $1 to be received in a given state divided by the probability is lower, the farther in the future the state and the smaller the cumulative return on the market over the period.

The pricing kernel formula reflects these assumptions. It indicates the price (present value) per chance (probability) of $1 for a given state in the future. To return to our matrix with 5,000 rows (scenarios) and 40 columns (years), imagine that there is present value of $1 in each cell, computed by multiplying the corresponding price per chance by the probability (in this case, 1/5,000). Then we can value each cell's payment by multiplying the amount paid by the present value. Moreover we can sum up the present values of all the payments made to various recipients. In effect, this is what the RIS software does (although it actually updates sums as it creates scenarios in order to avoid the need to store large matrices). In addition, the software provides similar computations for any fees paid to financial advisors or others (as if they were stored in a second matrix).

Here is a result (based on a standard RMD account with 1% fees), obtained by clicking the Present Values button under Analyses.

In this case the present value of the household's prospective incomes is roughly 70% of the total value of the future prospective payments. The present value of the estate's prospects are somewhat more than half of the remainder, with the present value of the financial advisor's prospects taking the rest.

I willl have more to say about these pie charts in future posts analyzing alternative strategies. But here are two generic observations.

First, only annuities offer reasonable income prospects for a household without providing valuable prospects for an estate. Moreover, different variations on a theme can substantially change the division of present value between the household and its estate, as can be seen by changing settings, then generating another set of scenarios.

Second, fees that may seem small may offer very valuable prospects for the person or organization charging the fees, with an associated decrease in the value of the prospects for the household and/or its estate. Note that in the figure above, a 1% fee is worth almost 1/8th of the total present value -- a cautionary tale.

Finally, some mechanics. 

Present value calculations using a limited number of scenarios are subject to some error. I'll explore this more in subsequent posts, but suffice it to say that it is important to have at least 5,000 scenarios before making such a calculation. You will see that if the scenario settings called for fewer than 5,000 scenarios when you produced multiple scenarios, the Present Values analysis will politely refuse to give you any results. The solution is, of course, to change the scenario settings, then produce a new set of multiple scenarios.

You will notice also that there are no numbers given for the percentages of the various recipients' prospects. This was not an oversight. I chose to show only the pie chart in deference to the imprecision of the estimates. You can test the dependence of the results on the scenarios generated by pressing the multiple scenarios button again, generating a new set of scenarios, then seeing the resulting present value pie chart. With luck, the variation may be relatively slight. And the larger the number of scenarios, the smaller it should be. 

Here's the bottom line. Any Monte Carlo analyses is at best an approximation of reality. It is a complex world out there. At best we can only hope that our models and simulations shed useful light on the relative prospects of different strategies.

Year/Year Incomes

The previous post ("Analyzing Multiple Scenarios") showed how to generate multiple scenarios using the RIS software, then view the yearly income ranges using the Analysis routine titled (appropriately enough) "Yearly Incomes". The resulting graph shows the probabilities of exceeding alternative levels of income in each of the future years in the selected range.

In economic analyses of multiple years' incomes, it is often assumed that people have "time separable utility functions." In principle, a person with such preferences can evaluate the probabilities of alternative incomes in any given year, then calculate the "expected utility" of that range of possibilities. If in this manner the results for each of the future years are expressed in terms of present-day happiness, the expected utilities for each of the years can be added together to get a single measure of the overall desirability of the prospective future incomes. An individual preferences of this type can in theory determine the desirability of a strategy by studying only a graph showing yearly incomes with the accompanying information about the probability of being alive in each year.  To be sure, even in this simplified setting, a couple would have to consider the probabilities of being alive and most people would want to include information about the possible amounts that could be left for an estate. But information on the changes in income from each year to the next would not be needed.

However, many people are in fact concerned with the extent to which their income might change from year to year. To accommodate them the RIS software also includes graphs showing at least some relevant information. On the Analysis page click the button labelled "Yr/Yr Incomes"; you will get a graph something like this:

In this case, each curve shows (on the vertical axis) the probability of exceeding a given value of a ratio of income divided by the prior year's income (X), shown on the horizontal axis. The ratios shown run from 0 to 2.0, with a ratio of 1.0 indicated by a green vertical line. As with the Yearly Income graph, the chance that one or both of you and your partner will be alive in a given year is shown in the upper right. Also, as with that graph, there are two variants. In the first, as in the case shown above, each graph starts at 100% and shows the ratios for cases in which at least one of the members of the household is alive. The other, shown below, plots the probabilities that the ratios will (1) exceed various levels and (2) that one or both of the members of the household will be alive.

The first item in the Analysis Settings determines which graph will be shown -- C(ontingent), as in our first example, or A(ctual), as in the second. As with the yearly income graph, the number of seconds of delay between years will be that indicated in the Analysis Settings.

It is important to understand how these graphs are constructed. With some exceptions, as each one of the multiple scenarios is generated, the ratio of each year's income to that of the prior year is calculated, with the result added to tables to be used to produce the graphs. Excluded are cases in which the recipient in a year differs from that in the prior year. Thus the final year in which any remaining savings goes to the estate is excluded as well as any year in which the household changes from two people to one (since some strategies call for changes in income in such instances that are not related to investment performance).

As with the yearly income graph, the year/year income graph should be interpreted as showing information available at the present time about the range of outcomes that could take place future years. For example, the range of possible year/year income ratios for years 19 and 20 is based on information available today. When year 19 actually arrives, the range of possible ratios of year 20 income to year 19 income will undoubtedly be very different, since many possible scenarios concerning the first 19 years will have failed to take place.

This distinction helps illustrate the difficulty of assessing the range of possible outcomes for any multi-year income strategy. There are simply too many income combinations to consider in detail. To emphasize the point, look at the graph below, showing the likelihoods of different combinations of income one and two years hence for a market-based strategy.

 It is possible that someone with training and patience could study such a graph, compare it with another showing the likelihoods of incomes from some other strategy, then choose the preferred strategy. But most retirement income strategies have many more dimensions (years of income) than two. It is simply impossible to portray and assess alternative 40-dimensional probability distributions of income.  One needs rather to concentrate on a manageable number of attributes. For some investors our yearly income graphs will suffice. Others may wish to also consider the year/year income graphs. More than that is beyond this project.

Analyzing Multiple Scenarios

Multiple Scenarios 

While it is very useful to see possible scenarios for future income and savings one at a time, there is merit in getting a view of the range of possible outcomes over many possible scenarios. Starting with RIS-20120120, my software on the Scratch site allows for the analyses of multiple scenarios. This post will describe the required procedures and show some examples, but will be short on analysis. Future posts will discuss the relevant economics and analyze alternative retirement income strategies.

In a previous post I discussed the scenario settings. There is now one additional setting that indicates the number of scenarios that you wish to create for multiple scenario analyses. It is the last one shown in the figure below. To obtain meaningful results you will need to analyze a great many scenarios so the setting is stated in thousands. The default is 5 thousand. I allow as few as one thousand but strongly recommend at least five thousand and more if you are willing to take the time.

 

 The Main Page

The main page now contains the ten buttons shown below. The new ones provide for  Analysis Settings, Multiple Scenarios and Analyses. 

 

The Multiple Scenarios Button

To create multiple scenarios, you need only click on the Multiple Scenarios button. After a short while you will be given an estimate of the time required to complete the process and asked whether you wish to proceed. If you say no, you will return to the main page but there will be no scenario statistics available to be used for any subsequent analyses. If you say yes, the desired number of scenarios will be generated and statistics gathered. You will see the progress on the screen  and it is very important that you do not interrupt the process. After it is completed, you may click the Analyses button at any time to see various properties of the scenarios.

Analysis Settings

The Analysis Settings button allows you to change the default settings for various analyses. At present there are only two settings, although more will be added. They are shown below.

The first setting indicates whether you want to see the actual probabilities of receiving income (A) or the contingent probabilities (C). As shown, the default is actual. I'll describe the two alternatives below.

The second setting indicates the time that the software should wait between years when plotting multiple yearly outcomes. The default is 0.25 seconds, which makes the plots come rather fast. You may want to use a larger value to slow down the display, although you can always stop temporarily by pressing and holding the 's' key (then resuming by pressing the space bar).

The Analysis Page

When you press the Analysis button on the main page, you will be transferred to the Analysis Page. The figure below shows its current state.

At present, only the first button is operative. The other two will be programmed in the future, and more may be added as well. Note that all analyses will use the multiple scenarios that you generated most recently. As usual, you may return to the main page by pressing the left arrow key.

Yearly Income Analyses

Now to the good part --what happens when you press the Yearly Income button.

I strongly suggest that you start by doing this using the software initial defaults settings (which include 5,000 multiple scenarios previously produced) to see the results in their full animated glory.

I'll start with graphs produced using all the default settings (based on an RMD account). In this case the  Analysis Setting calls for Actual probabilities. The figure below shows the graph after 18 years.

To produce this figure, I froze the display by pressing and holding the 's' key after the 18th year was shown. To produce the next figure, I simply pressed the space bar. (As usual, you can find context-sensitive help instructions by pressing the up arrow key to get a help message).

Let's look at the results. As shown at the top of the graph, the red curve is for year 18 (18 years in the future since the current year is year 0). The chance that you, your partner or both will be alive in that year is 89.7%. The horizontal axis shows levels of income from 0 to 80 $ thousand (the upper limit, taken from your scenario settings). There are twenty vertical grid lines, so in this case, each covers 4 ($thousand). Here the values shown on the horizontal axis are for real income, also taken from your scenario settings. You may change any of the Scenario Settings to produce different graphs, then producing a new set of multiple scenarios by pressing the Multiple Scenario button.

 

 As indicated, the vertical axis shows the chance that (1) income will exceed the value shown on the horizontal axis and (2) that one or both of you will be alive. Values range from 0% to 100% (or, for those of you who think in probability terms, from 0 to 1.0). Each horizontal grid line covers 5%, and the 50% (median) line is indicated as well. 

You can read this graph in either of two ways. 

You could pick a real income goal, say $40,000, find it on the horizontal axis, then go up to the curve and look over to the vertical axis to see your chances of doing that well or better -- in other words, beating that goal. Clearly, the better your chances, the happier you will be. Thus higher curves are better than lower ones.

Or you could pick a chance, say 50%, find it on the vertical axis, then go to the curve and look down to the horizontal axis to see the goal that you have a 50% chance of beating. The higher that goal, the happier you will be. Thus curves farther to the right are better than ones to the left.

If you have a statistical background, you may recognize this graph as similar to a cumulative probability distribution, but with one key difference. The typical statistical graph shows the probability of falling below the value on the horizontal axis, not the probability of exceeding it. I think this is not the way most human beings think about attaining goals and strongly prefer the approach I've employed in my prior research and incorporated in the RIS software. I'll probably have more to say about this "goal/chance" approach in future blogs.

Most of those who analyze retirement income strategies pick one, two or three probabilities (chances), then show the incomes associated with each of them in each future year. I feel that it is far better to show the entire ranges, as does the RIS software. I'll undoubtedly have more to say about this as well in the future.

To return to the figure, note that when income exceeds the maximum shown on the horizontal axis, the plot is just to the right of the vertical right edge of the graph box. Here, the actual income values are greater than 80 $thousand maximum plotted, but there is no way to tell how much greater they may be. If this is of concern you may want to change the Scenario Settings to provide higher maximum incomes, then run a new set of multiple scenarios, and analyze the results using the appropriate Analysis tools.

Now, back to the case at hand. The figure below shows the graph after all the years specified in the Scenario Settings have been shown. Not surprisingly, as time goes on, the chance of any income diminishes as mortality takes its toll. Moreover, there is a wide range of possible incomes in all but the initial year, and the range tends to be larger for later years. In future posts I'll discuss such matters at length when analyzing specific retirement income strategies.

I'll finish this post with graphs produced using the Contingent Probability setting in the Analysis Settings. (Happily, you do not have to run a new set of Multiple Scenarios to  change between Actual and Contingent probabilities).

The figure below shows years 0 through 14 in yellow and year 15 in red. For the first few years the graphs are virtually the same is in the previous case, since the probability that one or both of you will be alive in the near future close to or equal to 100%. The only difference is the heading for the vertical axis, which shows that the results indicate the chance that income will exceed the amount on the horizontal axis if one or both is alive. (In that sense, it is contingent). Note that this shows that for at least the next 15 years the median (50%) real income is larger in future years, the low-probability worst cases (90% and above) are somewhat worse, and the rosier low-probability cases (say, 25% and below) are considerably better.

The figure below completes the picture, including all the future years through year 49. As can be seen, the prospects for the very distant years become quite dismal. But of course the chances that anyone will be alive at the time are small.

Note also that in this case the curves for distant years are far from smooth. This reflects the fact that while the results were based on 5,000 scenarios, there are very few scenarios in later years in which anyone is alive, so the sample sizes are insufficient to provide good indications of the overall range of possible future outcomes. For example, in year 49 (shown in red), there were only 5 scenarios (0.1% of 5,000) -- far from sufficient to make well-informed estimates of the whole range of possibilities. Unfortunately, the only way to improve the reliability of distant forecasts is to take the (considerable) time required to run many more scenarios. But with at least 5,000 you should be able to get a rough idea of possible prospects.

There are profound differences between viewing future retirement income prospects using actual probabilities and conditional probabilities, as these figures show. Indeed, there is considerable debate about the extent to which people should weigh each of these two views when choosing among alternative retirement income strategies. I'll have more to say about this anon. Meanwhile, please do use the software to experiment with these new features.

RMD Accounts

As indicated in previous posts, an analysis using the RIS software can employ one or more accounts, each of which provides retirement income. Earlier I described the X% Rule account. This post covers another type, based on the Required Minimum Distribution requirements specified by the U.S. Internal Revenue Service for those older than 70 ½ holding investments in tax-deferred accounts such as Individual Retirement Accounts and 401(k)s.

The IRS rules are provided in IRS Publication 590. Required distributions each year are determined by dividing the value of an account by a life expectancy. Equivalently, the required distribution is equal to a percentage of the value of the account, with the percentage equal to the reciprocal of the life expectancy. For example, if the life expectancy is 20 years, the required withdrawal percentage is 1/20, or 5%. The required distribution amount each year must be moved from tax-deferred accounts and declared as income subject to regular income tax rates; otherwise a prohibitive tax is levied.


Life expectancies are given in three tables, each of which is applicable for taxpayers in a particular category. The simplest and most widely applicable is Table III, which is required for use by: “Unmarried Owners, Married Owners whose Spouses are Not More than 10 Years Younger, and Married Owners Whose Spouses are Not the Sole Beneficiaries of their IRAs” (IRS Publication 590, p. 109).


The first two columns of the table below are taken directly from publication 590. “Dist Period” is the Distribution Period (Life Expectancy). I have added the final column, which shows the percentage of an account that must be distributed at each age.

Age	Dist Period	Percent
70	27.4	3.65%
71	26.5	3.77%
72	25.6	3.91%
73	24.7	4.05%
74	23.8	4.20%
75	22.9	4.37%
76	22.0	4.55%
77	21.2	4.72%
78	20.3	4.93%
79	19.5	5.13%
80	18.7	5.35%
81	17.9	5.59%
82	17.1	5.85%
83	16.3	6.13%
84	15.5	6.45%
85	14.8	6.76%
86	14.1	7.09%
87	13.4	7.46%
88	12.7	7.87%
89	12.0	8.33%
90	11.4	8.77%
91	10.8	9.26%
92	10.2	9.80%
93	9.6	10.42%
94	9.1	10.99%
95	8.6	11.63%
96	8.1	12.35%
97	7.6	13.16%
98	7.1	14.08%
99	6.7	14.93%
100	6.3	15.87%
101	5.9	16.95%
102	5.5	18.18%
103	5.2	19.23%
104	4.9	20.41%
105	4.5	22.22%
106	4.2	23.81%
107	3.9	25.64%
108	3.7	27.03%
109	3.4	29.41%
110	3.1	32.26%
111	2.9	34.48%
112	2.6	38.46%
113	2.4	41.67%
114	2.1	47.62%
115 and over	1.9	52.63%

The calculations made by the IRS to generate this table are not specified. Presumably, mortality tables were utilized, with some sort of averaging across possible combinations of unmarried investors of both sexes and those married with spouses of both sexes and with differing ages.


There is no presumption that the owner of a tax-deferred account must spend the amount on which taxes must be paid. And many investors have additional sources of retirement income. This said, it has occurred to some investors and analysts that it might be desirable to adopt a retirement income strategy with a policy of spending the percentages of overall savings given in the final column above. Prominent studies of the efficacy of such an approach include:

Sun, Wei and Anthony Webb, 2012, “Should Households base Asset Decumulation Strategies onRequired Minimum Distribution Tables?” Center for Retirement Research at Boston College Working Paper (available here).

          
Blanchett, David, Maciej Kowara and Peng Chen, 2012, “Optimal Withdrawal Strategy for Retirement-Income Portfolios,” Retirement Management Journal, 2(3): 7-20

Blanchett, David M. 2013. “Simple Formulas to Implement Complex Withdrawal Strategies.” Journal of Financial Planning 26 (9): 40–48, available here .

In their 2012 paper, Sun and Webb concluded that the RMD approach was preferable to the 4% rule. In his 2013 paper, Blanchett found that “the RMD approach works well for periods less than 15 years...” and advocated the use of a more complex approach for subsequent years. I remain agnostic on the issue but feel that the approach is worthy of investigation.

Now, to the details of the RMD account.


To cover ages not shown in the IRS table, I have made the assumption that the life expectancy for any age younger than 70 will be (70 – age) years longer than that for age 70. Thus for a 65-year old the assumed life expectancy is 27.4 + 5, or 32.5 years. Moreover, when there are two participants (you and your partner), I base the withdrawal percentage each year on the age of the older participant in that year.


The settings for an RMD account are shown below.

The initial setting (line 2) is the usual one that determines whether or not the account is active. The second (line 4) indicates the initial balance – here, a million dollars (1,000 $ thousand). The next setting allows for a variation of the strategy in which the RMD longevity numbers are altered by adding or subtracting a constant number of years. For example, if the adjustment is 2, the life expectancy at age 70 will be 29.4 (27.4 + 2) years, and every other life expectancy will be adjusted by adding 2 years to the amount shown in the table. Lengthening the life expectancies in this manner will reduce the percentages of savings paid, lowering retirement income payments and increasing possible estate values. If desired, you may enter a negative number for this setting. This will reduce the life expectancies and increase the percentage payments. (Not to worry -- if this would result in any expectancies less than one, they are replaced with 1.0).


The remaining settings are the same as those for the X% Rule settings. The fee indicates the annual percentage of the account value charged as fees to third parties. The three settings for the investment strategy are, as for the X% Rule, the initial proportion of the account invested in the market portfolio, the number of years for any glide path, and the proportion of the account invested in the market portfolio at and after the end of the glide path period. As with the X% Rule, the default settings provide for a constant investment solely in the market portfolio in each year.

The RMD approach is a special case of a more general class that I have called Proportional Payout (PPO) strategies, in which a pre-specified proportion of an investment account is paid out to provide retirement income in each year. In previous research, I have used a set of proportions specified for the Fidelity Income Replacement 2042 Fund, which is designed to pay out all the assets in the portfolio by the end of 2042. For a detailed analysis, see my paper in the European Financial Management Journal, a version of which is here. While the Fidelity Funds are designed specifically for producing retirement income, they will pay out all assets by a target date no more than 30 years in the future. In contrast, the RMD approach as implemented in the RIS software will provide some income until the both participants are gone, leaving at least some funds for an estate. For this reason, and because it uses non-proprietary data, I chose to include the RMD method in the software. However, it would be a simple matter for a user to alter the longevity table used for the calculations to produce different results.

Do try the RMD account. Unlike the X% Rule, it conforms with two sensible criteria in each year:

      The amount you spend should depend on
              1. How much money you have, and
              2. How long you are likely to need it

This doesn't mean it is the best approach for you. But, combined with a sensible investment policy, it might provide a desirable component for your overall retirement income strategy.