Pedersen2013monte Carlo

Monte Carlo Simulation in Financial Valuation
By Magnus Erik Hvass Pedersen1 Hvass Laboratories Report HL-1302 First edition May 24, 2013 This revision June 4, 20132
Please ensure you have downloaded the latest revision of this paper from the internet: www.hvass-labs.org/people/magnus/publications/pedersen2013monte-carlo.pdf Source code and data files are also available from the internet: www.hvass-labs.org/people/magnus/publications/pedersen2013monte-carlo.zip Automatically updated figures are also available from the internet: www.value-yield.com
Abstract
This paper uses Monte Carlo simulation of a simple equity growth model with historical data to estimate the probability distribution of the future equity, earnings and payouts of companies, which are then used to estimate the probability distribution of the future return on the stock and stock options. The model is used on the S&P 500 stock market index and the Coca-Cola company. The relation between USA government bonds, the S&P 500 index and the Dow Jones Venture Capital index is also studied and it is found that there is no consistent and predictable risk premium between these three asset classes.
If you find any parts difficult to understand, or if you encounter any errors, whether logical, mathematical, spelling, grammatical or otherwise, please mark the document and e-mail it to: Magnus (at) Hvass-Labs (dot) Org 2 See last page for revision history.
Nomenclature
IID PDF CDF Independent and identically distributed stochastic variables. Probability Density Function. Cumulative Distribution Function (Empirical). Present value of future payouts and share-price. Annual growth rate used in valuation. Discount rate used in valuation. Kilo, a factor Million, a factor Billion, a factor Infinity. Price-to-Book ratio: Number of shares outstanding. Price per share. Market capitalization (also written market-cap): Capital supplied by shareholders as well as retained earnings, not per-share. Earnings available for payout to shareholders, not per-share. Dividend payout, pre-tax, not per-share. Amount used for share buyback. Amount from share issuance. Net payout from company: Ratio of earnings retained in company: Return on Assets: Return on Equity: means that equals means that is greater than or equal to means that is approximately equal to Implication: means that implies Bi-implication: means that and Multiplication: means that is multiplied by Summation: Multiplication: Maximum of or . Similarly for Logarithmic function with base (natural logarithm). Absolute value of (the sign of is removed). Probability of the stochastic variable being equal to . Expected (or mean) value of the stochastic variable . Variance of the stochastic variable . Standard deviation of the stochastic variable :
k m b
1.
Introduction
Present value calculations are ubiquitous in finance. But the average of stochastic inputs cannot be used to calculate the average present value because the formula is non-linear. This is due to the so-called Jensen inequality. It is therefore necessary to know the probability distribution of the inputs in order to calculate the mean present value. The probability distribution is also useful when the average present value is misleading because it is unlikely to occur. Monte Carlo simulation is the use of computers to simulate numerous outcomes of a mathematical model so as to estimate the probability distribution. This is useful when the model cannot be studied analytically. There are several problems in finance where Monte Carlo simulation is useful, see e.g. Glasserman [1]. This paper uses a simple equity growth model to simulate the future equity, earnings and payouts of companies, based on historical data for the return on equity and the fraction of earnings retained in the past. This is a reasonable model for companies whose earnings are related to their equity capital. The model can also estimate future share prices by multiplying the simulated equity with the historical distribution of the companys P/Book ratio (also called price-to-book or price-to-equity). The model can also estimate the annualized rate of return from buying shares at a given price and selling them after a given period or holding them for eternity. The model can also be used for valuing options. The equity growth model may also be used for stock market indices by considering the constituent companies of the index as if they were one big conglomerate. This is done for the Standard & Poors 500 (S&P 500) stock market index in section 6. Of particular interest may be the probability distribution of the future rate of return on the S&P 500 index which is given for an infinite holding period in section 6.7 and for finite holding periods in section 6.9. Options on the S&P 500 index are valued in section 6.10. The historical yield on USA government bonds is studied in section 7. The Dow Jones Venture Capital (DJVC) index is studied in section 8. It is found that there is no consistent and predictable risk premium between USA government bonds and the S&P 500 and DJVC indices. The Coca-Cola company is studied in section 9 including the probability distribution of its present value relative to the S&P 500 and DJVC indices.
1.1.

Paper Overview
Section 2 gives formulas for present value calculations and the equity growth model. Section 3 gives formulas for stochastic present value calculations. Section 4 gives formulas for stock option valuation. Section 5 gives algorithms for Monte Carlo simulation of the equity growth model. Section 6 studies the S&P 500 stock market index. Section 7 studies USA government bonds and their historical relation to the S&P 500 index. Section 8 studies the Dow Jones Venture Capital (DJVC) index. Section 9 studies the Coca-Cola company. Section 10 gives directions for obtaining the source-code and data files. Section 11 is the conclusion.
The paper is structured as follows:
2.
Present Value
The present value of the dividend for future year is the amount that would have to be invested today with an annual rate of return , also called the discount rate, so as to compound into becoming after years:
Eq. 2-1
An eternal shareholder is defined here as one who never sells the shares and thus derives value from the shares only through the receipt of dividends. In Williams theory of investment value [2], the value of a company to its eternal shareholders is defined as the present value of all future dividends. Let denote the present value of all future dividends prior to dividend tax and not per share. Assume the discount rate is constant forever. The present value is then:
Eq. 2-2
2.1.
Payout
Instead of making dividend payouts, companies may also buy back or issue shares, the net effect of which will be called payout and defined as:
Eq. 2-3
The term payout is a misnomer for share buybacks net of issuance as argued in [3] [4], because a share buyback merely reduces the number of shares outstanding which may have unexpected effects on the share-price and hence does not constitute an actual payout from the company to shareholders. A more accurate term for the combination of dividends and share buybacks is therefore needed but the term payout will be used here and the reader should keep this distinction in mind. The ratio of earnings retained in the company is:
Eq. 2-4
This means the payout can then also be written as:
Eq. 2-5
If the company has no excess cash, then the company will first have to generate earnings before making dividend payouts and share buybacks, so the present value is calculated with starting year :
Eq. 2-6
2.2.
Equity Growth Model
The companys equity is the capital supplied directly by shareholders as well as the accumulation of retained earnings. Earnings are retained for the purpose of investing in new assets that can increase future earnings. The Return on Equity (ROE) is defined as a years earnings divided by the equity at the beginning of the year. For year this is:
Eq. 2-7
This is equivalent to:
Eq. 2-8
Similarly the Return on Assets (ROA) is defined as a years earnings divided by the assets at the beginning of the year. For year this is:
Eq. 2-9
Let
be the fraction of earnings retained in year so that:
Eq. 2-10
This means that the equity grows by the rate
calculated as:
Eq. 2-11
The accumulation of equity is:
Eq. 2-12
Let be the amount used for dividend payouts and share buybacks net of issuance in year . Using Eq. 2-5 with Eq. 2-8 gives:
Eq. 2-13
The value to eternal shareholders is then:
Eq. 2-14
2.2.1. Normalized Equity The accumulation of equity in Eq. 2-12 can be normalized by setting so it is independent of the starting equity. This also normalizes the earnings calculated from Eq. 2-8 and the payouts calculated from Eq. 2-13, which allows for Monte Carlo simulation based solely on the probability distributions for ROE and Retain so the results can easily be used with a different starting equity. 2.2.2. Mean Growth Rate Consider and to be stochastic variables and their product growth rate, see Eq. 2-11. The mean equity growth rate is then:
is the equity
Eq. 2-15
The mean earnings also grow by this rate because earnings are derived from the equity. This follows from the definition of earnings from Eq. 2-8, the assumed independence of the stochastic variables and , and the identical distributions for and so their means are also identical:
Eq. 2-16
A similar derivation can be made from Eq. 2-13 to show that the mean payout growth rate is also
Independence If and
are independent stochastic variables then the mean growth rate is:
Eq. 2-17
2.3.
Equity Growth for Two Companies

