Abstract

After the dawn of the August 2007 financial crisis, banks became more aware of financial risk leading to the appearance of nonnegligible spreads between LIBOR and OIS rates and also between LIBOR of different tenors. This consequently led to the birth of multicurve models. This study establishes a new model; the multicurve cross-currency LIBOR market model (MCCCLMM). The model extends the initial LIBOR Market Model (LMM) from the single-curve cross-currency economy into the multicurve cross-currency economy. The model incorporates both the risk-free OIS rates and the risky forward LIBOR rates of two different currencies. The established model is suitable for pricing different quanto interest rate derivatives. A brief illustration is given on the application of the MCCCLMM on pricing quanto caplets and quanto floorlets using a Black-like formula derived from the MCCCLMM.

1. Introduction

Modeling of LIBOR rates has evolved since their inception in the late nineties. However, the most notable change occurred in August 2007 where a severe financial crisis caused a number of anomalies in the interest rate markets. Before the 2007 credit crunch, the spreads between the LIBOR rates and the overnight indexed swap (OIS) rates were negligible. In addition to this, the spreads between different LIBOR curves of different tenors were also considered to be negligible. Hence, a single interest rate curve was sufficient for both discounting and generating future cash flows.

However, after the 2007 financial crisis, the LIBOR-OIS spreads of different maturities began to evolve randomly over time. In addition to this, the spreads between different LIBOR curves of different tenors took the same fate. These spreads became substantially too large to be ignored making the negligibility assumption no longer hold. See images by [1] in Figures 1 and 2.

Hence, the possibility of using one curve for both discounting and generating future cash flows was greatly challenged. This led to the introduction of the multiple curve interest rate models. In this new framework, one curve is used for discounting (mostly the OIS (risk-free) curve) and the other curve(s) used for generating future cash flows (mostly different LIBOR curves of different tenors).

So far, there are various models proposed in different pieces of literature with regard to modeling LIBOR under the multiple curve framework: see [2ā€“12] and so on. However, to our knowledge, author of [13] was the first author to suggest this approach and this was ironically before the financial crisis occurred. It can be seen that all the proposed multicurve models can be categorized to fall either under the short rate approach such as [14ā€“16] or the Heath-Jarrow Morton (HJM) approach proposed in the early nineties by [17] or the LIBOR market model (LMM) approach first proposed by [18, 19] in the late nineties.

In addition to this, it was also noted that, under the multiple curve framework, one may choose to model the OIS and LIBOR rates directly. This is said to lead to tractable pricing formulas. However, one can not guarantee the positivity of the LIBOR-OIS spreads. Alternatively, one may choose to model the OIS and LIBOR-OIS spreads directly and then later infer the dynamics of the LIBOR. This ensures that the positivity of the LIBOR-OIS spreads is maintained. However, the pricing formulas derived under this approach have been found to be less tractable. Finally, one may choose to model the LIBOR and LIBOR-OIS spreads directly and later infer the dynamics of the OIS rates. However, the positivity of the OIS rates cannot be guaranteed. In this paper, we follow the first approach.

The aim of this paper is to construct a multicurve cross-currency LIBOR market model under the spot domestic risk neutral measure. The model involves taking into account both the risk-free and the risky interest rates. For this purpose, we adapt the lognormal volatility model for the LIBOR and OIS rates. We also adapt the geometric Brownian motion model for the spot foreign exchange rate. We practically aim to extend the works by [20, 21] from the single-curve cross-currency framework into the multicurve cross-currency framework so that it conforms to the modern practices.

2. Problem Formulation

In this paper, our approach closely resembles that of the Heath-Jarrow-Morton and the LIBOR market model (HJM-LMM) approach in coming up with the multicurve cross-currency LIBOR market model (MCCCLMM).

2.1. Model Notations

Let be the spot foreign exchange rate at time t quoted as the ratio of units of domestic currency to one unit of foreign currency. and denote the domestic and foreign markets, respectively. and denote the risk-free and risky curves, respectively. and are the domestic and foreign risk-free short rates of interest. and are the domestic and foreign risky short rates of interest. and are the corresponding risk-free domestic and foreign money market accounts. and are the spread between the domestic and foreign risky and risk-free short rates of interest, respectively. and are the corresponding risky domestic and foreign money market accounts. and denote the risk-free domestic and foreign zero-coupon bonds. and denote the risky domestic and foreign zero-coupon bonds. and are the domestic and foreign simply compounded t-forward rates associated with the domestic and foreign risk-free discount curves. Finally, and are the domestic and foreign simply compounded t-forward rates associated with the domestic and foreign risky fictitious zero-coupon bonds.

