Valuation of Quanto Caps and Floors in a Calibrated Multi-Curve Cross-Currency LIBOR Market Model

Show more

1. Introduction

Interest rate modeling has been a major interest amongst researchers. This is mostly because the interest rate markets have grown to dominate the financial world due to its vast number of traded financial products flooding the markets. According to [1], interest rate products form the largest portion of traded instruments in the financial markets. Initially, interest rates were modeled using the short rate models: [2] [3] [4] [5] [6] and many more. Later on, they were extended to model the instantaneous forward rates (see [7] ) which were not directly observable. In the late nineties, a new class of models better known as the market models were introduced into the financial markets by authors such as [8] [9] [10]. These models were quickly accepted as they made use of market observable rates such as LIBORs and swap rates. Examples of such market models are the LIBOR market model (LMM) and the swap market model (SMM).

The term LIBOR stands for the London Interbank Offered rate. It is the benchmark rate at which major banks borrow from each other on a short term basis. This rate is monitored by ICE and is published daily for five different currencies and seven maturities. LIBOR is used widely in financial markets as the underlying in standard interbank products, commercial products, hybrid products and consumer related products. LIBOR equivalents also accepted globally are the European Interbank Offered Rate (EURIBOR), Tokyo Interbank Offered Rate (TIBOR), Shanghai Interbank Offered Rate (SHIBOR), and Mumbai Interbank Offered Rate (MIBOR).

Before 2007, the spread between the different LIBOR tenors was almost negligible. Hence, a single curve was sufficient for both discounting and in generating future cashflows. However, after the 2007 financial crisis, [11] noticed that this assumption could no longer hold as the spreads started becoming larger and larger to an extent that they could no longer be considered as negligible. Giving rise to the birth of multi-curve LIBOR market models (see Figure 1 and Figure 2 plotted using data collected from [12] and [13] respectively).

Figure 1. USD LIBOR.

Figure 2. GBP LIBOR.

Under the multi-curve framework, one curve is used to generate future cashflows while the other one is used to discount the generated future cashflows. So far, many models have already been proposed in practice (see [11] [14] - [27] and so on).

Our main interest in this paper is in pricing quanto caps and floors using the multi-curve cross-currency LIBOR market model (MCCCLMM) dynamics introduced in [27]. The model parameters are first calibrated to exactly match the market observable cap prices. The estimated parameters from the calibration process are then used to essentially price quanto caps and floors in a Black MCCCLMM quanto caplet or floorlet formula derived in this paper. The quanto cap and floor prices are also estimated using Monte Carlo simulations and a comparison between the two models is done. The discretization scheme used to discretize the forward LIBORs is the Euler scheme.

However, the calibration problem, over the decades, has not been an easy one. It has seen researchers resort to both parametric and non-parametric techniques of calibrating the LIBOR market model (LMM). For instance, [28] developed a fast at the money (ATM) calibration of the LMM using Lagrange multipliers. He calibrated his model using ATM caps, swaptions and historical correlations. [29] came up with a numerical technique for calibrating financial models that essentially solves an inverse problem associated with some partial differential equations. [30] calibrated the LMM using cap and swaption price data collected on 16^{th} May 2010. They calibrated the cap volatilities using the Separable piecewise constant (SPC) parameterization technique and Linear-Exponential (L-E) formulation both under the general piecewise constant assumption. They estimated their model correlations using swaptions. [31] explains the different types of assumptions that can be made on the general piecewise constant technique for cap volatilities. [32] explains the concept behind the general piecewise constant. He goes further ahead to explain that the main problem in calibrating the LMM is in finding a volatility function that accurately reproduces a sample of market derivative prices e.g. cap and swaption prices.

In this paper, the instantaneous correlations are estimated from historical rates. The general piecewise constant assumption was assumed on the cap volatilities. The foreign exchange rate volatility was also estimated from historical data.

2. Materials and Methods

In this section, the relevant tools, models, methods and tests used are presented.

2.1. Data

The datasets considered in this study were: Historical overnight and 6 month USD and GBP LIBORs ranging from 2/1/2008 to 2/1/2018 collected from [12] and [13] respectively. The GBP/USD foreign exchange rate data ranging from 2/1/2008 to 2/1/2018 collected from [33]. UK gilts and US Treasury rates collected from [34] and [35] respectively on 2/1/2018. Finally, the GBP and USD cap prices collected from [36].

2.2. Data Analysis Tool

R open software version 3.1.2 was used in simulating and analysing all the data in this study. Useful packages considered were “MASS”, “sde”, and “lmtest”.

2.3. Model Notations

Let $X\left(t\right)$ be the spot foreign exchange rate at time t quoted as the ratio of units of United States Dollar (USD) to one unit of Pound Sterling (GBP). d be the US (domestic) economy and f be the UK (foreign) economy. D and L denote the risk-free and risky curves respectively. ${r}_{d}^{D}$ and ${r}_{f}^{D}$ be the risk-free short rates of interest associated with the US and UK economies respectively. In addition to this, ${P}_{d}^{D}\left(t\mathrm{,}T\right)$ and ${P}_{f}^{D}\left(t\mathrm{,}T\right)$ denote the risk-free zero coupon bonds in the domestic and foreign economies respectively. ${L}_{d}^{D}\left(t\mathrm{,}T\right)$ and ${L}_{f}^{D}\left(t\mathrm{,}T\right)$ denote the simply compounded overnight LIBORs associated with US and UK markets respectively. Last but not least, ${L}_{d}\left(t\mathrm{,}T\right)$ and ${L}_{f}\left(t\mathrm{,}T\right)$ denote the 6 month LIBOR associated with the US and UK markets.

2.4. The Multi-Curve Cross-Currency LIBOR Market Model (MCCCLMM)

The stochastic differential equations associated with the MCCCLMM dynamics that were considered in this study under the spot domestic risk neutral measure were given by:

$\{\begin{array}{l}\text{d}{P}_{d}^{D}\left(t,T\right)={r}_{d}^{D}\left(t\right){P}_{d}^{D}\left(t,T\right)\text{d}t-{\sigma}_{dD}\left(t,T\right){P}_{d}^{D}\left(t,T\right)\text{d}{W}_{dD}\left(t\right)\\ \text{d}{L}_{d}\left(t,T\right)={\rho}_{dD}^{dL}\ast {\sigma}_{dD}\left(t,T+\tau \right){\gamma}_{dL}\left(t,T\right){L}_{d}\left(t,T\right)\text{d}t\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+{\gamma}_{dL}\left(t,T\right){L}_{d}\left(t,T\right)\text{d}{W}_{dL}\left(t\right)\\ \text{d}{L}_{f}\left(t,T\right)=\left({\rho}_{fD}^{fL}{\sigma}_{fD}(t,T+\tau )-{\rho}_{xx}^{fD}{\sigma}_{X}\left(t\right)\right){\gamma}_{fL}\left(t,T\right){L}_{f}\left(t,T\right)\text{d}t\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}+{\gamma}_{fL}\left(t,T\right){L}_{f}\left(t,T\right)\text{d}{W}_{fL}\left(t\right)\\ \text{d}X\left(t\right)=\left({r}_{d}^{D}\left(t\right)-{r}_{f}^{D}\left(t\right)\right)X\left(t\right)\text{d}t+{\sigma}_{X}\left(t\right)\text{d}{W}_{XX}\left(t\right)\end{array}$ (1)

where

${\sigma}_{kD}\left(t,T\right)={\displaystyle \underset{j=1}{\overset{\lfloor {\tau}^{-1}\left(T-t\right)\rfloor}{\sum}}}\frac{\tau {L}_{k}^{D}\left(t,T-j\tau \right)}{1+\tau {L}_{k}^{D}\left(t,T-j\tau \right)}{\gamma}_{kD}\left(t,T-j\tau \right)$ (2)

${\sigma}_{kL}\left(t,T\right)={\displaystyle \underset{j=1}{\overset{\lfloor {\tau}^{-1}\left(T-t\right)\rfloor}{\sum}}}\frac{\tau {L}_{k}\left(t,T-j\tau \right)}{1+\tau {L}_{k}\left(t,T-j\tau \right)}{\gamma}_{kL}\left(t,T-j\tau \right)$ (3)

$\lfloor {\tau}^{-1}\left(T-t\right)\rfloor $ denotes the greatest integer that is less than ${\tau}^{-1}\left(T-t\right)$.

And

$W=\left[\begin{array}{c}{W}_{dD}\\ {W}_{dL}\\ {W}_{fD}\\ {W}_{fL}\\ {W}_{XX}\end{array}\right]$ (4)

is a correlated Wiener process with a correlation matrix given by:

