ファイナンス練習(2021年08月22日) - 「大人の教養・知識・気付き」を伸ばすブログ

　業務でC#を用いることになったので、最近勉強していなくて朧気になってきた知識をReviseする意味でも、以下の書籍を読みながらC#で実装してみる。今日はP.16まで。

2. グリークス

2. グリークス

　オプションは、理論価格があるために、入力パラメータが明示的にわかっている。あらためて挙げるとこのとおり：

原資産価格： $S$
行使価格： $K$
（初期時点を $0$ としたときの）将来時点(ここでは単位:年)： $T$
リスクフリーレート： $r$
原資産価格のボラティリティ（リスク・標準偏差）： $\sigma$

このため、この5つの変化がオプション価格に影響を与える。そこでリスク（価格変化）を管理するには、それらが変化したときにオプション価格がどの程度動くかを評価する。各パラメータに対する価格変化の指標をまとめてグリークス（Greeks）という。
　あるパラメータ $x$ （他のパラメータは固定）に対してオプション価格 $c$ を $x=x_{0}$ のまわりでTaylor展開してみると、 $\xi\in\mathbb{R}$ として

$\displaystyle{c(x)=c(x_{0})+\frac{c^{\prime}(x_{0})}{1!}(x-x_{0})+\cdots+\frac{c^{(n)}(\xi)}{n!}(x-x_{0})^n}$

が成り立つ。 $x-x_{0}$ がそこまで大きくなければ、1次項または2次項（1次微分もしくは2次微分）を評価すればリスク管理としては充分と言える。つまり、

$\displaystyle{\frac{\partial c}{\partial x}, \frac{\partial^{2}c}{\partial x^{2}}}$

に注目することになる。

2.1 デルタ

　原資産価格 $S$ の変化はもっとも重要なパラメータであるから、原資産価格の動きに対するオプション価格の変化がもっとも重要なリスク管理尺度である。
　原資産価格の変化に対するオプション価格の変化をデルタという。

$\begin{aligned} \Delta_{c}&=\displaystyle{\frac{\partial c}{\partial S}}\,\\ \Delta_{p}&=\displaystyle{\frac{\partial p}{\partial S}}=\Delta_{c}-1\ \end{aligned}$

なおプット・オプションのデルタは、プット・コール・パリティを用いれば*1

$\begin{aligned} \Delta_{p}&=\displaystyle{\frac{\partial}{\partial S}(c-S+Ke^{-rT})}\,\\ &=\Delta_{c}-1 \end{aligned}$

となることから分かる。
　さらに展開してみよう。

$\begin{aligned} \Delta_{c}&=\displaystyle{\frac{\partial c}{\partial S}}\,\\ &=\displaystyle{\frac{\partial}{\partial S}(S\Phi (d)-Ke^{-rT}\Phi(d-\sigma \sqrt{T}))}\,\\ &=\Phi(d)+S\phi(d)\frac{\partial d}{\partial S}-\frac{\partial}{\partial S}(Ke^{-rT}\Phi(d-\sigma\sqrt{T}))\ \end{aligned}$

である。ここに

$\begin{aligned} \frac{\partial d}{\partial S}&=\frac{\partial}{\partial S}(\frac{\log{(\frac{S}{K})}+(r+\frac{1}{2}\sigma_{2})T}{\sigma\sqrt{T}})\\ &=\frac{\frac{K}{S}\cdot\frac{1}{K}}{\sigma\sqrt{T}}\\ &=\frac{1}{S\sigma\sqrt{T}},\\ \frac{\partial (d-\sigma\sqrt{T})}{\partial S}&=\frac{1}{S\sigma\sqrt{T}} \end{aligned}$

を代入することで

$\begin{aligned} \frac{\partial c}{\partial S}&=\Phi(d)+S\phi(d)\frac{\partial d}{\partial S}-\frac{\partial}{\partial S}(Ke^{-rT}\Phi(d-\sigma\sqrt{T}))\\\ &=\Phi(d)+S\phi(d)\frac{1}{S\sigma\sqrt{T}}-\frac{1}{S\sigma\sqrt{T}}Ke^{-rT}\phi(d-\sigma\sqrt{T}) \end{aligned}$

さらに

