This section explains the usage of some basic framework classes and functions of DX Analytics. Mainly some helper functions, the discounting classes and the market environment class used to store market data and other parameters/data needed to model, value and risk manage derivative instruments.
In [1]:
from dx import *
In [2]:
np.set_printoptions(precision=3)
There are two helper functions used regulary:
get_year_deltas: get a list of year deltas (decimal fractions) relative to first value in time_listsn_random_numbers: get an array of standard normally distributed pseudo-random numbersSuppose we have a list object containing a number of datetime objects.
In [3]:
time_list = [dt.datetime(2015, 1, 1),
dt.datetime(2015, 4, 1),
dt.datetime(2015, 6, 15),
dt.datetime(2015, 10, 21)]
Passing this object to the get_year_deltas functions yields a list of year fractions representing the time intervals between the dates given. This is sometimes used e.g. for discounting purposes.
In [4]:
get_year_deltas(time_list)
Out[4]:
Monte Carlo simulation of course relies heavily an the use of random numbers. The function sn_random_numbers is a wrapper function around the pseudo-random number generator of the NumPy library. It implements antithetic variates and moment matching as generic variance reduction techniques. It also allows to fix the seed value for the random number generator. The shape parameter is a tuple object of three integers.
In [5]:
ran = sn_random_numbers((2, 3, 4), antithetic=True,
moment_matching=True, fixed_seed=False)
In [6]:
ran
Out[6]:
Using moment matching makes sure that the first and second moments match exactly 0 and 1, respectively.
In [7]:
ran.mean()
Out[7]:
In [8]:
ran.std()
Out[8]:
Setting the first value of the shape parameter to 1 yields a two-dimensional ndarray object.
In [9]:
ran = sn_random_numbers((1, 3, 4), antithetic=True,
moment_matching=True, fixed_seed=False)
In [10]:
ran
Out[10]:
In the risk-neutral valuation of derivative instrumente, discounting payoffs is a major task. The following discounting classes are implemented:
constant_short_rate: fixed short ratedeterministic_yield: deterministic yiels/term structureThe constant_short_rate class represents the most simple case for risk-neutral discounting. A discounting object is defined by instatiating the class and providing a name and a decimal short rate value only.
In [11]:
r = constant_short_rate('r', 0.05)
In [12]:
r.name
Out[12]:
In [13]:
r.short_rate
Out[13]:
The object has a method get_forward_rates to generate forward rates given, for instance, a list object of datetime objects.
In [14]:
r.get_forward_rates(time_list)
Out[14]:
Similarly, the method get_discount_factors returns discount factors for such a list object.
In [15]:
r.get_discount_factors(time_list)
Out[15]:
You can also pass, for instance, an ndarry object containing year fractions.
In [16]:
r.get_discount_factors(np.array([0., 1., 1.5, 2.]),
dtobjects=False)
Out[16]:
The deterministic_short_rate class allows to model an interest rate term structure. To this end, you need to pass a list object of datetime and yield pairs to the class.
In [17]:
yields = [(dt.datetime(2015, 1, 1), 0.02),
(dt.datetime(2015, 3, 1), 0.03),
(dt.datetime(2015, 10, 15), 0.035),
(dt.datetime(2015, 12, 31), 0.04)]
In [18]:
y = deterministic_short_rate('y', yields)
In [19]:
y.name
Out[19]:
In [20]:
y.yield_list
Out[20]:
The method get_interpolated_yields implements an interpolation of the yield data and returns the interpolated yields given a list object of datetime objects.
In [21]:
y.get_interpolated_yields(time_list)
Out[21]:
In similar fashion, the methods get_forward_rates and get_discount_factors return forward rates and discount factors, respcectively.
In [22]:
y.get_forward_rates(time_list)
Out[22]:
In [23]:
y.get_discount_factors(time_list)
Out[23]:
The market_environment class is used to collect relevant data for the modeling, valuation and risk management of single derivatives instruments and portfolios composed of such instruments. A market_environment object stores:
constants: e.g. maturity date of optionlists: e.g. list of datescurves: e.g. discounting objectsA market_environment object is instantiated by providing a name as a string object and the pricing date as a datetime object.
In [24]:
me = market_environment(name='me', pricing_date=dt.datetime(2014, 1, 1))
Constants are added via the add_constant method and providing a key and the value.
In [25]:
me.add_constant('initial_value', 100.)
In [26]:
me.add_constant('volatility', 0.25)
Lists of data are added via the add_list method.
In [27]:
me.add_list('dates', time_list)
The add_curve method does the same for curves.
In [28]:
me.add_curve('discount_curve_1', r)
In [29]:
me.add_curve('discount_curve_2', y)
The single data objects are stored in separate dictionary objects.
In [30]:
me.constants
Out[30]:
In [31]:
me.lists
Out[31]:
In [32]:
me.curves
Out[32]:
Data is retrieved from a market_environment object via the get_constant, get_list and get_curve methods and providing the respective key.
In [33]:
me.get_constant('volatility')
Out[33]:
In [34]:
me.get_list('dates')
Out[34]:
In [35]:
me.get_curve('discount_curve_1')
Out[35]:
Retrieving, for instance, a discounting object you can in one step retrieve it and call a method on it.
In [36]:
me.get_curve('discount_curve_2').get_discount_factors(time_list)
Out[36]:
Copyright, License & Disclaimer
© Dr. Yves J. Hilpisch | The Python Quants GmbH
DX Analytics (the "dx library" or "dx package") 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://tpq.io | dx@tpq.io | http://twitter.com/dyjh
Quant Platform | http://pqp.io
Python for Finance Training | http://training.tpq.io
Certificate in Computational Finance | http://compfinance.tpq.io
Derivatives Analytics with Python (Wiley Finance) | http://dawp.tpq.io
Python for Finance (2nd ed., O'Reilly) | http://py4fi.tpq.io