$\Sigma =\left(\begin{array}{ccccc}1& {\rho}_{dD}^{dL}& {\rho}_{dD}^{fD}& {\rho}_{dD}^{fL}& {\rho}_{xx}^{dD}\\ {\rho}_{dD}^{dL}& 1& {\rho}_{dL}^{fD}& {\rho}_{dL}^{fL}& {\rho}_{xx}^{dL}\\ {\rho}_{dD}^{fD}& {\rho}_{dL}^{fD}& 1& {\rho}_{fD}^{fL}& {\rho}_{xx}^{fD}\\ {\rho}_{dD}^{fL}& {\rho}_{dL}^{fL}& {\rho}_{fD}^{fL}& 1& {\rho}_{xx}^{fL}\\ {\rho}_{xx}^{dD}& {\rho}_{xx}^{dL}& {\rho}_{xx}^{fD}& {\rho}_{xx}^{fL}& 1\end{array}\right)$ (5)

where it is assumed that the Wiener processes are governed by constant correlation factors such that

${\rho}_{i}^{j}=(\begin{array}{l}{\rho}_{j}^{i}\text{\hspace{1em}}\text{if}\text{\hspace{0.17em}}j\ne k;\text{\hspace{1em}}j,k\in \left\{dD,dL,fD,fL,xx\right\},\text{\hspace{0.17em}}\left|{\rho}_{j}^{k}\right|\le 1\\ 1\text{\hspace{1em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{if}\text{\hspace{0.17em}}j=k\end{array}$ (6)

2.5. Calibration of MCCCLMM Parameters to Market Data

Calibration is the process of estimating the model parameters such that they match the market prices.

2.5.1. Estimation of MCCCLMM Correlations

The constant correlations ${\rho}_{i}^{j}$ can be estimated using the Pearson’s correlation formula given by:

$r=\frac{{\displaystyle {\sum}_{i=1}^{n}\left({X}_{i}-\stackrel{\xaf}{X}\right)\left({Y}_{i}-\stackrel{\xaf}{Y}\right)}}{\sqrt{{\displaystyle {\sum}_{i=1}^{n}{\left({X}_{i}-\stackrel{\xaf}{X}\right)}^{2}}{\displaystyle {\sum}_{i=1}^{n}{\left({Y}_{i}-\stackrel{\xaf}{Y}\right)}^{2}}}}$ (7)

2.5.2. Estimation of Foreign Exchange Rate Volatility from Historical Rates

Given that the SDE of the foreign exchange rate is given by:

$\frac{\text{d}{X}_{t}}{{X}_{t}}={\mu}_{X}\text{d}t+{\sigma}_{X}\text{d}{W}_{xx}(t)$

then the SDE can be discretised as follows:

$\begin{array}{c}\frac{\text{d}{X}_{t}}{{X}_{t}}=\text{d}\mathrm{ln}{X}_{t}\\ =\mathrm{ln}\left({X}_{t}\right)-\mathrm{ln}\left({X}_{t-1}\right)\\ =\mathrm{ln}\left(\frac{{X}_{t}}{{X}_{t-1}}\right)\\ ={\mu}_{X}\Delta t+{\sigma}_{X}\Delta {W}_{xx}\left(t\right)\end{array}$ (8)

and according to [37], for one to estimate the volatility parameter from historical rates, the following steps should be considered:

1) Calculate the logarithmic returns as follows:

${l}_{i}=\mathrm{ln}\left(\frac{{X}_{i}}{{X}_{i-1}}\right);\text{\hspace{0.17em}}\text{\hspace{0.17em}}i\in \left\{1,2,\cdots ,n\right\}$ (9)

2) Calculate the mean of the logarithmic returns as follows:

$\stackrel{\xaf}{l}={\displaystyle \underset{i=1}{\overset{n}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{l}_{i}$ (10)

3) Calculate the unbiased standard deviation of the logarithmic returns as follows:

${s}_{l}=\frac{1}{n-1}{\displaystyle \underset{i=1}{\overset{n}{\sum}}}{\left({l}_{i}-\stackrel{\xaf}{l}\right)}^{2}$ (11)

and now ${\sigma}_{X}$ can be estimated as

${\stackrel{^}{\sigma}}_{X}=\frac{{s}_{l}}{\Delta t}$ (12)

2.5.3. Parameterization of the Forward LIBOR Volatility

In this paper, the constant maturity-dependent volatility assumption under the general piecewise constant technique is considered. [31] explains in detail the possible assumptions that can be considered under this technique. It is assumed that

${\gamma}_{i}\left(t\right)={S}_{i};\text{\hspace{0.17em}}\text{\hspace{0.17em}}0<t\le {T}_{i-1}$ (13)

Calibrating the model to caplet amounts is equivalent to choosing deterministic LIBOR volatilities of forward rates ${\gamma}_{1}\mathrm{,}{\gamma}_{2}\mathrm{,}\cdots \mathrm{,}{\gamma}_{m}$ such that:

${v}_{{T}_{i}-cap}^{2}=\frac{1}{{T}_{i-1}}{\displaystyle {\int}_{0}^{{T}_{i-1}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\gamma}_{i}^{2}\left(t\right)\text{d}t;\text{\hspace{0.17em}}\text{\hspace{0.17em}}i=\left\{1,2,\cdots ,m\right\}$ (14)

The $\gamma $ values are as summarised in Table 1.

Under this assumption,

Table 1. Constant maturity volatility assumption.

$\begin{array}{l}{T}_{i-1}{v}_{{T}_{i}-cap}^{2}={T}_{i-1}{S}_{i}^{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}i\in \left\{1,2,\cdots ,m\right\}\\ {v}_{{T}_{i}-cap}^{2}={S}_{i}^{2}\end{array}$ (15)

The advantage of using this assumption is that the S parameters fit the market cap volatilities.

2.5.4. Calibration via the Black MCCCLMM Formula

It was assumed that in the U.S economy, the USD LIBOR is a domestic rate. In the same way, in the British economy, the GBP LIBOR is also a domestic rate. Hence under the domestic risk neutral measure, the MCCCLMM dynamics of the risky USD or GBP LIBOR in either domestic economies is given by:

$\frac{\text{d}{L}_{k}\left(t,T\right)}{{L}_{k}\left(t,T\right)}={\rho}_{kD}^{kL}\ast {\sigma}_{kD}\ast {\gamma}_{kL}\ast \text{d}t+{\gamma}_{kL}\ast \text{d}{W}_{kL}\left(t\right);\text{\hspace{0.17em}}\text{\hspace{0.17em}}k\in \left\{d,f\right\}$ (16)

The Black-like formula for calculating USD or GBP caplets or floorlets in the respective domestic currencies is given by:

${C}_{k}\left(t\right)=\omega \tau {P}_{k}^{D}\left(t,T\right)\left({L}_{k}\left(t,T\right)\ast {\text{e}}^{{\rho}_{kD}^{kL}\ast {\sigma}_{kD}\ast {\gamma}_{kL}\left(T-t\right)}\ast \Phi \left(\omega {d}_{11}\right)-{K}_{k}\ast \Phi \left(\omega {d}_{21}\right)\right)$ (17)

where

$\begin{array}{l}\omega =(\begin{array}{l}+1\text{\hspace{1em}}\text{if it is a caplet contract}\\ -1\text{\hspace{1em}}\text{if it is a floorlet contrac}\end{array}\\ {d}_{11}=\frac{\mathrm{ln}\left(\frac{{L}_{k}\left(T,T\right)}{{K}_{k}}\right)+\left({\rho}_{kD}^{kL}\ast {\sigma}_{kD}\ast {\gamma}_{kL}+0.5{\gamma}_{kL}^{2}\right)\left(T-t\right)}{{\gamma}_{kL}\ast \sqrt{T-t}}\\ {d}_{21}={d}_{11}-{\gamma}_{kL}\ast \sqrt{T-t}\end{array}$ (18)

The MCCCLMM parameters fitted in the Black-like formula described in Equation (16), were minimized to exactly match the cap prices collected from [36].

2.6. Yield Curve Bootstrapping

Bootstrapping in finance is the process of constructing a zero coupon yield curve from a set of coupon bearing instruments by filling in the missing yields. In this research, the Nelson Siegel Svenssons (NSS) method is used. The NSS method given by

$y\left(t\right)={\beta}_{0}+{\beta}_{1}\left(\frac{1-{\text{e}}^{\frac{-t}{{\lambda}_{1}}}}{\frac{t}{{\lambda}_{1}}}\right)+{\beta}_{2}\left(\frac{1-{\text{e}}^{\frac{-t}{{\lambda}_{1}}}}{\frac{t}{{\lambda}_{1}}}-{\text{e}}^{\frac{-t}{{\lambda}_{1}}}\right)+{\beta}_{3}\left(\frac{1-{\text{e}}^{\frac{-t}{{\lambda}_{2}}}}{\frac{t}{{\lambda}_{2}}}-{\text{e}}^{\frac{-t}{{\lambda}_{2}}}\right)$ (19)

is used to bootstrap the yield curve from the USD treasury and UK gilt rates. The NSS model parameters ${\beta}_{0}\mathrm{,}{\beta}_{1}\mathrm{,}{\beta}_{2}\mathrm{,}{\beta}_{3}\mathrm{,}{\lambda}_{1}$ and ${\lambda}_{2}$ are estimated by minimizing the sum of squared errors (SSE):

$\mathrm{min}{\displaystyle \sum}{\left({y}^{A}-y\right)}^{2}$ (20)

where ${y}^{A}$ are the actual market rates and y are the rates estimated via the NSS method. The zero coupon bond prices are then estimated using the formula:

$Z\left(t\right)={\left(1+y\left(t\right)\right)}^{-t}$ (21)

2.7. Mean Error Analysis

The option prices estimated using the Black’s formula were compared with the mean Monte Carlo simulated option prices using the mean absolute percentage error (MAPE) given by:

$\text{MAPE}=\frac{1}{n}{\displaystyle \underset{i=1}{\overset{n}{\sum}}}\frac{\left|B{P}_{i}-M{P}_{i}\right|}{B{P}_{i}}\ast 100\%$ (22)

where $B{P}_{i}$ is Black price and $M{P}_{i}$ is the Monte Carlo price.

According to [38], if $\text{MAPE}\le \mathrm{10\%}$, then the model is considered to be highly accurate. If $10\%<\text{MAPE}\le 20\%$, then the model is considered to be a good model. If $20\%<\text{MAPE}\le 50\%$, then the model is considered to be reasonable, and if $\text{MAPE}>50\%$ then the model is taken to be inaccurate.

2.8. MCCCLMM Simulation

In this section, a brief introduction of the discretization scheme used in our simulations is done, how the correlated Wiener processes were generated, and finally how the cap and floor prices were simulated is given.

2.8.1. Discretization Scheme

The Euler discretization scheme [39] was used to discretize our calibrated model dynamics so as to enable us to simulate the risky forward LIBORs in discrete time.

Consider an Itô process ${\left\{{X}_{t}\right\}}_{0\le t\le T}$ with a stochastic differential Equation (SDE) given by:

$\text{d}{X}_{t}=\mu \left(t,{X}_{t}\right)\text{d}t+\sigma \left(t,{X}_{t}\right)\text{d}{W}_{t}$ (23)

and an initial deterministic value of ${X}_{0}={x}_{0}$. Then the Euler approximation of ${X}_{t}$ for the interval $t\in \left[\mathrm{0,}T\right]$ ; is a process Y such that $0={t}_{0}<{t}_{1}<\cdots <{t}_{N}=T$ satisfies the iterative scheme:

${Y}_{i+1}={Y}_{i}+\mu \left({t}_{i},{Y}_{i}\right)\left({t}_{i+1}-{t}_{i}\right)+\sigma \left({t}_{i},{Y}_{i}\right)\left({W}_{i+1}-{W}_{i}\right);\text{\hspace{0.17em}}\text{\hspace{0.17em}}i=0,1,2,\cdots ,N-1$ (24)

where ${Y}_{0}={X}_{0}$.

2.8.2. Generation of Correlated Wiener Processes

According to [40], the correlated Wiener process, W can be simulated by applying Cholesky decomposition as follows:

$W=AZ$ (25)

In our case, this is given by:

$\left(\begin{array}{c}\text{d}{W}_{dD}\\ \text{d}{W}_{dL}\\ \text{d}{W}_{fD}\\ \text{d}{W}_{fL}\\ \text{d}{W}_{xx}\end{array}\right)=\left(\begin{array}{ccccc}1& 0& 0& 0& 0\\ a& b& 0& 0& 0\\ c& d& e& 0& 0\\ f& g& h& i& 0\\ j& k& l& m& n\end{array}\right)\left(\begin{array}{c}\text{d}{Z}_{1}\\ \text{d}{Z}_{2}\\ \text{d}{Z}_{3}\\ \text{d}{Z}_{4}\\ \text{d}{Z}_{5}\end{array}\right)$ (26)

where:

${Z}_{i}\mathrm{;}\text{\hspace{0.17em}}i\in \left\{\mathrm{1,2,3,4,5}\right\}$ are independent standard normal variables.

$a={\rho}_{dD}^{dL};\text{\hspace{1em}}b=\sqrt{1-{a}^{2}};\text{\hspace{1em}}c={\rho}_{dD}^{fD}$

$d=\frac{{\rho}_{dD}^{fD}-ac}{b};\text{\hspace{1em}}e=\sqrt{1-{c}^{2}-{d}^{2}};\text{\hspace{1em}}f=\frac{{\rho}_{dD}^{fL}-ac}{b}$

$g=\frac{{\rho}_{dL}^{fL}-af}{b};\text{\hspace{1em}}h=\frac{{\rho}_{fD}^{fL}-cf-dg}{e};\text{\hspace{1em}}i=\sqrt{1-{f}^{2}-{g}^{2}-{h}^{2}}$

$j={\rho}_{xx}^{dD};\text{\hspace{1em}}k=\frac{{\rho}_{xx}^{dL}-aj}{b};\text{\hspace{1em}}l=\frac{{\rho}_{xx}^{fD}-cj-dk}{e}$

$m=\frac{{\rho}_{xx}^{fL}-fj-gk-hl}{i};\text{\hspace{1em}}n=\sqrt{1-{j}^{2}-{k}^{2}-{l}^{2}-{m}^{2}}$

Remark: See proof in [27].

2.8.3. Simulation of USD Cap and Floor Prices

To price USD interest rate options, it was assumed that the USD LIBOR is a domestic rate in the US. The dynamics used to simulate the USD Forward risky LIBOR in the US economy were given by:

$(\begin{array}{l}{a}_{1}={\rho}_{dD}^{dL}\\ \text{d}{W}_{dL}={a}_{1}\ast {Z}_{1}+{Z}_{2}\sqrt{1-{a}_{1}^{2}}\\ \mathrm{ln}{L}_{d}\left({T}_{i}\right)=\mathrm{ln}{L}_{d}\left({T}_{i-1}\right)+{\rho}_{dD}^{dL}\ast {\sigma}_{dD}\ast {\gamma}_{dL}\ast \Delta t+{\gamma}_{dL}\ast \sqrt{\Delta t}\ast \text{d}{W}_{dL}\end{array}$ (27)

The formula expressed in Equation (28) was used to simulate USD cap or floor prices at time t.

$U\left(t\right)=\Delta t\ast {N}_{d}\ast {\displaystyle \underset{i=1}{\overset{n}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{P}_{d}^{D}\left(t,{T}_{i}\right)\ast \mathrm{max}\left(\omega \left({L}_{d}\left({T}_{i}\right)-{K}_{d}\right),0\right)$ (28)

where ${N}_{d}$ is the notional principal of the USD interest rate option.

2.8.4. Simulation of GBP Cap and Floor Prices.

To price GBP interest rate options, it was first assumed that the GBP LIBOR is a domestic rate in the U.K. The dynamics used to simulate the GBP Forward risky LIBOR in the U.K economy were given by:

$(\begin{array}{l}{a}_{2}={\rho}_{fD}^{fL}\\ \text{d}{W}_{fL}={a}_{2}\ast {Z}_{1}+{Z}_{2}\sqrt{1-{a}_{2}^{2}}\\ \mathrm{ln}{L}_{f}\left({T}_{i}\right)=\mathrm{ln}{L}_{f}\left({T}_{i-1}\right)+{\rho}_{fD}^{fL}\ast {\sigma}_{fD}\ast {\gamma}_{fL}\ast \Delta t+{\gamma}_{fL}\ast \sqrt{\Delta t}\ast \text{d}{W}_{fL}\end{array}$ (29)

The formula expressed in Equation (30) was used to simulate GBP cap or floor prices at time, t.

$G\left(t\right)=\Delta t\ast {N}_{f}\ast {\displaystyle \underset{i=1}{\overset{n}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{P}_{f}^{D}\left(t,{T}_{i}\right)\ast \mathrm{max}\left(\omega \left({L}_{f}\left({T}_{i}\right)-{K}_{f}\right),0\right)$ (30)

where ${N}_{f}$ is the notional principal of the GBP interest rate option.

2.8.5. Simulation of Quanto Cap and Floor Prices

To price GBP options in the US that remits payments in GBP, it is first assumed that the GBP LIBOR is a foreign rate in the US. The dynamics used to simulate the foreign Forward risky LIBORs under the spot domestic risk neutral measure, ${\mathbb{Q}}^{D}$, in the domestic economy was given by;

$\{\begin{array}{l}\text{d}{W}_{fL}=f\ast {Z}_{1}+g\ast {Z}_{2}+h\ast {Z}_{3}+i\ast {Z}_{4}\\ \mathrm{ln}{L}_{f}\left({T}_{i}\right)=\mathrm{ln}{L}_{f}\left({T}_{i-1}\right)+\left({\rho}_{fD}^{fL}\ast {\sigma}_{fD}-{\rho}_{xx}^{fL}\ast {\sigma}_{x}\right)\ast {\gamma}_{fL}\ast \Delta t+{\gamma}_{fL}\ast \sqrt{\Delta t}\ast \text{d}{W}_{fL}\end{array}$ (31)

The formula expressed in Equation (32) was used to simulate Quanto cap or floor prices at time, t.

$Q\left(t\right)=\Delta t\ast {X}^{0}\ast {N}_{f}\ast {\displaystyle \underset{i=1}{\overset{n}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{P}_{d}^{D}\left(t,{T}_{i}\right)\ast \mathrm{max}\left(\omega \left({L}_{f}\left({T}_{i}\right)-{K}_{f}\right),0\right)$ (32)

where ${N}_{f}$ is the notional principal of the GBP interest rate option.

2.9. Black MCCCLMM Quanto Caplet or Floorlet Formula

A quanto is an interest rate derivative that allows the holder to receive payment in a currency different from that of the underlying. A caplet is a call optional derivative that offers payment to the holder whenever the interest rate exceeds the cap price at maturity and a cap is a series of caplets. A floorlet is the opposite of a caplet. It is a put optional derivative that offers payment to the holder whenever the put rate exceeds the interest rate at maturity and in the same way, a floor is a series of floorlets. In this section, it was assumed that there exists a domestic investor interested in hedging against foreign interest rate risk. It was also assumed that the investor prefers using their domestic currency in trading as opposed to using the foreign currency. Hence, the underlying in the quanto options are considered to be struck in foreign currency and payments converted into domestic currency using a fixed exchange rate.

Theorem 1 (Black Scholes MCCCLMM Quanto-Caplet/Floorlet Formula). Assuming that the dynamics of the multi-curve cross currency LIBOR market model under the spot domestic martingale measure ${\mathbb{Q}}^{D}$ is as defined in Equation (1) and further assuming the underlying ${L}_{f}\left(\mathrm{.}\right)$, is struck in foreign currency. Then the payoff at maturity of a quanto caplet or floorlet expressed in domestic currency is given by:

$\tau {N}_{f}{X}^{0}{\left(\omega \ast \left({L}_{f}\left(T\mathrm{,}T\right)-{K}_{f}\right)\right)}^{+}$ (33)

and the arbitrage free price of the quanto caplet or floorlet at time $t\le T$ is given by:

$\begin{array}{l}Q\left(t,T,{K}_{f},{X}^{0},{N}_{f}\right)\\ =\tau \omega {N}_{f}{X}^{0}{P}_{d}^{D}\left(t,T\right)\left[{L}_{f}\left(t,T\right){\text{e}}^{\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}\ast \left(T-t\right)}\Phi \left(\omega {d}_{1}\right)-{K}_{f}\ast \Phi \left(\omega {d}_{2}\right)\right]\end{array}$

where

$\omega =(\begin{array}{l}+1\text{\hspace{1em}}\text{if it is a caplet contract}\\ -1\text{\hspace{1em}}\text{if it is a floorlet contract}\end{array}$

${d}_{1}=\frac{\mathrm{ln}\left(\frac{{L}_{f}\left(t,T\right)}{{K}_{f}}\right)+\left(\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}+\frac{1}{2}\ast {\gamma}_{fL}^{2}\right)\left(T-t\right)}{{\gamma}_{fL}\ast \sqrt{T-t}}$

${d}_{2}={d}_{1}-{\gamma}_{fL}\ast \sqrt{T-t}$

$\Phi \left(\mathrm{.}\right)$ is the cdf of a standard normal distribution.

Remark See derivation of Black Scholes MCCCLMM Quanto Caplet formula in Appendix.

3. Numerical Results

The MCCCLMM model parameters were calibrated to real world data using ATM cap prices and historical rates. The parameters were then used to price quanto caps and floors under the MCCCLMM.

3.1. Data Description

Six months spaced ten year historical data was used in this study to estimate the model correlations and foreign exchange rate volatility. The data was taken for the period beginning from 2nd January 2008 to 2nd January 2018. The data consisted of the Overnight and 6 month GBP and USD LIBOR term structures obtained from [12], [13], and the GBP/USD foreign exchange rate obtained from [33]. The descriptive statistics of the data was as illustrated in Table 2 below.

Table 2. Descriptive statistics of the data.

3.2. Calibration of the MCCCLMM Parameters

Our MCCCLMM model was calibrated to match the cap market prices provided by [36]. The parameters estimated were as summarised in the subsections of this section. For presentability purposes, in this paper, the parameter estimates stated were rounded off to the nearest 5 decimal places. However, the values used while pricing were not rounded off.

3.2.1. Estimation of the MCCCLMM Correlation Parameters

The correlation parameters were estimated using Equation (7) using data collected from [12] [13]. These parameter estimates were as summarized in Table 3.

3.2.2. Estimation of the Foreign Exchange Rate Volatility Parameter

The foreign exchange rate volatility parameter was estimated using the formula expressed in Equation (12) from historical rates collected from [33]. The ${\sigma}_{X}$ estimate was as summarised in Table 4. We assumed that $\Delta t=\mathrm{0.5.}$

3.2.3. Forward LIBOR Volatility Estimates

Just as described in Section 2.5.3, we assumed that the S parameters fit the market cap volatilities. The S parameters for the GBP and USD markets are as summarized in Table 5 and Table 6 respectively. The GBP LIBOR market volatility parameter ${\gamma}_{fL}$ was extracted directly from actual GBP ATM cap prices dataset struck on 2/1/2018 and summarized in Table 5.

Table 3. MCCCLMM correlation parameter estimates from historical rates.

Table 4. ${\sigma}_{X}$ parameter estimate from historical rates.

Table 5. ${\gamma}_{fL}$ estimate.

Table 6. ${\gamma}_{dL}$ estimate (%).

The USD LIBOR volatility parameter ${\gamma}_{dL}$ was also extracted directly from actual USD ATM cap prices dataset struck on 2/1/2018 and summarized in Table 6.

3.2.4. Calibrating the ${\sigma}_{dD}$ Parameter

The ${\sigma}_{dD}$ parameter was calibrated using Equation (17) to exactly match the actual USD ATM cap prices. The optimal values of ${\sigma}_{dD}$ were as summarized in Table 7.

3.2.5. Calibrating the ${\sigma}_{fD}$ Parameters

The ${\sigma}_{fD}$ parameter was calibrated using Equation (17) to exactly match the actual GBP ATM cap prices. The optimal values of ${\sigma}_{fD}$ were as summarized in Table 8.

3.3. The Discount Curve

The US treasury and UK gilt yields collected from [35] and [34] respectively were assumed to be the risk free rates. The USD and GBP yield curves were then estimated using the NSS method.

3.3.1. Bootstrapping of the USD and GBP Yield Curves

The USD and GBP yield curves extracted from [35] and [34] respectively were bootstrapped using the NSS method described in Section 2.6. The estimated NSS model parameters of the two curves were as summarised in Table 9.

Figure 3 and Figure 4 below show how the NSS yield curve matches the actual yields.

3.3.2. Estimation of the Discount Curve

The USD or GBP discount factors were estimated using the formula in Equation (21) where y(t) was taken to be the USD or GBP NSS yields. The estimated USD and GBP discount curves were then plotted in Figure 5 and Figure 6.

3.4. Valuation of Caps and Floors in the Black MCCCLMM Formula

Using the calibrated model parameters, the GBP and USD caps and floors were priced using the Black like formula expressed in Section 2.5.4.

Figure 3. USD NSS yield curve.

Figure 4. GBP NSS yield curve.

Figure 5. USD discount curve.

Figure 6. GBP discount curve.

Table 7. ${\sigma}_{dD}$ estimate.

Table 8. ${\sigma}_{fD}$ estimate.

Table 9. NSS parameter estimates.

3.4.1. USD Cap and Floor Prices

USD ATM Cap and Floor prices for options struck on 2/1/2018, with a range of maturities, were calculated using the calibrated Black-like formula expressed in Equation (17). Table 10 gives the results of the USD ATM Cap and Floor prices in basis points (bp).

Table 10. Black USD ATM CAP and FLOOR prices in bp struck on 2/1/2018.

3.4.2. GBP Cap and Floor Prices

GBP ATM Cap and Floor prices for options struck on 2/1/2018, with a range of maturities, were calculated using the calibrated Black-like formula expressed in Equation (17). Table 11 gives the results of the GBP ATM Cap and Floor prices in basis points (bp).

3.5. Valuation of Quanto Caps and Floors in the Black MCCCLMM Formula

ATM Quanto Cap and Floor prices for options struck on 2/1/2018 with a range of maturities were calculated using the calibrated Black MCCCLMM Quanto Caplet or Floorlet Formula expressed in Theorem 1. Table 12 gives the results of the Quanto ATM Cap and Floor prices in basis points (b.p).

Where ${X}^{0}$ was taken to be the GBP/USD foreign exchange rate closing price on 2/1/2018 given by 1.3588. The time evolution of the quanto option prices is as shown in Figure 7.

From Figure 7 it can be seen that for the 10 year period, the quanto caps are expected to be sold at a higher price compared to the quanto floors. However, it was noted that the rate of appreciation of the floor prices is quite high such that it is expected that they will eventually overtake the cap prices as time goes by.

3.6. Simulation of Cap and Floor Prices

The domestic and foreign LIBOR model dynamics were first discretized using the Euler scheme.

3.6.1. Generation of the Correlated Wiener Processes

The correlation matrix of our observed data was found to be:

$\Sigma =\left(\begin{array}{ccccc}1.00000& 0.95207& 0.85917& 0.79783& 0.50421\\ 0.95207& 1.00000& 0.89546& 0.86769& 0.44993\\ 0.85917& 0.89546& 1.00000& 0.98662& 0.74092\\ 0.79783& 0.86769& 0.98662& 1.00000& 0.74244\\ 0.50421& 0.44993& 0.74092& 0.74244& 1.00000\end{array}\right)$ (34)

Figure 7. Quanto option prices evolution.

Table 11. Black GBP ATM CAP and FLOOR prices in bp struck on 2/1/2018.

Table 12. Black quanto ATM CAP and FLOOR prices struck on 2/1/2018.

The lower triangular matrix defined in Equation (25) was then calculated and found to be:

$A=\left(\begin{array}{ccccc}1.00000& 0.00000& 0.00000& 0.00000& 0.00000\\ 0.95207& 0.30587& 0.00000& 0.00000& 0.00000\\ 0.85917& 0.25325& 0.44462& 0.00000& 0.00000\\ 0.79783& 0.35339& 0.47601& 0.10949& 0.00000\\ 0.50421& -0.09844& 0.74816& 0.17190& 0.38312\end{array}\right)$ (35)

The correlated Wiener process was then estimated as:

$\left(\begin{array}{c}\text{d}{W}_{dD}\left(t\right)\\ \text{d}{W}_{dL}\left(t\right)\\ \text{d}{W}_{fD}\left(t\right)\\ \text{d}{W}_{fL}\left(t\right)\\ \text{d}{W}_{xx}\left(t\right)\end{array}\right)=\sqrt{0.5}\left(\begin{array}{ccccc}1.00000& 0.00000& 0.00000& 0.00000& 0.00000\\ 0.95207& 0.30587& 0.00000& 0.00000& 0.00000\\ 0.85917& 0.25325& 0.44462& 0.00000& 0.00000\\ 0.79783& 0.35339& 0.47601& 0.10949& 0.00000\\ 0.50421& -0.09844& 0.74816& 0.17190& 0.38312\end{array}\right)\left(\begin{array}{c}{Z}_{1}\\ {Z}_{2}\\ {Z}_{3}\\ {Z}_{4}\\ {Z}_{5}\end{array}\right)$ (36)

where ${Z}_{i};i\in \left\{1,2,3,4,5\right\}$ are independent standard normal random variables.

3.6.2. Simulation of USD Cap and Floor Prices

The formula expressed in Equation (28) was used to simulate USD cap or floor prices at time, t. ${N}_{d}=\mathrm{\$10000}$. ${P}_{d}^{D}\left(t\mathrm{,.}\right)$ is the risk-free discount factor associated with the USD discount curve. The mean simulated USD at the money (ATM) Cap and Floor prices stuck on 2/1/2018 were as summarised in Table 13.

3.6.3. Simulation of GBP Cap and Floor Prices

The formula expressed in Equation (30) was used to simulate GBP cap or floor prices at time, t. ${N}_{f}=\pounds 10000$, ${P}_{f}^{D}\left(t\mathrm{,.}\right)$ is the GBP risk-free discount factor associated with the U.K discount curve, $\Delta t=0.5$. The mean simulated GBP at the money (ATM) Cap and Floor prices stuck on 2/1/2018 were as summarised in Table 14.

Table 13. Simulated USD ATM CAP and FLOOR prices struck on 2/1/2018.

Table 14. Simulated GBP ATM CAP and FLOOR prices stuck on 2/1/2018.

3.6.4. Simulation of Quanto Cap and Floor Prices

The formula expressed in Equation (32) was used to simulate quanto cap or floor prices at time, t. ${N}_{f}=\pounds 10000$, ${P}_{d}^{D}\left(t\mathrm{,.}\right)$ is the USD risk-free discount factor associated with the US discount curve, $\Delta t=0.5$. The mean simulated at the money (ATM) quanto Cap and Floor prices stuck on 2/1/2018 were as summarised in Table 15.

3.6.5. Comparison between the Black MCCCLMM Formulas and the Monte Carlo Simulations

The Black MCCCLMM cap and floor prices were compared to those estimated via the Monte Carlo simulation technique so as to ascertain the pricing accuracy of the Black MCCCLMM formulas using the mean absolute percentage error (MAPE) technique. The MAPE results were as given in Table 16. According to [38], if $\text{MAPE}\le \mathrm{10\%}$, then the model is considered to be highly accurate. If $10\%<\text{MAPE}\le 20\%$, then the model is considered to be a good model. If $20\%<\text{MAPE}\le 50\%$, then the model is considered to be reasonable, and if $\text{MAPE}>50\%$ then the model is taken to be inaccurate. Hence from Table 16, we can see that the simulated results can be considered to be highly accurate when compared to the Black MCCCLMM formula as the MAPEs of all the options priced lie below 1%. Hence it was concluded that the two models are almost similar which is a good result. Hence the derived Black MCCCLMM model is reliable.

3.7. Cap and Floor Volatility Surfaces

Using the calibrated parameters, different USD, GBP and Quanto Cap and Floor prices were estimated under different strike price assumptions. The cap volatility

Table 15. Simulated quanto ATM CAP and FLOOR prices struck on 2/1/2018.

Table 16. MAPE results of the quanto cap and floor prices.

surfaces of the USD, GBP and GBP/USD Quanto options were as shown in Figure 8 below: From Figure 8, we noticed that the cap volatility surfaces of the 3 different options had roughly the same shape. It was also noted that in general over the three options, when the strike prices increase, the CAP prices reduced. In addition to this, cap prices increased with maturity.

The Floor volatility surfaces of the USD, GBP and GBP/USD Quanto options were as shown in Figure 9 below: From Figure 9, we noticed that the floor volatility surfaces of the 3 different options also had roughly the same shape. It was also noted that in general, when the strike prices increase, the floor prices also increase. In addition to this, floor prices increased with maturity.

4. Conclusions

This study aimed at pricing quanto caps and floors using the multi-curve cross-currency LIBOR market model dynamics. It focused on calibrating the model to real world USD and GBP cap market prices. The MCCCLMM correlation and foreign exchange rate volatility parameters were estimated from 10 years historical market rates. The generated payoffs were later on discounted using USD treasury rates. The estimated MCCCLMM parameters were then used to price quanto cap and floor options both under the derived Black MCCCLMM quanto caplet and floorlet formula and also via 100,000 Monte Carlo simulations. The two methods were then compared using the mean absolute percentage error (MAPE) to ascertain the pricing accuracy of the derived model. The MAPEs between the two models were found to lie below 1% hence

Figure 8. Cap volatility surfaces.

Figure 9. Floor volatility surfaces.

deeming the two models as similar implying that the derived MCCCLMM quanto caplet or floorlet formula is highly accurate. The volatility surfaces of the quanto options were similar to those of the GBP and USD options.

It should be noted that due to data availability constraints, the MCCCLMM parameters in this paper were neither calibrated using swaptions nor discounted using overnight indexed swap (OIS) rates. Hence in this paper we were constrained on pricing our caps and floors on a single tenor (6 months). Also, the model was discounted using treasury rates as they are also proxies for the risk-free rate. In the future, it is advisable for the model to be discounted using the OIS rates and compare between the two methods.

Acknowledgements

We thank the Editor and the reviewers for their comments. The first author acknowledges financial support from Pan African University Institute of Basic and Applied Sciences Scholarship. This support is greatly appreciated.

Appendix

A. Derivation of the Black MCCCLMM Quanto Caplet Formula

Assuming the parameters are constant such that the dynamics of the foreign risky LIBOR under the spot domestic risk-neutral measure can be expressed as follows:

$\frac{\text{d}{L}_{f}\left(t,T\right)}{{L}_{f}\left(t,T\right)}=\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}\ast \text{d}t+{\gamma}_{fL}\ast \text{d}{W}_{fL}\left(t\right)$ (37)

Then the price at time t of a quanto caplet struck on this LIBOR is given by:

$\begin{array}{l}QCplt\left(t\mathrm{,}T\mathrm{,}{K}_{f}\mathrm{,}{X}^{0}\mathrm{,}{N}_{f}\right)\\ =\tau {N}_{f}{X}^{0}{P}_{d}^{D}\left(t\mathrm{,}T\right){E}^{{\mathbb{Q}}^{D}}\left[{\left({L}_{f}\left(T\mathrm{,}T\right)-{K}_{f}\right)}^{+}\mathrm{|}{\mathcal{F}}_{t}\right]\\ =\tau {N}_{f}{X}^{0}{P}_{d}^{D}\left(t\mathrm{,}T\right){E}^{{\mathbb{Q}}^{D}}\left[\left({L}_{f}\left(T\mathrm{,}T\right)-{K}_{f}\right){I}_{A}\mathrm{|}{\mathcal{F}}_{t}\right]\\ =\tau {N}_{f}{X}^{0}{P}_{d}^{D}\left(t\mathrm{,}T\right)\underset{M}{\underset{\ufe38}{{E}^{{\mathbb{Q}}^{D}}\left[{L}_{f}\left(T\mathrm{,}T\right){I}_{A}\mathrm{|}{\mathcal{F}}_{t}\right]}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}-\tau {N}_{f}{X}^{0}{P}_{d}^{D}\left(t\mathrm{,}T\right)\underset{N}{\underset{\ufe38}{{E}^{{\mathbb{Q}}^{D}}\left[{K}_{f}{I}_{A}\mathrm{|}{\mathcal{F}}_{t}\right]}}\end{array}$ (38)

where ${I}_{A}$ is an indicator function such that

${I}_{A}=(\begin{array}{l}1\text{\hspace{1em}}\text{if}\text{\hspace{0.17em}}{L}_{f}\left(T,T\right)-{K}_{f}>0\\ 0\text{\hspace{1em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{elsewhere}\end{array}$ (39)

M and N were then solved separately and the results inserted back into Equation (38) at the end.

The risky foreign LIBOR defined in Equation (37) under the spot domestic risk neutral measure ${\mathbb{Q}}^{D}$ has a solution given by:

${L}_{f}\left(T,T\right)={L}_{f}\left(t,T\right){\text{e}}^{\left(\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}-\frac{1}{2}\ast {\gamma}_{fL}^{2}\right)\left(T-t\right)+{\gamma}_{fL}\left({W}_{fL}\left(T\right)-{W}_{fL}\left(t\right)\right)}$ (40)

Solving for M: We saw that

$\begin{array}{c}M={\mathbb{E}}^{{\mathbb{Q}}^{D}}\left[{L}_{f}\left(T\mathrm{,}T\right){I}_{A}\mathrm{|}{\mathcal{F}}_{t}\right]\\ ={\mathbb{E}}^{{\mathbb{Q}}^{D}}\left[{L}_{f}\left(t\mathrm{,}T\right){\text{e}}^{\left(\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}-\frac{1}{2}\ast {\gamma}_{fL}^{2}\right)\left(T-t\right)+{\gamma}_{fL}\left({W}_{fL}\left(T\right)-{W}_{fL}\left(t\right)\right)}{I}_{A}\mathrm{|}{\mathcal{F}}_{t}\right]\\ ={L}_{f}\left(t\mathrm{,}T\right){\text{e}}^{\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}\ast \left(T-t\right)}\ast {\mathbb{E}}^{{\mathbb{Q}}^{D}}\left[{\text{e}}^{-\frac{1}{2}\ast {\gamma}_{fL}^{2}\ast \left(T-t\right)+{\gamma}_{fL}\left({W}_{fL}\left(T\right)-{W}_{fL}\left(t\right)\right)}{I}_{A}\mathrm{|}{\mathcal{F}}_{t}\right]\end{array}$ (41)

Defining a new probability measure V given by the Radon Nikodym derivative

$\frac{\text{d}V}{\text{d}{Q}^{D}}={\text{e}}^{-\frac{1}{2}\ast {\gamma}_{fL}^{2}\ast \left(T-t\right)+{\gamma}_{fL}({W}_{fL}\left(T\right)-{W}_{fL}(t)}$

such that

$\text{d}{W}_{fL}^{V}\left(t\right)=\text{d}{W}_{fL}\left(t\right)-{\gamma}_{fL}\ast \text{d}t$ (42)

then

$\begin{array}{c}M={L}_{f}\left(t\mathrm{,}T\right){\text{e}}^{\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}\ast \left(T-t\right)}\mathrm{*}{\mathbb{E}}^{{\mathbb{Q}}^{D}}\left[{\text{e}}^{\frac{1}{2}\ast {\gamma}_{fL}^{2}\ast \left(T-t\right)+{\gamma}_{fL}\left({W}_{fL}\left(T\right)-{W}_{fL}\left(t\right)\right)}{I}_{A}\mathrm{|}{\mathcal{F}}_{t}\right]\\ ={L}_{f}\left(t\mathrm{,}T\right){\text{e}}^{\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}\ast \left(T-t\right)}{\mathbb{E}}^{V}\left[{I}_{A}\mathrm{|}{\mathcal{F}}_{t}\right]\\ ={L}_{f}\left(t\mathrm{,}T\right){\text{e}}^{\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\mathrm{*}{\gamma}_{fL}\mathrm{*(}T-t\mathrm{)}}\ast 0\ast P{r}^{\mathbb{V}}\left({L}_{f}\left(T\mathrm{,}T\right)-{K}_{f}\le \mathrm{0|}{\mathcal{F}}_{t}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}+{L}_{f}\left(t\mathrm{,}T\right){\text{e}}^{\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}\ast \left(T-t\right)}\ast 1\ast P{r}^{\mathbb{V}}\left({L}_{f}\left(T\mathrm{,}T\right)-{K}_{f}>\mathrm{0|}{\mathcal{F}}_{t}\right)\\ ={L}_{f}\left(t\mathrm{,}T\right){\text{e}}^{\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}\ast \left(T-t\right)}\ast 1\ast P{r}^{\mathbb{V}}\left({L}_{f}\left(T\mathrm{,}T\right)-{K}_{f}>\mathrm{0|}{\mathcal{F}}_{t}\right)\end{array}$ (43)

Now solving for $P{r}^{\mathbb{V}}\left({L}_{f}\left(T\mathrm{,}T\right)-{K}_{f}>\mathrm{0|}{\mathcal{F}}_{t}\right)$

Under the V-measure, the dynamics of ${L}_{f}\left(t\mathrm{,}T\right)$ is given by:

$\frac{\text{d}{L}_{f}\left(t\mathrm{,}T\right)}{{L}_{f}\mathrm{(}t\mathrm{,}T\mathrm{)}}=\left(\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}+{\gamma}_{fL}^{2}\right)\ast \text{d}t+{\gamma}_{fL}\ast \text{d}{W}_{fL}^{V}\left(t\right)$ (44)

Solving this using Itô formula, we get:

$\begin{array}{c}{L}_{f}\left(T\mathrm{,}T\right)={L}_{f}\left(t,T\right){\text{e}}^{\left(\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}+\frac{1}{2}\ast {\gamma}_{fL}^{2}\right)\left(T-t\right)+{\gamma}_{fL}\left({W}_{fL}^{V}\left(T\right)-{W}_{fL}^{V}\left(t\right)\right)}\\ ={L}_{f}\left(t,T\right){\text{e}}^{R}\end{array}$ (45)

where

$R=\left(\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}+\frac{1}{2}\ast {\gamma}_{fL}^{2}\right)\left(T-t\right)+{\gamma}_{fL}\left({W}_{fL}^{V}\left(T\right)-{W}_{fL}^{V}\right)(t)$

$\mathbb{E}\left[R\right]={\mu}_{R}=\left(\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}+\frac{1}{2}\ast {\gamma}_{fL}^{2}\right)\left(T-t\right)$

$\mathbb{V}\left(R\right)={\sigma}_{R}^{2}={\gamma}_{fL}^{2}\ast \left(T-t\right)$

such that $R~N\left({\mu}_{R}\mathrm{,}{\sigma}_{R}^{2}\right)$. Therefore:

$\begin{array}{l}P{r}^{\mathbb{V}}\left({L}_{f}\left(T\mathrm{,}T\right)-{K}_{f}>\mathrm{0|}{\mathcal{F}}_{t}\right)\\ =Pr\left({L}_{f}\left(t\mathrm{,}T\right){\text{e}}^{R}-{K}_{f}>\mathrm{0|}{\mathcal{F}}_{t}\right)=Pr\left({\text{e}}^{R}>\frac{{K}_{f}}{{L}_{f}\left(t\mathrm{,}T\right)}\mathrm{|}{\mathcal{F}}_{t}\right)\\ \mathrm{=}Pr\left(R>\mathrm{ln}\left(\frac{{K}_{f}}{{L}_{f}\left(t,T\right)}\right)|{\mathcal{F}}_{t}\right)={\displaystyle {\int}_{{L}^{\mathrm{*}}}^{\infty}}\frac{1}{{\sigma}_{R}\sqrt{2\text{\pi}}}{\text{e}}^{-\frac{1}{2}{\left(\frac{R-{\mu}_{R}}{{\sigma}_{R}}\right)}^{2}}\text{d}R\\ ={\displaystyle {\int}_{{Z}_{1}^{\mathrm{*}}}^{\infty}}\frac{1}{\sqrt{2\text{\pi}}}{\text{e}}^{-\frac{1}{2}{Z}_{1}^{2}}\text{d}{Z}_{1}=\underset{\Phi (d1)}{\underset{\ufe38}{{\displaystyle {\int}_{-\infty}^{{d}_{1}}}\frac{1}{\sqrt{\text{2\pi}}}{\text{e}}^{-\frac{1}{2}{Z}_{1}^{2}}\text{d}{Z}_{1}}}\end{array}$

where

${L}^{*}=\mathrm{ln}\left(\frac{{K}_{f}}{{L}_{f}\left(t,T\right)}\right)$

${Z}_{1}=\frac{R-{\mu}_{R}}{{\sigma}_{R}}{Z}_{1}^{*}=\frac{1}{{\sigma}_{R}}\left(\mathrm{ln}\left(\frac{{K}_{f}}{{L}_{f}\left(t,T\right)}\right)-{\mu}_{R}\right)$

${d}_{1}=-{Z}_{1}^{*}=\frac{1}{{\sigma}_{R}}\left({\mu}_{R}+\mathrm{ln}\left(\frac{{L}_{f}\left(t,T\right)}{{K}_{f}}\right)\right)$

Hence M becomes:

$M={L}_{f}\left(t,T\right){\text{e}}^{\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}\ast \left(T-t\right)}\Phi \left({d}_{1}\right)$ (46)

where

${d}_{1}=\frac{\mathrm{ln}\left(\frac{{L}_{f}\left(t,T\right)}{{K}_{f}}\right)+\left(\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}+\frac{1}{2}\ast {\gamma}_{fL}^{2}\right)\left(T-t\right)}{{\gamma}_{fL}\ast \sqrt{T-t}}$ (47)

Solving for N:

$\begin{array}{c}N={\mathbb{E}}^{{\mathbb{Q}}^{D}}\left[{K}_{f}{I}_{A}\mathrm{|}{\mathcal{F}}_{t}\right]\\ ={K}_{f}\mathrm{*}{\mathbb{E}}^{{\mathbb{Q}}^{D}}\left[{I}_{A}\mathrm{|}{\mathcal{F}}_{t}\right]\\ ={K}_{f}\ast 0\ast P{r}^{{\mathbb{Q}}^{D}}\left[{L}_{f}\left(T\mathrm{,}T\right)-{K}_{f}\le \mathrm{0|}{\mathcal{F}}_{t}\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}+{K}_{f}\ast 1\ast P{r}^{{\mathbb{Q}}^{D}}\left[{L}_{f}\left(T\mathrm{,}T\right)-{K}_{f}\le \mathrm{0|}{\mathcal{F}}_{t}\right]\\ ={K}_{f}\ast P{r}^{{\mathbb{Q}}^{D}}\left[{L}_{f}\left[T\mathrm{,}T\right]-{K}_{f}\le \mathrm{0|}{\mathcal{F}}_{t}\right]\end{array}$ (48)

Solving for $P{r}^{{\mathbb{Q}}^{D}}\left[{L}_{f}\left(T\mathrm{,}T\right)-{K}_{f}\le \mathrm{0|}{\mathcal{F}}_{t}\right]$

We saw that under the ${\mathbb{Q}}^{D}$ -measure,

$\begin{array}{c}{L}_{f}\left(T\mathrm{,}T\right)={L}_{f}\left(t\mathrm{,}T\right){\text{e}}^{\left(\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}-\frac{1}{2}\ast {\gamma}_{fL}^{2}\right)\left(T-t\right)+{\gamma}_{fL}\left({W}_{fL}\left(T\right)-{W}_{fL}\left(t\right)\right)}\\ ={L}_{f}\left(t\mathrm{,}T\right){\text{e}}^{S}\end{array}$ (49)

where

$S=\left(\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}-\frac{1}{2}\ast {\gamma}_{fL}^{2}\right)\left(T-t\right)+{\gamma}_{fL}\left({W}_{fL}\left(T\right)-{W}_{fL}\left(t\right)\right)$ (50)