and , the
Consider two companies whose starting equities are and with whose annual equity growth rates are and with . Because of the higher growth rate lower grows to exceed after years, which is calculated as:
Eq. 2-18
2.4.
Share Issuance & Buyback
Share issuance and buyback may result in considerable changes to the present value per share depending on the price per share at the time of the share issuance or buyback, see [3] [4]. But for simplicity, it will be assumed here that future share issuances and buybacks will always be made at the equilibrium where the share price equals the present value per share and hence does not affect the value to eternal shareholders.
2.5.
Debt Change
Changes in the level of debt relative to equity are taken implicitly into account by the use of historical data for ROE and Retain which are affected by the change in debt levels. This is deemed acceptable because the changes are small and not permanent in these case studies. Large and permanent changes should be taken into account in the valuation models.
2.6.
Market Capitalization
Let be the number of shares outstanding and let be the market-price per share. The market-cap (or market capitalization, or market value) is the total price for all shares outstanding:
Eq. 2-19
The market-cap is frequently considered relative to the equity, which is also known as the price-to-bookvalue or P/Book ratio:
Eq. 2-20
Eq. 2-21
Because the starting equity is normalized to one in these Monte Carlo simulations as described in section 2.2.1, it is often convenient to express the formulas involving the market-cap in terms of the P/Book ratio instead.
2.7.
Value Yield
The value yield is defined as the discount rate which makes the market-cap equal to the present value:3
Eq. 2-22
This may be easier to understand if the notation makes clear that the present value is a function of the discount rate by writing the present value as . The value yield is then the choice of discount rate that causes the present value to equal the market-cap:
Eq. 2-23
The value yield for an eternal shareholder must satisfy the equation:
Eq. 2-24
For a shareholder who owns the shares and receives payouts for years after which the shares are sold at a price of , the value yield must satisfy the equation:
Eq. 2-25
2.7.1. Interpretation as Rate of Return The value yield is the annualized rate of return on an investment over its life, given the current market price of that investment. This follows from the duality of the definition of the present value from section 2, in which the present value may be considered as the discounting of a future payout using a discount rate , or equivalently the future payout may be considered the result of exponential growth of the present value using as the growth rate. The choice of that makes the present value equal to the market-cap is called the value yield. This interpretation also extends to multiple future payouts that spread over a number of
3
The value yield is also called the Internal Rate of Return (IRR) but that may be confused with the concept of the Return on Equity (ROE) and is therefore not used here.
years or perhaps continuing for eternity, where the present value is merely the sum of all those future payouts discounted at the same rate. Note that the value yield is not the rate of return on reinvestment of future payouts, which will depend on the market price of the financial security at the time of such future payouts. The value yield may be identical for different holding periods and different timing of payouts as shown in the examples below. Individual investors may have different preferences for when they would like to receive payouts even though the annualized rate of return is identical. This is known as a utility function and is ignored here because it differs amongst investors. 2.7.2. Use as Discount Rate Because the value yield is the rate of return that can be obtained from an investment over a given holding period, the value yield can be used as the discount rate in calculating the present value of another investment so as to find their relative value. If the value yield changes with the holding period, as is the case in section 6.9, then the value yield must be chosen to match the period of the investment whose present value is being calculated. An example of such a present value calculation is given in section 9.6. 2.7.3. Payout Growth If the payout grows by the rate each year for eternity so that value yield is calculated using Eq. 12-2:
then the
Eq. 2-26
If the payout growth first starts next year so that
then Eq. 12-3 gives:
Eq. 2-27
In both cases the lower bound of the value yield is the growth rate
as
approaches infinity.
2.7.4. Examples Assume all payouts are zero and that the market-cap grows to become after years. This corresponds to a bond paying all of its interest when the principal is returned (also called zero coupon). It follows from Eq. 2-25 that the value yield equals the interest rate :
Eq. 2-28
Now assume and which corresponds to a bond paying annual interest (also called coupons) of and the principal is returned after years. It follows from Eq. 2-25 and Eq. 12-5 that the value yield again equals the interest rate :
Eq. 2-29
Now assume the payout grows by the rate each year so that the selling price after years grows to become then calculated using Eq. 2-25 and Eq. 12-5:
and . The value yield is
Eq. 2-30
This equals the value yield for eternal growth in Eq. 2-26 with Now assume the payout first starts growing next year so that the selling price after years is still
. but . The value yield is then:
Eq. 2-31
This equals the value yield for eternal growth in Eq. 2-27 with
2.7.5. Existence of Value Yield The so-called Cauchy-Hadamard theorem states that the power series used in the calculation of present value, converges if the payout growth does not exceed exponential growth indefinitely, which is a reasonable assumption for payout growth in the real world. According to that theorem, the present value exists for some continuous range of discount rates. In case the power series is finite, meaning all payouts are zero after some year, then the present value exists for all choices of discount rate except . In case all payouts are non-negative then the value yield exists, that is, the discount rate can be chosen so as to make the present value equal the market-cap, because the present value decreases continuously and monotonically from a limit of infinity when towards zero as the discount rate increases. The value yield may not be well-defined when some of the payouts are negative. For example, if positive payouts are followed by negative payouts (corresponding to raising new capital) and then zero payouts for eternity, which may be the pattern of a company going bankrupt, then the present value does not change monotonically with a changing discount rate because the positive payouts may dominate the present value when the discount rate is sufficiently high and the negative payouts may dominate the present value when the discount rate is sufficiently low. This means there is a maximum present value when varying the
10
discount rate. If the market-cap is higher than this maximum present value then the value yield is not welldefined. The value yield is well-defined in the Monte Carlo simulations in this paper and is found using a numerical root-finding method which is commonly available in mathematical software packages.
2.8.
Terminal Value
The present value of the equity growth model in Eq. 2-14 is defined from an infinite number of iterations, but the Monte Carlo simulation must terminate after a finite number of iterations. Estimating the present value is therefore done by separating into which is the present value of the payout in the years that have been Monte Carlo simulated, and which is an approximation to the remaining value if the Monte Carlo simulation had been allowed to continue for an infinite number of iterations:
Eq. 2-32
2.8.1. Iterations Required The number of iterations required to achieve a given accuracy of is now derived. To estimate the present value assume the starting payout is and the annual payouts grow by each for eternity so that . Let denote the partial present value where the year is between and :
Eq. 2-33
Note that:
Eq. 2-34
The geometric series
is reduced using Eq. 12-2:
Eq. 2-35
The terminal value consisting of the payouts in year
and onwards is calculated using Eq. 12-4:
Eq. 2-36
11
The limit of the terminal value is zero as
approaches infinity provided
Eq. 2-37
The limit of
approaches
as
approaches infinity:
Eq. 2-38
The number of iterations
required in order for
to be within relative error
of the value
is:
Eq. 2-39
In the derivation of this, the inequality is reversed because assumed here.
when
which is
If the payout growth rate is stochastic then the required number of iterations can be estimated using Eq. 2-39 with the mean growth rate. For example, in the equity growth model from section 2.2, the mean growth rate from Eq. 2-15 may be used. If the discount rate is the unknown variable that is sought determined, as is the case in the Monte Carlo simulation in section 5.2, then the discount rate is unknown when is calculated in Eq. 2-39. The discount rate may then be set to the expected discount rate or slightly above the mean growth-rate. Furthermore, may be increased by a margin because the Monte Carlo simulation may produce a value that is significantly different from the average used to derive here.
2.8.2. Mean Terminal Value Let be the stochastic payout in year which starts growing next year by the stochastic rate independently of . Note that a single growth rate is used for all years but the choice of that growth rate is stochastic. Let be a deterministic discount rate which is assumed to be greater than the growth rate . The mean terminal value is calculated using Eq. 3-13 with additional discounting to calculate the present value today:
Eq. 2-40
This would have been an equality if a new stochastic growth rate had been used each year and the growth rates were independent and identically distributed (IID), see section 3.3.1.
12
For the equity growth model in section 2.2, the mean growth rate
would be calculated from Eq. 2-15:
When the equity at the start of year
is known, the expected payout can be calculated from Eq. 2-13:
Eq. 2-41
No Growth If the equity at the beginning of year is known and all earnings are being paid out so the equity and payouts do not grow, then the mean payout is:
Eq. 2-42
As the mean payout is constant for eternity, the mean terminal value is:
Eq. 2-43
This approaches infinity as the discount rate approaches zero, while the mean terminal value in Eq. 2-40 approaches infinity as the discount rate approaches the mean growth rate .
3.
Stochastic Present Value
There may be uncertainty about future payouts in the present value calculation. There may also be uncertainty about the discount rate, which is the rate of return that can be obtained from the alternative investment used as the benchmark in calculating the present value. In such cases, the present value may be considered a stochastic variable taking on values according to some probability distribution . Using the mean payout growth and discount rates in a present value calculation does not result in the mean present value. This is because the present value formula is non-linear so Jensens inequality applies.
3.1.
Jensens Inequality
be a stochastic variable and let be a convex function, then Jensens inequality holds:
In general, let
Eq. 3-1
This becomes a strict inequality if
is strictly convex and
13
Let be the reciprocal function which is strictly convex for according to Jensens inequality:
and assume
, then
Eq. 3-2
Let
which is strictly convex for
, then according to Jensens inequality:
Eq. 3-3
Let
which is strictly convex for
and
, then according to Jensens inequality:
Eq. 3-4
Let
be another stochastic variable which is independent of , then Eq. 3-2 gives:
Eq. 3-5
This inequality also holds for powers of
Eq. 3-6
Jensens inequality has been studied previously for present value calculations, see e.g. Newell and Pizer [5].
3.2.
Stochastic Discount Rate
Let be a stochastic variable for the payout in year . Let be a stochastic variable for the discount rate assumed to be independent of . Note that a single discount rate is used for all years but the choice of that discount rate is stochastic. The present value is the derived stochastic variable:
Eq. 3-7
14
Using linearity of expectation, which applies to geometric series that converge, and Jensens inequality from Eq. 3-4, gives:
Eq. 3-8
That is, the present value may be underestimated when calculated from the mean discount rate. 3.2.1. Time-Varying Discount Rate Let be a stochastic variable for the discount factor in year so the present value is:
Eq. 3-9
Using Jensens inequality in Eq. 3-2 gives:
Eq. 3-10
That is, the present value may be underestimated when it is calculated from the mean discount rates.
3.3.
Stochastic Growth Rate
Let be a stochastic variable for the starting payout and let be a stochastic variable for the payout growth which starts next year. Note that a single growth rate is used for all years but the choice of that growth rate is stochastic. Assume the stochastic variables and are independent. The payout in year is the derived stochastic variable:
Eq. 3-11
Let be a deterministic discount rate greater than . The present value calculated using Eq. 12-3:
is a derived stochastic variable
Eq. 3-12
15
The mean present value is calculated using Jensens inequality in Eq. 3-2:
Eq. 3-13
That is, the mean present value may be underestimated when calculated from the mean growth rate. 3.3.1. Time-Varying Growth Rate Let be the payout growth rate in year which is assumed to start next year and be independent of the stochastic variable for the starting payout. The payout for year is the derived stochastic variable:
Eq. 3-14
Where the accumulated growth is denoted:
Eq. 3-15
Further assume the growth rates are independent and identically distributed (IID) so they have identical means . The mean accumulated growth can then be reduced to:
Eq. 3-16
The mean payout in year is then:
Eq. 3-17
Using this with Eq. 12-3 and linearity of expectation, which applies to geometric series that converge, gives the mean present value:
Eq. 3-18
Had the growth rates not been IID then the mean present value could not have been reduced in this manner. Also note that Eq. 3-13 is an inequality because a single stochastic growth rate is used for all years and the exponentiation of that growth rate invokes Jensens inequality.
16
3.4.
Stochastic Discount & Growth Rates
If the discount rate and the payout growth rate are independent stochastic variables then it follows from Eq. 3-8 (for a time-constant discount rate) and Eq. 3-10 (for a time-varying discount rate) that the mean present value may be underestimated when calculated from the mean discount and growth rates. This is because Jensens inequality is invoked for the stochastic discount rate regardless of whether the payout growth rate is stochastic or not.
3.5.
Stochastic Value Yield
Let be a stochastic variable for the starting payout and let be a stochastic variable for the payout growth which starts next year. Note that a single growth rate is used for all years but the choice of that growth rate is stochastic. The payout in year is a derived stochastic variable:
Eq. 3-19
Let be a stochastic discount rate greater than . The present value calculated using Eq. 12-3:
is a derived stochastic variable
Eq. 3-20
The value yield is a stochastic variable derived from the market-cap:
and
which makes the present value equal to
Eq. 3-21
Because
is a constant, the mean value yield is:
Eq. 3-22
Note that the mean value yield approaches
as the market-cap approaches infinity.
3.5.1. Time-Varying Growth Rate Let be the stochastic payout growth rate in year and assume the accumulated growth in year as defined in Eq. 3-15.
s are IID. Let
be the
17
The value yield is the discount rate that makes the present value equal to the market-cap, assuming it exists and is well-defined (see section 2.7.5):
3-23
This geometric series is assumed to converge, so we can select a single stochastic growth rate to construct a geometric series that converges to the same value according to Eq. 12-3:
and use it
3-24
This can be rearranged to derive the value yield from the stochastic variables
and
The mean value yield is then:
Eq. 3-25
The probability distribution of is unknown and may have another distribution than , for example, may be discrete while is continuous. The mean payout growth rate is and it is known from Eq. 3-17 that the mean payout growth rate also equals so the mean value yield is:
Eq. 3-26
This simple formula for the mean value yield is possible because the payout growth rates are assumed to be IID. Also note that the mean value yield approaches as the market-cap approaches infinity. 3.5.2. Equity Growth Model In the equity growth model from section 2.2, the mean starting payout is calculated using Eq. 2-13 assuming the stochastic variables and are independent of the starting equity:
Eq. 3-27
18
Using this in Eq. 3-26 with the mean payout growth rate from Eq. 2-15 gives the mean value yield:
Eq. 3-28
If instead becomes:
and
are independent stochastic variables, then the mean payout in Eq. 3-27
Eq. 3-29
Using this and the mean payout growth rate from Eq. 2-17 with Eq. 3-26 gives the mean value yield:
Eq. 3-30
3.5.3. Mean Value Yield as Discount Rate The value yield of one investment can be used as the discount rate for calculating the present value of another investment, thus giving the value of one investment relative to the other. For example, the value yield of the S&P 500 stock market index is used as the discount rate in calculating the present value of the Coca-Cola company in section 9.6.1. If the mean present value is wanted but not its probability distribution, then it would be desirable to use the mean value yield as the discount rate because it can be calculated from the simple Eq. 3-28 under certain assumptions. But using the mean value yield as the discount rate underestimates the mean present value due to Jensens inequality as shown in section 3.2. So the probability distribution of the value yield is still needed.
3.6.
Comparing Stochastic Value Yields
Comparing stochastic value yields is useful in estimating the probability that the value yield of one investment is greater than that of another investment with some risk premium. Let and denote the stochastic value yields of investments in a company and a stock market index, respectively. The probability of the companys value yield being the greatest is denoted:
Eq. 3-31
19
Eq. 3-32
In this form, the probability can be calculated directly when the value yields result from Monte Carlo simulations, by simply counting the number of value yields that satisfy the condition and dividing by the total number of simulations. This way of calculating the probability also has the advantage of working for value yields that are dependent, for example if the company is itself a part of the stock market index or if the value yields are dependent in some other way. 3.6.1. Alternative Calculation There is another way of calculating the probability of one value yield being greater than another. Let denote a value yield taken on by the stochastic variable so Eq. 3-31 becomes:
Eq. 3-33
If
and
are independent then this probability equals:
Eq. 3-34
The CDF for
is denoted
and defined as:
Eq. 3-35
This means:
Eq. 3-36
Using this with Eq. 3-34 gives:
Eq. 3-37
20
This formula is only valid when and are independent otherwise the joint probability distribution must be used. The significantly simpler Eq. 3-32 works in both dependent and independent cases and is therefore recommended.
4.
Option Valuation
An option is a financial instrument which gives the holder of the option the right but not the obligation to either buy (in case of a call option) or sell (in case of a put option) an underlying financial security at a predetermined exercise price (or strike price). Options have expiration dates after which they cannot be exercised. The period until the expiration date is known as the options life or maturity period. A call option will not get exercised if the share price is lower than the exercise price because the option holder could instead buy the share directly in the market at a lower price, so the exercise value is:
Eq. 4-1
Similarly, a put option will not get exercised if the share price is higher than the exercise price because the option holder could instead sell the share directly in the market at a higher price, so the exercise value is:
Eq. 4-2
The so-called intrinsic value of an option is the exercise value calculated using the current share price rather than the share price at the time of exercise. The difference between the intrinsic value and the option price is called the time value and is the premium that must be paid for the possibility that the share price will change sufficiently over time so as to increase the exercise value of the option. In some cases the time value is almost zero as shown in Figure 38. The profit (or loss) from exercising an option is the exercise value minus the options price:
Eq. 4-3
This is the profit from the perspective of an option buyer which may be made explicit by writing instead. The profit for the option seller is described below. When considering Monte Carlo simulated exercise values, the probability of a (positive) profit is calculated by counting the number of simulated profits that are positive and dividing them by the total number of simulated profits. The probability of profit is denoted:
Eq. 4-4
21
The present value of the exercise value is calculated using the discount rate (possibly not an integer) until the option is exercised:
and the number of years
Eq. 4-5
The value yield is the discount rate which makes the present value equal to the options price:
Eq. 4-6
If the exercise value is zero then the value yield is -1, which represents an annual rate of return that is a total loss. Call options are worthless when the market price for the share is below the exercise price, and put options are worthless when the share price is above the exercise price, see Eq. 4-1 and Eq. 4-2.
4.1.
Sellers Profit
The formulas above are from the perspective of an option buyer. From the perspective of an option seller (also known as an option writer), the profit (or loss) is the option price minus the exercise value:
Eq. 4-7
Using this with the definition of a buyers profit from Eq. 4-3 shows that the buyers profit is the sellers loss, and vice versa:
Eq. 4-8
If the exercise value is zero then the option price is all profit for the seller. The exercise value of a call option does not have an upper bound so the profit to a call option seller is upper bounded by the option price but not lower bounded. The exercise value of a put option is both lower and upper bounded because the price of the underlying stock cannot be lower than zero, so the put option sellers profit is lower bounded by the option price minus the exercise price and upper bounded by the option price. The probability of (positive) profit to the seller is defined from Eq. 4-7:
Eq. 4-9
Using this with Eq. 4-4 and assuming
gives:
Eq. 4-10
22
4.2.
Value Estimation
The exercise value of an option depends on the stock price at the time of exercise which is inherently unpredictable. Several models exist for estimating future stock prices and hence the value of options, with a popular one being the so-called Black-Scholes formula [6] [7] which assumes stock prices follow a random walk, a so-called geometric Brownian motion (GBM) with constant volatility, that stocks and options are correctly priced so there is no possibility for increasing expected returns without also increasing the risk, that there is unlimited borrowing ability, no transaction costs, interest rates are constant and known in advance, etc. Although several of these assumptions were criticized as being unrealistic by one of the original authors [8], the fundamental assumption of the Black-Scholes formula was not critiqued, namely that stock prices are assumed to follow a random walk in which changes to stock prices are independent log-normal distributed stochastic variables. Figure 1 shows the daily stock price returns for the S&P 500 stock market index for the period 1984-2011 which are clearly not log-normal distributed and section 6.8 shows that successive price changes are dependent. This means the Black-Scholes model is oversimplified and becomes increasingly inaccurate with the option duration, as also noted by Buffett [9]. The so-called binomial model by Cox et al. [10] uses a different model to forecast stock prices but it is still based on assumptions similar to the Black-Scholes model, namely that stock prices follow a random walk with a simple probability distribution, capital markets are perfect, etc., and the model therefore exhibits problems that are similar to those of Black-Scholes. Monte Carlo simulation can also be used to estimate the value of options by simulating future stock prices. This is typically done by simulating a stochastic process as a GBM variant, see e.g. Glasserman [1]. This paper takes another approach by first simulating the underlying economics of a company using probability distributions from historical data and then simulating the stock price as a function of these economics.
5.
Monte Carlo Simulation
Monte Carlo simulation is computer simulation of a stochastic model repeated numerous times so as to estimate the probability distribution of the outcome of the stochastic model. This is useful when the probability distribution is not possible to derive analytically, either because it is too complex or because the stochastic variables of the model are not from simple, well-behaved probability distributions. Monte Carlo simulation allows for arbitrary probability distributions so that very rare events can also be modelled. It was shown in section 3 that using the mean growth and discount rates will underestimate the present value. Monte Carlo simulation allows better estimation of the present value by repeating the calculation with the growth and discount rates selected at random according to their probability distribution.
23
5.1.
1. 2. 3. 4. 5.
Equity, Earnings & Payout
A single Monte Carlo simulation of future equity, earnings and payout, consists of these steps: Load historical data for and and determine probability distributions and dependencies. Determine the required number of Monte Carlo iterations using Eq. 2-39. Randomly generate and for to from historical data and a suitable model. Set . Calculate from Eq. 2-8, from Eq. 2-13 and from Eq. 2-10, for to .
The probability distribution is found by repeating steps 3-5 and recording the resulting values.
5.2.
Value Yield
A single Monte Carlo simulation of the value yield extends the simulation in section 5.1 by adding steps 6-7: 6. Calculate the terminal value from Eq. 2-40. 7. Use a numerical optimization method to find the value yield in Eq. 2-23. The probability distribution is found by repeating steps 3-7 and recording the resulting value yields.
5.3.
Present Value
Calculating the present value consists of first Monte Carlo simulating the payouts in section 5.1. The discount rate is the rate of return that can be obtained from an alternative investment, thus giving the value of the payouts relative to that alternative investment. The discount rate may be sampled from the probability distribution of the value yield that can be obtained from the alternative investment, see section 5.2. Monte Carlo simulations can also create a sequence of compounded returns from an alternative investment which can be used as the discount rate. The payouts and discount rate are used with Eq. 2-6 to calculate the present value of one Monte Carlo simulation and this is repeated to get the probability distribution.
5.4.
P/Book
A single Monte Carlo simulation of future P/Book ratios consists of these steps: 1. Determine the probability distribution for the daily P/Book change. 2. Initialize with some value. 3. For to first select from the P/Book change distribution and then calculate:
The probability distribution is found by repeating steps 2-3 and recording the resulting P/Book ratios.
5.5.
Price
The future equity and P/Book ratios resulting from the Monte Carlo simulations in sections 5.1 and 5.4 can be multiplied so as to simulate the share price, see Eq. 2-21. This in turn may be used to value options.
24
6.
S&P 500
The Standard & Poors 500 (S&P 500) stock market index consists of 500 large companies traded on the stock markets in USA and operating in a wide variety of industries including energy and utility, financial services, health care, information technology, heavy industry, manufacturers of consumer products, etc. The S&P 500 index may be used as a proxy for the entire USA stock market as it covers about 75% of that market.4
6.1.
Stock Market Forecasting
Estimating the stock markets future rate of return is important in many aspects of financial valuation. 6.1.1. Extrapolation of Total Return A simple approach to forecasting future returns of the stock market is to extrapolate its historical total return, that is, the annual rate of return from capital gains and reinvestment of dividends. But this approach is flawed because it ignores the current stock market price relative to the future earnings of the underlying companies, hence implicitly assuming that the stock market is correctly priced relative to its future earnings. 6.1.2. Equity Risk Premium Another approach to forecasting future stock market returns is to use an Equity Risk Premium (ERP) derived from historical observations of stock market returns relative to low-risk government bonds. Then adding the historical average ERP to the current yield on government bonds gives the expected future rate of return of the stock market. This is an attractive technique for estimating future stock market returns because it is simple and uses two known numbers to estimate an unknown number. Several variations to ERP estimation have been proposed for taking into account macro-economic indicators, supply and demand, etc., see e.g. Hammond et al. [11] and Damodaran [12]. Unfortunately, the fundamental idea of ERP is flawed, firstly because it is incorrect to use averages in valuations with non-linear calculations, see section 3, and secondly because the ERP changes unpredictably over time, see section 7.1. 6.1.3. Value Yield This paper takes another approach to estimating future stock market returns. Instead of trying to forecast the stock markets future rate of return from its historical total return or relation to government bonds, the stock markets accumulation of equity capital and its earnings resulting from that capital is modelled and Monte Carlo simulated. Combined with the stock markets historical P/Book ratio it results in probability distributions for the stock markets rate of return, or value yield, depending on the current stock price and the holding period of the stock. This circumvents the problem of determining the relationship between government bonds and the stock market, and through its use of historical data implicitly takes into account a wide range of potential changes in the macro economy, government bond yields, ERP, etc.
6.2.
Historical Data
Table 1 shows the historical Return on Assets (ROE), Return on Equity (ROE), earnings retaining ratio (Retain) and equity-to-assets for the S&P 500 index during the period 1984-2011. The summary statistics are shown in Table 2 and suggest that the earnings of the S&P 500 index are related to its assets which in
4
S&P 500 Fact Sheet, retrieved April 11, 2013: www.standardandpoors.com/indices/articles/en/us/?articleType=PDF&assetID=1221190434733
25
turn are related to its equity, and this means the equity growth model from section 2.2 can be used to simulate the future equity, earnings and payouts. Appendix 12.3 gives details about the compiling and limitations of this data, most notably that the calculations using e.g. equity and earnings of the S&P 500 companies are un-weighted because the weights used in the official S&P 500 index are proprietary, so the equity and earnings are merely the sum of the equity and earnings of all the companies in the S&P 500 index. This is not a problem as it still provides a useful estimate of the rate of return to be expected from investing in a broad stock market index, it is just slightly different from the official S&P 500 index. 6.2.1. Sampling Statistical dependencies of the data in Table 1 should be modelled because the data is being repeatedly sampled in the Monte Carlo simulations. The relations between ROE and Retain when selected from the same year or time-shifted one year are shown in the scatter-plots of Figure 2. When ROE and Retain are selected from the same year they tend to increase together, while no such tendency is apparent when ROE and Retain are time-shifted from each other. Although the data set is small and this tendency may not be statistically significant, it seems more reasonable to select ROE and Retain from the same year when sampling the data in the Monte Carlo simulations. The autocorrelations for ROE and Retain are shown in Figure 3 and Figure 4, which suggest that there is no statistically significant linear correlation over time. This means ROE and Retain can be randomly sampled in the Monte Carlo simulations. 6.2.2. Mean Growth Rate Because ROE and Retain are sampled from the same year, the mean equity growth rate calculated using Eq. 2-15 with the ROE and Retain pairs from Table 1:
can be
Eq. 6-1
The growth rate of the mean earnings and payouts is the same because they are derived from the equity. If ROE and Retain were instead sampled independently then the mean growth rate would be calculated using Eq. 2-17 with the mean ROE and Retain from Table 2:
6.3.
Curve Fittings
The results of the Monte Carlo simulations below are curve fitted so as to provide simple formulas for estimating various probability distributions from input variables such as the P/Book ratio. No theoretical justification is given for selecting those particular functions for the curve fittings. The guiding principle has merely been to select functions that are simple and provide a good fit to the observed data.
26
6.4.
Earnings
The probability distribution for the future earnings of the S&P 500 stock market index can be estimated using the Monte Carlo simulation from section 5.1 which is repeated 1,000 times here. 6.4.1. Mean & Standard Deviation Figure 5 and Figure 6 show histograms of the earnings in different future years resulting from these Monte Carlo simulations. The earnings are evidently not from a simple probability distribution. The earnings mean and standard deviation increase over time as shown in Figure 7 where the fitted functions are:
Eq. 6-2
Eq. 6-3
The growth rates for the earnings mean and standard deviation are:
Eq. 6-4
Eq. 6-5
That is, the mean earnings increase about 4.46% per year and the standard deviation about 4.71%. The growth rate for the mean earnings is the same as for the mean equity, see section 6.2.2. Adjusting for Normalized Equity Because the equity starts at 1 in these Monte Carlo simulations, the stochastic variable for the earnings must be multiplied by the actual starting equity which is denoted here. The adjusted mean earnings are derived from Eq. 6-2 using linearity of the expectation operator:
Eq. 6-6
The standard deviation of the adjusted earnings is derived using Eq. 6-3 with a well-known property of the standard deviation and assuming :
Eq. 6-7
27
For example, at the time of this writing the last known equity for the S&P 500 index was USD 661.93 per share on September 28, 2012,5 so the expected earnings per share are calculated using Eq. 6-6:
The standard deviation is calculated using Eq. 6-7:
This assumes there is no change in the number of shares.
6.5.
Payout
The probability distribution for the future payout of the S&P 500 stock market index can be estimated using the Monte Carlo simulation from section 5.1 which is repeated 1,000 times here. 6.5.1. Mean & Standard Deviation Figure 8 and Figure 9 show Q-Q log-normal plots of the payouts in different future years resulting from these Monte Carlo simulations. In the first 50 years the payouts are evidently not from a simple probability distribution but as the years increase the Q-Q plots become increasingly linear which suggests the payouts approach a log-normal distribution. The payout mean and standard deviation increase over time as shown in Figure 10 where the fitted functions are:
Eq. 6-8
Eq. 6-9
The growth rates for the payout mean and standard deviation are:
Eq. 6-10
Eq. 6-11
That is, the mean payouts increase about 4.46% per year and the standard deviation about 4.73%. The growth rate for the mean payouts is the same as for the mean equity, see section 6.2.2.
S&P Dow Jones Indices (retrieved early March 2013): http://us.spindices.com/documents/additional-material/sp-500-eps-est.xls
28
Adjusting for Normalized Equity Because the equity starts at 1 in these Monte Carlo simulations, the stochastic variable for the payout must be multiplied by the actual starting equity which is denoted here. The adjusted mean payout is derived from Eq. 6-8:
Eq. 6-12
The standard deviation of the adjusted payout is derived using Eq. 6-9:
Eq. 6-13
For example, at the time of this writing the last known equity for the S&P 500 index was USD 661.93 per share on September 28, 2012,6 so the expected payouts per share are calculated using Eq. 6-12:
This assumes there is no change in the number of shares. 6.5.2. Payout Sum The cumulative sum of payouts for the years 1 to is:
Eq. 6-14
Figure 11 shows histograms for the payout sums resulting from the Monte Carlo simulations described above. Also shown are the fitted log-normal distributions which indicate that the probability distributions approximate the log-normal distribution as the years increase. Assume 12 shows the parameters functions: for year is a log-normal distributed stochastic variable. Figure when the year ranges between 1 and 290. Also shown are the fitted
and
Eq. 6-15
S&P Dow Jones Indices (retrieved early March 2013): http://us.spindices.com/documents/additional-material/sp-500-eps-est.xls
29
Eq. 6-16
Adjusting for Normalized Equity The stochastic variable must be multiplied by the actual starting equity which is denoted without a subscript, because the equity in these Monte Carlo simulations starts at one for normalization purposes. Using a well-known property of the log-normal distribution, the stochastic variable is also log-normal distributed with parameter derived from Eq. 6-15:
Eq. 6-17
The parameter
is the same as in Eq. 6-16.
6.6.
Equity
The probability distribution for the future equity of the S&P 500 stock market index can be estimated using the Monte Carlo simulation from section 5.1 which is repeated 1,000 times here. 6.6.1. Log-Normal Parameters Figure 13 and Figure 14 show Q-Q log-normal plots for the equity in different years resulting from these Monte Carlo simulations. As the years increase the Q-Q plots approach straight lines, thus indicating the equity is approximately log-normal distributed. Let the starting equity be and assume for year is a log-normal distributed stochastic variable. Figure 15 shows the parameters and resulting from the Monte Carlo simulations when the year ranges between 1 and 290. Also shown are the fitted functions:
Eq. 6-18
Eq. 6-19
For example, in year parameters:
the equity is estimated to be a log-normal distributed stochastic variable with
Figure 16 compares this distribution to the results of the Monte Carlo simulation. Adjusting for Normalized Equity The stochastic variable must be multiplied by the actual starting equity which is denoted without a subscript, because was set to equal one in these Monte Carlo simulations for
30
normalization purposes. Using a well-known property of the log-normal distribution, the stochastic variable is also log-normal distributed with parameter derived from Eq. 6-18:
Eq. 6-20
The parameter
is the same as in Eq. 6-19. is:
For example, the last known equity at the time of this writing is USD 662 per share so
This assumes there is no change in the number of shares. 6.6.2. Mean & Standard Deviation Using the log-normal parameters and gives:
for
from Eq. 6-18 and Eq. 6-19 with Eq. 12-6 and Eq. 12-7
Eq. 6-21
Eq. 6-22
The mean equity grows by the rate:
Eq. 6-23
This almost equals the mean growth rate calculated directly from the ROE and Retain data, see section 6.2.2, with the small difference likely due to sampling and rounding error. Adjusting for Normalized Equity Because the equity starts at 1 in these Monte Carlo simulations, the stochastic variable must be multiplied by the actual starting equity which is denoted here. The adjusted mean equity is derived from Eq. 6-21:
Eq. 6-24
The standard deviation of the adjusted equity is derived from Eq. 6-22:
Eq. 6-25
31
For example, the last known equity for the S&P 500 index at the time of this writing is USD 662 per share so the expected equity per share is calculated using Eq. 6-24:
This assumes there is no change in the number of shares.
6.7.
Value Yield
The value yield of the S&P 500 stock market index is the annual rate of return an investor would obtain from buying shares at the given price and owning the shares for eternity. The value yield can be used as the discount rate when calculating the present value of another investment, which then values that investment relative to the return that could be obtained from investing in the S&P 500 index. The actual value yield of the S&P 500 index is not known in advance because it depends on future earnings and payouts of the constituent companies of the S&P 500 index which are not known in advance. The probability distribution of the value yield is therefore estimated using the Monte Carlo simulation from section 5.2. A Monte Carlo run consists of 1,000 simulations with the number of iterations (years) per simulation set to 300 as calculated from Eq. 2-39 for the S&P 500 data-set. 6.7.1. Log-Normal Parameters Figure 17 shows the value yield histograms resulting from Monte Carlo runs with different P/Book ratios. Also shown are the fitted log-normal distributions. Figure 18 shows the Q-Q log-normal plots which appear to be approximately linear except for a few outliers that seem to be caused by error in the numerical optimization of the value yield and are therefore removed in subsequent experiments. The value yields appear to be approximately log-normal distributed. Assume is a log-normal distributed stochastic variable with parameters and . Figure 20 shows the parameters and resulting from a series of Monte Carlo runs when P/Book ranges between 0.2 and 10. Also shown are the fitted functions:
Eq. 6-26
Eq. 6-27
For example, the P/Book ratio is about 2.3 at the time of this writing in late January 2013, which gives:
32
That is, an investor who bought shares in the S&P 500 stock market index at a P/Book ratio of 2.3 could expect the value yield to be log-normal distributed with parameters and . 6.7.2. Mean & Standard Deviation Figure 19 shows the Q-Q normal plots (not log-normal) which appear to be approximately linear except for a few outliers that seem to be caused by error in the numerical optimization of the value yield and are therefore removed in subsequent experiments. The value yields appear to be approximately normal distributed. Figure 21 shows the value yield mean and standard deviation resulting from a series of Monte Carlo runs when P/Book ranges between 0.2 and 10. Also shown are the fitted functions:
Eq. 6-28
Eq. 6-29
For example, when
the value yield mean and standard deviation is:
That is, an investor who bought shares in the S&P 500 stock market index at a P/Book ratio of 2.3 could expect a value yield around 8.4% with standard deviation about 0.4%. Table 3 shows the minimum and maximum and for each year in the period 1984-2011, calculated using Eq. 6-28 and Eq. 6-29 with each years minimum and maximum P/Book ratios from Figure 22. During the period 1984-2011 the average was about 8.3% and the average was about 0.4%. 6.7.3. Comparison to Analytic Mean In the Monte Carlo simulations above, ROE and Retain were sampled in pairs from the same year because they are believed to be dependent in that manner. Using the ROE and Retain data pairs from Table 1 with Eq. 3-28 gives the mean value yield:
Eq. 6-30
33
This is slightly higher than the mean value yield in Eq. 6-28 derived from the Monte Carlo simulations. The difference is likely due to sampling and rounding error.
6.8.
Price
The probability distribution for the price of the S&P 500 stock market index at a given future time can be estimated by multiplying the P/Book ratio with the equity at that time. Recall that is just another notation for so multiplying this ratio with the equity results in the market-cap which is the price of all the shares:
Eq. 6-31
This may be divided by the number of shares to give the price per share. These terms and notations will be used interchangeably in the following. 6.8.1. P/Book Data In the data-set for the S&P 500 index used here, the book-values are only known at the last trading day of each year during 1983-2011 and the daily book-values are interpolated from these. It is reasonable to use interpolated book-values because the actual book-values are not very volatile during a year. The volatility in the P/Book ratio is caused by the volatility of the share-price. The daily price is available from another dataset7. Dividing the daily prices with the interpolated book-values estimates the P/Book ratio on each trading day from the beginning of 1984 to the end of 2011 as shown in Figure 22. 6.8.2. Constant P/Book Assuming the P/Book ratio is constant, the future price of the S&P 500 index is estimated from Eq. 6-31 by multiplying the constant P/Book ratio with the future equity at that time which is estimated in section 6.6. 6.8.3. Historical P/Book Distribution Figure 23 shows the historical distribution of P/Book ratios for the S&P 500 index during the period 19842011, which is clearly not log-normal or any other simple distribution. If a future P/Book ratio is assumed to be chosen randomly from the historical distribution then the price of the S&P 500 index is estimated from Eq. 6-31 using samples of the historical P/Book ratio multiplied by samples of the equity distribution, which are either taken from Monte Carlo simulations or estimated from the formulas in section 6.6. Figure 24 shows histograms for the price distributions in different years, where each distribution results from multiplying one million samples of the historical P/Book ratios with samples of the Monte Carlo simulated equity in that year. The price distribution can be seen to approach a log-normal distribution as the years increase. These price distributions assume that the equity starts at 1 and must therefore be multiplied by the actual equity which is USD 662 per share at the time of this writing. 6.8.4. Simulated P/Book Distribution Figure 25 shows the autocorrelation of the P/Book ratio for the S&P 500 index in the period 1984-2011. This shows that successive P/Book ratios are strongly correlated which should be taken into account when estimating future P/Book ratios. Figure 26 shows the distribution of daily P/Book changes, that is, the
7
http://research.stlouisfed.org/fred2/series/SP500/
34
P/Book of one day divided by the P/Book of the previous day, which is a narrow bell-shaped distribution with mean 1.000113 and standard deviation 0.01175446, but it is neither normal nor log-normal distributed (Figure 26 has a log-normal fitting). Also shown in that figure is a scatter-plot of the P/Book versus P/Book change for each day in the period 1984-2011, which exhibits dense clusters around certain P/Book ratios such as 1.8, 2.8 and 4.8, etc. Also note the outlier with a daily change of 0.8, that is, a decrease in the price of the S&P 500 index of about 20% which occurred on October 19, 1987. A simple way of modelling the dependency amongst successive P/Book ratios is to perform a Monte Carlo simulation as outlined in section 5.4, where the P/Book change is sampled from the distribution in Figure 26. The distribution of P/Book changes is found by first separating the P/Book ratios into bins similar to the histogram in Figure 23 and for each bin assemble the corresponding P/Book changes. The Monte Carlo simulation starts with a preset P/Book ratio and samples from the appropriate bin to get the P/Book change which is then multiplied with the P/Book ratio to get the P/Book ratio for the next day; this is repeated for each day being simulated. Figure 27 shows examples of the P/Book time-series resulting from these Monte Carlo simulations when the starting P/Book ratio is set to 2.3. Compare these to the historical time-series in Figure 22 which indicates that the simulated P/Book time-series are considerably more volatile. Figure 28 shows the autocorrelations for the P/Book time-series resulting from the Monte Carlo simulations. Compare these to the historical time-series in Figure 25 which again indicates that the simulated P/Book time-series are more volatile. Figure 29 shows the histograms for the simulated P/Book ratios in future years. Figure 30 shows the corresponding CDF. Figure 31 shows Q-Q plots comparing the distributions for the simulated P/Book ratios from Figure 29 to the historical P/Book ratios from Figure 23, which indicates that the simulated distribution approaches the historical distribution as the years increase. In these Monte Carlo simulations the starting P/Book ratio was set to 2.3.Had the starting P/Book ratio been different then the simulated P/Book ratios would have been affected in the short-term but in the long-term the P/Book distributions are unaffected by the starting P/Book ratio. This is because the simulation model is memory-less, meaning that only the current P/Book is considered when sampling the next P/Book, so the time-series would merely start at another P/Book ratio but otherwise exhibit similar long-term behaviour. Figure 32 shows the price distribution in future years where each distribution is estimated from one million samples of the P/Book ratios resulting from the above Monte Carlo simulation and the equity resulting from the Monte Carlo simulations described in section 6.6. Compare these to the price distributions in Figure 24 where the P/Book ratio is instead sampled from its historical distribution. In the near future the distributions depend on the starting P/Book ratio for the Monte Carlo simulations, but as the years increase both price distributions approach the log-normal distribution with similar parameters, because the P/Book distribution resulting from Monte Carlo simulations approaches the historical P/Book distribution. Figure 33 shows the CDFs for the price of the S&P 500 index in future years where the near future again depends on the starting P/Book ratio. Figure 34 shows the price mean and standard deviation for different years and starting P/Book ratios, and subsets are shown in Table 5 and Table 6. Because the starting equity was assumed to be one, the prices in these plots and tables must be multiplied with the actual equity which is
35
USD 662 per share at the time of this writing in early 2013. For example, Figure 33 shows that in five years from now the probability is about 0.07 that the price is below 2 and the probability is about 0.95 that the price is below 4, that is, the probability is (or 88%) that the price is between 2 and 4. This should be multiplied with the equity of USD 662 per share, so the probability is 0.88 that the price per share of the S&P 500 index is between USD 1,324 and USD 2,648 five years from now, that is, in the beginning of year 2018. This assumes there will be no change in the number of shares outstanding.
6.9.
Payouts & Price
Consider the scenario where shares of the S&P 500 index are bought and held for years after which they are sold. The shareholder will receive payouts during the holding period and will receive the proceeds from the sale of the shares. Let denote the current price of the shares and let be the price after years. Let be the payout in year . The value yield is the discount rate that makes the present value of the payouts and final selling price equal to the current share price, see Eq. 2-25:
Eq. 6-32
This is dominated by the present value of the payouts when the holding period is large and it is dominated by the present value of the selling price when the holding period is small. The price after years can be calculated from the P/Book ratio multiplied by the equity at that time:
Eq. 6-33
Note that in these calculations, the prices and payouts are all divided by the starting equity normalization purposes. This means the starting price is actually the starting P/Book ratio instead.
for
6.9.1. Historical P/Book Figure 35 shows the value yield distributions for different choices of starting P/Book ratios and years . The shares are assumed to be sold after years and the selling price is calculated by sampling the historical P/Book ratios of the S&P 500 index from Figure 22 and multiplying the sampled P/Book ratio with the equity for the th year resulting from the Monte Carlo simulations described in section 5.1. As can be seen from the plots in Figure 35, when the starting P/Book is 1 and the shares are sold after two years then the mean value yield is more than 70%. That is, the mean annualized rate of return a shareholder would get from buying S&P 500 shares at a P/Book ratio of 1 and selling the shares after two years, is more than 70%. After 5 years the mean value yield is more than 30% and after 10 years it is more than 20%. If the starting P/Book ratio had instead been 3 then the mean value yield for a 2 year holding period would be around zero. For a 5 year holding period the mean value yield would be about 4% and for a 10 year holding period it would be about 6%. If instead the starting P/Book ratio had been 5 then the mean value yield for a 2 year holding period would be negative by almost -25%. After 5 years the mean value yield would be about -6% and after 10 years it would be about zero. This shows that the mean value yield depends on the starting P/Book ratio and the duration of the holding period.
36
Regardless of the starting P/Book ratio, the standard deviation for the value yield is high if the shares are sold after only a few years, which is because the price of the S&P 500 index is more volatile than the annual payouts and the selling price dominates the value yield as calculated in Eq. 6-32 when the holding period is short. The standard deviation for the value yield decreases as the holding period increases. 6.9.2. Simulated P/Book As shown in section 6.8.4, successive P/Book ratios of the S&P 500 index are strongly correlated. Figure 36 shows the value yield distributions when the selling price is calculated using the P/Book ratios resulting from the Monte Carlo simulations in section 6.8.4. As the holding period increases, the value yield distributions approach those from Figure 35 which used random samples from the historical P/Book distributions. This is because the simulated P/Book distribution approaches the historical distribution as the number of years increase. For short holding periods the value yield distributions are quite different. Figure 36 shows for a starting P/Book ratio of 1 and a holding period of 2 years, the mean value yield is about 40% compared to more than 70% in Figure 35. For a starting P/Book ratio of 3 and a 2 year holding period the mean value yield is about 4% compared to zero in Figure 35. For a starting P/Book ratio of 5 and a 2 year holding period the mean value yield is about -3% compared to -25% in Figure 35. Although the value yield distributions in Figure 36 are obviously not normal, it is useful to plot the mean and standard deviations for varying P/Book ratios and holding periods. This is done in Figure 37 and subsets of the numbers are shown in Table 7 and Table 8. The standard deviation for the value yield has a clear tendency to decrease as the number of years increase and is almost independent of the starting P/Book ratio. The mean value yield is more complicated. If the starting P/Book ratio is lower than about 3 then the mean value yield decreases as the number of years increase. But if the starting P/Book ratio is around 3 or higher than 4 then the mean value yield increases as the number of years increase. Although there are also exceptions to this, for example as shown in Table 7, if the starting P/Book ratio is 4 then the mean value yield will decrease until the holding period is about 4 years after which it increases.
6.10. Option Valuation

