Without doubt, the Markowitz (1952) mean-variance portfolio theory is a cornerstone of modern financial theory. This section illustrates the use of the mean_variance_portfolio class to implement this approach.
In [1]:
from dx import *
We start by instantiating a market environment object which in particular contains a list of ticker symbols in which we are interested in.
In [2]:
ma = market_environment('ma', dt.date(2010, 1, 1))
ma.add_list('symbols', ['AAPL', 'GOOG', 'MSFT', 'FB'])
ma.add_constant('source', 'google')
ma.add_constant('final date', dt.date(2014, 3, 1))
Using pandas under the hood, the class retrieves historial stock price data from either Yahoo! Finance of Google.
In [3]:
%%time
port = mean_variance_portfolio('am_tech_stocks', ma)
# instantiates the portfolio class
# and retrieves all the time series data needed
Since no portfolio weights have been provided, the class defaults to equal weights.
In [4]:
port.get_weights()
# defaults to equal weights
Out[4]:
Given these weights you can calculate the portfolio return via the method get_portfolio_return.
In [5]:
port.get_portfolio_return()
# expected (= historical mean) return
Out[5]:
Analogously, you can call get_portfolio_variance to get the historical portfolio variance.
In [6]:
port.get_portfolio_variance()
# expected (= historical) variance
Out[6]:
The class also has a neatly printable string representation.
In [7]:
print port
# ret. con. is "return contribution"
# given the mean return and the weight
# of the security
Via the method set_weights the weights of the single portfolio components can be adjusted.
In [8]:
port.set_weights([0.6, 0.2, 0.1, 0.1])
In [9]:
print port
You cal also easily check results for different weights with changing the attribute values of an object.
In [10]:
port.test_weights([0.6, 0.2, 0.1, 0.1])
# returns av. return + vol + Sharp ratio
# without setting new weights
Out[10]:
Let us implement a Monte Carlo simulation over potential portfolio weights.
In [11]:
# Monte Carlo simulation of portfolio compositions
rets = []
vols = []
for w in range(500):
weights = np.random.random(4)
weights /= sum(weights)
r, v, sr = port.test_weights(weights)
rets.append(r)
vols.append(v)
rets = np.array(rets)
vols = np.array(vols)
And the simulation results visualized.
In [12]:
import matplotlib.pyplot as plt
%matplotlib inline
plt.scatter(vols, rets, c=rets / vols, marker='o')
plt.grid(True)
plt.xlabel('expected volatility')
plt.ylabel('expected return')
plt.colorbar(label='Sharpe ratio')
Out[12]:
One of the major application areas of the mean-variance portfolio theory and therewith of this DX Analytics class it the optimization of the portfolio composition. Different target functions can be used to this end.
The first target function might be the portfolio return.
In [13]:
port.optimize('Return')
# maximizes expected return of portfolio
# no volatility constraint
In [14]:
print port
Instead of maximizing the portfolio return without any constraints, you can also set a (sensible/possible) maximum target volatility level as a constraint. Both, in an exact sense ("equality constraint") ...
In [15]:
port.optimize('Return', constraint=0.225, constraint_type='Exact')
# interpretes volatility constraint as equality
In [16]:
print port
... or just a an upper bound ("inequality constraint").
In [17]:
port.optimize('Return', constraint=0.4, constraint_type='Bound')
# interpretes volatility constraint as inequality (upper bound)
In [18]:
print port
The class also allows you to minimize portfolio risk.
In [19]:
port.optimize('Vol')
# minimizes expected volatility of portfolio
# no return constraint
In [20]:
print port
And, as before, to set constraints (in this case) for the target return level.
In [21]:
port.optimize('Vol', constraint=0.175, constraint_type='Exact')
# interpretes return constraint as equality
In [22]:
print port
In [23]:
port.optimize('Vol', constraint=0.20, constraint_type='Bound')
# interpretes return constraint as inequality (upper bound)
In [24]:
print port
Often, the target of the portfolio optimization efforts is the so called Sharpe ratio. The mean_variance_portfolio class of DX Analytics assumes a risk-free rate of zero in this context.
In [25]:
port.optimize('Sharpe')
# maximize Sharpe ratio
In [26]:
print port
Another application area is to derive the efficient frontier in the mean-variance space. These are all these portfolios for which there is no portfolio with both lower risk and higher return. The method get_efficient_frontier yields the desired results.
In [27]:
%%time
evols, erets = port.get_efficient_frontier(100)
# 100 points of the effient frontier
The plot with the random and efficient portfolios.
In [28]:
plt.scatter(vols, rets, c=rets / vols, marker='o')
plt.scatter(evols, erets, c=erets / evols, marker='x')
plt.grid(True)
plt.xlabel('expected volatility')
plt.ylabel('expected return')
plt.colorbar(label='Sharpe ratio')
Out[28]:
The capital market line is another key element of the mean-variance portfolio approach representing all those risk-return combinations (in mean-variance space) that are possible to form from a risk-less money market account and the market portfolio (or another appropriate substitute efficient portfolio).
In [29]:
%%time
cml, optv, optr = port.get_capital_market_line(riskless_asset=0.05)
# capital market line for effiecient frontier and risk-less short rate
In [30]:
cml # lambda function for capital market line
Out[30]:
The following plot illustrates that the capital market line has an ordinate value equal to the risk-free rate (the safe return of the money market account) and is tangent to the efficient frontier.
In [31]:
plt.figure(figsize=(8, 4))
plt.plot(evols, erets, lw=2.0, label='efficient frontier')
plt.plot((0, 0.4), (cml(0), cml(0.4)), lw=2.0, label='capital market line')
plt.plot(optv, optr, 'r*', markersize=10, label='optimal portfolio')
plt.legend(loc=0)
plt.grid(True)
plt.ylim(0)
plt.xlabel('expected volatility')
plt.ylabel('expected return')
Out[31]:
Portfolio return and risk of the efficient portfolio used are:
In [32]:
optr
Out[32]:
In [33]:
optv
Out[33]:
The portfolio composition can be derived as follows.
In [34]:
port.optimize('Vol', constraint=optr, constraint_type='Exact')
In [35]:
print port
Or also in this way.
In [36]:
port.optimize('Return', constraint=optv, constraint_type='Exact')
In [37]:
print port
As a larger, more realistic example, consider all symbols of the Dow Jones Industrial Average 30 index.
In [38]:
symbols = ['AXP', 'BA', 'CAT', 'CSCO', 'CVX', 'DD', 'DIS', 'GE',
'GS', 'HD', 'IBM', 'INTC', 'JNJ', 'JPM', 'KO', 'MCD', 'MMM',
'MRK', 'MSFT', 'NKE', 'PFE', 'PG', 'T', 'TRV', 'UNH', 'UTX',
'V', 'VZ','WMT', 'XOM']
# all DJIA 30 symbols
In [39]:
ma = market_environment('ma', dt.date(2010, 1, 1))
ma.add_list('symbols', symbols)
ma.add_constant('source', 'google')
ma.add_constant('final date', dt.date(2014, 3, 1))
Data retrieval in this case takes a bit.
In [40]:
%%time
djia = mean_variance_portfolio('djia', ma)
# defining the portfolio and retrieving the data
In [41]:
%%time
djia.optimize('Vol')
print djia.variance, djia.variance ** 0.5
# minimium variance & volatility in decimals
Given the larger data set now used, efficient frontier ...
In [42]:
%%time
evols, erets = djia.get_efficient_frontier(25)
# efficient frontier of DJIA
... and capital market line derivations take also longer.
In [43]:
%%time
cml, optv, optr = djia.get_capital_market_line(riskless_asset=0.01)
# capital market line and optimal (tangent) portfolio
In [44]:
plt.figure(figsize=(8, 4))
plt.plot(evols, erets, lw=2.0, label='efficient frontier')
plt.plot((0, 0.4), (cml(0), cml(0.4)), lw=2.0, label='capital market line')
plt.plot(optv, optr, 'r*', markersize=10, label='optimal portfolio')
plt.legend(loc=0)
plt.grid(True)
plt.ylim(0)
plt.xlabel('expected volatility')
plt.ylabel('expected return')
Out[44]:
Copyright, License & Disclaimer
© Dr. Yves J. Hilpisch | The Python Quants GmbH
DX Analytics (the "dx library") is licensed under the GNU Affero General Public License version 3 or later (see http://www.gnu.org/licenses/).
DX Analytics comes with no representations or warranties, to the extent permitted by applicable law.
http://www.pythonquants.com | analytics@pythonquants.com | http://twitter.com/dyjh
Python Quant Platform | http://quant-platform.com
Derivatives Analytics with Python (Wiley Finance) | http://derivatives-analytics-with-python.com
Python for Finance (O'Reilly) | http://shop.oreilly.com/product/0636920032441.do