$\mathbb{E}\left[S\right]={\mu}_{S}=\left(\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}-\frac{1}{2}\ast {\gamma}_{fL}^{2}\right)\left(T-t\right)$ (51)

$\mathbb{V}\left(S\right)={\sigma}_{S}^{2}={\gamma}_{fL}^{2}\ast \left(T-t\right)$ (52)

such that $S~N\left({\mu}_{S}\mathrm{,}{\sigma}_{S}^{2}\right)$. Therefore:

$\begin{array}{l}P{r}^{{\mathbb{Q}}^{D}}\left({L}_{f}\left(T\mathrm{,}T\right)-{K}_{f}>\mathrm{0|}{\mathcal{F}}_{t}\right)\\ =Pr\left({L}_{f}\left(t\mathrm{,}T\right){\text{e}}^{S}-{K}_{f}>0\right)=Pr\left({\text{e}}^{S}>\left(\frac{{K}_{f}}{{L}_{f}\left(t,T\right)}\right)\right)\\ =Pr\left(S>\mathrm{ln}\left(\frac{{K}_{f}}{{L}_{f}\left(t,T\right)}\right)\right)={\displaystyle {\int}_{{L}^{\mathrm{*}}}^{\infty}}\frac{1}{{\sigma}_{S}\sqrt{2\text{\pi}}}{\text{e}}^{-\frac{1}{2}{\left(\frac{S-{\mu}_{S}}{{\sigma}_{S}}\right)}^{2}}\text{d}S\\ ={\displaystyle {\int}_{{Z}_{2}^{\mathrm{*}}}^{\infty}}\frac{1}{\sqrt{2\text{\pi}}}{\text{e}}^{-\frac{1}{2}{Z}_{2}^{2}}\text{d}{Z}_{2}=\underset{\Phi (d2)}{\underset{\ufe38}{{\displaystyle {\int}_{-\infty}^{{d}_{2}}}\frac{1}{\sqrt{2\text{\pi}}}{\text{e}}^{-\frac{1}{2}{Z}_{2}^{2}}\text{d}{Z}_{2}}}\end{array}$

