Konvexita

Oceňování a analýza dluhopisů v Pythonu

Joshua Mayhew

Options Trader

Predikované vs. skutečné ceny – graf

import numpy as np
import numpy_financial as npf
import pandas as pd
import matplotlib.pyplot as plt

price = -npf.pv(rate=0.05, nper=20, pmt=5, fv=100)
price_up = -npf.pv(rate=0.06, nper=20, pmt=5, fv=100)
price_down = -npf.pv(rate=0.04, nper=20, pmt=5, fv=100)
duration = (price_down - price_up) / (2 * price * 0.01)
dollar_duration = duration * price * 0.01
print("Bond Price (USD): ", price)
print("Dollar Duration (USD): ", dollar_duration)

Bond Price (USD): 100.00
Dollar Duration (USD): 12.53

Predikované vs. skutečné ceny – graf

bond_yields = np.arange(0, 10, 0.1)
bond = pd.DataFrame(bond_yields, columns=['bond_yield'])
bond['price'] = -npf.pv(rate=bond['bond_yield'] / 100, nper=20, pmt=5, fv=100)

    bond_yield       price
0          0.0  200.000000
1          0.1  196.978503
..         ...         ...
98         9.8   58.570780
99         9.9   57.997210

[100 rows x 2 columns]

Predikované vs. skutečné ceny – graf

bond['yield_change'] = bond['bond_yield'] - 5

    bond_yield       price  yield_change
0          0.0  200.000000        -5.0
1          0.1  196.978503        -4.9
2          0.2  194.013231        -4.8
..         ...         ...           ...
97         9.7   59.153044         4.7
98         9.8   58.570780         4.8
99         9.9   57.997210         4.9

[100 rows x 3 columns]

Predikované vs. skutečné ceny – graf

$ \text{Price Change} = -100 \times \text{Dollar Duration} \times \Delta y$

bond['price_change'] = -100 * dollar_duration * bond['yield_change'] / 100

    bond_yield       price  yield_change  price_change
0          0.0  200.000000          -5.0     62.650619
1          0.1  196.978503          -4.9     61.397607
..         ...         ...           ...           ...
98         9.8   58.570780           4.8    -60.144594
99         9.9   57.997210           4.9    -61.397607

[100 rows x 4 columns]

Predikované vs. skutečné ceny – graf

bond['predicted_price'] = price + bond['price_change']

   bond_yield       price  yield_change  price_change  predicted_price
0          0.0  200.000000          -5.0     62.650619       162.650619
1          0.1  196.978503          -4.9     61.397607       161.397607
2          0.2  194.013231          -4.8     60.144594       160.144594
..         ...         ...           ...           ...              ...
97         9.7   59.153044           4.7    -58.891582        41.108418
98         9.8   58.570780           4.8    -60.144594        39.855406
99         9.9   57.997210           4.9    -61.397607        38.602393

[100 rows x 5 columns]

Predikované vs. skutečné ceny – graf

plt.plot(bond['bond_yield'], bond['price'])
plt.plot(bond['bond_yield'], bond['predicted_price'])
plt.xlabel('Yield (%)')
plt.ylabel('Price (USD)')
plt.title("Actual Bond Prices vs. 
    Predicted Prices Using Duration")
plt.legend(["Actual Price", "Predicted Price"])
plt.show()

Omezení durace

Durace je lineární míra
Ceny dluhopisů jsou konvexní
Durace je přesná pouze pro malé změny výnosu

Co je konvexita?

Měří zakřivení cen dluhopisů
Zlepšuje predikci cen a měření rizika
Vyšší konvexita = větší zakřivení vztahu ceny a výnosu

Vzorec konvexity

Použijeme zjednodušený vzorec pro konvexitu:

${\large Convexity = \frac{P_{down}\ +\ P_{up} \ - \ 2\times P}{P\ \times\ (\Delta y)^2}}$

$P_{down}$ = Cena dluhopisu při výnosu o 1 % nižším
$P_{up}$ = Cena dluhopisu při výnosu o 1 % vyšším
$2 \times P$ = Dvojnásobek ceny dluhopisu při aktuálním výnosu
$(\Delta y)^2$ = Druhá mocnina změny výnosu (použijeme 1 % ^ 2)

Příklad konvexity

10letý dluhopis, 5% roční kupón, 4% výnos do splatnosti – jaká je jeho konvexita?

${ Convexity = \frac{P_{down}\ +\ P_{up} \ - \ 2\times P}{P\ \times\ (\Delta y)^2}}$

price = -npf.pv(rate=0.04, nper=10, pmt=5, fv=100)

price_up = -npf.pv(rate=0.05, nper=10, pmt=5, fv=100)
price_down = -npf.pv(rate=0.03, nper=10, pmt=5, fv=100)

convexity = (price_down + price_up - 2 * price) / (price * 0.01 ** 2)
print("Convexity: ", convexity)

Convexity:  77.56

Shrnutí

Durace je lineární míra
Ceny dluhopisů jsou zakřivené
Durace je přesná pouze pro malé změny výnosu
Konvexita toto zakřivení měří

${ Convexity = \frac{P_{down}\ +\ P_{up} \ - \ 2\times P}{P\ \times\ (\Delta y)^2}}$

Pojďme si procvičit!

Oceňování a analýza dluhopisů v Pythonu