Option valuation for the S&P 500 stock market index is done by considering the options for the SPY Exchange Traded Fund (ETF) which mimics the S&P 500 index. The share price of the SPY ETF is 1/10 of the S&P 500 index and its financial data such as equity per share is therefore also 1/10 of that of the S&P 500 index. Figure 38 shows the so-called option chains for SPY call and put options expiring on December 18, 2015 as quoted on March 15, 2013. These are the options with the farthest expiration date that are publicly available at the time of this writing. For simplicity, transaction costs are ignored in this study and the option prices are the averages of the bid and ask quotes. In practise, an option buyer should use the asking prices in valuation and an option seller should use the bid prices. The difference between bid and ask prices (or the spread) is sometimes large. The options will be valued on their expiration date only (also known as European options) but the valuation method can be extended to cover the entire life of the option (known as American options), or an average of the last period of the options life (known as Asian options), etc. Figure 39 shows the probability distribution of SPY share prices that result from Monte Carlo simulations similar to those described in section 6.8.4 for the S&P 500 index. The simulation period starts March 15, 2013 and ends December 18, 2015 which is the expiration date of the options being considered. This period has an estimated 696 trading days (about 252 trading days per calendar year and about 2.76 years in the
37
period considered). The last known equity for the S&P 500 index was USD 661.93 per share on September 28, 2012. According to Eq. 6-23, the mean annual growth rate for the equity of the S&P 500 index is 4.46% so the equity on March 15, 2013, which is about 6 months after the latest known equity, is expected to grow about 2.2% to become USD 676.49 per share, or USD 67.65 per share of the SPY ETF. The closing share price for SPY on March 14, 2013 was USD 156.32 so the P/Book ratio was about 2.31. This is the starting P/Book ratio for the Monte Carlo simulations of future P/Book ratios. The equity is simulated from its last known date on September 28, 2012 until September 28, 2015 and there is then an additional 81 days until the expiration date of the options on December 18, 2015. The simulated equity on September 28, 2015 is therefore multiplied by the expected growth rate for the remaining 81 days, which is about 1%. Figure 40 shows the CDF for the profits on the SPY call options with different exercise prices, and Figure 41 shows the same for the put options. The profits are calculated for the option expiration date only, using Eq. 4-3 with the share prices from the Monte Carlo simulations described above and the option prices from Figure 38. Figure 42 shows the profit means and standard deviations for the SPY call options with various exercise prices, and Figure 43 shows the same for the put options. For the call options, the maximum profit mean occurs when the exercise price is at or below USD 100 per share where the profit mean is almost USD 25 per share and the standard deviation is about USD 40 per share. The minimum profit mean is about USD 1 per share which occurs when the exercise price is USD 250 per share. For the put options, the profit means are all negative which means the put options are all expected to result in a loss. The lowest mean loss on the put options occurs when the exercise price is USD 10 per share where the profit mean is just below zero and the standard deviation is zero, this is because the option price is almost zero, see Figure 38. The biggest mean loss on the put options is almost USD -30 per share which occurs when the exercise price is USD 250 per share. Figure 44 shows the probability of the profit being positive, that is, for the call and put options. For example, for the SPY call option with an exercise price of USD 150 per share the probability of (positive) profit is about 0.55 (or 55%). For a call option with an exercise price at or below USD 110 the probability of profit is almost 0.68 (or 68%). For a call option with an exercise price of USD 250 the probability of profit is about 0.05 (or 5%). For a put option with an exercise price of USD 150 the probability of profit is about 0.075 (or 7.5%). For a put option with an exercise price at or below USD 100 the probability of profit is zero because the share prices in the Monte Carlo simulations were never less than USD 100, see Figure 39. For a put option with an exercise price of USD 250 the probability of profit is about 0.25 (or 25%). Figure 45 shows the CDFs for the value yields of SPY call options for different exercise prices, and Figure 46 shows the same for the put options. The value yield is the annual rate of return that would be obtained from buying an option at the given price and holding it until the expiration date at which point it is exercised if it has positive value, otherwise it expires worthless. The value yield is calculated for the expiration date only, using Eq. 4-6 with the share prices resulting from the Monte Carlo simulations described above. In many of these simulations the call options become worthless because the simulated share prices are below the exercise price and the option price is then a loss to the buyer of the call option. This makes the value yield calculated from Eq. 4-6 equal to -1. A put option also has value yield -1 if it expires worthless, which happens when the share price is higher than the exercise price. Figure 47 shows the value yield means and standard deviations for the call options and Figure 48 shows them for the put
38
options. For the call options, the maximum value yield mean occurs when the exercise price is around USD 110 per share where the value yield mean is about 11% and the standard deviation is about 25%. For the put options, all the value yield means are negative because the put options are expected to cause a loss on average. The maximum value yield mean for the put options occurs when the exercise price is USD 250 per share where the value yield mean is about -17% and the standard deviation is around 25%. 6.10.1. Comparison of Puts & Calls According to Figure 44 the call options with exercise prices at or below USD 165 have probabilities of profit above 0.5 (or 50%), while the put options all have profit probabilities below 0.25 (or 25%) and most put options have much lower or even zero profit probability. Figure 42 shows that the mean profit is always positive for the call options and Figure 43 shows that the mean profit is always negative for the put options. The reason is that the Monte Carlo simulated share prices used in these option valuations, as shown in Figure 39, are often higher than the combined exercise and option prices. If the equity growth model used in these Monte Carlo simulations is deemed a reasonable way of estimating the probability distribution of the future equity, earnings and payout for the S&P 500 index (and hence the SPY ETF), and if the simulation of future P/Book ratios from their historical probability distribution is also deemed reasonable, then it means the options must be mispriced; the call options are underpriced and the put options are overpriced. 6.10.2. Trading Strategies An option buyer would want to maximize both the profit probability and mean profit while simultaneously minimizing the option price. An option seller would want to do the opposite. Both buyers and sellers would want to minimize the standard deviation of profits. This is known as a multi-objective optimization problem because several objectives must be optimized simultaneously. Figure 49 shows the call option price versus profit mean and probability which both increase as the option price increases until the option price is slightly over USD 50, which corresponds to an exercise price of about USD 105, see Figure 38. When the option price is above USD 50 the profit mean and probability do not change (actually, the mean profit decreases slightly), which means that the call option would be more expensive but it would not lead to a higher mean profit or profit probability. An investor would want to invest as little capital as possible for the same profit mean and probability so the call options priced above USD 50 should not be bought. Figure 47 is an alternative way of reaching this conclusion by considering the mean value yield (annualized rate of return) that could be obtained from buying the different call options at their respective prices. An investor would want the highest rate of return on the capital used to purchase the call option and Figure 47 shows that the maximum value yield mean for a call option is when the exercise price is USD 105 where the value yield mean is about 10% and the standard deviation about 24%. This may be compared to the value yield estimates for buying the underlying stock as shown in Figure 37 for the S&P 500 index (which the SPY shares mimic) with various starting P/Book ratios and holding periods. In the option valuations here, the starting P/Book ratio is about 2.3 and the period before the options expire is 2.76 years. Interpolating the nearby values from Table 7 and Table 8 gives for the underlying stock a value yield mean about 10% and standard deviation about 8%. The value yield mean is about the same for the call option and the underlying stock in this case but their standard deviations differ greatly with the stock having about 8% and the call option about 24%, which means the value yield of the call option may be
39
significantly larger or smaller than the value yield of the stock. This is because an option is effectively a gearing of the stock. Figure 50 shows put option price versus profit mean and probability. In this example the profit means are all negative and the profit probability is below 0.25 (or 25%) for all the put options. This means a loss is expected from buying these put options and they will give a corresponding profit to the option seller, as shown in section 4.1. A seller of put options would want to maximize the option price while minimizing the buyers profit mean and probability. Figure 50 shows that the put option buyers mean profit is inversely related to the profit probability, that is, the mean profit decreases when the profit probability increases. This makes it impossible to simultaneously minimize both the mean profit and probability. The seller of put options must therefore decide whether to sell a put option at a high price and with a comparatively high probability of loss, or sell a put option at a low price and with a comparatively small probability of loss. The profit probability for the put seller is almost 100% when the option price is less than about USD 5 which is when the exercise price is less than USD 105, see Figure 38. So it is better to sell the put option with exercise price USD 105 than the options with lower exercise prices because the seller will receive more in sales proceeds with the same profit probability. Because the profit probability and mean are inversely related, the choice of which put option to sell will depend on the individual seller and what the sales proceeds will be used for. If the seller cannot reinvest the sales proceeds elsewhere, then the seller may want to minimize the probability of loss upon exercising of the options and should therefore sell the options with exercise price USD 105 and option price USD 5 which have almost zero probability of a loss to the option seller. But if the seller can invest the sales proceeds elsewhere at a sufficiently high rate of return and the put option can only be exercised on its expiration date, then the put option with an exercise price of USD 250 will give the option seller USD 100 in sales proceeds to reinvest elsewhere. The probability of a loss to the option seller is about 26% upon exercising of such a put option, see Figure 50. But this probability as well as the magnitude of any loss both decrease if the sales proceeds are reinvested elsewhere at a sufficiently high rate of return until the put option is exercised.
7.
USA Government Bonds
Figure 51 shows the historical yield on USA government bonds for the periods 1798-2012 (average annual yields) and 1962-2013 (average daily yields). During this period the bonds have had varying maturity period, terms and taxation. Table 9 shows basic statistics and a mean yield about 4.5-5% with standard deviation about 2%. Figure 52 shows the histograms of these historical bond yields and Figure 53 shows the autocorrelation plots with significant correlation between successive bond yields.
7.1.
Equity Risk Premium
Section 6.1.2 discussed the so-called Equity Risk Premium (ERP), which is a simple way of forecasting future stock market returns from the historical average premium of stock market returns relative to low-risk government bond yields. There are different ways of calculating the ERP, a few of which are now considered.
40
Figure 54 shows for the period 1957-2013 the ERP calculated as the difference between the monthly returns of the S&P 500 stock-market index and the monthly yield on USA government bonds with 10 year maturity periods, which have been converted from annual to monthly yields using the formula:
Eq. 7-1
Figure 54 also shows the histogram which is not normal distributed. Figure 55 shows the scatter-plot and cross-correlation for this data which show no significant correlations. Figure 56 shows the historical government bond yields for the period 1984-2011 from Figure 51 along with the mean value yield for the S&P 500 stock market index as estimated from Eq. 6-28 with the P/Book data from Figure 22. Their difference is an estimate of the ERP for eternal shareholders which is also shown in the plot and ranges from about -2% around year 1985 to almost +8% in years 2009 and 2011. The arithmetic mean ERP is about 2.1%, geometric mean slightly less, harmonic mean about 4.8%, and standard deviation about 1.8%. The ERP is not normal distributed as shown in the histogram in Figure 57. Figure 58 shows the scatter-plot and cross-correlation of the S&P 500 mean value yields and the USA government bond yields during the period 1984-2011, which shows significant cross-correlation but not in a simple pattern that is predictable from this data alone. Recall that Figure 37 and section 6.9 showed that the value yield of the S&P 500 index is affected by its current P/Book ratio and the holding period, so the ERP is also affected by this and hence cannot be a single number. All this corroborates the argument in section 6.1.2 that future stock market returns cannot be reliably estimated from merely adding the historical ERP to the current yield on long-term government bonds. For example, at the time of this writing in March 2013 the yield on USA government bonds with 10 year maturity is 1.94%. If the historical mean ERP of about 2.1% is merely added to the bond yield then it would suggest a mean value yield for the S&P 500 index about 4%, while in fact the P/Book ratio for the S&P 500 index is about 2.3 thus giving a mean value yield around 8.4% according to Eq. 6-28 which gives an ERP of 6.5%. There is currently much speculation in the news-media from even the most experienced and respected people in finance that if the interest rate increases then stock prices would decrease. But the current government bond yields around 2% are historically low according to Table 9 and the S&P 500 P/Book ratio and hence its mean value yield is also lower than its historical average. If both of these are to revert to their historical averages in the future then it would mean that the government bond yields should increase (bond prices decrease) and stock prices should increase (value yields decrease). But as shown with the many statistical studies in this paper it is not certain that this will happen.
8.
Dow Jones Venture Capital Index
The Dow Jones Venture Capital (DJVC) / Sand Hill index shown in Figure 59 covers approximately 18,000 venture companies during the period 1991-2010. This is a so-called total return index which measures the return, before taxes and fees for the venture capital funds through which investments are typically made,
41
from investing in all the venture companies of the index in proportion to their market-cap, and re-investing all payouts from public offerings and acquisitions into the remaining venture companies.
8.1.
Monthly Returns
The monthly return to the DJVC index is calculated as:
Eq. 8-1
The monthly returns are shown in Figure 59, from which it can be seen that the volatility of the DJVC index increased during the years 1999-2003 which is known as the Dot-Com bubble where the market price increased greatly and then collapsed for both venture and established public companies. The index was also more volatile during a severe financial crisis in the years 2008-2009. Figure 60 shows a histogram of the monthly returns which have arithmetic mean 1.65%, geometric mean 1.24% and a comparatively high standard deviation of 9%. The compounded annual return has geometric mean about 16%. Also shown in Figure 60 is the autocorrelation of monthly returns which shows no significant correlation between successive monthly returns.
8.2.
Risk Premium
The monthly returns of the DJVC index are highly volatile and future returns to the index appear to be unpredictable from its prior returns. It would therefore be desirable to forecast future returns from the historical risk premium on another investment whose future return is known. 8.2.1. USA Government Bonds Figure 61 shows the scatter plot of monthly returns of the DJVC index and the yield on USA government bonds with 10 year maturity periods, which have been converted from annual to monthly yields using Eq. 7-1. This shows the bond yields are uncorrelated with the return experienced on the DJVC index. Figure 61 also shows the cross-correlation between the DJVC index and government bond yields, that is, the correlation of the two time-series with varying time-lags. Although there are some statistically significant correlations, the correlation coefficient is at most 0.2 which is not enough to suggest a clear predictive relationship between the two time-series. This means the historical difference in return between the DJVC index and government bonds cannot be used as a predictable risk premium in forecasting future returns of the DJVC index. 8.2.2. S&P 500 Figure 62 shows the scatter plot of monthly returns of the DJVC venture index and the S&P 500 stockmarket index. Note that DJVC is a total return index while the S&P 500 index data used here is not because the monthly total return data for the S&P 500 index was not publicly available at the time of this writing. However, the scatter plot still shows a clear correlation with coefficient 0.87 and the linear fit has coefficient of determination and this should not be materially affected if the total return was used for the S&P 500 index. Also shown in Figure 62 is the cross-correlation which shows that there is only strong correlation between the DJVC and S&P 500 indices when the time-lag is zero.
42
Although a strong correlation was found between the monthly changes in the DJVC and S&P 500 indices, it is questionable whether this can be used to reliably predict future returns of the DJVC index. Figure 63 shows the highly volatile difference between the monthly changes of the DJVC and S&P 500 indices with arithmetic mean about 1.1% (percentage points) and standard deviation about 5.7% (percentage points). Furthermore, the correlation is merely a statistical relationship for an observation period of about 18 years and there is no intrinsic argument or proof that this will continue in the future.
8.3.
Simulating Future Returns
As shown above, consecutive monthly changes of the DJVC index are uncorrelated so the simulation of future changes is done by randomly sampling the historical monthly changes of the DJVC index and multiplying the samples to calculate the compounded returns. Figure 64 shows examples of this.
8.4.
Discount Rate
Valuing an investment relative to the return that could be obtained from the DJVC index is done by simulating future returns of the DJVC index as described above and using the compounded return as the discount rate in the present value calculation. An example of this is given in section 9.6.
9.
Coca-Cola
The company Coca-Cola (abbreviated as Coke) was incorporated in 1919 in USA but its origin is significantly older. The company produces beverages that are sold worldwide. In 2012 the company employed 150,000 people and had revenue of about USD 48b with net income of about USD 9b.8 At the time of this writing in late March 2013, Cokes market-cap is about USD 180b and the latest known equity was almost USD 33b in December 2012, which means the P/Book ratio is currently about 5.5.
9.1.
Simulation Settings
The companys future equity, earnings and payouts are Monte Carlo simulated using the equity growth model from section 2.2 which was also used to simulate the S&P 500 index in section 6. Although the same aspects could be studied here as was studied in section 6, this section focuses on the value yield of the Coca-Cola company and its relation to the value yield of the S&P 500 index, as well as the present value. The ROE and Retain data from Table 10 is used in the Monte Carlo simulation described in section 5.2. Figure 65 shows that there is no significant correlation between the ROE and Retain data so they are sampled in pairs from the same year. The Monte Carlo simulation is repeated 1,000 times and each simulation is for 320 years, which is calculated using Eq. 2-39 so as to achieve sufficient numerical accuracy.
9.2.
Value Yield
The value yield is the annualized rate of return a shareholder would get from buying Coke shares at the given price and holding the shares for eternity. Although the payouts for the first year are only expected to be about 5.3% of the market-cap, the retained earnings will cause future payouts to increase exponentially. For example, after 10 years the mean payout is 16% of the initial share price, after 50 years the mean payout is 22 times (2200%) the initial share price, and after 100 years the mean payout is 58k (5,800,000%)
8
Form 10-K annual report for 2012 filed with US SEC: www.sec.gov/Archives/edgar/data/21344/000002134413000007/a2012123110-k.htm
43
times the initial share price. The exponentially growing payouts cause the value yield to be much greater than the first years rate of return of about 5.3%. Figure 66 shows the probability distribution of the value yield which has mean 18.3% and standard deviation 1.5%. However, this rate of return is extremely high and unrealistic as shown next.
9.3.
Equity Growth
On September 28, 2012, the equity of the S&P 500 stock market index was about USD 5,918b9 and the equity of Coke was about USD 33b.10 The mean equity growth rate of the S&P 500 index for the period 1984-2011 was 4.46% as calculated in Eq. 6-1. The mean equity growth rate for Coke is calculated using Eq. 2-15 with the data from Table 10 for the period 1993-2012:
Eq. 9-1
If these mean equity growth rates are assumed to continue in the future, then the equity of Coke is expected to exceed the equity of the S&P 500 index after years, which is calculated using Eq. 2-18:
That is, after 65 years the equity of Coke will exceed the equity of the S&P 500 index under these assumptions. This calculation ignored the fact that Coke is itself a part of the S&P 500 index which means the two cannot have different growth rates forever. If Coke grows to dominate the S&P 500 index then either the growth rate of Coke will decrease or the growth rate of the S&P 500 index will increase. However, it seems unrealistic that a beverage producer will grow to dominate the USA stock markets so the growth assumptions must be wrong.
9.4.
Decreasing ROE & Retain Ratios
In modelling the future growth of a companys equity, earnings and payouts, a number of factors can be taken into account, such as the expected change in the overall market size for the companys products, the companys expected market share, decreasing return on equity if competition is believed to increase, etc. Such adjustments rely at least partially on qualitative assessments. A simple way of reducing the future growth of Coke is used here for demonstration purposes and the reader is encouraged to experiment with other models and assumptions using the computer source-code and data-files described in section 10. The mean equity growth rate for Coke was 13.2%, see Eq. 9-1, while it was 4.46% for the S&P 500 index, see Eq. 6-1. The future equity growth rate of Coke can be made to converge to that of the S&P 500 index by exponentially decreasing the ROE and Retain ratios in the Monte Carlo simulations. The total factor of reduction must converge to so the reduction in ROE can be set to converge to e.g.
9
Equity of the S&P 500 index is estimated as the equity per share (USD 661.93) multiplied by the divisor (8940.96). Data was retrieved early March 2013 from S&P Dow Jones Indices: http://us.spindices.com/documents/additional-material/sp-500-eps-est.xls 10 Form 10-Q for Coca-Cola filed with US SEC for the quarter ending September 28, 2012: www.sec.gov/Archives/edgar/data/21344/000002134412000051/a2012092810q.htm
44
0.67 and the reduction in the Retain ratio is then made to converge to 0.5 because . The choice of how quickly ROE and Retain should be reduced as a function of the number of years is somewhat arbitrarily set to a factor 0.95 per year for ROE and 0.9 for the Retain ratio. The Retain ratio scale is then:
The ROE scale is:
The Retain scale is 0.95 in year , in year 10 it is about 0.67, in year 50 it has almost converged to 0.5. The ROE scale is 0.9835 in year 1, in year 10 it is about 0.87, in year 50 it is about 0.70, and it converges to 0.67 as increases. The last ROE scale is also used in the estimation of the terminal value in Eq. 2-43. Applying these scales to the sampled ROE and Retain ratios in the Monte Carlo simulations results in Cokes mean payout for the first year to be about 5.4% of the market-cap, which is about the same as in the nonscaled version in section 9.2 because the scale is almost 1. After 10 years the mean payout is about 12% of the initial share price, after 50 years the mean payout is about 88% of the initial share price, and after 100 years the mean payout is about 743% of the initial share price. Figure 67 shows the probability distribution of the value yield which has mean 12.6% and standard deviation 1.4%.
9.5.
Comparison to S&P 500
The value yield of Coke can be compared to the value yield of the S&P 500 stock market index plus some risk premium so as to assess the probability that Coke offers a higher risk-adjusted rate of return. This probability is denoted:
Figure 68 shows the difference where Cokes value yields result from the Monte Carlo simulations described in section 9.4 and the value yields for the S&P 500 index result from Monte Carlo simulations similar to those described in section 6.7. Coke is a part of the S&P 500 index so they are dependent in at least this way, but they have been assumed here to be independent for the sake of simplicity. Two plots are made in Figure 68 for different P/Book ratios. The left plot uses the P/Book ratios at the time of this writing in late March 2013 which were about 5.5 for Coke and about 2.3 for the S&P 500 index. The mean value yield difference is then about 4.1% with standard deviation 1.4%. The right plot uses a hypothetical case where Coke is somewhat cheaper with a P/Book ratio of 5 and the S&P 500 index is more expensive with a P/Book ratio of 3.5. The mean value yield difference is then 6.2% with standard deviation 1.5%. Figure 69 shows the probability of Cokes value yield being greater than that of the S&P 500 index with varying risk premiums and for different P/Book ratios. Note how the curve shifts on the x-axis when the P/Book ratios change. For example, using the P/Book ratios at the time of this writing (left plot) gives a
45
probability of almost 30% of Cokes value yield being greater than that of the S&P 500 index plus a risk premium of 5%. In the hypothetical case (right plot) where Coke is cheaper and S&P 500 is more expensive, the probability is more than 60% of Cokes value yield being greater than that of the S&P 500 index plus a risk premium of 5%.
9.6.
Present Value
The present value of Coke can be calculated by Monte Carlo simulating Cokes future payouts for years, estimating a terminal value, and discounting all these with the compounded rate of return that can be obtained from alternative investments. Let the discount-factor for year be denoted . The present value calculation from Eq. 2-6 then becomes:
Eq. 9-2
Cokes payouts are simulated as described in section 9.4. The number of simulation years is 320 as calculated from the formulas in section 2.8.1 so as to achieve sufficient numerical precision. 9.6.1. Present Value Relative to S&P 500 Cokes future payouts can be valued relative to the return on an investment in the S&P 500 index by using the value yield for the S&P 500 index as the discount rate in the present value calculation. This is done by sampling a discount rate from the log-normal distribution for S&P 500 value yields described in section 6.7.1 and setting the discount factor in Eq. 9-2 to . This present value calculation is repeated several times to estimate the probability distribution. Figure 70 shows the resulting probability distributions for Cokes present value with different choices of P/Book ratios for the S&P 500 index and risk premiums demanded in excess of the rate of return on the S&P 500 index. The CDFs for these probability distributions are shown in Figure 71. Note that these plots have been normalized by dividing the present values with Cokes current equity (denoted PV/Equity or equivalently PV/Book) and should therefore be compared to the companys P/Book ratio which is currently about 5.5. For example, if the risk premium demanded in excess of the return from the S&P 500 index is 2% (percentage points) and the P/Book ratio for the S&P 500 index is 2.3 as it currently is, then Figure 71 shows that there is a probability of about 0.07 (or 7%) that the present value of Coke is less than its current market-cap to someone who can own Coke shares for eternity. The same plot shows that there is a probability of about 0.05 (or 5%) that the present value of Coke is greater than 11 which is twice the current market price of Coke (11 equals two times its current P/Book ratio of 5.5). This means there is a probability of about 0.88 (or 88%) that the present value of Coke lies between its current market price and twice that price. Using present value calculations such as these it is also possible to make a plot showing the probability of the present value being greater than the market-cap for varying risk premiums. Such a plot would be
46
identical to those in Figure 69 that are made with value yields which by definition is equivalent to comparing the present value to the market-cap. 9.6.2. Present Value Relative to DJVC Cokes future payouts can also be valued relative to an investment in the DJVC index by using the compounded returns resulting from the Monte Carlo simulations described in section 8.4 as the discount factor in Eq. 9-2. Although it is also possible to add a risk premium, it is not done here. Figure 72 shows the resulting probability distribution for the present value of Coke. The distribution has a single present value above 800 but the plots have been limited to only show the present values below 50 to improve legibility. As above, the present values are normalized by dividing with Cokes equity and should therefore be compared to the P/Book ratio which is currently 5.5. The CDF plot shows that the probability is more than 0.6 (or 60%) that the present value of Coke is less than its current market-cap, and there is a probability of about 0.15 (or 15%) that the present value is greater than twice its current market-cap. Now consider the annualized rate of return which is the exponential growth rate that causes the starting index level to become in years and is calculated as:
Eq. 9-3
The annualized rate of return gives a smooth exponential increase from compounded return gives a volatile increase, see Figure 74.
to
while the
Figure 73 shows the probability distribution for the present value of Coke when the discount rate is the annualized rate of return of the DJVC index. The CDF plot shows that the probability is now about 0.9 (or 90%) that the present value of Coke is less than its current market-cap, and there is about zero probability that the present value is greater than twice its current market-cap. Figure 75 shows a large difference in the probability distributions of Cokes present value when the discount rate is either the compounded or annualized returns of the DJVC index. The reason is that present value calculations are affected by changes in the discount rate between years so that using the volatile compounded return as the discount rate gives a different present value than using the constant annualized return. Ideally, the compounded return should be used as the discount rate because it more accurately accounts for the timing in the present value calculation, but this also requires dependencies between Cokes payouts and the DJVC index to be modelled which is not done here. The same is true when using the value yield of the S&P 500 as the discount rate. 9.6.3. Comparison of S&P 500 and DJVC Discount Rates The discount rates that were sampled from the S&P 500 value yield distribution ranged between 6.9% and 10.6% with arithmetic mean 8.5% and standard deviation 0.5%, which resulted in the present values of Coke (that is, PV/Equity) ranging between 6.1 and 28 with arithmetic mean 12.6 and standard deviation 3.5. No risk premium was used and the P/Book ratio of the S&P 500 index was 2.3. If a risk premium of 2% (percentage points) is used then the mean discount rate is 10.5% and the present value ranges between 3.8 and 18.4 with arithmetic mean 7.9 and standard deviation 1.9. 47
Using the simulated compounded returns of the DJVC index as the discount rate, resulted in the present values ranging between 0.6 and 809 with arithmetic mean 9.5 and standard deviation 36.9. Over the entire simulation period of 320 years the annualized rate of return of the DJVC index ranged between 6.2% and 21.3% with arithmetic mean 16% and standard deviation 2%. Using this annualized rate of return as the discount rate resulted in a present value between 1.7 and 22.9 with arithmetic mean 3.9 and standard deviation 1.3. The DJVC returns were highly volatile as shown in Figure 76, which means that some simulations result in the returns of DJVC being much lower than the returns of S&P 500 for many years. In such cases the present values calculated using the compounded returns of DJVC as the discount rate are much greater than if the S&P 500 returns had been used as the discount rate, because the present value is dominated by the payouts that are nearest in time. This also explains why the mean present value of Coke is higher when calculated with the DJVC returns even though the mean compounded return is higher for the DJVC index than for the S&P 500 index. This shows a severe weakness in using only the arithmetic mean of present values when evaluating investments and the probability distribution should be used instead. For example, although the mean present value of Coke is 7.9 when calculated with the S&P 500 returns as discount rate and a risk premium of 2%, Figure 71 shows that there is a probability of about 0.93 (or 93%) that the present value is greater than the current market-cap. Conversely, the mean present value is 9.5 when calculated relative to the DJVC returns, but Figure 72 shows that there is a probability of less than 0.4 (or 40%) that the present value is greater than the current market-cap. 9.6.4. Jensens Inequality Another reason for using the probability distribution in present value calculations is due to Jensens inequality, which states that if the discount rate is a stochastic variable then the present value may be underestimated by using the mean discount rate in the present value calculation, see section 3.2. Figure 77 shows the difference in present value when calculated using as discount rate samples of the S&P 500 value yield, the DJVC compounded returns, and their means. This shows the present value is significantly affected from using mean discount rates. The S&P 500 value yield from section 9.6.1 is used here and its mean is about 10.5% including a 2% risk premium. Using this as the discount rate results in Cokes mean present value being about 7.75. Instead using samples of the S&P 500 value yield as the discount rate, gives a mean present value of about 7.81. This confirms Jensens inequality in Eq. 3-8 that using a mean discount rate underestimates the mean present value. The small difference in present value is due to the comparatively low standard deviation of about 0.5% for the S&P 500 value yield. The mean annualized rate of return of the DJVC index from section 9.6.2 is about 15.9%. Using this as the discount rate results in Cokes mean present value being about 3.70. Instead using samples of the DJVC annualized rate of return as the discount rate gives a mean present value of about 8.35. This again confirms Jensens inequality that using a mean discount rate may underestimate the mean present value. The large difference here is due to the comparatively large standard deviation of about 2% for the annualized rate of return of the DJVC index.
48
10.
Source Code & Data
The experiments in this paper have been implemented in the statistical programming language R which is freely available from the internet. The source-code and data-files are available at: www.hvass-labs.org/people/magnus/publications/pedersen2013monte-carlo.zip
10.1. Time Usage