2.2. Model Assumptions

We model a frictionless market, free of arbitrage opportunities. In this market, we assume that trading takes place continuously for a given time interval , where is some positive final date. The market uncertainty is assumed to be modeled by the filtered probability space, . In our model, we consider only one tenor, . We also assume that there exist two markets, domestic (d) and foreign (f), which can be linked with the foreign exchange rate markets. In addition to this, we assume that there exist both risk-free and risky rates in the two markets. The OIS rate is taken to be the most preferable risk-free rate and it is the one that is used to construct the discount curve. On the other hand, we take the LIBOR of a single tenor to be the risky rate and this is the rate that is used to generate the future cash flows.

We assume that there exists, at any time t, a risk-free (domestic or foreign) zero-coupon bond ; such that for all . Assuming that the mapping at any time is differentiable, then the simply compounded risk-free forward interest rate is given as where the risk-free domestic or foreign zero-coupon bond price in [6] is defined as We further assume that there exists a risk-free (domestic or foreign) money market account whose price process is given by where is the (domestic or foreign) risk-free short rate of interest.

We assume that there also exists, at any time t, a risky fictitious (domestic or foreign) zero-coupon bond ; . Assuming that the mapping at any time is differentiable, then the fictitious domestic or foreign bond in [22] is defined as We further assume that there exists a risky (domestic or foreign) money market account whose price process is given by where is the (domestic or foreign) risky short rate of interest and is the k LIBOR-OIS short rate spread.

Finally, given the filtered probability space , it is assumed that the dynamics of the entire economy under the objective measure, , associated with the real world probabilities is given by

2.3. Conditions for No Arbitrage

We see that we have five sources of risk: the domestic risk-free markets, domestic risky markets, foreign risk-free markets, foreign risky markets, and the foreign exchange rate markets and according to Meta Theorem [23], if A is the number of underlying traded assets excluding the risk-free asset and R is the number of random sources of risk, then we have the following.(i)The model is free of arbitrage-free if and only if .(ii)The model is complete if and only if .(iii)The model is complete and arbitrage-free if and only if .

Hence assuming that all these markets are correlated, then by Metaā€™s theorem, a 5-dimensional correlated Wiener process must be used to ensure that the model is complete and arbitrage-free. The correlated 5-dimensional Wiener process is defined as such that

Also, according to the First Fundamental Theorem [23], a model is arbitrage-free if and only if there exists an equivalent (local) martingale measure . Hence, assuming that the risk-free assets in both the foreign and domestic markets are basic traded instruments and according to [4, 6], the risky LIBOR bonds in the two markets are fictitious (i.e., nontraded) assets. Assuming further that there exists a usual domestic risk neutral probability measure , then to ensure that there is no arbitrage, we have the following.(i)Under , all domestically traded assets with a price process of say must have the domestic risk-free short rate of interest, , as its local rate of return [23]. That is, its dynamics will be of the formwhere the volatility vector is the same as the one under the objective measure .(ii)Under , all the domestic fictitious assets with a price process of say must have the domestic risky short rate of interest, , as the local rate of interest. That is, its dynamics will be of the formwhere the volatility vector is the same as the one under the objective measure .(iii)Under , all normalized asset price processes of all domestically traded assets with a price process of say discounted using the domestic risk-free money market account, , as the numeraire, must be martingales [21, 23]. That is, the normalized price processis a martingale.(iv)Under , all normalized asset price processes of all domestic fictitious assets with a price process of say discounted using the domestic risky money market account, , as the numeraire, must be martingales [6]. That is, the normalized price processis a martingale.(v)All normalized asset price processes of all domestically traded assets discounted using the domestic risk-free zero-coupon bond, , are martingales. That is, the normalized price process

2.4. The Multicurve Cross-Currency LIBOR Market Model

In this section, a brief description of how the dynamics of the MCCCLMM is derived under the domestic risk neutral measure, , satisfying the no arbitrage assumptions described in Section 2.3 above, is given. It is clear from our model assumptions that and are domestically traded assets. In addition to this, and are also assumed to be domestic fictitious assets. Hence from the no arbitrage conditions described in Section 2.3 requirements (i) and (ii) and from the dynamics of the entire economy described in (6), we get the dynamics of the domestic assets under the spot domestic risk neutral measures to be given by