where

${L}^{*}=\mathrm{ln}\left(\frac{{K}_{f}}{{L}_{f}\left(t,T\right)}\right);\text{\hspace{0.17em}}\text{\hspace{0.17em}}{Z}_{2}=\frac{S-{\mu}_{S}}{{\sigma}_{S}}$

${Z}_{2}^{*}=\frac{1}{{\sigma}_{S}}\left(\mathrm{ln}\left(\frac{{K}_{f}}{{L}_{f}\left(t,T\right)}\right)-{\mu}_{S}\right)$

${d}_{2}=-{Z}_{2}^{*}=\frac{1}{{\sigma}_{S}}\left({\mu}_{S}+\mathrm{ln}\left(\frac{{L}_{f}\left(t,T\right)}{{K}_{f}}\right)\right)$

Hence N becomes:

$N={K}_{f}\ast \Phi \left({d}_{2}\right)$ (53)

where

$\begin{array}{c}{d}_{2}=\frac{\mathrm{ln}\left(\frac{{L}_{f}\left(t,T\right)}{{K}_{f}}\right)+\left(\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}-\frac{1}{2}\ast {\gamma}_{fL}^{2}\right)\left(T-t\right)}{{\gamma}_{fL}\ast \sqrt{T-t}}\\ ={d}_{1}-{\gamma}_{fL}\ast \sqrt{T-t}\end{array}$ (54)