Executing this implementation on a consumer-level computer from the year 2011 typically requires only seconds or minutes for an experiment consisting of a thousand Monte Carlo simulations. When a larger number of simulations is required, e.g. if the simulations in section 5.2 are repeated for numerous P/Book ratios, then the time usage may extend to an hour. Implementing a parallelized version or using another programming language might significantly decrease the time usage.
11.
Conclusion
This paper used a simple equity growth model combined with historical financial data in Monte Carlo simulation of the future equity, earnings and payouts of the S&P 500 stock market index and the Coca-Cola company. This estimated the probability distribution of the equity, earnings and payouts, which was in turn used to estimate the annualized rate of return with different holding periods, as well as stock option values. The paper also studied the returns from the S&P 500 index in relation to USA government bonds and the Dow Jones Venture Capital (DJVC) index, where it was found that there is no consistent and predictable risk premium between these asset classes. Several assumptions were made in this study and the results should be interpreted with some caution. The computer source-code and data files are provided on the internet so the reader can experiment with other assumptions.
49
12.
Appendix
12.1. Geometric Series

The following is a well-known property of geometric series:
Starting the summation at
gives:
So for a growth rate
and discount rate
with
we have:
Eq. 12-1
Eq. 12-2
Eq. 12-3
Instead starting the summation at
gives:
Eq. 12-4
Summing from
to
gives:
Eq. 12-5
50
12.2. Properties of Log-Normal Distribution