Lemma 1. Assuming that the market is arbitrage-free, then the possibility of investing in a certain foreign asset at a foreign risk-free short rate of interest should be equivalent to investing in a domestic asset with a price process (see [23]) where

Therefore, from the dynamics of the entire economy described in (6), it is trivial that the dynamics of the price process is given byHowever, since is also a domestically traded asset, then it should also satisfy the no arbitrage condition under Section 2.3 (i). Hence, under the usual domestic risk neutral measure, ,and from (14), it is clear that the spot foreign exchange rate can be expressed as and hence using the dynamics expressed in (13) and (16), we see that the dynamics of the spot foreign exchange rate under the domestic risk neutral measure, , will be given by

Lemma 2. Assuming that the market is free of arbitrage, then the possibility of investing in a foreign risk-free zero-coupon bond should be equivalent to investing in a general risk-free domestic zero-coupon bond with a price process (see [23]) whereTherefore, using the dynamics of the entire economy described in (6), the dynamics of is derived to be

It should also be noted that is a domestically traded asset; hence, it should satisfy the no arbitrage condition discussed under Section 2.3 (i). Therefore, under the domestic risk neutral measure, ,From (19), it is clear that the foreign risk-free zero-coupon bond can be expressed asand hence using the dynamics expressed in (18) and (21), we find the dynamics of the foreign risk-free zero-coupon bond to be

Remark 3. In practice, we notice that and are the time t domestic prices of the floating leg part of the FRA rates associated with the risk-free simply compounded domestic and foreign forward rates: and , respectively, where is the tenor. In addition to this, and are the time t domestic prices of the floating leg part of the FRA rates associated with the risky simply compounded domestic and foreign LIBOR forward rates: and , respectively. Hence, as we can see, they are also domestically traded assets. Therefore, they should satisfy the no arbitrage condition defined under Section 2.3 (i).

Hence as a consequence of the above remark and in using ItĆ“ā€™s Lemma, the dynamics of the different simply compounded forward rates of interest under the usual domestic risk neutral measure are found to be given by

Remark 4. All the domestically traded assets were checked and it was confirmed that they all satisfy the no arbitrage condition defined in Section 2.3 (iii). This means that all the normalized domestically traded assets discounted using the domestic risk-free money market account as the numeraire were all found to be martingales.

Proof. See Proof in Appendix A.

Theorem 5 (the multicurve cross-currency LIBOR market model). Given the filtered probability space and assuming that the 5-dimensional Wiener process, , is a vector of correlated Wiener processes such thatthen the dynamics of the multicurve cross-currency LIBOR market model under the spot domestic risk neutral martingale measure is given bywherewhere denotes the greatest integer that is less than .

Remark 6. Assuming there exists an equivalent risk neutral measure known as the -forward measure, then, for there to be no arbitrage (as explained in Section 2.3 (v)), all the normalized domestically traded assets discounted using the domestic risk-free zero-coupon bond, as the numeraire, must be martingales. Hence, as consequence,

Proof. See proof in Appendix B.

Remark 7. According to [6], the normalized asset processis a martingale under the measure if and only if the risky domestic simply compounded LIBOR satisfies the relationshipHowever, for in (26) to be a martingale, thenand hence to ensure that there are no arbitrage opportunities, this result must be equivalent to the one in (25). This means thatand it is logical to assume that there will always be a positive relationship between the risk-free domestic (or foreign) market and risky domestic (or foreign) market, respectively. HenceImplying that, for the condition under (29) to hold, then,which is a very strict condition. That is why, in this paper, we assume that the equality in (27) is not necessarily true. Hence on (26) need not be a martingale.

Lemma 8. Let be a 5-dimensional standard Wiener process defined on and Girsanovā€™s kernel, , be a 5-dimensional adapted column vector process given bysuch that the domestic market price of risk is given bywhereand then for a fixed time , we shall define a process on as the process given bysuch that . Therefore, a new probability measure on is now defined by the processsuch thatwhere denotes the transpose and is the 5-dimensional standard Wiener process defined on .