$\begin{aligned} \phi(d-\sigma\sqrt{T})&=\frac{d}{dd}\Phi(d-\sigma\sqrt{T})\\ &=\frac{1}{\sqrt{2\pi}}\exp(-\frac{1}{2}(d-\sigma\sqrt{T})^{2})\\ &=\frac{1}{\sqrt{2\pi}}\exp{(-\frac{1}{2}d^{2})}\exp{(d\sigma\sqrt{T})}\exp(-\frac{1}{2}\sigma^{2}T) \end{aligned}$

に $d$ の定義から得られる

$\begin{aligned} \frac{1}{2}\sigma^{2}T=d\sigma\sqrt{T}-\log(\frac{S}{K})-rT \end{aligned}$

を代入して得た、

$\begin{aligned} \phi(d-\sigma\sqrt{T})&=\frac{1}{\sqrt{2\pi}}\exp{(-\frac{1}{2}d^{2})}\exp{(d\sigma\sqrt{T})}\exp(-\frac{1}{2}\sigma^{2}T)\\ &=\frac{1}{\sqrt{2\pi}}\exp{(-\frac{1}{2}d^{2})}\exp{\log{(S/K)}}\exp(rT)\\ &=\phi(d)(\frac{S}{K})e^{rT} \end{aligned}$

を代入することで

$\begin{aligned} \frac{\partial c}{\partial S}&=\Phi(d)+S\phi(d)\frac{1}{S\sigma\sqrt{T}}-\frac{1}{S\sigma\sqrt{T}}Ke^{-rT}\phi(d)(\frac{S}{K})e^{rT}\\ &=\Phi(d)+\frac{1}{\sigma\sqrt{T}}\phi(d)-\frac{1}{\sigma\sqrt{T}}\phi(d)\\ &=\Phi(d) \end{aligned}$

を得る。
　なおプット・オプションについては

$\begin{aligned} \Delta_{p}&=\Delta_{c}-1\\ &=\Phi(d)-1\\ &=-\Phi(-d) \end{aligned}$

まとめよう：

コール・オプションのデルタ $\Delta_{c}=\Phi(d)$

プット・オプションのデルタ $\Delta_{p}=-\Phi(-d)$

2.2 ガンマ

　原資産価格 $S$ の変化がもっとも重要であることは既に述べた。それだけでなく原資産価格 $S$ はそれなりに大きく動くことがあるため、その変化による価格変化としてデルタのみで捉えようとすると誤差が大きくなってしまう可能性がある。そこで二次微分まで考える。
　原資産価格 $S$ の変化に対するデルタの変化、すなわち原資産価格 $S$ の変化に対するオプション価格の二次微分としてガンマを以下で定義する：

$\begin{aligned} \Gamma_{c}&=\displaystyle{\frac{\partial \Delta_{c}}{\partial S}}=\displaystyle{\frac{\partial^{2}c}{\partial^{2} S}}\\ \Gamma_{p}&=\displaystyle{\frac{\partial \Delta_{c}}{\partial S}}=\displaystyle{\frac{\partial^{2}p}{\partial^{2} S}} \end{aligned}$

　これを展開してみよう。 $\phi(d)$ の対称性、すなわち $\phi(d)=\displaystyle{\frac{1}{\sqrt{2\pi}}e^{\frac{1}{2}d^{2}}}=\phi(-d)$ に注意すれば

$\begin{aligned} \Gamma_{c}&=\displaystyle{\frac{\partial \Delta_{c}}{\partial S}}\\ &=\frac{\partial \Phi(d)}{\partial S}\\ &=\phi(d)\cdot\frac{\partial d}{\partial S}\\ &=\frac{1}{S\sigma\sqrt{T}}\phi(d)\\ \Gamma_{p}&=\displaystyle{\frac{\partial \Delta_{p}}{\partial S}}\\ &=-\frac{\partial \Phi(-d)}{\partial S}\\ &=-\phi(-d)\cdot\frac{\partial -d}{\partial S}\\ &=\frac{1}{S\sigma\sqrt{T}}\phi(-d)\\ &=\frac{1}{S\sigma\sqrt{T}}\phi(d)=\Gamma_{c} \end{aligned}$

が得られる。

まとめよう：

コール・オプションのガンマ $\Gamma_{c}=\frac{1}{S\sigma\sqrt{T}}\phi(d)$

プット・オプションのガンマ $\Gamma_{p}=\frac{1}{S\sigma\sqrt{T}}\phi(d)$

2.3 C#での実装