The mean of a log-normal distributed stochastic variable is known to be:
Eq. 12-6
The standard deviation is known to be:
Eq. 12-7
12.3. S&P 500 Data

Data for the S&P 500 stock market index has been collected by the staff at the Customized Research Department of Compustat. 11 The data was then used to calculate the various financial ratios. This appendix details the compromises made in collecting and making calculations on the data. The S&P 500 stock market index consists of 500 companies that are weighted according to certain changes and events that affect their capitalization. The weights are proprietary and could not be obtained. Instead, the data for the individual companies in the S&P 500 index has merely been aggregated (summed). In the period 1983-2011 the S&P 500 index consisted mostly of companies reporting their financial statements in USD currency. Of the 500 companies in the index an average of 497 companies used USD currency each year. To avoid distorting effects of non-USD currencies those would either have to be converted into USD or removed from the data. For simplicity and because so few companies reported in non-USD currency they were removed from the data. The S&P 500 index is being studied here as if it was one big conglomerate. This means the financial data such as assets, equity, earnings, and market-cap are aggregated. This differs from the accounting used in actual conglomerates where the consolidated financial statements would adjust for inter-company, intraconglomerate dependencies such as revenue and liabilities. Making such consolidated financial statements is a complex process requiring access to financial details of the companies in question which is only available to those companies and their auditors. The sum of financial data is deemed sufficiently accurate for this study. The Compustat database contains an item named MKVALT, which is the market price for the common stock of a company, or market-cap as it is referred to here. To find the market-cap for the S&P 500 conglomerate, the MKVALT items should be summed for all companies in the S&P 500 index. However, prior to 1998 the MKVALT item does not exist in the Compustat database so the research staff at Compustat had to make a formula for estimating this by taking several factors into account, such as multiple share classes. The formula gave estimates of MKVALT that were reasonably close to the existing values in the Compustat database for the period 1998-2011 so the MKVALT estimates for the period 1983-2011 are used with sufficient confidence of their accuracy, which is deemed to be within the precision required for this study.
11
The Compustat database was accessed through the facilities of the Collaborative Research Center 649 on Economic Risk at the Humboldt University of Berlin, Germany.
51
The equity per share of the S&P 500 index is calculated for the last day of each year from the known price per share, the aggregated equity and MKVALT estimates for the entire index (book-value is another name for the equity):
This calculation may have some error because the S&P 500 share price is for weighted constituent companies but MKVALT is un-weighted. The error is expected to be negligible. The equity per share is estimated for each day of a year by linearly interpolating between the known equity values at year end. This is also used to estimate the P/Book ratios for all days of a year by dividing a known price per share for each day with the estimated equity per share for that day. ROE is calculated using the reported net income available to common shareholders for a given year, divided by the equity at the end of the prior year. The constituent companies of the S&P 500 index change each year so the equity for the prior year may be for different companies than the net income of the current year. In the period 1983-2011, covering the data of this study, the number of changes to the constituent companies of the S&P 500 index was 24 companies per year on average with standard deviation 11. That is, on average less than 5% of the constituent companies of the S&P 500 index were changed each year, whose impact on the financial data is likely negligible as large companies typically remain in the S&P 500 index, and large companies dominate the aggregated financial data. So this way of calculating ROE for the S&P 500 conglomerate is considered satisfactory in terms of numerical precision. Similarly for ROA. Compustat provides the data items PRSTKC and SSTK for the amount of share buyback and issuance, respectively, where the preferred and common stocks are combined. However, the other data items being considered in this study are for the common stock alone, which means comparisons and calculations made using these numbers for share buyback and issuance contain an error as the preferred stock is included. But the error is negligible because the preferred equity is only 4% of the common equity on average for the period 1983-2011. Cash flows associated with tax benefits of stock options have been ignored as they are negligible.
52
Daily Returns 1984-2011 (S&P 500)
Daily Returns 1984-2011 (S&P 500)

50
0.10
Density
30
Sample Quantiles
20
0.80
0.85
0.90
0.95 Daily Return
1.00
1.05
1.10
-0.20
-4
-0.15
10
-0.10
-0.05
0.00
0.05
40
-2
0 Theoretical Quantiles
Figure 1: Daily price returns for the S&P 500 stock market index during the period 1984-2011, calculated as . The histogram (left) has a fitted log-normal PDF and the Q-Q plot (right) also compares to the log-normal PDF. This clearly shows that the distribution of daily price returns is not log-normal. The minimum price return is about 0.795 (corresponding to a change of -20.5%), the maximum is about 1.116 (a change of +11.6%), the mean is 1 (a change of 0%) and the standard deviation is 0.0118 (or 1.18%).
53
Year 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011
ROA 4.0% 3.0% 2.8% 3.2% 4.1% 3.2% 2.6% 1.7% 0.8% 1.8% 2.9% 3.1% 3.3% 3.2% 3.1% 3.3% 3.0% 1.2% 0.5% 2.4% 2.6% 2.9% 3.2% 2.4% 0.9% 1.9% 3.0% 3.2%
ROE Retain Equity/Assets 13.7% 49.3% 28% 10.9% 29.2% 26% 10.9% 30.5% 25% 12.6% 27.4% 25% 16.7% 33.2% 21% 14.8% 41.6% 20% 12.6% 28.5% 20% 8.3% 28.9% 20% 4.0% (49.3)% 18% 9.7% 29.3% 17% 16.8% 50.6% 18% 17.2% 35.9% 18% 18.7% 47.9% 18% 17.6% 35.9% 18% 17.9% 32.1% 17% 20.0% 45.7% 17% 17.5% 50.4% 18% 6.7% 5.9% 18% 2.9% (140.8)% 16% 15.5% 49.8% 16% 15.7% 41.1% 17% 17.1% 26.6% 17% 18.7% 21.4% 17% 13.9% (16.8)% 16% 5.3% 42.0% 16% 12.0% 47.5% 19% 15.8% 44.0% 21% 15.7% 28.3% 21%
Table 1: Return on Assets (ROA), Return on Equity (ROE), earnings retaining ratio (Retain) and equity-to-assets for the aggregated companies in the S&P 500 stock market index, calculated using Eq. 2-9, Eq. 2-7 and Eq. 2-4. Statistics are shown in Table 2. This data is available for download from the internet as described in section 6.9.
Mean Stdev Min Max ROA 2.6% 0.9% 0.5% 4.1% ROE 13.5% 4.7% 2.9% 20.0% Retain 26.4% 39.0% (140.8)% 69.4% Equity/Assets 19.6% 3.6% 15.6% 28.8%
Table 2: Statistics for Table 1.
54

Retain vs. ROE (S&P 500, 1984-2011)
0.20 ROEt 0.05 0.10 0.15
-1.0
-0.5 Retaint
0.0
0.5
Retain vs. ROE (S&P 500, 1984-2011)

0.20
Retain vs. ROE (S&P 500, 1984-2011)
0.15
ROEt
0.10
ROEt 0.05
-1.0 -0.5 Retain t
1
0.0
0.5
0.05
0.10
0.15
0.20
-1.0
-0.5 Retain t
1
0.0
0.5
Figure 2: Relations between and for the S&P 500 stock market index using the data from Table 1. The upper plot is for , the lower left plot is for and the lower right plot is for , that is, the lower plots have shifted the series one year backwards and forwards, respectively. All three plots show a cluster of around 25-50% and around 10-20%, but the outliers behave differently in these plots. For the two lower plots where the series is time-shifted, the outliers do not show a clear tendency e.g. of Retain and ROE increasing together, as is evident in the upper plot for which is therefore used when sampling in Monte Carlo simulations.
55
ROE (S&P 500)

1.0 1.0
Retain (S&P 500)
0.8
0.6
Autocorrelation
Autocorrelation
0.4
0.2
0.0
-0.2
-0.4
6 Lag
10
12
14
-0.4
-0.2
0.0
0.2
0.4
0.6
0.8
6 Lag
10
12
14
Figure 3: Autocorrelation of and with varying time lag for the S&P 500 stock market index using the data from Table 1. The left plot shows the autocorrelation of with the dotted line being the level of statistical significance, which shows that is obviously perfectly correlated with itself (that is, when ), and there appears to be a positive correlation between and which is explored further in Figure 4, and an almost significant negative correlation between and and . The right plot shows the autocorrelation of the series whose largest autocorrelation occurs with a 10 year lag.
ROE (S&P 500)

0.20 0.5
Retain (S&P 500)
0.15
Retain t
0.05 0.10 ROE 0.15 0.20
ROE t
0.10
0.05
-1.0
-0.5
0.0
-1.0
-0.5 Retaint
0.0
0.5
t Figure 4: Autocorrelation of and with a one-year time lag, for the S&P 500 stock market index using the data from Table 1. The left plot shows the relation between and whose linear fit has a low coefficient of determination and hence does not appear to be significantly linearly correlated. The right plot shows the relation between and which has a cluster around 25-50% but appears to be uncorrelated otherwise.
56
Year 2
14
Year 3
12
Year 4
15
12
10
10
Density
Density
Density
0.05 0.10 Earnings 0.15 0.20
0.05
0.10 Earnings
0.15
0.20
10
0.05
0.10
0.15 Earnings
0.20
0.25
Year 5
12
Year 6
Year 7
10
10
Density
Density
Density
0.05 0.10 0.15 Earnings 0.20 0.25
0.05
0.10
0.15 Earnings
0.20
0.25
0.05
0.10
0.15 Earnings
0.20
0.25
0.30
Year 8
Year 9
Year 10
10
Density
Density
Density
0.05 0.10 0.15 0.20 Earnings 0.25 0.30
0.05
0.10
0.15
0.20 Earnings
0.25
0.30
0.35
0.05
0.10
0.15
0.20 Earnings
0.25
0.30
0.35
Figure 5: Histograms of the earnings of the S&P 500 stock market index at different years, resulting from the Monte Carlo simulations described in section 6.4.
57
Year 20
Year 55
0.14
Year 90
0.8
0.6
Density
Density
0.4
Density
0.5 1.0 1.5 2.0 Earnings 2.5 3.0 3.5
0.2
0.1
0.2
0.3 Earnings
0.4
0.5
0.6
0.00
0
0.0
0.02
0.04
0.06
0.08
0.10
0.12
10 Earnings
15
20
Year 125
Year 160
0.0014
Year 195
0.030
0.006
0.025
0.005
0.020
0.004
Density
Density
0.015
0.003
Density
0 100 200 300 Earnings 400 500
0.010
0.002
0.005
0.001
20
40
60 Earnings
80
100
0.0000
0
0.000
0.000
0.0002
0.0004
0.0006
0.0008
0.0010
0.0012
500
1000
1500 Earnings
2000
2500
Year 230
0.00025
Year 265
1.2e-05
Year 300
5e-05
0.00020
4e-05
0.00015
Density
Density
0.00010
Density
0 20000 40000 Earnings 60000 80000
0.00005
5000
10000 Earnings
15000
20000
0.0e+00
0e+00
0.00000
0e+00
2.0e-06
1e-05
4.0e-06
2e-05
6.0e-06
3e-05
8.0e-06
1.0e-05
1e+05
2e+05
3e+05 Earnings
4e+05
5e+05
Figure 6: Histograms of the earnings of the S&P 500 stock market index in different years, resulting from the Monte Carlo simulations described in section 6.4.
58
Earnings Mean (S&P 500)
Earnings Stdev (S&P 500)
1e+03
Earnings Stdev
Earnings Mean
1e+01
1e-01
50
100
150 Year
200
250
300
1e-01
1e+01
1e+03
50
100
150 Year
200
250
300
Figure 7: Earnings mean and standard deviation for the S&P 500 stock market index, resulting from a series of Monte Carlo simulations as described in section 6.4. The small circles are the results of the Monte Carlo simulations and the fitted lines are given in Eq. 6-2 and Eq. 6-3.
59
Year 2
Year 3
Year 4
-2.0
-2.0
Sample Quantiles
Sample Quantiles
Sample Quantiles
-2.5
-2.5
-3.0
-3.0
-3.5
-3.5
-3
-2
-1
-3
-2
-1
-3.5
-3.0
-2.5
-2.0
-3
-2
-1
Theoretical Quantiles
Year 5
Year 6
Year 7
-1.5
-1.5
-2.0
-2.0
Sample Quantiles
Sample Quantiles
Sample Quantiles
-2.5
-2.5
-3.0
-3.0
-3.5
-3
-2
-1
-3.5
-3
-2
-1
-3.0
-2.5
-2.0
-1.5
-3
-2
-1
Year 8
Year 9
Year 10
-1.5
-1.5
-2.0
-2.0
Sample Quantiles
Sample Quantiles
Sample Quantiles
-2.5
-2.5
-3.0
-3.0
-3
-2
-1
-3
-2
-1
-3.0
-2.5
-2.0
-1.5
-3
-2
-1
Figure 8: Q-Q log-normal plots of the payouts of the S&P 500 stock market index at different years, resulting from the Monte Carlo simulations described in section 6.5. The x-axes show theoretical quantiles and the y-axes show sample quantiles.
60
Year 20
Year 55
Year 90
-1.0
1.0
0.5
Sample Quantiles
-1.5
Sample Quantiles
Sample Quantiles
-3 -2 -1 0 1 2 3
-2.0
0.0
-0.5
-2.5
-1.0
-3.0
-3
-2
-1
Year 125
-1.5
0.0
0.5
1.0
1.5
2.0
2.5
Year 160
-3
-2
-1
Year 195
6.0
4.0
5.5
3.5
5.0
Sample Quantiles
Sample Quantiles
Sample Quantiles
-3 -2 -1 0 1 2 3
3.0
4.5
4.0
2.5
3.5
2.0
3.0
1.5
2.5
-3
-2
-1
Year 230
4.0
4.5
5.0
5.5
6.0
6.5
7.0
7.5
Year 265
-3
-2
-1
Year 300
11
10
Sample Quantiles
Sample Quantiles
Sample Quantiles
-3 -2 -1 0 1 2 3
-3
-2
-1
10
11
12
-3
-2
-1
Figure 9: Q-Q log-normal plots of the payouts of the S&P 500 stock market index in different years, resulting from the Monte Carlo simulations described in section 6.5. The x-axes show theoretical quantiles and the y-axes show sample quantiles.
61
Payout Mean (S&P 500)
Payout Stdev (S&P 500)
1e+03
Payout Stdev
Payout Mean
1e+01
1e-01
50
100
150 Year
200
250
300
1e-01
1e+01
1e+03
50
100
150 Year
200
250
300
Figure 10: Payout mean and standard deviation for the S&P 500 stock market index, resulting from a series of Monte Carlo simulations as described in section 6.5. The small circles are the results of the Monte Carlo simulations and the fitted lines are given in Eq. 6-18 and Eq. 6-19.
62
Year 2
8
Year 5
3.5
Year 10
Density
Density
Density
0.0
0.5
1.0
1.5
2.0
10
2.5
3.0
15
0.10
0.15
0.20 PayoutSum
0.25
0.30
0.3
0.4
0.5 PayoutSum
0.6
0.7
0.8
1.0
1.2 PayoutSum
1.4
1.6
Year 15
1.5
Year 20
0.7
Year 30
2.5
2.0
1.0
Density
1.5
Density
Density 0.5 0.0
1.0
0.5
0.0
1.4
1.6
1.8 PayoutSum
2.0
2.2
2.4
2.0
2.5
3.0 PayoutSum
3.5
0.0
0.1
0.2
0.3
0.4
0.5
0.6
6 PayoutSum
Year 50
Year 75
Year 100
0.04
0.15
0.10
0.03
Density
Density
0.02
Density
0.05
0.01
0.000
0.00
0.00
0.002
0.004
0.006
0.008
0.010
10
15
20 PayoutSum
25
20
40
60 PayoutSum
80
100
50
100
150
200 PayoutSum
250
300
350
Figure 11: Histograms with fitted log-normal PDFs of the cumulative payout sums for the S&P 500 stock market index in different years, resulting from the Monte Carlo simulations described in section 6.5.2.
63
PayoutSum
(S&P 500)
PayoutSum
(S&P 500)
10
PayoutSum
PayoutSum
0.2
0.3
0.4
0.5
50
100
150 Year
200
250
300
50
100
150 Year
200
250
300
Figure 12: Log-normal probability distribution parameters and for the payout sum of the S&P 500 stock market index, resulting from a series of Monte Carlo runs as described in section 6.5.2. The small circles are the results of the Monte Carlo simulations and the fitted curves are given in Eq. 6-15 and Eq. 6-16.
64
Year 2
Year 3
Year 4
0.08
0.15
0.06
Sample Quantiles
0.04
Sample Quantiles
0.10
Sample Quantiles
0.00
0.02
0.00
0.05
-0.02
-0.04
-0.05
-3
-2
-1
-3
-2
-1
-0.1
0.0
0.1
0.2
-3
-2
-1
Theoretical Quantiles Year 5

0.4
0.3
0.2
0.3
Sample Quantiles
Sample Quantiles
0.2
Sample Quantiles
0.1
0.1
0.0
0.0
-3
-2
-1
-0.1
-3
-2
-1
0.0
0.1
0.2
0.3
0.4
-3
-2
-1