Theorem 9 (the multicurve cross-currency LIBOR market model). Given the filtered probability space and assuming that the 5-dimensional Wiener process, , is a vector of correlated Wiener processes such thatthen the dynamics of the multicurve cross-currency LIBOR market model under the forward domestic risk neutral martingale measure is given by

3. Methodology

In this section, the relevant tools, models, methods, and tests used to achieve the objectives of this study are presented.

3.1. Data

Five datasets were considered in this study. The datasets consisted of the British pound (GBP) overnight and 3-month LIBOR, the overnight and 3-month United States dollar (USD) LIBOR, and the GBP/USD foreign exchange rate. The test datasets consisted of 421 daily trading daysā€™ historical data collected from [24ā€“26] for the period between 3/1/2017 and 3/9/2018. The train dataset consisted of 23 daily trading days for the period between 4/9/2018 and 4/10/2018. These datasets only consist of rates recorded on working days in both the United States and the British economy. Hence, it excludes weekends and any public holiday in either economy. We assumed that, in a year, there was an average of 245 working days in both economies.

3.2. Data Analysis Tool

R open source software version 3.1.2 was used in analyzing all the data in this study. Useful packages used were ā€œmassā€, ā€œstats4ā€, ā€œxtsā€, and ā€œlmtestā€.

3.3. Parameter Estimation

It was noted that the dynamics of the MCCCLMM were actually forms of the Geometric Brownian motion; hence a brief description of how the model looks like and how the model parameters were estimated is given in this section.

3.3.1. The Geometric Brownian Motion

The Geometric Brownian motion (GBM) solves the stochastic differential equation given bywhere and . The solution of a GBM is given byThe conditional density function of a GBM model is log normal with a mean ofand a variance of Hence

3.3.2. Maximum Likelihood Estimation Method

The parameters of the MCCCLMM were estimated using the maximum likelihood estimation method. The method tends to maximize the likelihood function. The maximum likelihood estimate is given bywhere is the likelihood function. In R, this is achieved using package ā€œstats4".

3.4. Test of Significance

Individual significance tests performed on the estimated parameters ensure that the fitted parameters are significant. The test hypotheses are where is the parameter estimate under consideration. The null hypothesis is rejected at level of significance when the p-value is less than .

3.5. Model Simulation

According to [20], the correlated Wiener process, , can be simulated by applying Cholesky decomposition as follows:where are independent standard normal variables.

The Euler Maruyama discretization scheme [27] was then used to simulate the fitted MCCCLMM model.

4. Numerical Results

The MCCCLMM model was fitted to real world data and the results were as recorded in this section.

4.1. Data Description

The data used in this study consisted of 421 trading daysā€™ data from 3rd January 2017 to 3rd September 2018. The data consisted of the overnight and 3-month GBP and USD LIBOR term structures obtained from [24, 25] and the GBP/USD foreign exchange rate obtained from [26]. The descriptive statistics of the various sets of the data were recorded in Table 1.

4.2. Parameter Estimation

The parameter estimates of the various dynamics of the multicurve cross-currency LIBOR market model were estimated from the 421 points of the test data via the maximum likelihood estimation (MLE) method using ā€œstats4ā€ package in R. The parameter estimates, their standard errors, and p values were rounded off to the nearest 5 decimal places. The results were as recorded in Table 2.

From Table 2, it can be seen from the respective p-values that all the parameter estimates are significant at level of significance. Also, according to [28], if the standard errors of our estimates are less than 5%, then it is said that the calibration procedure was successful.

4.3. The Simulation Process of the MCCCLMM

The correlation matrix of our observed data was found to beThe correlated Wiener process was then estimated aswhere are independent standard normal random variables and . The overnight and the 3M USD LIBOR path, GBP LIBOR path, and the GBP/USD foreign exchange rate path were simulated and their results were as recorded in Sections 4.3.1ā€“4.3.4.

4.3.1. Simulation of the Overnight USD LIBOR

1000 simulations of the overnight USD LIBOR term structure were done on the train dataset and the descriptive statistics of the simulated path was summarized in Table 3.

The visual plot of the actual path, the 1000 simulated paths, the mean simulated path, and the 95% confidence intervals of the simulated paths were as seen in Figure 3.

It is observed that all the data points lie within the 95% confidence interval of the simulated paths.

4.3.2. Simulation of the 3-Month USD LIBOR