　コール（プット）・オプションのデルタを実装した。なお理論値と数値微分による数値解の2つを実装した。またデルタの数値微分は前回実装したオプション価格を出力するメソッド*2を数値的に二次微分している。

using System;
using MathNet.Numerics;

/// <summary>
/// デルタの理論値を返す。
/// </summary>
/// <param name="S">原資産価格</param>
/// <param name="K">行使価格</param>
/// <param name="r">リスクフリーレート</param>
/// <param name="sigma">原資産価格のボラティリティ</param>
/// <param name="T">満期までの期間</param>
/// <param name="CallPutFlg">コール・プットを表すフラグ。0ならばコール・オプション、1ならばプット・オプションの価格を返す。</param>
/// <returns>コール・オプションまたはプット・オプションのデルタ。</returns>
public static double Delta(double S, double K, double r, double sigma, double T, int CallPutFlg)
{
    double d = Math.Log(S/K) + (r + 0.5 * Math.Pow(sigma, 2.0)) * T;
    d /= (sigma * Math.Sqrt(T));

    double Sign = -2.0 * (double)CallPutFlg + 1.0;

    double delta = Sign * MathNet.Numerics.Distributions.Normal.CDF(0.0, 1.0, Sign * d);

    return delta;
}

/// <summary>
/// デルタの数値微分による値を返す。
/// </summary>
/// <param name="S">原資産価格</param>
/// <param name="K">行使価格</param>
/// <param name="r">リスクフリーレート</param>
/// <param name="sigma">原資産価格のボラティリティ</param>
/// <param name="T">満期までの期間</param>
/// <param name="CallPutFlg">コール・プットを表すフラグ。0ならばコール・オプション、1ならばプット・オプションの価格を返す。</param>
/// <returns>コール・オプションまたはプット・オプションのデルタ。</returns>
public static double DeltaByNumericalDifferential(double S, double K, double r, double sigma, double T, int CallPutFlg)
{
    double Bin = 0.00001;

    double delta = OptionPricing(S + Bin,K,r,sigma,T,CallPutFlg);
    delta -= OptionPricing(S - Bin,K,r,sigma,T,CallPutFlg);
    delta /= (2.0 * Bin);

    return delta;
}

/// <summary>
/// ガンマの理論値を返す。
/// </summary>
/// <param name="S">原資産価格</param>
/// <param name="K">行使価格</param>
/// <param name="r">リスクフリーレート</param>
/// <param name="sigma">原資産価格のボラティリティ</param>
/// <param name="T">満期までの期間</param>
/// <param name="CallPutFlg">コール・プットを表すフラグ。0ならばコール・オプション、1ならばプット・オプションの価格を返す。</param>
/// <returns>コール・オプションまたはプット・オプションのガンマ。</returns>
public static double Gamma(double S, double K, double r, double sigma, double T, int CallPutFlg)
{
    double d = Math.Log(S / K) + (r + 0.5 * Math.Pow(sigma, 2.0)) * T;
    d /= (sigma * Math.Sqrt(T));

    double Gamma = Math.Exp(-0.5 * Math.Pow(d, 2.0));
    Gamma /= (S * sigma * Math.Sqrt(T) * Math.Sqrt(2 *Math.PI));

    return Gamma;
}

/// <summary>
/// ガンマの数値微分による値を返す。
/// </summary>
/// <param name="S">原資産価格</param>
/// <param name="K">行使価格</param>
/// <param name="r">リスクフリーレート</param>
/// <param name="sigma">原資産価格のボラティリティ</param>
/// <param name="T">満期までの期間</param>
/// <param name="CallPutFlg">コール・プットを表すフラグ。0ならばコール・オプション、1ならばプット・オプションの価格を返す。</param>
/// <returns>コール・オプションまたはプット・オプションのガンマ。</returns>
public static double GammaByNumericalDifferential(double S, double K, double r, double sigma, double T, int CallPutFlg)
{
    double Bin = 0.001;

    double P_BinPlus = OptionPricing(S + Bin, K, r, sigma, T, CallPutFlg);
    double P = OptionPricing(S, K, r, sigma, T, CallPutFlg);
    double P_BinMinus = OptionPricing(S - Bin, K, r, sigma, T, CallPutFlg);

    double Gamma = P_BinPlus - 2.0 * P + P_BinMinus;
    Gamma /= Math.Pow(Bin, 2.0);

    return Gamma;
}

*1:https://power-of-awareness.com/entry/2021/08/21/010000参照。

*2:同上参照。