0.6
Theoretical Year Quantiles 10
0.5
0.4
0.5
Sample Quantiles
Sample Quantiles
0.4
Sample Quantiles
0.3
0.3
0.2
0.2
0.1
0.1
0.0
-3
-2
-1
0.0
-3
-2
-1
0.1
0.2
0.3
0.4
0.5
0.6
-3
-2
-1
Figure 13: Q-Q log-normal plots of the equity of the S&P 500 stock market index in different years, resulting from the Monte Carlo simulations described in section 6.6. The x-axes show theoretical quantiles and the y-axes show sample quantiles.
65
Year 20
1.2
Year 55
Year 90
3.0
Sample Quantiles
1.0
Sample Quantiles
Sample Quantiles
0.8
2.5
0.6
0.4
2.0
-3
-2
-1
-3
-2
-1
3.0
3.5
4.0
4.5
-3
-2
-1

6.5 8.0
6.0
7.5
Sample Quantiles
Sample Quantiles
Sample Quantiles
7.0
5.5
5.0
6.5
4.5
6.0
-3
-2
-1
5.5
-3
-2
-1
7.0
7.5
8.0
8.5
9.0
9.5
-3
-2
-1

15
11
13
Sample Quantiles
Sample Quantiles
12
Sample Quantiles
10
11
10
-3
-2
-1
-3
-2
-1
11
12
13
14
-3
-2
-1
Figure 14: Q-Q log-normal plots of the equity of the S&P 500 stock market index in different years, resulting from the Monte Carlo simulations described in section 6.6. The x-axes show theoretical quantiles and the y-axes show sample quantiles.
66
Equity (S&P 500)
Equity (S&P 500)
12
10
Equity
Equity
50
100
150 Year
200
250
300
0.0
0.1
0.2
0.3
0.4
0.5
50
100
150 Year
200
250
300
Figure 15: Log-normal probability distribution parameters and for the equity of the S&P 500 stock market index, resulting from a series of Monte Carlo runs as described in section 6.6.1. The small circles are the results of the Monte Carlo simulations and the fitted curves are given in Eq. 6-18 and Eq. 6-19.
Equity at Year 8 (S&P 500)
Density
0.0
0.5
1.0
1.5
2.0
2.5
3.0
1.0
1.2 Equity
1.4
1.6
Figure 16: Distribution of equity for the S&P 500 stock market index in year 8. Also shown is the fitted log-normal PDF. The Monte Carlo simulation resulting in this plot is described in section 6.6. The Q-Q plot is shown in Figure 13.
67
P/Book=1 (S&P 500)
P/Book=1.625 (S&P 500)
P/Book=2.25 (S&P 500)
50
60
40
50
Density
30
Density
40
30
Density
20
20
10
10
0.10
0.11
0.12
0.13 Value Yield
0.14
0.15
0.16
0.08
0.09
0.10 Value Yield
0.11
0.12
20
40
60
80
0.070
0.075
0.080
0.085 Value Yield
0.090
0.095
0.100
P/Book=2.875 (S&P 500)

120
P/Book=3.5 (S&P 500)

120
P/Book=4.125 (S&P 500)
100
80
Density
60
Density
80
Density
40
60
20
40
20
0.065
0.070
0.075
0.080 Value Yield
0.085
0.090
0.060
0.065
0.070 Value Yield
0.075
0.080
20
40
60
80
100
0.055
0.060
0.065
0.070
0.075
0.080
Value Yield
P/Book=4.75 (S&P 500)
P/Book=5.375 (S&P 500)
P/Book=6 (S&P 500)
100
120
150
Density
80
Density
100
60
Density 0 50
20
40
0.055
0.060
0.065 Value Yield
0.070
0.075
0.040
0.045
0.050
0.055 Value Yield
0.060
0.065
0.070
50
100
150
0.035
0.040
0.045
0.050 Value Yield
0.055
0.060
0.065
Figure 17: Value yield distributions for the S&P 500 stock market index with different P/Book ratios. Also shown are the fitted log-normal PDFs. The Monte Carlo simulations resulting in these plots are described in section 6.7. The Q-Q plots are shown in Figure 18 and Figure 19.
68
P/Book=1 (S&P 500)

-1.8
P/Book=1.625 (S&P 500)

-2.3
P/Book=2.25 (S&P 500)
-1.9
-2.2
Sample Quantiles
-2.0
Sample Quantiles
Sample Quantiles
-3 -2 -1 0 1 2 3
-2.3
-2.1
-2.2
-2.4
-2.3
-3
-2
-1
-2.5
-2.6
-2.5
-2.4
-3
-2
-1
Theoretical Quantiles P/Book=2.875 (S&P 500)
-2.4
-2.55
-2.50
Sample Quantiles
-2.5
Sample Quantiles
-2.60
Sample Quantiles
-3 -2 -1 0 1 2 3
-2.65
-2.6
-2.70
-2.7
-2.80
-2.75
-3
-2
-1
-2.9
-2.8
-2.7
-2.6
-3
-2
-1
Theoretical(S&P Quantiles P/Book=6 500)
-2.60
-2.65
-2.7
-2.70
-2.8
Sample Quantiles
Sample Quantiles
Sample Quantiles
-3 -2 -1 0 1 2 3
-2.75
-2.9
-2.80
-3.0
-2.85
-3.1
-2.90
-3.2
-3.3
-2.95
-3
-2
-1
-3.4
-3.3
-3.2
-3.1
-3.0
-2.9
-2.8
-2.7
-3
-2
-1
Figure 18: Q-Q log-normal plots for the value yield of the S&P 500 stock market index with different P/Book ratios. The Monte Carlo simulations resulting in these plots are described in section 6.7. The x-axes are theoretical quantiles and the y-axes are sample quantiles.
69
P/Book=1 (S&P 500)

0.12
P/Book=1.625 (S&P 500)

0.100
P/Book=2.25 (S&P 500)
0.16
0.15
0.11
Sample Quantiles
0.14
Sample Quantiles
Sample Quantiles
-3 -2 -1 0 1 2 3
0.13
0.12
0.09
0.11
0.10
-3
-2
-1
0.08
0.070
0.075
0.080
0.085
0.10
0.090
0.095
-3
-2
-1
0.090
0.080
0.085
0.075
Sample Quantiles
Sample Quantiles
Sample Quantiles
-3 -2 -1 0 1 2 3
0.080
0.070
0.075
0.065
0.070
0.065
0.060
-3
-2
-1
0.055
0.060
0.065
0.070
0.075
0.080
-3
-2
-1
Theoretical(S&P Quantiles P/Book=6 500)
0.075
0.070
0.070
0.065
0.060
Sample Quantiles
Sample Quantiles
Sample Quantiles
-3 -2 -1 0 1 2 3
0.065
0.055
0.060
0.050
0.045
0.055
0.040
-3
-2
-1
0.035
0.040
0.045
0.050
0.055
0.060
0.065
-3
-2
-1
Figure 19: Q-Q normal plots for the value yield of the S&P 500 stock market index with different P/Book ratios. The Monte Carlo simulations resulting in these plots are described in section 6.7. The x-axes are theoretical quantiles and the y-axes are sample quantiles.
70
Value Yield
(S&P 500)
Value Yield
(S&P 500)
-1.0
Value Yield
-1.5
Value Yield
0 2 4 P/Book 6 8 10
-2.5
-2.0
-3.0
0.06
0
0.08
0.10
0.12
4 P/Book
10
Figure 20: Log-normal probability distribution parameters and for the value yield of the S&P 500 stock market index, resulting from a series of Monte Carlo runs with different P/Book ratios as described in section 6.7. The small circles are the results of the Monte Carlo simulations and the fitted curves are given in Eq. 6-26 and Eq. 6-27.
Value Yield Mean (S&P 500)

0.5
Value Yield Stdev (S&P 500)
0.4
Value Yield Mean
Value Yield Stdev

0 2 4 P/Book 6 8 10
0.3
0.2
0.1
0.01
0
0.02
0.03
0.04
0.05
0.06
4 P/Book
10
Figure 21: Mean and standard deviation for the value yield of the S&P 500 stock market index, resulting from a series of Monte Carlo runs with different P/Book ratios as described in section 6.7. The small circles are the results of the Monte Carlo simulations and the fitted curves are given in Eq. 6-28 and Eq. 6-29.
71
Value Yield
Mean Stdev
Year Min Max Min Max
1983 10.9% 12.5% 0.7% 0.8% Year Min Max Min Max Year Min Max Min Max
1984 11.3% 12.3% 0.7% 0.8% 1994 8.0% 8.4% 0.5% 0.5%
1985 10.2% 11.9% 0.6% 0.8% 1995 7.5% 8.6% 0.4% 0.5%
1986 9.5% 10.7% 0.6% 0.7% 1996 7.2% 7.9% 0.4% 0.4%
1987 8.4% 10.4% 0.5% 0.6% 1997 6.7% 7.5% 0.4% 0.4%
1988 9.6% 10.5% 0.6% 0.7% 1998 6.3% 6.9% 0.3% 0.4%
1989 8.8% 10.2% 0.5% 0.6% 1999 6.1% 6.5% 0.3% 0.3%
1990 9.0% 10.1% 0.5% 0.6%
1991 8.5% 9.9% 0.5% 0.6%
1992 8.2% 8.7% 0.5% 0.5%
Value Yield
Mean Stdev
1993 7.9% 8.3% 0.4% 0.5%
2000 6.2% 6.6% 0.3% 0.4%
2001 6.6% 7.5% 0.4% 0.4%
2002 6.9% 8.2% 0.4% 0.5%
Value Yield
Mean Stdev
2003 7.2% 8.3% 0.4% 0.5%
2004 7.3% 7.7% 0.4% 0.4%
2005 7.5% 7.9% 0.4% 0.4%
2006 7.4% 7.9% 0.4% 0.4%
2007 7.4% 7.8% 0.4% 0.4%
2008 7.5% 10.3% 0.4% 0.6%
2009 8.2% 10.8% 0.5% 0.7%
2010 8.3% 9.2% 0.5% 0.5%
2011 8.4% 9.3% 0.5% 0.6%
Table 3: Value yield mean and standard deviation for the S&P 500 stock market index during the period 1984-2011, calculated using Eq. 6-28 and Eq. 6-29 with each years minimum and maximum P/Book ratios from Figure 22.
72
P/Book 1984-2011 (S&P 500)

600
Book-Value Per Share 1984-2011 (S&P 500)
Book-Value Per Share
P/Book
200
300
400
500
1985
1990
1995 Year
2000
2005
2010
1985
1990
1995 Year
2000
2005
2010
Figure 22: P/Book ratio (left) and book-value per share (right) for the S&P 500 stock market index during the period 1984-2011. The price observations are available daily while the book-value is interpolated from year-end observations which are depicted as circles in the plot to the right. As can be seen the book-value is increasing comparatively steadily and the volatility of the P/Book ratio is caused by the volatile share-price.
P/Book 1984-2011 (S&P 500)
1.0
P/Book CDF 1984-2011 (S&P 500)

1.0
0.8
Density
0.6
0.4
F(x)
0.2
0.0
3 P/Book
0.0
0.2
0.4
0.6
0.8
3 P/Book
Figure 23: Histogram with fitted log-normal PDF (left) and CDF (right) for the daily P/Book ratios of the S&P 500 stock market index as shown in Figure 22.
Arithmetic Geometric Harmonic Stdev Min Max Mean Mean Mean 2.56 2.45 2.32 0.88 1.16 5.07
Table 4: Statistics for the P/Book ratios of the S&P 500 stock market index during the period 1984-2011 as shown in Figure 22.
73
Year 2
0.5 0.7
Year 4
0.4
Year 7
0.6
0.4
0.5
Density
0.4
Density
0.3
Density
0.3
0.2
0.2
0.0
0.1
0.0
0.1
3 Price
4 Price
0.0
0.1
0.2
0.3
4 Price
Year 10
0.35
Year 15
Year 20
0.30
0.25
0.20
0.25
0.20
Density
Density
0.15
Density
0.15
0.10
0.10
0.00
0.05
0.00
0.05
4 Price
10
6 Price
10
12
0.00
0.05
0.10
0.15
0.20
10 Price
15
Year 30
0.14
Year 50
Year 100
0.12
0.08
Density
Density
0.06
Density
0.04
0.02
10
15 Price
20
25
20
40 Price
60
80
0.000
0.00
0.00
0.001
0.01
0.002
0.02
0.003
0.03
0.004
0.10
0.04
0.005
0.05
200
400
600 Price
800
1000
Figure 24: Price distribution for the S&P 500 stock market index when the price is calculated by sampling the P/Book ratio from its historical distribution shown in Figure 23 and multiplying the P/Book ratio with the equity sampled from the probability distribution resulting from the Monte Carlo simulations described in section 6.6. Also shown are the fitted log-normal PDFs which show that the log-normal distribution is approached as the number of years increase.
74
P/Book Autocorrelation 1984-2011 (S&P 500)

1.0 Autocorrelation -0.5 0.0 0.5
1000
2000
3000
4000
5000
6000
7000
Lag/Days
Figure 25: Autocorrelation of the P/Book ratios of the S&P 500 stock market index for the period 1984-2011.
P/Book Daily Change 1984-2011 (S&P 500)

50
P/Book vs. P/Book Change 1984-2011 (S&P 500)
40
30
Density
P/Book Change
20
10
0.80
0.85
0.90
0.95 Change
1.00
1.05
1.10
0.80
0.85
0.90
0.95
1.00
1.05
1.10
3 P/Book
Figure 26: Histogram with fitted log-normal PDF (left) for the daily P/Book changes of the S&P 500 stock market index in the period 1984-2011, and scatter-plot (right) for the daily P/Book and P/Book changes. Note that the histogram is similar to Figure 1 for the daily price changes because the daily change in P/Book ratio is dominated by the price change as the equity (or bookvalue) is almost constant on successive days.
75
P/Book Simulated 1 (S&P 500)

4.5
3.0
4.0
3.5
2.5
P/Book
P/Book
P/Book 2.0 1.5
3.0
2.5
2.0
1.5
2000
4000 Days
6000
2000
4000 Days
6000
1.5
2.0
2.5
3.0
3.5
2000
4000 Days
6000

3.0

5.0
3.0
2.5
2.5
P/Book
P/Book
P/Book 2.0 1.5
2.0
1.5
2000
4000 Days
6000
2000
4000 Days
6000
1.5
2.0
2.5
3.0
3.5
4.0
4.5
2000
4000 Days
6000
4.0
3.5
3.0
P/Book
3.0
P/Book
2.5
2.5
P/Book
0 2000 4000 Days 6000
1.5
2.0
1.5
2.0
2000
4000 Days
6000
1.5
0
2.0
2.5
3.0
2000
4000 Days
6000
Figure 27: Simulated P/Book for the S&P 500 stock market index 30 years into the future (about 252 trading days per year), as described in section 6.8.4. The starting P/Book is set to 2.3.
76
P/Book ACF 1 (S&P 500)

1.0 1.0

1.0
0.8
0.6
0.6
0.8
Autocorrelation
Autocorrelation
0.2
0.2
Autocorrelation
0.4
0.4
0.0
-0.2
-0.2
0.0
2000
4000 Lag/Days
6000
2000
4000 Lag/Days
6000
-0.5
0.0
0.5
2000
4000 Lag/Days
6000

1.0 1.0

1.0
0.8
0.8
Autocorrelation
0.6
Autocorrelation
0.6
Autocorrelation
0.4
0.4
0.2
0.2
0.0
0.0
-0.2
2000
4000 Lag/Days
6000
-0.2
2000
4000 Lag/Days
6000
-0.2
0.0
0.2
0.4
0.6
0.8
2000
4000 Lag/Days
6000

1.0 1.0

1.0
0.6
0.8
Autocorrelation
Autocorrelation
0.5
0.2
Autocorrelation 0.0 -0.5
0.4
-0.2
0.0
2000
4000 Lag/Days
6000
2000
4000 Lag/Days
6000
-0.5
0.0
0.5
2000
4000 Lag/Days
6000
Figure 28: Autocorrelation for the simulated P/Book ratios from Figure 27 for the S&P 500 stock market index.
77
Simulated P/Book Year 2 (S&P 500)

1.4

1.4
1.2
1.2
1.0
1.0
0.8
Density
0.8
Density
0.6
Density
0.6
0.4
0.4
0.2
0.0
0.2
0.0
1.5
2.0
2.5
3.0 P/Book
3.5
4.0
4.5
5.0
1.5
2.0
2.5
3.0 P/Book
3.5
4.0
4.5
0.0
0.2
0.4
0.6
0.8
1.0
1.2
3 P/Book

1.5 1.5
1.0
1.0
Density
Density
Density 0.5 0.0
0.5
0.0
3 P/Book
3 P/Book
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.5
2.0
2.5
3.0
3.5
4.0
4.5
5.0
P/Book

1.4

1.0
1.0
1.2
0.8
1.0
Density
0.6
Density
0.8
0.6
Density
0.4
0.4
0.2
0.0
0.0
0.2
3 P/Book
3 P/Book
0.0
0.2
0.4
0.6
0.8
3 P/Book
Figure 29: Distribution of P/Book ratios for the S&P 500 stock market index in different years that resulted from the Monte Carlo simulations described in section 6.8.4. Also shown are the fitted log-normal PDFs. The CDFs are shown in Figure 30.
78
Simulated P/Book CDF Year 2 (S&P 500)

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
3 P/Book
3 P/Book
0.0
0.2
0.4
0.6
0.8
3 P/Book

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
3 P/Book
3 P/Book
0.0
0.2
0.4
0.6
0.8
3 P/Book

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
3 P/Book
3 P/Book
0.0
0.2
0.4
0.6
0.8
3 P/Book
Figure 30: CDF of P/Book ratios for the S&P 500 stock market index in different years that resulted from the Monte Carlo simulations described in section 6.8.4. The histograms are shown in Figure 29.
79
Year 2
5.0
Year 3
5
Year 5
4.5
4.0
4.0
4.5
Simulated P/Book
Simulated P/Book
3.5
Simulated P/Book
3.5
3.0
3.0
2.5
2.0
2.0
2.5
1.5
3 Historical P/Book
1.5
3 Historical P/Book
3 Historical P/Book
Year 7
5
Year 10
5.0
Year 15
Simulated P/Book
Simulated P/Book
Simulated P/Book
3 Historical P/Book
3 Historical P/Book
1.5
2.0
2.5
3.0
3.5
4.0
4.5
3 Historical P/Book
Year 20
5.0 5
Year 25
5
Year 30
4.0
4.5
Simulated P/Book
Simulated P/Book
Simulated P/Book
3.0
3.5
2.0
2.5
3 Historical P/Book
1.5
3 Historical P/Book
3 Historical P/Book
Figure 31: Q-Q plots comparing the distributions of the simulated P/Book ratios from Figure 29 to the historical P/Book ratios from Figure 23.
80
Year 2 (S&P 500)

0.8
Year 3 (S&P 500)
Year 5 (S&P 500)
1.0
0.8
0.6
Density
0.6
Density
Density
0.4
0.4
0.2
0.2
0.0
0.0
3 Price
3 Price
0.0
0.1
0.2
0.3
0.4
0.5
0.6
4 Price
Year 7 (S&P 500)
Year 10 (S&P 500)
Year 15 (S&P 500)
0.5
0.4
0.4
0.3
Density
0.3
Density
Density
0.2
0.2
0.1
0.1
4 Price
4 Price
0.00
0.0
0.0
0.05
0.10
0.15
0.20
0.25
0.30
6 Price
10
12
Year 20 (S&P 500)

0.20
Year 25 (S&P 500)
Year 30 (S&P 500)
0.25
0.20
0.15
0.15
Density
Density
0.10
Density 0.00 0.05
0.00
0.05
0.00
0.05
0.10
0.10
0.15
8 Price
10
12
14
16
10 Price
15
20
10
15 Price
20
25
Figure 32: Price distribution for the S&P 500 stock market index when the price is calculated by Monte Carlo simulating the P/Book ratio with a starting value of 2.3 as described in section 6.8.4, and multiplying with Monte Carlo simulated equity as described in section 6.6. Also shown are the fitted log-normal PDFs which show that log-normal distributions are approached as the number of years increase.
81
Year 2 (S&P 500)

1.0 1.0
Year 3 (S&P 500)

1.0
Year 5 (S&P 500)
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
3 Price
3 Price
0.0
0.2
0.4
0.6
0.8
4 Price
Year 7 (S&P 500)

1.0 1.0
Year 10 (S&P 500)

1.0
Year 15 (S&P 500)
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
4 Price
4 Price
0.0
0.2
0.4
0.6
0.8
6 Price
10
Year 20 (S&P 500)

1.0 1.0
Year 25 (S&P 500)