1000 simulations of the 3-month USD LIBOR term structure were done on the train dataset and the descriptive statistics of the simulated path was summarized in Table 4.

The visual plot of the actual path, the 1000 simulated paths, the mean simulated path, and the 95% confidence intervals of the simulated paths were as seen in Figure 4.

It is observed that all the data points lie within the 95% confidence interval of the simulated paths.

4.3.3. Simulation of the Overnight GBP LIBOR

1000 simulations of the GBP overnight LIBOR term structure were done on the train dataset and the descriptive statistics of the simulated path was summarized in Table 5.

The visual plot of the actual path, the 1000 simulated paths, the mean simulated path, and the 95% confidence intervals of the simulated paths were as seen in Figure 5.

It is observed that all the data points lie within the 95% confidence interval of the simulated paths.

4.3.4. Simulation of the 3 Month GBP LIBOR

1000 simulations of the GBP 3-month LIBOR term structure were done on the train dataset and the descriptive statistics of the simulated path was summarized in Table 6.

The visual plot of the actual path, the 1000 simulated paths, the mean simulated path, and the 95% confidence intervals of the simulated paths were as seen in Figure 6.

It is observed that all the data points lie within the 95% confidence interval of the simulated paths.

4.3.5. Simulation of the GBP/USD Foreign Exchange Rate

1000 simulations of the GBP/USD Foreign exchange rate term structure were done on the train dataset and the descriptive statistics of the simulated path was summarized in Table 7.

The visual plot of the actual path, the 1000 simulated paths, the mean simulated path, and the 95% confidence intervals of the simulated paths were as seen in Figure 7.

It is observed that all the data points lie within the 95% confidence interval of the simulated paths.

4.4. Testing the Fitted MCCCLMM Dynamics

The mean absolute percentage error (MAPE) was performed on the difference between the actual and the 1000 simulated paths of the overnight and 3-month USD LIBOR, GBP LIBOR, and the GBP/USD foreign exchange rate and the results were as recorded in Table 8.

According to [29], if , then the forecasts are highly accurate. If , then the forecasts are good. If , then the forecasts are reasonable. However, if , then the forecasts are inaccurate. Hence, we can conclude that the MCCCLMM produces highly accurate forecasts for the 23 working daysā€™ period.

5. Application to Pricing Quanto Caplets and Floorlets

In this section, we illustrate briefly how the dynamics derived under the multicurve cross-currency LIBOR market model can be used to price quanto interest rate derivatives such as quanto caplets and quanto floorlets. A caplet is a call optional type of interest rate derivative where the investor or hedger receives payments if the interest rate exceeds the agreed strike price at maturity. In the same way, a floorlet is a put optional type of interest rate derivative where the hedger receives payment if the interest rate falls below the agreed strike price at maturity. In addition to this, a quanto is a type of derivative whereby the underlying instrument is denominated in one currency but settled in another currency.

We start by assuming that there exists a domestic investor who is interested in hedging against a foreign interest rate risk. We also assume that the investor is more comfortable in using his/her domestic currency in trading and hence all the quanto caplets and floorlet described hereby are priced in terms of the domestic currency. We shall consider a scenario whereby a foreign caplet or floorlet is struck in foreign currency but has to be converted into domestic currency using either a fixed exchange rate or a floating exchange rate.

5.1. Fixed Exchange Rate

Consider a domestic investor who wishes to buy a caplet or floorlet contract struck in foreign currency at a foreign strike price of where the underlying is the risky foreign LIBOR forward rate, . If the seller of the option fixes the exchange rate at an agreed value say at the inception of the contract, then the payoff of this contract at time expressed in domestic currency will be given bywhere is the pricipal of the option expressed in units of foreign currency, is the fixed exchange rate agreed upon at the inception of the contract expressed as the ratio of domestic currency to one unit of foreign currency, and is a binary operator such thatLet the value of the foreign caplet or floorlet at time t expressed in domestic currency be denoted by . Therefore, since a fixed exchange rate is considered then it should be noted that the dynamics of the foreign risky forward LIBOR, , in this case will be valued in the domestic economy. Hence under the domestic risk neutral forward measure , Hence the pricing boundary value problem on used to price such a contract will be given by

The value at time of the foreign caplet or floorlet struck in foreign currency but expressed in terms of domestic currency to one unit of foreign currency given that a fixed exchange rate is considered will therefore be given bywhere and is the cumulative standard normal distribution.