Hence, the quanto caplet price at time t, with payments given at maturity in domestic currency whenever the foreign risky LIBOR exceeds the foreign strike price, ${K}_{f}$ is given by:

$\begin{array}{l}QCplt\left(t\mathrm{,}T\mathrm{,}{K}_{f}\mathrm{,}{X}^{0}\mathrm{,}{N}_{f}\right)\\ =\tau {N}_{f}{X}^{0}{P}_{d}^{D}\left(t\mathrm{,}T\right)\ast \left[{L}_{f}\left(t\mathrm{,}T\right){\text{e}}^{\left({\rho}_{fD}^{fL}\ast {\sigma}_{dD}-{\rho}_{xx}^{fL}\ast {\sigma}_{X}\right)\ast {\gamma}_{fL}\ast \left(T-t\right)}\Phi \left({d}_{1}\right)-{K}_{f}\ast \Phi \left({d}_{2}\right)\right]\end{array}$ (55)

where ${d}_{1}$ and ${d}_{2}$ is as expressed in Equations (47) and (54) respectively.

Remark: The Black MCCCLMM quanto floor let formula can be derived in the same way.

B. ATM CAP Market Price Data

Cap prices struck on 2/1/2018 on the 6 month USD and GBP LIBORs with maturities ranging from 1 year to 4 years were used to calibrate some of the MCCCLMM parameters to market. The USD and GBP cap price datasets provided by [36] were as summarised in Table B1 and Table B2 below.