1.0
Year 30 (S&P 500)
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
5 Price
10
15
10 Price
15
20
0.0
0.2
0.4
0.6
0.8
10 Price
15
20
25
Figure 33: CDF of the prices for the S&P 500 stock market index in different years where the price is calculated by Monte Carlo simulating the P/Book ratio as described in section 6.8.4 and multiplying with Monte Carlo simulated equity as described in section 6.6.
82
3.0 8 2.5
Price Stdev
30 25 5 4 3 2 5 1 10 20 15
Price Mean
2.0
1.5
4 1.0
0.5 25 5 4 3 2 5 1 10 20 15
30
Year
Year
P/Book
P/Book
Figure 34: Mean and standard deviation for the future price of the S&P 500 stock market index with varying starting P/Book ratios. The price is calculated by Monte Carlo simulating the P/Book ratio as described in section 6.8.4 and multiplying with Monte Carlo simulated equity as described in section 6.6. Subsets are shown in Table 5 and Table 6.
1 1 2 4 6 8 10 12 14 16 18 20
1.5
2.5
P/Book 3 3.5
4.5
Year
1.42 1.73 2.30 2.74 3.15 3.59 4.00 4.49 5.00 5.47 6.05
1.99 2.21 2.56 2.94 3.34 3.73 4.14 4.54 5.02 5.57 6.09
2.12 2.29 2.63 3.00 3.37 3.75 4.15 4.63 5.09 5.62 6.14
2.71 2.82 3.08 3.36 3.65 3.98 4.35 4.77 5.15 5.68 6.26
2.99 3.08 3.32 3.58 3.84 4.16 4.50 4.90 5.31 5.80 6.32
3.70 3.75 3.90 4.05 4.25 4.53 4.83 5.19 5.61 6.04 6.48
4.13 4.15 4.20 4.28 4.43 4.66 4.93 5.25 5.63 6.08 6.57
4.53 4.49 4.50 4.57 4.66 4.84 5.10 5.36 5.71 6.13 6.56
4.62 4.57 4.56 4.62 4.72 4.88 5.11 5.43 5.77 6.11 6.52
Table 5: Subset of the price means from Figure 34.
1 1 2 4 6 8 10 12 14 16 18 20
1.5
2.5
P/Book 3 3.5
4.5
Year
0.15 0.34 0.48 0.60 0.75 0.92 1.07 1.30 1.53 1.69 1.89
0.25 0.35 0.53 0.70 0.89 1.01 1.16 1.31 1.53 1.74 1.93
0.27 0.40 0.56 0.72 0.88 1.06 1.18 1.38 1.54 1.72 1.96
0.37 0.51 0.72 0.91 1.02 1.16 1.32 1.48 1.62 1.81 2.02
0.41 0.58 0.83 1.02 1.14 1.26 1.46 1.57 1.68 1.90 2.10
0.72 0.90 1.09 1.22 1.32 1.45 1.55 1.68 1.84 2.02 2.22
0.73 0.90 1.11 1.26 1.37 1.49 1.61 1.75 1.83 2.01 2.20
0.56 0.75 1.06 1.25 1.37 1.48 1.65 1.74 1.90 2.03 2.14
0.52 0.75 1.04 1.25 1.39 1.54 1.65 1.83 1.97 2.10 2.22
Table 6: Subset of the price standard deviations from Figure 34.
83
P/Book=1, Year 2
2.0
P/Book=1, Year 5
P/Book=1, Year 10
1.5
Density
Density
Density
1.0
0.5
0.0
0.2
0.4
0.6
0.8 Value Yield
1.0
1.2
1.4
0.2
0.3 Value Yield
0.4
0.5
10
12
0.15
0.20 Value Yield
0.25
0.30
P/Book=3, Year 2
3.5 8
P/Book=3, Year 5
P/Book=3, Year 10
3.0
2.5
Density
2.0
Density
1.5
Density
1.0
0.5
0.0
-0.2
0.0 Value Yield
0.2
0.4
10
12
-0.10
-0.05
0.00
0.05 Value Yield
0.10
0.15
0.20
0.00
0.05 Value Yield
0.10
0.15
P/Book=5, Year 2
P/Book=5, Year 5
15
P/Book=5, Year 10
Density
Density
Density
-0.5
-0.4
-0.3
-0.2 Value Yield
-0.1
0.0
0.1
-0.20
-0.15
-0.10
-0.05 Value Yield
0.00
0.05
10
-0.05
0.00 Value Yield
0.05
Figure 35: Value yield distributions for the S&P 500 stock market index for different starting P/Book ratios. The shares are assumed to be sold after the given number of years and the selling price is calculated from Monte Carlo simulated equity multiplied by a P/Book ratio sampled from the historical distribution, as described in section 6.9.1. Also shown are the fitted normal PDFs.
84
P/Book=1, Year 2
6
P/Book=1, Year 5
P/Book=1, Year 10
10
Density
Density
Density 0 2
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.15
0.20
0.25
0.30 Value Yield
0.35
0.40
0.45
10
15
0.15
0.20 Value Yield
0.25
0.30
Value Yield
P/Book=3, Year 2
P/Book=3, Year 5
P/Book=3, Year 10
12
10
Density
Density
Density
-0.2
-0.1
0.0
0.1 Value Yield
0.2
0.3
0.4
-0.10
-0.05
0.00
0.05 Value Yield
0.10
0.15
0.20
10
15
0.00
0.05 Value Yield
0.10
0.15
P/Book=5, Year 2
P/Book=5, Year 5
P/Book=5, Year 10
10
Density
Density
Density
-0.2
-0.1 Value Yield
0.0
0.1
-0.15
-0.10
-0.05 Value Yield
0.00
0.05
0.10
10
15
-0.04
-0.02
0.00
0.02 Value Yield
0.04
0.06
0.08
Figure 36: Value yield distributions for the S&P 500 stock market index for different starting P/Book ratios. The shares are assumed to be sold after the given number of years and the selling price is calculated from the equity and P/Book ratios resulting from Monte Carlo simulations, as described in section 6.9.2. Also shown are the fitted normal PDFs.
85
Year
20 0.5 15 10 5 1 2
P/Book
3 4
Year
20
5
P/Book
5 1 2 3 4 5
15
10
0.20
0.4
Value Yield Stdev
Value Yield Mean
0.15
0.3
0.2
0.10
0.1
0.05
0.0
Figure 37: Value yield mean and standard deviation for the S&P 500 stock market index. The shares are bought at different P/Book ratios and are assumed to be sold after the given number of years. The selling price is calculated from the equity and P/Book ratios resulting from Monte Carlo simulations as described in section 6.9.2. Subsets are shown in Table 7 and Table 8.
Year
1 2 4 6 8 10 12 14 16 18 20
1 51% 39% 30% 25% 22% 20% 19% 18% 17% 16% 16%
1.5 39% 26% 19% 17% 15% 14% 14% 13% 13% 12% 12%
2 11% 11% 11% 11% 11% 10% 10% 10% 9.9% 9.8% 9.7%
2.5 13% 9.8% 8.5% 8.3% 8.1% 8.1% 8.2% 8.2% 8.2% 8.1% 8.1%
P/Book 3 3.5 2.4% 7.5% 3.8% 5.3% 4.9% 4.4% 5.4% 4.7% 6.0% 4.7% 6.1% 4.9% 6.3% 5.2% 6.5% 5.4% 6.6% 5.6% 6.7% 5.7% 6.8% 5.9%
4 4.5% 3.2% 3.0% 3.3% 3.4% 3.8% 4.1% 4.4% 4.6% 4.8% 5.0%
4.5 2.5% 0.9% 1.4% 1.9% 2.3% 2.8% 3.2% 3.5% 3.8% 4.0% 4.3%
5 (6.4)% (2.7)% (0.6)% 0.4% 1.2% 1.8% 2.3% 2.7% 3.0% 3.4% 3.6%
Table 7: Subset of the value yield means from Figure 37.
Year
1 2 4 6 8 10 12 14 16 18 20
1 16% 12% 5.6% 3.7% 2.9% 2.4% 2.1% 1.8% 1.6% 1.4% 1.3%
1.5 17% 9.1% 5.2% 3.8% 2.9% 2.4% 2.1% 1.8% 1.6% 1.4% 1.3%
2 15% 8.9% 5.5% 3.9% 3.0% 2.4% 2.0% 1.7% 1.5% 1.4% 1.3%
2.5 16% 9.7% 5.9% 4.2% 3.2% 2.5% 2.1% 1.9% 1.7% 1.5% 1.3%
P/Book 3 14% 8.8% 5.9% 4.4% 3.3% 2.6% 2.2% 1.9% 1.6% 1.5% 1.3%
3.5 20% 13% 6.8% 4.8% 3.6% 2.8% 2.3% 2.0% 1.7% 1.5% 1.3%
4 18% 11% 6.6% 4.7% 3.5% 2.9% 2.3% 2.0% 1.7% 1.5% 1.3%
4.5 13% 9.1% 6.0% 4.5% 3.6% 2.9% 2.4% 2.0% 1.8% 1.5% 1.4%
5 11% 8.1% 5.7% 4.3% 3.4% 2.8% 2.3% 2.0% 1.7% 1.5% 1.4%
Table 8: Subset of the value yield standard deviations from Figure 37.
86
Call Option Price (SPY ETF, S&P 500)

150 100
Put Option Price (SPY ETF, S&P 500)
Option Price (Bid Ask Avg.)
100
Option Price (Bid Ask Avg.)
50
50
100
150
200
250
20
40
60
80
50
100
150
200
250
Exercise Price
Exercise Price
Figure 38: Option chains for the SPY ETF (S&P 500) for call options (left) and put options (right) expiring on December 18, 2015 as quoted on March 15, 2013. The option quotes are shown with circles that are connected with lines for clarity. The dashed lines 12 show the closing share price on the previous trading day which was USD 156.32. The bold lines show the intrinsic values.
MC Simulated Share Price (SPY ETF, S&P 500)

0.012
MC Simulated Share Price (SPY ETF, S&P 500)

1.0
0.010
Density
0.008
0.006
0.004
F(x)
0.002
0.000
100
150
200
250 Share Price
300
350
0.0
0.2
0.4
0.6
0.8
100
150
200
250
300
350
Share Price
Figure 39: Histogram with fitted log-normal PDF (left) and CDF (right) of Monte Carlo simulated share prices for the SPY ETF (S&P 500) when the starting P/Book is 2.31 and the simulation period is 696 trading days, corresponding to the period between March 15, 2013 and December 18, 2015. The share price of SPY is 1/10 of the S&P 500 index. The closing share price on March 14, 2013 is shown as dashed lines.
12
Data source for the SPY options: http://quote.morningstar.com/Option/Options.aspx?ticker=SPY
87
Exercise Price 80 (SPY, S&P 500)

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
-50
50 Profit
100
150
200
-50
50 Profit
100
150
200
0.0
0.2
0.4
0.6
0.8
-50
50
100 Profit
150
200

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
-50
50
100 Profit
150
200
50
100 Profit
150
200
0.0
0.2
0.4
0.6
0.8
50 Profit
100
150

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
50 Profit
100
150
50 Profit
100
150
0.0
0.2
0.4
0.6
0.8
20
40
60 Profit
80
100
120
Figure 40: CDF for the profits on call options of SPY ETF (S&P 500) expiring on December 18, 2015 as quoted on March 15, 2013. The profits are calculated for the expiration date only, using Eq. 4-3 with the share prices resulting from Monte Carlo simulations as described in section 6.10.
88

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
-2.15
-2.10
-2.05 Profit
-2.00
-1.95
-1.90
-5
-4
-3
-2 Profit
-1
0.0
0.2
0.4
0.6
0.8
-10
-5
5 Profit
10
15

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
-20
-10
10 Profit
20
30
-20
0 Profit
20
40
0.0
0.2
0.4
0.6
0.8
-40
-20
0 Profit
20
40

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
-60
-40
-20
0 Profit
20
40
60
-50 Profit
50
0.0
0.2
0.4
0.6
0.8
-100
-50 Profit
50
Figure 41: CDF for the profits on put options of SPY ETF (S&P 500) expiring on December 18, 2015 as quoted on March 15, 2013. The profits are calculated for the expiration date only, using Eq. 4-3 with the share prices resulting from Monte Carlo simulations as described in section 6.10.
89
Call Option Profit Mean (SPY ETF, S&P 500)
Call Option Profit Stdev (SPY ETF, S&P 500)
25
20
Profit Stdev
50 100 150 200 250
Profit Mean
15
10
10
15
20
25
30
35
40
50
100
150
200
250
Exercise Price
Exercise Price
Figure 42: Profit mean and standard deviation for call options of SPY ETF (S&P 500) expiring on December 18, 2015 as quoted on March 15, 2013. The profits are calculated for the expiration date only, using Eq. 4-3 with the share prices resulting from Monte Carlo simulations as described in section 6.10. The options are shown with circles that are connected with lines for clarity. The dashed lines show the closing share price on the previous trading day.
Put Option Profit Mean (SPY ETF, S&P 500)

0
Put Option Profit Stdev (SPY ETF, S&P 500)
-5
-10
Profit Stdev
50 100 150 200 250
Profit Mean
-20
-15
-30
-25
10
20
30
50
100
150
200
250
Exercise Price
Exercise Price
Figure 43: Profit mean and standard deviation for put options of SPY ETF (S&P 500) expiring on December 18, 2015 as quoted on March 15, 2013. The profits are calculated for the expiration date only, using Eq. 4-3 with the share prices resulting from Monte Carlo simulations as described in section 6.10. The options are shown with circles that are connected with lines for clarity. The dashed lines show the closing share price on the previous trading day.
90
Call Option Profit Probability (SPY ETF, S&P 500)
Put Option Profit Probability (SPY ETF, S&P 500)
0.7
0.6
Profit Probability
Profit Probability
50 100 150 200 250
0.4
0.3
0.1
0.00
0.05
0.2
0.10
0.15
0.5
0.20
50
100
150
200
250
Exercise Price
Exercise Price
Figure 44: Profit probability for call options (left) and put options (right) of SPY ETF (S&P 500) expiring on December 18, 2015 as quoted on March 15, 2013. The profits are calculated for the expiration date only, using share prices resulting from Monte Carlo simulations as described in section 6.10. The options are shown with circles that are connected with lines for clarity. The dashed lines show the closing share price on the previous trading day. Note that the scales of the y-axes (profit probability) are different in the two plots.
91

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
-0.4
-0.2
0.0
0.2
0.4
0.6
-1.0
-0.5
0.0 Value Yield
0.5
0.0
0.2
0.4
0.6
0.8
-1.0
-0.5
0.0 Value Yield
0.5
1.0
Value Yield

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
-1.0
-0.5
0.0 Value Yield
0.5
1.0
1.5
-1.0
-0.5
0.0
0.5 Value Yield
1.0
1.5
2.0
0.0
0.2
0.4
0.6
0.8
-1
1 Value Yield

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
-1
1 Value Yield
-1
2 Value Yield
0.0
0.2
0.4
0.6
0.8
2 Value Yield
Figure 45: CDF for the value yields on call options of SPY ETF (S&P 500) expiring on December 18, 2015 as quoted on March 15, 2013. The value yields are calculated for the expiration date only, using Eq. 4-6 with the share prices resulting from Monte Carlo simulations as described in section 6.10.
92

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
-1.06
-1.04
-1.02
-1.00 Value Yield
-0.98
-0.96
-0.94
-1.0
-0.5 Value Yield
0.0
0.0
0.2
0.4
0.6
0.8
-1.0
-0.5 Value Yield
0.0
0.5

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
-1.0
-0.5 Value Yield
0.0
0.5
-1.0
-0.5 Value Yield
0.0
0.5
0.0
0.2
0.4
0.6
0.8
-1.0
-0.5 Value Yield
0.0
0.5

1.0 1.0

1.0
0.8
0.8
0.6
0.6
F(x)
F(x)
0.4
0.4
F(x) 0.2 0.0
0.2
0.0
-1.0
-0.5 Value Yield
0.0
-1.0
-0.8
-0.6
-0.4 Value Yield
-0.2
0.0
0.2
0.0
0.2
0.4
0.6
0.8
-1.0
-0.8
-0.6
-0.4 Value Yield
-0.2
0.0
0.2
Figure 46: CDF for the value yields on put options of SPY ETF (S&P 500) expiring on December 18, 2015 as quoted on March 15, 2013. The value yields are calculated for the expiration date only, using Eq. 4-6 with the share prices resulting from Monte Carlo simulations as described in section 6.10.
93
Call Option, Value Yield Mean (SPY ETF, S&P 500)

1.2
Call Option, Value Yield Stdev (SPY ETF, S&P 500)
0.0
-0.2
Value Yield Stdev

50 100 150 200 250
Value Yield Mean
-0.4
-0.6
0.2
0.4
0.6
0.8
1.0
50
100
150
200
250
Exercise Price
Exercise Price
Figure 47: Value yield mean and standard deviation for call options of SPY ETF (S&P 500) expiring on December 18, 2015 as quoted on March 15, 2013. The value yields are calculated for the expiration date only, using Eq. 4-6 and the share prices resulting from Monte Carlo simulations as described in section 6.10. The options are shown with circles that are connected with lines for clarity. The dashed lines show the closing share price on the previous trading day.
Put Option, Value Yield Mean (SPY ETF, S&P 500)

-0.2
Put Option, Value Yield Stdev (SPY ETF, S&P 500)
-0.4
Value Yield Stdev

50 100 150 200 250
Value Yield Mean
-0.6
-0.8
-1.0
0.0
0.1
0.2
0.3
0.4
50
100
150
200
250
Exercise Price
Exercise Price
Figure 48: Value yield mean and standard deviation for put options of SPY ETF (S&P 500) expiring on December 18, 2015 as quoted on March 15, 2013. The value yields are calculated for the expiration date only, using Eq. 4-6 and the share prices resulting from Monte Carlo simulations as described in section 6.10. The options are shown with circles that are connected with lines for clarity. The dashed lines show the closing share price on the previous trading day.
94
Call Option Price vs. Profit Mean (SPY ETF, S&P 500)
Call Option Price vs. Profit Probability (SPY ETF, S&P 500)
25
20
Profit Probability (Maximize)

0 50 100 150
Profit Mean (Maximize)
15
10
0.1
0
0.2
0.3
0.4
0.5
0.6
0.7
50
100
150
Option Price (Minimize)
Figure 49: Call option price versus profit mean (left) and profit probability (right) for options of SPY ETF (S&P 500) expiring on December 18, 2015 as quoted on March 15, 2013. The profits are calculated from Monte Carlo simulations as described in section 6.10. The buyer of a call option would want to maximize the profit mean and probability (y-axes) and minimize the option price (x-axes). The seller would want to do the opposite.
Put Option Price vs. Profit Mean (SPY ETF, S&P 500)
0
Put Option Price vs. Profit Probability (SPY ETF, S&P 500)
-5
Profit Probability (Maximize)

0 20 40 60 80 100
-10
Profit Mean (Maximize)
-15
-20
-30
-25
0.00
0
0.05
0.10
0.15
0.20
20
40
60
80
100
Figure 50: Put option price versus profit mean (left) and profit probability (right) for options of SPY ETF (S&P 500) expiring on December 18, 2015 as quoted on March 15, 2013. The profits are calculated from Monte Carlo simulations as described in section 6.10. The buyer of a call option would want to maximize the profit mean and probability (y-axes) and minimize the option price (x-axes). The seller would want to do the opposite.
95
USA Gov. Bond Yield Long Maturity / %
USA Gov. Bond Yield 10-Year Maturity / %
15
10
Yield / %
Yield / % 5 5
10
15
1800
1850
1900 Year
1950
2000
1960
1970
1980
1990 Date
2000
2010
Figure 51: USA government bond yields for the period 1798-2012 (left plot, annual average, bonds have varying maturity, terms and taxation during this period) and for the period 1962-2013 (right plot, daily average, bonds have 10 year maturity period). 13 Because of the annual averaging in the left plot the extreme bond yields are smoothened somewhat.
USA Gov. Bond Yields (1798-2012)
USA Gov. Bond Yields (1962-2013)
Frequency
Frequency
8 Yield / %
10
12
200
400
600
5 Yield / %
10
15
Figure 52: Histograms for the bond yields in Figure 51.
13
Bond yields for the period 1798-2012 are from Homer & Sylla [13] tables 38, 46, 48, 51, 87, as well as from the Federal Reserve. Bond yields for the period 1962-2013 are also from the Federal Reserve:
www.federalreserve.gov/releases/h15/data.htm 96
Arithmetic Geometric Harmonic Stdev Min Max Mean Mean Mean 1798-2012 4.6% 4.2% 3.8% 2.0% 1.8% 12.9% 1962-2013 6.6% 6.1% 5.5% 2.7% 1.4% 15.8% Period
Table 9: Statistics for the USA government bond yields in Figure 51. Because the bond yields for the period 1798-2012 are annual averages, the minimum and maximum yields for that period are inaccurate.
USA Gov. Bond Yield (1798-2012)