The advantage of such a contract is that it will shield the domestic investor from the risk of exposure to exchange rate risk.

5.2. Floating Exchange Rate

Consider a domestic investor who wishes to buy a caplet or floorlet contract struck in foreign currency at a foreign strike price of where the underlying is the risky foreign LIBOR forward rate, . If the seller of the option assumes that the exchange rate that will be considered will be the spot exchange rate at maturity, , then the payoff of this contract at time expressed in domestic currency will be given bywhere is the pricipal of the option expressed in units of foreign currency, is the spot foreign exchange rate at maturity of the contract expressed as the ratio of domestic currency to one unit of foreign currency, and is a binary operator such that Let the value of the foreign caplet or floorlet at time t expressed in domestic currency be denoted by . Therefore, since a floating exchange rate is considered, then it should be noted that the position is unhedged meaning that the contract is exposed to foreign exchange rate risk. This means that the dynamics of the foreign risky forward LIBOR, , in this case will rely directly on the foreign economy. Hence under the foreign risk neutral forward measure ,Hence the pricing boundary value problem on used to price such a contract will be given by

The value at time of the foreign caplet or floorlet struck in foreign currency but expressed in terms of domestic currency to one unit of foreign currency given that a fixed exchange rate is considered will therefore be given bywhere and is the cumulative standard normal distribution. It should be noted that this type of contract is directly affected by exchange rate movements. The advantage of such a contract is that if the exchange rate moves upwards (or downwards for the floorlet option), then the domestic investor is set to make a profit. However, if the exchange rate moves downwards (or upwards for the floorlet option), then the payout from the option is set to reduce.

6. Conclusion

This study extends the concepts by [20, 21] into the multiple curve cross-currency setting. It focused on constructing a model that can be used to model the simply compounded forward rates of both the domestic and foreign markets under the domestic risk neutral probability measure. The study assumed that interest rates are strictly positive, hence the choice of lognormal type of models. However, in recent times, negative IBOR rates have been recorded. Therefore, this study can be extended to include models that incorporate negative interest rates.

Model parameters were estimated using 421 trading daysā€™ datasets of the GBP overnight and 3-month LIBOR, USD overnight and 3-month LIBOR, and the GBP/USD foreign exchange rate data for the period between 3rd January, 2017, and 3rd September, 2018. The estimated parameters were then used to simulate 1,000 sample paths of 23 trading days out of sample estimates from 4th September, 2018, to 4th October, 2018. From the mean absolute percentage errors (MAPE) calculated, it was seen that the MCCCLMM produces accurate results for this period.

An illustration of how quanto optional interest rate derivatives such as the quanto caplets and floorlets can be valued under the MCCCLMM was also done. However, this should not limit the research on valuing other types of derivatives as the derived MCCCLMM is a robust model that can be used to price numerous interest rate derivatives.

In addition to this, the derived MCCCLMM model is a relatively new model yet to be tested on various term structures or interest rate derivatives and hence pricing performance tests can still be done on it to ascertain its accuracy when applied to more real world data. This will help illustrate how the model can be effectively used in pricing various interest rate derivatives including quantos.

Appendix

A. Proof That the No Arbitrage Condition in Section 2.3 (iii) Is Satisfied

(1)is a martingale.

Proof. which is a martingale under the measure.

(2)is a martingale.

Proof. which is a martingale under the measure.

(3) is a martingale.

Proof. which is a martingale under the measure.

(4)is a martingale.

Proof. which is a martingale under the measure.

(5)is a martingale.

Proof. which is a martingale under the measure.

(6)is a martingale.

Proof. which is a martingale under the measure.

B. Proof That the No Arbitrage Condition in Section 2.3 (v) Is Satisfied

(1)is a martingale.

Proof. which is a martingale under the measure, where

(2)is a martingale.

Proof. which is a martingale under the measure, where

(3)is a martingale.

Proof. which is a martingale under the measure.

(4)is a martingale.

Proof. which is a martingale under the measure.

(5)is a martingale.

Proof. which is a martingale under the measure, where

(6) is a martingale.

Proof. which is a martingale under the measure,

where

Data Availability

The data used is freely available on the websites as at the date last accessed. See [24ā€“26].

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The first author acknowledges financial support from Pan African University Institute of Basic and Applied Sciences Scholarship.