Table B1. Actual GBP ATM CAP prices struck on 2/1/2018.

Table B2. Actual USD ATM CAP prices struck on 2/1/2018.

Data Availability

The historical datasets used are freely available on the websites [11], [12], [33]. [34], [35] as at the date last accessed. The USD and GBP ATM market cap prices used (provided by [36] ) are given in Appendix B.

References

[1] Papapantoleon, A. (2010) Old and New Approaches to Libor Modeling. Statistica Neerlandica, 64, 257-275.

https://doi.org/10.1111/j.1467-9574.2010.00458.x

[2] Vasicek, O. (1977) An Equilibrium Characterization of the Term Structure. Journal of Financial Economics, 5, 177-188.

https://doi.org/10.1016/0304-405X(77)90016-2

[3] Hull, J. and White, A. (1990) Pricing Interest-Rate-Derivative Securities. The Review of Financial Studies, 3, 573-592.

https://doi.org/10.1093/rfs/3.4.573

[4] Chan, K.C., Karolyi, G.A., Longstaff, F.A. and Sanders, A.B. (1992) An Empirical Comparison of Alternative Models of the Short-Term Interest Rate. The Journal of Finance, 47, 1209-1227. https://doi.org/10.1111/j.1540-6261.1992.tb04011.x

[5] Andersen, T.G. and Lund, J. (1997) Estimating Continuous-Time Stochastic Volatility Models of the Short-Term Interest Rate. Journal of Econometrics, 77, 343-377.