1.0 1.0
USA Gov. Bond Yield (1962-2013)
0.8
Autocorrelation
0.6
Autocorrelation
0.4
0.2
0.0
-0.2
-0.4
50
100 Lag/Years
150
200
-0.4
0
-0.2
0.0
0.2
0.4
0.6
0.8
2000
4000
6000 Lag/Days
8000
10000
12000
Figure 53: Autocorrelation for the USA government bond yields from Figure 51 which show significant correlation between successive bond yields. The reason the shapes of these plots differ is perhaps because the left plot is for annual averages and the right is for daily averages.
97
Equity Risk Premium / %

80
10 Equity Risk Premium / %
Frequency
-10
-20
20
40
60
1960
1970
1980 Date
1990
2000
2010
-20
-10
10
Figure 54: Equity risk premium for the period 1957-2013, calculated as the difference between the monthly yields of USA government bonds with 10 year maturity and the monthly returns of the S&P 500 stock market index. Left plot shows the time series and right plot shows the histogram.
S&P 500 Monthly Return vs. USA Gov. Bond Monthly Yield
1.2
USA Gov. Bond Monthly Yield / %
0.8
1.0
Cross-correlation
0.4
0.6
0.2
-0.05
0.00
0.05
-20
-10
10
-600
-400
-200
0 Lag/Months
200
400
600
S&P 500 Monthly Return / %
Figure 55: Scatter-plot (left) of the monthly returns of the S&P 500 stock market index versus the monthly yields on USA government bonds with 10 year maturity. Data is for the period 1957-2013 from Figure 54. Right plot shows cross-correlation.
98
USA Gov. Bond Yield 10-Year Maturity / %

14 14
S&P 500 Mean Value Yield / %

8
12
12
10
10
Value Yield / %
Yield / %
1990 Date
2000
2010
1985
1990
1995 Date
2000
2005
2010
-2
1985
1990
1995 Date
2000
2005
2010
Figure 56: USA government bond yields with 10 year maturity (left plot), mean value yield for the S&P 500 stock market index (middle plot) estimated from Eq. 6-28 with the P/Book data from Figure 22. Their difference is the Equity Risk Premium (ERP) (right plot). Equity Risk Premium (1984-2011)
300 Frequency 0
-2
50
100
150
200
250
Figure 57: Histogram of the Equity Risk Premium (ERP) in Figure 56.
14
12
USA Gov. Bond Yields / %
10
Cross-correlation
10
11
12
-0.4
-0.2
0.0
0.2
0.4
0.6
-6000
-4000
-2000
0 Lag/Days
2000
4000
6000
S&P 500 Mean Value Yield / %
Figure 58: Scatter-plot (left) of the mean value yield for the S&P 500 stock market index versus USA government bond yields with 10 year maturity. Data is for the period 1984-2011 from Figure 56. Cross-correlation of the data is shown in the right plot.
99
Dow Jones Venture Capital Index (DJVC)

0.3
3000
2500
2000
Monthly Return
Index Level
1000
1500
500
1995
2000 Date
2005
2010
-0.2
-0.1
0.0
0.1
0.2
1995
2000 Date
2005
2010
Figure 59: Dow Jones Venture Capital (DJVC) index total return (left) and monthly returns (right).
14

30
1.0
25
Frequency
20
Autocorrelation
15
10
-0.3
-0.2
-0.1
0.0 Monthly Return
0.1
0.2
0.3
0.0
0
0.2
0.4
0.6
0.8
50
100 Lag/Months
150
200
Figure 60: Dow Jones Venture Capital (DJVC) monthly returns histogram (left) and autocorrelation (right).
14
Data retrieved from Sand Hill Econometrics in April 2013: www.sandhillecon.com/downloads/DJIndexOfVentureCapital-2010Q2.xls
100

DJVC Monthly Return vs. USA Gov. Bond Yield
0.20 0.006 USA Gov. Bond Monthly Yield 0.005 Cross-correlation 0.004 0.003 0.002 -0.10 -0.05 0.00 0.05 0.10 0.15
-0.2
-0.1
0.0
0.1
0.2
0.3
-200
-100
0 Lag/Months
100
200
DJVC Monthly Total Return
Figure 61: Scatter-plot (left) and cross-correlation (right) for the monthly returns on the Dow Jones Venture Capital (DJVC) index versus the yields on USA government bonds with 10 year maturity, during the period December 1991 to March 2010.
DJVC Monthly Total Return vs. S&P 500 Monthly Return
0.10
0.05
S&P 500 Monthly Return
0.00
Cross-correlation
-0.2 -0.1 0.0 0.1 0.2 0.3
-0.05
-0.10
-0.15
0.0
0.2
0.4
0.6
0.8
-200
-100
0 Lag/Months
100
200
DJVC Monthly Total Return
Figure 62: Scatter-plot (left) and cross-correlation (right) for the monthly returns on the Dow Jones Venture Capital (DJVC) index versus the monthly percentage returns on the S&P 500 stock market index, during the period December 1991 to March 2010.
101

Monthly Return Difference, DJVC - S&P 500
40
Monthly Return Difference, DJVC - S&P 500
0.2
0.1
Frequency
Return Difference
0.0
-0.1
-0.2
10
20
30
1995
2000 Date
2005
2010
-0.2
-0.1
0.0 Return Difference
0.1
0.2
Figure 63: Difference between monthly returns of the Dow Jones Venture Capital (DJVC) index and the S&P 500 stock market index. Left plot shows the difference over time and right plot shows the histogram.
102

Simulated Compounded Return (DJVC)
1e+08
1e+06
1e+06
1e+04
1e+04
1e+02
1e+02
1e+00
20
40 Year
60
80
100
1e+00
20
40 Year
60
80
100
1e+00
1e+02
1e+04
1e+06
20
40 Year
60
80
100
10000
10000
100
20
40 Year
60
80
100
20
40 Year
60
80
100
1e+00
1e+02
1e+04
100
1e+06
1e+08
20
40 Year
60
80
100
10000
1000
1e+06
100
1e+04
10
1e+02
20
40 Year
60
80
100
1e+00
20
40 Year
60
80
100
1e+00
1e+02
1e+04
1e+06
20
40 Year
60
80
100
Figure 64: Compounded returns of the Dow Jones Venture Capital (DJVC) index which are simulated by randomly sampling the historical monthly returns of the DJVC index and compounding (multiplying) the samples. The y-axes are logarithmic.
103
Year 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 ROE 56.0% 55.7% 57.0% 64.8% 67.1% 48.6% 28.9% 22.9% 42.6% 26.8% Retain 34.8% 16.6% 5.6% 24.3% 39.5% 22.4% 41.3% 31.7% 52.0% 15.7% Year 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 ROE 65.8% 54.2% 30.6% 31.1% 35.3% 26.7% 33.3% 47.6% 27.7% 28.5% Retain 54.8% 48.0% 7.6% (2%) 43.7% 30.9% 31.8% 54.6% 15.6% 15.0%
Table 10: Return on Equity (ROE) and earnings retaining ratio of Coca-Cola during the period 1993-2012. Statistics are shown in 15 Table 11.
Mean Geo. Mean Stdev ROE 42.6% 41.8% 15.0% Retain 29.2% 28.1% 17.0%
Table 11: Statistics for Table 10.
ROE vs. Retain 1993-2012 (Coca-Cola)
ROE vs. Retain 1993-2012 (Coca-Cola)
0.5
0.4
Cross-correlation
Retain
0.3
0.2
0.1
0.0
-0.4
-0.2
0.0
0.2
0.4
0.3
0.4 ROE
0.5
0.6
-10
-5
0 Lag / Years
10
Figure 65: Scatter-plot (left) and cross-correlation (right) for the ROE and Retain ratios of Coca-Cola during the period 1993-2012. Data is shown in Table 10.
15
Form 10-K annual reports filed with US SEC: www.sec.gov/cgi-bin/browse-edgar?action=getcompany&CIK=0000021344&type=10-K
104

Value Yield (Coca-Cola)
Value Yield (Coca-Cola)

30
1.0
25
Density
20
15
F(x)
10
0.14
0.16
0.18 Value Yield
0.20
0.22
0.0
0.2
0.4
0.6
0.8
0.14
0.16
0.18 Value Yield
0.20
0.22
Figure 66: Histogram with fitted log-normal PDF (left) and CDF (right) for the value yield of the Coca Cola company when its P/Book ratio is 5.5 and its future payouts result from Monte Carlo simulations as described in section 9.2. The distribution converges to a log-normal PDF as the number of simulations increase.
Value Yield, Scaled MCS (Coca-Cola)
Value Yield, Scaled MCS (Coca-Cola)

1.0
30
Density
20
F(x)
10
0.10
0.12
0.14 Value Yield
0.16
0.18
0.0
0.2
0.4
0.6
0.8
0.08
0.10
0.12 Value Yield
0.14
0.16
Figure 67: Histogram with fitted log-normal PDF (left) and CDF (right) for the value yield of the Coca Cola company when its P/Book ratio is 5.5 and its future payouts result from Monte Carlo simulations where the ROE and Retain ratios are decreased exponentially over time so as to make the equity growth rate converge to that of the S&P 500 index, as described in section 9.4. The distribution converges to a log-normal PDF as the number of simulations increase.
105
Coca-Cola P/Book=5.5, S&P 500 P/Book=2.3

40
Coca-Cola P/Book=5, S&P 500 P/Book=3.5
25
30
Density
20
Density
10
15
0.00
0.02
0.04
0.06
0.08
0.10
10
20
30
0.02
0.04
0.06
0.08
0.10
Coca-Cola VY - S&P 500 VY
Coca-Cola VY - S&P 500 VY
Figure 68: Difference between the value yield of Coca-Cola and the value yield of the S&P 500 stock market index for different P/Book ratios. The Monte Carlo simulations are described in section 9.5. Also shown are the fitted normal PDFs.
Coca-Cola P/Book=5.5, S&P 500 P/Book=2.3

10 10
Coca-Cola P/Book=5, S&P 500 P/Book=3.5
Probability
Probability 2 0
-0.05 0.00 0.05 Risk Premium 0.10 0.15
0
-0.05
0.00
0.05 Risk Premium
0.10
0.15
Figure 69: Probability of the value yield of Coca Cola being greater than the value yield of the S&P 500 stock market index for different P/Book ratios and with varying risk premiums: . The Monte Carlo simulations are described in section 9.5.
106
Coca-Cola PV/Equity (S&P 500 P/Book=2.3, Risk Premium=0%)

0.15
0.10
Density
Density 0.05 0.00
10
15
20
25
30
35
0.00
0.05
0.10
0.15
0.20
0.25
10 Present Value / Equity
15
Present Value / Equity
0.08
Density
0.06
Density
0.04
0.00
0.02
10
20
40
50
0.00
0.05
0.10
0.15
0.20
10
20
Figure 70: Probability distributions for the present value of Coca-Cola when the future payouts are Monte Carlo simulated as described in section 9.4, and the discount rates used in the present value calculations are sampled from the log-normal distribution of the value yield for the S&P 500 index described in section 6.7.1. The present values are normalized by dividing with the equity and should therefore be compared to the companys P/Book ratio which is currently about 5.5 and shown as dashed lines in these plots. Also shown are the fitted log-normal PDFs. The CDFs are shown in Figure 71.
107

1.0 1.0 F(x)
5 10 15 20 25 30 35
0.8
0.6
F(x)
0.4
0.2
0.0
0.0
0.2
0.4
0.6
0.8
15

1.0 1.0 F(x)
10 20 30 40 50
0.8
0.6
F(x)
0.4
0.2
0.0
0.0
0.2
0.4
0.6
0.8
10
20
Figure 71: CDF for the present value of Coca-Cola when the future payouts are Monte Carlo simulated as described in section 9.4, and the discount rates used in the present value calculations are sampled from the log-normal distribution of the value yield for the S&P 500 index described in section 6.7.1. The present values are normalized by dividing with the equity and should therefore be compared to the companys P/Book ratio which is currently about 5.5 and shown as dashed lines in these plots. The histograms are shown in Figure 70.
108
Coca-Cola PV/Book (DJVC Compounded)

1.0
Coca-Cola PV/Equity (DJVC Compounded)
0.15
Density
0.10
F(x)
0 10 20 30 40 50
0.05
0.00
0.0
0.2
0.4
0.6
0.8
10
20
30
40
50
Present Value
Figure 72: Histogram with fitted log-normal PDF (left) and CDF (right) for the present value of Coca-Cola when the future payouts are Monte Carlo simulated as described in section 9.4, and the discount rates used in the present value calculations are simulations of the compounded returns on the Dow Jones Venture Capital (DJVC) index as described in section 9.6.2. No risk premium is used here. The present values are normalized by dividing with the equity and should therefore be compared to the companys P/Book ratio which is currently about 5.5 and shown as dashed lines in these plots. There is a single present value above 800 in this distribution but the plots have been limited to only show present values below 50 to improve legibility.
Coca-Cola PV/Book (DJVC Annualized)

1.0
Coca-Cola PV/Equity (DJVC Annualized)
0.4
0.3
Density
0.2
F(x) 0.1 0.0 0.0

0
0.2
0.4
0.6
0.8
10 Present Value
15
20
10
15
20
25
Figure 73: Histogram with fitted log-normal PDF (left) and CDF (right) for the present value of Coca-Cola when the future payouts are Monte Carlo simulated as described in section 9.4, and the discount rates used in the present value calculations are simulations of the annualized returns on the Dow Jones Venture Capital (DJVC) index as described in section 9.6.2. No risk premium is used here. The present values are normalized by dividing with the equity and should therefore be compared to the companys P/Book ratio which is currently about 5.5 and shown as dashed lines in these plots.
109
Compounded vs. Annualized Growth (DJVC)

1e+20 1e+00
0
1e+05
1e+10
Index
1e+15
50
100
150 Year
200
250
300
Figure 74: Comparison of compounded (volatile line) and annualized (straight line) growth which show identical results for the last year. The y-axis is logarithmic.
Coca-Cola PV (DJVC Compounded vs. Annualized)
Annualized 5
0
10
15
20
200
400 Compounded
600
800
Figure 75: QQ-plot comparing the probability distributions for the present value of Coca-Cola when the discount rate is the compounded (x-axis) and annualized (y-axis) returns of the DJVC index from Figure 72 and Figure 73.
110
Minimum Compounded Returns (DJVC vs. S&P 500)

1e+15 1e+24
Maximum Compounded Returns (DJVC vs. S&P 500)
1e+11
Compounded Return
Compounded Return
0 50 100 150 Year 200 250 300
1e+07
1e+03
1e-01
1e+00
0
1e+06
1e+12
1e+18
50
100
150 Year
200
250
300
Figure 76: Simulated compounded returns for the S&P 500 index (straight lines) and the DJVC index (irregular lines). Note that the y-axes are logarithmic. The left plot shows the lowest compounded returns and the right plot shows the highest compounded returns for the simulations in section 9.6, which are used as the discount rate in the calculation of the present value of Coca-Cola.
Coca-Cola Present Value Difference (S&P 500)

0.7
Coca-Cola Present Value Difference (DJVC)
0.6
0.4
0.5
Density
0.2
0.3
Density
0.1
-3
-2
-1
0.00
0.0
0.05
0.10
0.15
10
20
30
40
50
Present Value Difference
Present Value Difference
Figure 77: Present value of Coca-Cola with the discount rate being sampled from the stochastic value yield of the S&P 500 stock market index, minus the present value with the discount rate being the mean S&P 500 value yield (left plot). Similarly when the discount rate is the annualized return of the DJVC venture index (right plot, x-axis is limited at 50 for clarity). This shows the present value is significantly affected when calculated using the mean discount rate.
111
13.
Bibliography
[1] P. Glasserman, Monte Carlo Methods in Financial Engineering.: Springer, 2003. [2] J.B. Williams, The Theory of Investment Value.: Harvard University Press, 1938. [3] M.E.H. Pedersen, "Introduction to Share Buyback Valuation," Hvass Laboratories, Report HL-1301, 2013. [Online]. http://www.hvass-labs.org/people/magnus/publications/pedersen2013share-buybackintro.pdf [4] M.E.H. Pedersen, "The Value of Share Buybacks," Hvass Laboratories, Report HL-1201, 2012. [Online]. http://www.hvass-labs.org/people/magnus/publications/pedersen2013share-buyback-intro.pdf [5] R.G. Newell and W.A. Pizer, "Discounting the Distant Future: How Much Do Uncertain Rates Increase Valuations?," Journal of Environmental Economics and Management, vol. 46, pp. 52-71, 2003. [6] F. Black and M. Scholes, "The Pricing of Options and Corporate Liabilities," Journal of Political Economy, vol. 81, no. 3, pp. 637-654, 1973. [7] R.C. Merton, "Theory of Rational Option Pricing," Bell Journal of Economics and Management Science, vol. 4, no. 1, pp. 141-183, 1973. [8] F. Black, "How to Use the Holes in Black-Scholes," Journal of Applied Corporate Finance, vol. 1, no. 4, pp. 67-73, 1989. [9] W.E. Buffett, "Letter to Shareholders of Berkshire Hathaway," February 27, 2009. [10] J.C. Cox, S.A. Ross, and M. Rubinstein, "Option Pricing: A Simplified Approach," Journal of Financial Economics, vol. 7, no. 3, pp. 229-263, 1979. [11] P.B. Hammond, M.L. Leibowitz, and L.B. Siegel, Eds., Rethinking the Equity Risk Premium.: Research Foundation of CFA Institute, 2011. [12] A. Damodaran, "Equity Risk Premiums (ERP): Determinants, Estimation and Implications," Stern School of Business, Working Paper 2012. [13] S. Homer and R. Sylla, A History of Interest Rates, 4th ed.: John Wiley & Sons, 2005.
112
14.

Revision History
2013, June 4: Extended sections 7.1, 9.6. Added Table 5, Table 6, Figure 54, Figure 55, Figure 73, Figure 74, Figure 75. Minor clarifications and corrections. (113 Pages) 2013, May 24: First edition. (109 pages)
113

Pedersen2013monte Carlo

Uploaded by

Copyright:

Available Formats

Pedersen2013monte Carlo

Uploaded by

Document Information

Original Description:

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Pedersen2013monte Carlo

Uploaded by

Copyright:

Available Formats

Monte Carlo Simulation in Financial Valuation

Monte Carlo Simulation in Financial Valuation

Monte Carlo Simulation in Financial Valuation

The paper is structured as follows:

Monte Carlo Simulation in Financial Valuation

This means the payout can then also be written as:

Monte Carlo Simulation in Financial Valuation

Equity Growth Model

This is equivalent to:

be the fraction of earnings retained in year so that:

This means that the equity grows by the rate

Monte Carlo Simulation in Financial Valuation

The accumulation of equity is:

The value to eternal shareholders is then:

Monte Carlo Simulation in Financial Valuation

Equity Growth for Two Companies

Share Issuance & Buyback

Monte Carlo Simulation in Financial Valuation

This is equivalent to:

Monte Carlo Simulation in Financial Valuation

If the payout growth first starts next year so that

then Eq. 12-3 gives:

Monte Carlo Simulation in Financial Valuation

and . The value yield is

. but . The value yield is then:

Monte Carlo Simulation in Financial Valuation

The geometric series

is reduced using Eq. 12-2:

The terminal value consisting of the payouts in year

and onwards is calculated using Eq. 12-4:

Monte Carlo Simulation in Financial Valuation

The limit of the terminal value is zero as

approaches infinity provided

The number of iterations

required in order for

to be within relative error

In the derivation of this, the inequality is reversed because assumed here.

Monte Carlo Simulation in Financial Valuation

would be calculated from Eq. 2-15:

When the equity at the start of year

is known, the expected payout can be calculated from Eq. 2-13:

Stochastic Present Value

This becomes a strict inequality if

is strictly convex and

Monte Carlo Simulation in Financial Valuation

which is strictly convex for

, then according to Jensens inequality:

which is strictly convex for

, then according to Jensens inequality:

be another stochastic variable which is independent of , then Eq. 3-2 gives:

This inequality also holds for powers of

Stochastic Discount Rate

Monte Carlo Simulation in Financial Valuation

Using Jensens inequality in Eq. 3-2 gives:

Stochastic Growth Rate

is a derived stochastic variable

Monte Carlo Simulation in Financial Valuation

Where the accumulated growth is denoted:

The mean payout in year is then: