We apply Markowitz's MVO model to the problem of constructing a portfolio of US stocks, bonds and cash. We use historical data for the returns of these three asset classes: The S&P 500 index for the returns on stocks, the 10-year Treasury bond index for the returns on bonds, and we assume that the cash is invested in a money market account whose return is the 1-day federal fund rate. The times series for the "Total Return" are given below for each asset between 1960 and 2003.
| Year | Stocks | Bonds | MM |
| 1960 | 20.2553 | 262.935 | 100.00 |
| 1961 | 25.6860 | 268.730 | 102.33 |
| 1962 | 23.4297 | 284.090 | 105.33 |
| 1963 | 28.7463 | 289.162 | 108.89 |
| 1964 | 33.4484 | 299.894 | 113.08 |
| 1965 | 37.5813 | 302.695 | 117.97 |
| 1966 | 33.7839 | 318.197 | 124.34 |
| 1967 | 41.8725 | 309.103 | 129.94 |
| 1968 | 46.4795 | 316.051 | 137.77 |
| 1969 | 42.5448 | 298.249 | 150.12 |
| 1970 | 44.2212 | 354.671 | 157.48 |
| 1971 | 50.5451 | 394.532 | 164.00 |
| 1972 | 60.1461 | 403.942 | 172.74 |
| 1973 | 51.3114 | 417.252 | 189.93 |
| 1974 | 37.7306 | 433.927 | 206.13 |
| 1975 | 51.7772 | 457.885 | 216.85 |
| 1976 | 64.1659 | 529.141 | 226.93 |
| 1977 | 59.5739 | 531.144 | 241.82 |
| 1978 | 63.4884 | 524.435 | 266.07 |
| 1979 | 75.3032 | 531.040 | 302.74 |
| 1980 | 99.7795 | 517.860 | 359.96 |
| 1981 | 94.8671 | 538.769 | 404.48 |
| 1982 | 115.308 | 777.332 | 440.68 |
| 1983 | 141.316 | 787.357 | 482.42 |
| 1984 | 150.181 | 907.712 | 522.84 |
| 1985 | 197.829 | 1200.63 | 566.08 |
| 1986 | 234.755 | 1469.45 | 605.20 |
| 1987 | 247.755 | 1424.91 | 646.17 |
| 1988 | 288.116 | 1522.40 | 702.77 |
| 1989 | 379.409 | 1804.63 | 762.16 |
| 1990 | 367.636 | 1944.25 | 817.87 |
| 1991 | 479.633 | 2320.64 | 854.10 |
| 1992 | 516.178 | 2490.97 | 879.04 |
| 1993 | 568.202 | 2816.40 | 905.06 |
| 1994 | 575.705 | 2610.12 | 954.39 |
| 1995 | 792.042 | 3287.27 | 1007.84 |
| 1996 | 973.897 | 3291.58 | 1061.15 |
| 1997 | 1298.82 | 3687.33 | 1119.51 |
| 1998 | 1670.01 | 4220.24 | 1171.91 |
| 1999 | 2021.40 | 3903.32 | 1234.02 |
| 2000 | 1837.36 | 4575,33 | 1313.00 |
| 2001 | 1618.98 | 4827.26 | 1336.89 |
| 2002 | 1261.18 | 5558.40 | 1353.47 |
| 2003 | 1622.94 | 5588.19 | 1366.73 |
Let I
it denote the above "Total Return" for asset i = 1,2,3 and t = 0, ... T, where t = 0 corresponds to
1960 and t = T to 2003. For each asset i, we can convert the raw data I
it, t=0,...,T, into rates of returns r
it, t=1,...,T, using the formula.
r
it = I
i,t - I
i, t-1 / I
i, t-1| Year | Stocks | Bonds | MM |
| 1961 | 26.81 | 2.20 | 2.33 |
| 1962 | -8.78 | 5.72 | 2.93 |
| 1963 | 22.69 | 1.79 | 3.38 |
| 1964 | 16.36 | 3.71 | 3.85 |
| 1965 | 12.36 | 0.93 | 4.32 |
| 1966 | -10.10 | 5.12 | 5.40 |
| 1967 | 23.94 | -2.86 | 4.51 |
| 1968 | 11.00 | 2.25 | 6.02 |
| 1969 | -8.47 | -5.63 | 8.97 |
| 1970 | 3.94 | 18.92 | 4.90 |
| 1971 | 14.30 | 11.24 | 4.14 |
| 1972 | 18.99 | 2.39 | 5.33 |
| 1973 | -14.69 | 3.29 | 9.95 |
| 1974 | -26.47 | 4.00 | 8.53 |
| 1975 | 37.23 | 5.52 | 5.20 |
| 1976 | 23.93 | 15.56 | 4.65 |
| 1977 | -7.16 | 0.38 | 6.65 |
| 1978 | 6.57 | -1.26 | 10.03 |
| 1979 | 18.61 | -1.26 | 13.78 |
| 1980 | 32.50 | -2.48 | 18.90 |
| 1981 | -4.92 | 4.04 | 12.37 |
| 1982 | 21.55 | 44.28 | 8.95 |
| 1983 | 22.56 | 1.29 | 9.47 |
| 1984 | 6.27 | 15.29 | 8.38 |
| 1985 | 31.17 | 32.27 | 8.27 |
| 1986 | 18.67 | 22.39 | 6.91 |
| 1987 | 5.25 | -3.03 | 6.77 |
| 1988 | 16.61 | 6.84 | 8.76 |
| 1989 | 31.69 | 18.54 | 8.45 |
| 1990 | -3.10 | 7.74 | 7.31 |
| 1991 | 30.46 | 19.36 | 4.43 |
| 1992 | 7.62 | 7.34 | 2.92 |
| 1993 | 10.08 | 13.06 | 2.96 |
| 1994 | 1.32 | -7.32 | 5.45 |
| 1995 | 37.58 | 25.94 | 5.60 |
| 1996 | 22.96 | 0.13 | 5.29 |
| 1997 | 33.36 | 12.02 | 5.50 |
| 1998 | 28.58 | 14.45 | 4.68 |
| 1999 | 21.04 | -7.51 | 5.30 |
| 2000 | -9.10 | 17.22 | 6.40 |
| 2001 | -11.89 | 5.51 | 1.82 |
| 2002 | -22.10 | 15.15 | 1.24 |
| 2003 | 28.68 | 0.54 | 0.98 |
Let R
i denote the random rate of return of asset i. From the above historical data, we can compute the arithmetric mean rate of return for each asset:
_ T
r
i = 1/T SUM r
it t=1
which gives
| | Stocks | Bonds | MM |
| Arithmetric mean | 12.06% | 7.85% | 6.32% |
Because
the rates of return are multiplicative over time, we prefer t use the geometric mean instead of the arithmetric mean. The geometric mean is the constant yearly rate of return that needs to be applied in years t = 0 through t = T-1 in order to get the compounded Total Return I
iT, starting from I
i0. The formula for the geometric mean is:
ui = ( Π (1+r
it) )
1/T -1
t=1...T
We get the following results.
| | Stocks | Bonds | MM |
| Geometric mean | 10.73% | 7.37% | 6.27% |
We also compute the covariance matrix:
_ _
cov(Ri, Rj) = (1/T) SUM (r
it - r)(r
jt - r).
t=1...T
| Covariance | Stocks | Bonds | MM |
| Stocks | 0.02778 | 0.00387 | 0.00021 |
| Bonds | 0.00387 | 0.01112 | -0.00020 |
| MM | 0.00021 | -0.00020 | 0.00115 |
Although not needed to slolve the Markowitz model, it is interesting to compute the volatility of the rate of return on each asset σ
i =sqrt(cov(Ri, R,)):
| | Stocks | Bonds | MM |
| Volatility | 16.67% | 10.55% | 3.40% |
and the correlation matrix p
ij = cov(R
i, R
j) / σ
iσj
| Correlation | Stocks | Bonds | MM |
| Stocks | 1 | 0.2199 | 0.0366 |
| Bonds | 0.2199 | 1 | -0.0545 |
| MM | 0.0366 | -0.0545 | 1 |
Setting up the QP for portfolio optimization
min 0.02778x
2 S + 2×0.00387x
Sx
B + 2×0.00387x
Sx
M + 0.01112x
2B - 2×0.00020x
Bx
M + 0.00115x
2M + 0.1073x
S + 0.0737x
B + 0.0627x
M0.1073x
S + 0.0737x
B + 0.0627x
M ≧ R
x
S + x
B + x
M = 1
x
S, x
B, x
M ≧ 0
and solving it for R=6.5% to R=10.5% with increments of 0.5% we get the optimal portfolios shown below and the corresponding variance. The optimal allocations on the efficient frontier are also depicted in the graph below.
| Rate of Return R | Variance | Stocks | Bonds | MM |
| 0.065 | 0.0010 | 0.03 | 0.10 | 0.87 |
| 0.070 | 0.0014 | 0.13 | 0.12 | 0.75 |
| 0.075 | 0.0026 | 0.24 | 0.14 | 0.62 |
| 0.080 | 0.0044 | 0.35 | 0.16 | 0.49 |
| 0.085 | 0.0070 | 0.45 | 0.18 | 0.37 |
| 0.090 | 0.0102 | 0.56 | 0.20 | 0.24 |
| 0.095 | 0.0142 | 0.67 | 0.22 | 0.11 |
| 0.100 | 0.0189 | 0.78 | 0.22 | 0 |
| 0.105 | 0.0246 | 0.93 | 0.07 | 0 |

Based on the first two columns of table above, the graph plots the maximum expected rate of return R of a portfolio as a function of its volatility (standard deviation). This curve is called the effecient frontier. Every possible portfolio of Stocks/Bonds/MM is represented by a point lying on or below the efficient frontier in the expected return/standard deviation plane.
[Quote]:Optimization Methods in Finance