https://doi.org/10.1016/S0304-4076(96)01819-2

[6] Cox, J.C., Ingersoll Jr., J.E. and Ross, S.A. (1985) A Theory of the Term Structure of Interest Rates. Econometrica, 53, 321-346.

[7] Heath, D., Jarrow, R. and Morton, A. (1992) Bond Pricing and the Term Structure of Interest Rates: A New Methodology for Contingent Claims Valuation. Econometrica: Journal of the Econometric Society, 60, 77-105.

https://doi.org/10.2307/2951677

[8] Brace, A., Musiela, M. and Ģatarek, D. (1997) The Market Model of Interest Rate Dynamics. Mathematical Finance, 7, 127-155.

https://doi.org/10.1111/1467-9965.00028

[9] Miltersen, K.R., Sandmann, K. and Sondermann, D. (1997) Closed Form Solutions for Term Structure Derivatives with Log-Normal Interest Rates. The Journal of Finance, 52, 409-430.

https://doi.org/10.1111/j.1540-6261.1997.tb03823.x

[10] Jamshidian, F. (1997) LIBOR and Swap Market Models and Measures. Finance and Stochastics, 1, 293-330.

https://doi.org/10.1007/s007800050026

[11] Henrard, M. (2007) The Irony in the Derivatives Discounting. SSRN Electronic Journal.

https://doi.org/10.2139/ssrn.970509

[12] IBORate (2018) USD LIBOR History, 2018.

http://www.iborate.com/usd-libor/

[13] IBORate (2018) GBP LIBOR History, 2018.

http://www.iborate.com/gbp-libor/

[14] Ametrano, F.M. and Bianchetti, M. (2013) Everything You Always Wanted to Know about Multiple Interest Rate Curve Bootstrapping but Were Afraid to Ask.

https://doi.org/10.2139/ssrn.2219548

[15] Bianchetti, M. (2012) Two Curves, One Price: Pricing and Hedging Interest Rate Derivatives Decoupling Forwarding and Discounting Yield Curves.

arXiv: 0905.2770.

[16] Bianchetti, M. and Carlicchi, M. (2011) Interest Rates after the Credit Crunch: Multiple Curve Vanilla Derivatives and SABR. Quantitative Finance, 1-26.

https://doi.org/10.2139/ssrn.1783070

[17] Crépey, S., Grbac, Z. and Nguyen, H.-N. (2012) A Multiple-Curve HJM Model of Interbank Risk. Mathematics and Financial Economics, 6, 155-190.

https://doi.org/10.1007/s11579-012-0083-4

[18] Cuchiero, C., Fontana, C. and Gnoatto, A. (2016) A General HJM Framework for Multiple Yield Curve Modeling. Finance and Stochastics, 20, 267-320.

https://doi.org/10.1007/s00780-016-0291-5

[19] Grasselli, M. and Miglietta, G. (2016) A Flexible Spot Multiple-Curve Model. Quantitative Finance, 16, 1465-1477.

https://doi.org/10.1080/14697688.2015.1108521

[20] Grbac, Z., Papapantoleon, A., Schoenmakers, J. and Skovmand, D. (2015) Affine LIBOR Models with Multiple Curves: Theory, Examples and Calibration. SIAM Journal on Financial Mathematics, 6, 984-1025.

https://doi.org/10.1137/15M1011731

[21] Henrard, M. (2010) The Irony in Derivatives Discounting Part II: The Crisis. Wilmott Journal, 2, 301-316.

https://doi.org/10.1002/wilj.39

[22] Henrard, M. (2014) Interest Rate Modelling in the Multi-Curve Framework: Foundations, Evolution and Implementation. Springer, New York.

https://doi.org/10.1057/9781137374660

[23] Hess, M. (2016) An Arithmetic Pure-Jump Multi-Curve Interest Rate Model. Browser Download This Paper.

https://doi.org/10.2139/ssrn.2817816

[24] Mercurio, F. (2009) Interest Rates and the Credit Crunch: New Formulas and Market Models. Bloomberg Portfolio Research Paper.

[25] Mercurio, F. (2010) Libor Market Models with Stochastic Basis. Bloomberg Education and Quantitative Research Paper.

https://doi.org/10.2139/ssrn.1563685

[26] Mercurio, F. (2017) The Present of Futures: Valuing Eurodollar-Futures Convexity Adjustments in a Multi-Curve World. SSRN Electronic Journal.

https://doi.org/10.2139/ssrn.2987832

[27] Wamwea, C., Ngare, P. and Le Doux Mbele Bidima, M. (2019) A Multicurve Cross-Currency Libor Market Model. Journal of Applied Mathematics, 2019, Article ID: 8246578.

https://doi.org/10.1155/2019/8246578

[28] Wu, L. (2002) Fast at-the-Money Calibration of the Libor Market Model Using Lagrange Multipliers. Journal of Computational Finance, 6, 39-77.

https://doi.org/10.21314/JCF.2002.093

[29] Lagnado, R. and Osher, S. (1997) A Technique for Calibrating Derivative Security Pricing Models: Numerical Solution of an Inverse Problem. Journal of Computational Finance, 1, 14-25.

https://doi.org/10.21314/JCF.1997.002

[30] Brigo, D., Capitani, C. and Mercurio, F. (2001) On the Joint Calibration of the LIBOR Market Model to Caps and Swaptions Market Volatilities.

[31] Riga, C. (2011) The Libor Market Model: From Theory to Calibration. PhD Thesis.

[32] Sidenius, J. (2000) Libor Market Models in Practice. Journal of Computational Finance, 3, 5-26.

https://doi.org/10.21314/JCF.2000.044

[33] Fusion Media Limited (2018) GBP/USD-British Pound US Dollar.

https://www.investing.com/currencies/gbp-usd

[34] Fusion Media Limited (2018) United Kingdom-Government Bonds.

https://www.investing.com/rates-bonds/uk-government-bonds?maturity_from=40maturity_to=310

[35] Fusion Media Limited (2018) United States-Government Bonds.

https://www.investing.com/rates-bonds/usa-government-bonds?maturity_from=40maturity_to=290

[36] TPICAP (2018) Tullett Prebon Information, Interest Rate Options.

https://www.tpicap.com/

[37] Dmouj, A. (2006) Stock Price Modelling: Theory and Practice. Master’s Thesis, Vrije Universiteit, Amsterdam.

[38] Abidin, S.N.Z. and Jaffar, M.M. (2014) Forecasting Share Prices of Small Size Companies in Bursa Malaysia Using Geometric Brownian Motion. Applied Mathematics & Information Sciences, 8, 107-112.

https://doi.org/10.12785/amis/080112

[39] Kloeden, P. and Platen, E. (1992) Numerical Solution of Stochastic Differential Equations. Springer, New York.

https://doi.org/10.1007/978-3-662-12616-5

[40] Benner, W., Zyapkov, L. and Jortzik, S. (2009) A Multi-Factor Cross-Currency Libor Market Model. The Journal of Derivatives, 16, 53-71.

https://doi.org/10.3905/JOD.2009.16.4.053