Abstract
The concern of this article is to derive a regime switching model that can be utilized to price European call options for a financial market that exhibits structural changes with time. The model is formulated based on the fact that the underlying asset process is described by a geometric Brownian motion that is modulated by a continuous-time Markov chain with two regimes. Moreover, by an application of the change of measure technique, an option price is derived under the risk neutral valuation and the model parameter estimates is performed by use of the maximum likelihood estimation. The model implementation is carried out by utilizing the Russell 2000 and Facebook in dices data sets. The model results are compared with that of the Black-Scholes model in order to establish the model with better results in terms of predicting the European call option prices. In general, the data sets have common characteristics of financial time series across the regimes and the volatility process spends longer time in regime 2 than it stays in regime 1. The predicted call option prices from both models are more or less similar across the market indices; however, the results of the Black-Scholes model are a bit closer to the market prices than that of the regime-switching model across the two markets. Therefore, the Black-Scholes model slightly gives better results for the Russell 2000 and Facebook indices data sets as compared with the RS model.
Public Interest Statement
In view of the existing literature, there is limited evidence that regime-switching (RS) model performs better than the famous Black-Scholes model in terms of pricing European options. This study develops a RS model for pricing European options in the event that the underlying asset dynamics depend on market regimes. The model developed here, when implemented to compute the European options, can help investors in making decision on when to trade their assets in order to get positive returns.
1 Introduction
The current literature portrays the Black-Scholes formula as the most broadly utilized in option pricing. The model utilizes the Geometric Brownian motion in capturing the dynamics of the underlying asset and it is assumed that the two model parameters, that is, the expected rate of return and the volatility, are deterministic, see Yao et al. (Citation2006). Furthermore, the model assumes that volatility is steady until the lapse date, however, volatility is known to fluctuate with time, that is, it is a well-known fact that the financial market is stochastic in nature, see Mitsui and Satoyoshi (Citation2010). It is thus evident that the Black-Scholes model cannot capture both the stochastic variability of the market dynamics and the volatility smile (the implied volatility of the underlying security). To address these shortfalls of the B-S model in modeling the financial market dynamics, the RS model of Hamilton (Citation1989) has received an extensive application in the context of pricing derivatives. This model allows its parameters to flexibly change over time in accordance with the underlying state process that follow a Markov chain. Furthermore, the model has the ability to reproduce the common stylized facts of financial time series such as fat tails, volatility clustering, among others. The models are also utilized to price and hedge long-dated options in addition to allowing recovery of volatility smiles exhibited by empirical option prices; see Yao et al. (Citation2006) and Godin and Trottier (Citation2021). In support of RS models exhibiting autonomous shifts in mean and variances, Bollen et al. (Citation2000) and Hardy (Citation2001), proposed a RS log-normal (RSLN) model in which log-returns follow a normal distribution with mean and variance that depend on regime variable. In a study to price both American and European options using RS model, Bollen (Citation1998) utilizes lattice method and simulation, whereas Buffington and Elliott (Citation2002) use risk-neutral pricing and derive a set of partial differential equations for option price.
The evaluation of European options has previously been done using Black and Scholes model. It is thus very important to formulate the variation in volatility and evaluate option prices. In an effort to understand the variation in volatility, we refer to Engle (Citation1982) who proposed the ARCH model that formulates the volatility at each time as the linear function of the square of the past unexpected shock. Additionally, Bollerslev (Citation1986), added the past volatility values to the explanatory variables and extended the GARCH (generalized ARCH) model. Empirical studies of options utilizing the GARCH model based on local risk neutrality, which was proposed by Duan (Citation1995), have been conducted by Mitsui (Citation2000), Duan and Zhang (Citation2001), Bauwens and Lubrano (Citation1998), Mitsui and Watanabe (Citation2003). Incidentally, it is known that in RS model the persistence of shock on volatility is extremely high. However, Dieobold (Citation1986) and Lamoureux and Lastrapes (Citation1990) reported that such persistence is considered to be caused by the structural volatility change. It is therefore based on this fact, Hamilton and Susmel (Citation1994) and Cai (Citation1994) proposed the Markov-Switching ARCH (MS-ARCH) model using a state variable that follows the Markov process in the formulation of the ARCH model so as to take into account the structural change. Moreover, Gray (Citation1996) proposed the Markov Switching GARCH model by taking into account the structural change in the GARCH model. Siu et al. (Citation2004) proposed an option pricing model based on GARCH assumption for underlying assets in the context of dynamic version of Gerber-Shiu’s option-pricing model and utilized conditional Esscher Transforms to develop a martingale measure in the context of an incomplete market situation. Their results were consistent with that of Duan (Citation1995) under the assumption of conditional normality for the stock innovation they their results are justified based on the dynamic framework of utility maximization problems as per the Gerber-Shiu’s option-pricing model. Barone-Adesi et al. (Citation2008) in their research proposed option pricing method that was based on GARCH models with filtered historical innovations. Their model outperformed other GARCH pricing models and Black—Scholes models empirically for S&P 500 index options. The model did explain implied volatility smiles by the negative asymmetry of the filtered historical innovations. Empirical evidence and the extent of deterioration of the delta hedging in the presence of large volatility shock were revealed by the study.
Recently, studies on option pricing have been carried out, for instance, Biswas et al. (Citation2018) used a RS stochastic model with a semi-Markov modulated square root mean reverting process to study European option pricing. Their study discovered the local risk-minimizing European-type vanilla options price and it was proven that the price function adequately meet a non-local degenerate parabolic PDE which is the general case of the Heston PDE model. A study by Deelstra et al. (Citation2020) considers risk-neutral pricing of Vanilla, digital and down-and-out call options with the price of the underlying asset evolving according to the exponential of a Markov-modulated Brownian motion (MMBM) with two-sided phase-type jumps. Further, he documents that such options are closely related to the MMBM’s first passage properties which he analyzed by randomizing the time horizon using Erlang distribution. Moreover, Nasri et al. (Citation2020) used a RS Copula model to option prices where serial independence of error terms for time series is presented which was subjected to several time series analysis. Moreover, Kirkby and Nguyen (Citation2020) utilized duality and FFT-based density projection implementation to establish a novel and efficient transform approach for pricing Asian options for general asset dynamics, such as RS Levy processes and stochastic volatility models with jumps. In his research, Lin and He (Citation2020) addressed the problem of pricing the European options using a RS finite moment log-stable model. Their model has the ability to represent the key features of asset returns as well as the impact of regime transition, and it is compatible with market findings. Furthermore, they state that option prices are driven by a coupled FPDE (Fractional Partial Differential Equation) system in their model. They note that their model can mimic option prices emanating from other existing models with particular parameter settings while also bringing into significance pricing differences relative to current models by changing parameter values, indicating that it has the potential to be used in practice. Citation2021) in his study reveals that Esscher transform approach for obtaining risk neutrality is famous in pricing options under RS models. However, this approach creates path-dependence in the dynamics of option price since under the physical measure; the underlying asset is integrated in a Markov process. To address the path-dependence they develop a novel and intuitive risk-neutral measures that integrates risk-aversion in regimes and consequently removing the path dependence side effects. Later on, the study of Citation2021) was complimented by Godin and Trottier (Citation2019) who developed and utilized extended Girsanov principle which also removes the path-dependence effect. This model is easy to interpret in terms of consistence with hedging agents locally minimizing their risk-adjusted discounted squared hedging errors.
In view of the available literature, this study is certainly not the first one to deal with RS dynamics in financial markets or pricing options under RS. However, in terms of pricing European options there seem to be limited evidence that RS model gives better pricing results than the famous Black-Scholes models and thus further research in this area of European option pricing remains necessary. Therefore, the key objective of this study is to develop a RS model for pricing options when the underlying asset dynamics depend on the market regimes. The formulation of this model is based on a geometric Brownian motion that is governed by a continuous-time Markov chain with two states and applying a change of measure, an option price is derived under the risk-neutral valuation. Moreover, the model parameters are estimated using the maximum likelihood estimation method and finally, an implementation of the model is done by computing the European call option prices for some chosen stock market indices. The results of this model are compared with those of the Black-Scholes model with an aim to establish whether the RS model is a better and reliable model compared with the Black-Scholes model.
The organization of the rest of the paper is as follows. Section 1 outlines the derivation of the European option pricing RS model and estimation of parameters. Section 2 is the empirical analysis of the data to implement the model and lastly, section 3 concludes the paper.
2 Methodology
2.1 Regime switching market dynamics
Let be a stochastic process on a discrete time set and be a probability space, where is a physical probability measure. Further let, be Markov chain with states or regimes and transition probability . We consider a risk-free asset, that is assumed to be continuously compounded in value at a constant risk-free rate, (but not necessarily non-negative) over the trading time interval [t,T]. Taking , the price of the asset at time , is given by . Consider a second market which is risky and whose price is defined as
where is the initial price of the asset. The log-returns for the asset are given by
The constants represents the mean and standard deviation of log-returns under each regime, and is a sequence of i.i.d random variables with zero mean and unit variance under measure and independent from the Markov chain . The filtration on the probability space contains all the information generated by all the underlying asset prices and realized regimes upto time . However, the regimes are latent variables hence this filtration characterizes only a partial information to investors. Therefore, we assume that for
where is the transition probability from state of the Markov chain
2.2 Regime-switching model
We consider a risk-free asset and a risky asset as discussed earlier. These assets are tradable continuously over time in a finite time horizon where . In order to describe uncertainty, we consider a complete probability space where is a real world probability measure. According to Yao et al. (Citation2006), in the RS world, one typically modulates the rate of return and the volatility by a finite-state Markov chain , which represents the market regime. Let be continuous, finite-state, Markov process on with state space and that is the stock price at time satisfying the stochastic differential equation defined by
where is the initial price, is a standard Brownian motion independent of .
The parameters and denote the expected rate of return on the asset and volatility of the asset price, respectively, and are assumed to be constant and distinct for each regime. In this way, is regarded as a variable that chooses one of the states in at time of the market. EquationEquation (4)(4) (4) is solved using Itô lemma to give
The stock price process in EquationEquation (5)(5) (5) is assumed to exhibit regime switches. Let the price process undergo discrete shifts between regimes and that it follows a first-order Markov chain, then the transition probability from state at time to state at time is denoted by
where and . The transition matrix of the Markov chain is given by
In particular, for a very small time interval ,
That is, in a short interval of time , the probability of leaving state is approximately . This means, is the transition rate out of state and can formally be expressed as
Since is the probability to move from state to state , the quantity is the transition rate from state to state .
Definition 2.1 (Generator matrix)
Define a generator matrix, for a continuous-time Markov chain. The element of the transition matrix is given by
The generator matrix of the Markov chain is thus given by
and, for each state, the elements of the generator matrix satisfy the equation . For small time interval and the probability to be in state at time given that the random variable was in state at time is given by
where is the probability from state to state in more than one step but .
The model developed here postulates that the regimes are unobservable and thus state transitions until maturity are considered. Denote by , the variable representing the regime in which the process was at time . When the change in log price is in regime , it is presumed to have been drawn from a normal distribution with mean, and variance, . Hence, the probability density function of conditional on taking on the value is given by
for and is a vector of parameters , that is, Furthermore, it is presumed that the unobserved regime have been generated by a probability distribution for which the conditional probability that given the information up to time is denoted by , that is
Thus, the vector also includes , that is, . Next is to find the probability of the joint event that and falls within some time interval [t,T]. This is determined by integrating
over all values of between and . The EquationEquation (13)(13) (13) is the joint probability density function of and and by utilizing Equations (11) and (12), it is expressed as
The unconditional probability density function of is determined by summing EquationEquation (14)(14) (14) over all values for :
Again, based on the definition of conditional probability,
Note that the magnitude in EquationEquation (16)(16) (16) for each observation in the sample can be calculated if the knowledge of the parameters is known by use of EquationEquations (11(11) (11) -Equation14(14) (14) ). EquationEquation (16)(16) (16) represent probability, given the observed data, that the unobserved regime responsible for observation was regime . The probability of given the information up to time , that is, can be calculated by first calculating
with the following equation
where is the transition probability as earlier defined. Next, adding the data at the time leads to the equation
Note that the probability, for can be obtained by repeating the calculations of EquationEquations (18(18) (18) ,Equation19(19) (19) ), and the results substituted into EquationEquation (14)(14) (14) .
In order to price options, it is key to derive a pricing model that produces no arbitrage (a situation where investors can make a guaranteed profit without incurring risk). To avoid arbitrage, the options are priced under risk neutral measure or martingale measure. Under the physical probability measure , we wish to have to be a martingale. Here, is the risk-free interest rate and it is assumed to be constant across the regimes. Since denotes the price of a stock at time which satisfies EquationEquation (4)(4) (4) , we assume that are mutually independent; and , for all . Suppose that denote the sigma field generated by . We note that all(local) martingales concerned are with respect to the filtration. It is clear that and , are both martingales. Define an equivalent measure under which the discounted stock process is a martingale. Let the risk-free rate be denoted by . Then for , let
where . By Girsanov’s theorem, the process
is a -Brownian motion. Moreover, applying Ito’s rule results to where is a local martingale, with . An equivalent measure is defined via the Radon-Nikodym derivative, . Combining Equations (4) and(21) gives
This model has two types of random sources, and the inclusion of makes the underlying market incomplete.
Let denote the total time spent in regime in the interval in trials, given that at time , the state is . Denote the probability by for . For simplicity, we restrict ourselves to two regimes, i.e, , hence the transition matrix, discussed earlier, reduces to
The fraction of times spent by the Markov chain in each regime, as the numbers of transitions become large, can be calculated using a time-average (invariant) distribution of the Markov chain, that is , where is the transition matrix and is the average fraction of the time spent in state over steps as n approaches infinity. Let , then
This results into two equations as both of which simplify to . Again, since is a valid probability distribution, and solving we get
Since, is the trading period, the total time spent in regimes 1 and 2 can now be calculated as , respectively. In view of the research by Duan et al. (Citation2002) and Hardy (Citation2001) the distribution of log returns, , conditional on the total time spent in regime , can be developed such that there exist a normal density function with mean and variance , that is,
where . Since is the probability function for ,
where is the standard normal probability distribution function. This implies that the probability density function for is
where is the standard normal density function.
Now, define , the European call option (under RS world) with strike price that matures after time , and is valued at at an initial time . Since in a RS market, the parameter switches regimes, we can define a parameter conditional on knowing , the total time spent in regime . This implies that, the asset price has a log normal distribution with parameters that depend on , that is, the parameters are as defined earlier. Now, to derive a RS option pricing model, the Black-Scholes formula is considered, where the parameter is replaced with to give the desired model as below;
2.2.1 Model parameter estimates
EquationEquation (15)(15) (15) best describes the actually observed data since the regime is unobserved. If the regime variable is i.i.d across different dates, , then the log likelihood for the observed data can be calculated from EquationEquation (15)(15) (15) as
The maximum likelihood of is obtained by maximizing EquationEquation (29)(29) (29) subject to the constraints that and for . To obtain the maximum likelihood estimates(MLEs) of EquationEquation (29)(29) (29) we form the Lagrangian
and set the derivative w.r.t equal to zero. From EquationEquation (29)(29) (29) , the derivative of the log likelihood is given by
Note that from EquationEquation (15)(15) (15)
Thus
Recalling EquationEquation (16)(16) (16) the derivative in EquationEquation (33)(33) (33) can be written as
From Equation (33) the derivative of Equation (30) with respect to (w.r.t) is given by
This implies that
Summing EquationEquation (36)(36) (36) over produces
or . Replacing with T in EquationEquation (36)(36) (36) produces
Next is to find the MLE of . From EquationEquation (15)(15) (15) it follows that
hence
Applying EquationEquation (16)(16) (16) we have
Setting the derivative of EquationEquation (30)(30) (30) w.r.t equal to zero implies that
Hence
The estimate of follows as below. From EquationEquation (15)(15) (15) we have
hence
Equating Equation (45) to zero leads us to finding the MLE of ,that is,
This implies that
If we restrict the transition probability only by the conditions and for all and , then Hamilton (Citation1990) showed that the MLEs for the transition probability is given by
This implies that the estimated is essentially the number of times state seems to have been followed by state divided by the number of times the process was in state .
2.3 Root Mean Square Error(RMSE)
A Root Mean Square Error (RMSE) is computed in order to compare the two models in terms of prediction of the option prices. It is computed by utilizing the model predicted option prices and the observed option prices. That is, the RMSE of a prediction model with respect to the observed option prices is defined as the square root of the mean squared error.
where is the predicted option prices, the observed option prices and is the total number of the observations. Note that, lower values of the RMSE is an indication of a model with a better prediction.
3 Results and discussion
3.1 Empirical data
This study utilized the daily closing prices data as reported in the Russell 2000 (RUT) and Facebook (FB) indices for the period 4 January 2016 to 29 October 2021. The daily indices returns, , are continuously compounded as where and are the daily closing indices prices at day and , respectively. Also, the call options data from the two markets expiring in 34 and 60 days for RUT and FB indices, respectively, is utilized in pricing the European options.
3.2 Descriptive statistics
Table presents the basic statistics for daily RUT and FB indices returns. The mean daily returns are 0.000497 & 0.000785 for the Russell 2000 and Facebook indices, respectively, whereas the standard deviations are 0.015011 and 0.020547 for RUT and FB indices, respectively. The positive mean is an indication that investors in these markets have realized positive returns on the investment. Both indices have negative skewness and an excess kurtosis greater than 3, the value for normal distribution. This is an indication that the distribution of the indices returns has a thicker tail and is left skewed as compared with the normal distribution and that this tail thickness may be as a result of the temporal fluctuations in volatility. According to Jarque and Bera (Citation1980) if a sample data is drawn from a normal distribution, then its expected skewness and excess kurtosis are zero. It is thus clear that this data is not drawn from a normal distribution as it can be inferred from the descriptive statistics and which is further confirmed by the Jarque-Bera(JB) test statistics in Table which are significant at 1% level of significance. That is, the Jarque-Bera(JB) test statistics rejects the null hypothesis of normality for both indices returns. Moreover, the correlation coefficient for the two indices return series is 0.4845 which is an indication of a moderate positive relationship between the two series and this means that results emanating from an analysis of the two return series can be compared.
The time series data is further described by plotting the empirical density versus the normal distribution as well as the quantile-quantile(qq) plots for the two indices returns as shown in Figures Clearly, the plots show that the empirical densities are sharp peaked and are left skewed, that is, they are different from the normal distribution and hence it can be concluded that the data are not normally distributed. This finding is in an agreement with the report by the values of skewness and excess kurtosis from Table which are significantly not equal to those of normal distribution of zero and three.
3.3 Empirical findings and discussion
The plots for the RUT and FB stock indices prices and indices returns are presented in Figures , respectively. It is evidenced by Figure that the stock price series have trend implying that the mean and variances of the data are non-constant and as a consequence the data is not stationary.
The returns series plots on Figure show the common properties of time series data such as the volatility clustering, leverage effects, leptokurtic distribution and heavy tails. In this connection, the Facebook index returns are characterized by pronounced alternating long spikes compared to the Russell 2000 index returns and these alternating long and short spikes is evidence of volatility clustering in the underlying assets. The non-stationarity of the data calls for the series differencing in order to remove the series dependence before further analysis of the data. The data are therefore differenced once and subjected to statistical tests among them the Ljung-Box, Langrange Multiplier, ARCH effects, Jarque-Bera(JB) and ADF tests in order to confirm that it is stationary and to attain desirable results. The results of these Tests are presented in Table . The Ljung-Box test value is significant at 1% significance level for both series which implies that the null hypothesis that the series are not stationary is rejected and hence concludes that the return series are stationary. Moreover, the Ljung-Box and ADF tests imply that the return series have no auto correlations and ARCH effects as well as no unit roots.
Table reports the maximum likelihood parameter estimates for the two indices which are computed according to EquationEquation 28(28) (28) and it is clear that the parameters are more or less similar across the two return series. The mean and standard deviations are different across the regimes for the two indices, for instance, the mean in regime 1 is 0.00104 and 0.00163 for Russell 2000 and Facebook indices, respectively. Moreover, the mean in high-volatility regime is negative for both indices whereas it is positive in the low-volatility regime, that is is positive. It is evident that in regime 1, the assets has higher returns and low risk compared to lower returns and high risk in regime 2. On the other hand, the probability of switching from the low-volatility regime to the high-volatility regime, estimated at 0.0651 for Russell 2000 and at 0.2087 for Facebook, is significantly different for the two indices. Furthermore, the estimates indicate that the volatility persistence in the low-volatility regime is less in Facebook index than in Russell 2000 index. That is, the probability of switching regime from low-volatility to high-volatility regime is high in Facebook index as compared to Russell 2000 index. In fact, the duration of stay in each regime is computed using the expression and thus once in regime 1, the process stays there for approximately 15.35 days for Russell 2000 index compared to 4.79 days for Facebook index. In both indices, the process has the tendency of staying in high-volatility regime longer than in low-volatility regime, that is, the process approximately stays in high volatility regime for 58.43 and 16.83 days for Russell 2000 and Facebook indices, respectively.
The model log likelihood is higher for Russell 2000 index returns than that of Facebook index returns (estimated at 4404.76 and 3862.082 for RUT and FB returns, respectively), whereas the model AIC and BIC are small for RUT and FB indices returns, respectively. This implies that the RS) model provides a reasonable fit to the Russell 2000 index returns compared with the Facebook index returns. Moreover, Figures show the smoothed probabilities for the Russell 2000 and Facebook indices returns and clearly the volatility process stays in regime 2 longer most of the periods than in regime 1. That is, although the two regimes are highly persistent, regime two is more persistent than regime one in the two market indices and is expected to continue for approximately 58 and 17 days for RUT and FB indices, respectively, before exiting to regime 1. Furthermore, a comparison of the two markets shows that the regimes dynamics in RUT and FB indices appears rather different.
The results presented in Table for the parameter estimates of the Facebook and Russell 2000 indices returns are implemented in EquationEquation 28(28) (28) to calculate the European call option prices for some strike prices obtained from the Facebook and Russell 2000 stock indices. Table presents the calculated option prices for 34 and 60 days call options for Russell 2000 and Facebook indices, respectively. This research assumes a risk free-free interest rate of 6% per annum and the initial stock price is 2277.1 and 330.64 for Russell 2000 and Facebook indices, respectively. The calculated call option prices from both the Black-Scholes and RS models are more or less similar for the two market indices. However, the calculated prices from the Black-Scholes (BS) model are slightly closest to the market prices across the two markets. This implies that the BS model gives better results for the Russell 2000 and Facebook indices data sets as compared with the RS model. In addition, the plot of the strike prices versus the estimated call options and the observed option prices in the market indices is presented in Figure . The plots indicates similar curve patterns for the two stock indices with the estimated price curve from Black-Scholes model being closer to the observed market price curve than the curve for the RS model. This may imply that the Black-Scholes model gives a better estimated European call options than the RS model for the utilized data sets.
Furthermore, the estimated European call option prices from the Black-Scholes and RS models and the observed market option prices are used to compute the Root Mean Square Error (RMSE) in order to establish which model gives better results for pricing the European options. The results of the Root Mean Square Error (RMSE) are presented in Table and it can be seen that the Black-Scholes model has smaller values than the RS model across the two data sets. It is thus implied that the Black-Scholes model gives a better option price prediction than the RS model.
4 Conclusion
A RS model for pricing European call options in a financial market that exhibits structural changes with time is derived in this paper. The model is formulated based on the basis that the underlying asset dynamics are described by a geometric Brownian motion that is modulated by a continuous-time Markov chain with switching regimes. That is, this model allows the volatility to follow a RS process that is governed by a Markov chain process and it is implemented by considering two regimes, that is, the low and high volatility regimes. The model parameters are estimated by utilizing data from the Russell 2000 and Facebook indices and these model parameter estimates are then utilized to compute the European call option prices. The European call option prices reported by the RS model are compared with those that are reported by the Black-Scholes model. Moreover, a comparison of the two models is done in order to establish the model that gives better results in terms of predicting the European call option prices by computing a Root Mean Square Error for each model. This study finds out that across the two market indices, regime 1 is characterized by low risk and high returns whereas regime 2 is characterized by high risk and low returns. The transition probabilities imply that the probability of switching from the low-volatility regime to the high-volatility regime is significantly different across the two market indices and it is higher than the probability of reverting back. That is, the volatility process has the tendency of staying in high-volatility regime longer than it stays in low-volatility regime across the two markets. However, the volatility persistence in the low-volatility regime is less in Facebook index than in Russell 2000 index, that is, the probability of switching regimes from low-volatility to high-volatility regime and vice versa is high in Facebook index as compared to Russell 2000 index. The calculated European call option prices from both the Black-Scholes and RS models are more or less similar for the two market indices. However, the price values reported by the Black-Scholes (BS) model are a bit closer to the market price than that of RS model across the two markets. This implies that the BS model slightly gives better results for the Russell 2000 and Facebook indices data sets as compared with the RS model.
This study focused on developing a RS model for pricing options when the underlying asset dynamics depend on the market regimes. The model formulation is based on a geometric Brownian motion that is governed by a continuous-time Markov chain with two states, and the implementation of the model is done by computing the European call option prices for some chosen stock market indices. Future studies can focus on implementing this model for a case of more than two regimes. Moreover, future studies can consider the application of this model in the American options pricing and draw its performance comparison with that of the Black-Scholes model for both long and short-dated data.
Disclosure statement
No potential conflict of interest was reported by the authors.
Additional information
Notes on contributors
Sebastian Kaweto Kalovwe
Sebastian Kaweto Kalovwe, is a Doctor of Philosophy in Financial Mathematics from the University of Nairobi-Kenya. He also holds Msc in Mathematics (Statistics) from the Catholic University of Eastern Africa. Previously, he has worked as part-time lecturer (for five years) in the school of pure and applied sciences, South Eastern Kenya University (SEKU). Currently, he is a teacher of Mathematics and Physics at Our Lady of Mercy-Mutomo Girls Secondary School, and a part-time lecturer at Strathmore University. His research interest include; applications of stochastic time series models in modeling financial market data, regime switching modeling and Options pricing under regime switching model. Dr. Kalovwe has vast experience in teaching mathematics and research, and has published articles in different peer-reviewed and reputed journals.
References
- Barone-Adesi, G., Engle, R. F., & Mancini, L. (2008). A garch option pricing model with filtered historical simulation. The Review of Financial Studies, 21(3), 1223–19. https://doi.org/10.1093/rfs/hhn031
- Bauwens, L., & Lubrano, M. (1998). Bayesian inference on garch models using the gibbs sampler. The Econometrics Journal, 1(1), C23–46. https://doi.org/10.1111/1368-423X.11003
- Biswas, A., Goswami, A., & Overbeck, L. (2018). Option pricing in a regime switching stochastic volatility model. Statistics & Probability Letters, 138, 116–126. https://doi.org/10.1016/j.spl.2018.02.056
- Bollen, N. P. (1998). Valuing options in regime-switching models. Journal of Derivatives, 6(1), 38–50. https://doi.org/10.3905/jod.1998.408011
- Bollen, N. P., Gray, S. F., & Whaley, R. E. (2000). Regime switching in foreign exchange rates: Evidence from currency option prices. Journal of Econometrics, 94(1–2), 239–276. https://doi.org/10.1016/S0304-4076(99)00022-6
- Bollerslev, T. (1986). Generalized autoregressive conditional heteroskedasticity. Journal of Econometrics, 31(3), 307–327. https://doi.org/10.1016/0304-4076(86)90063-1
- Buffington, J., & Elliott, R. J. (2002). American options with regime switching. International Journal of Theoretical and Applied Finance, 5(05), 497–514. https://doi.org/10.1142/S0219024902001523
- Cai, J. (1994). A markov model of switching-regime arch. Journal of Business & Economic Statistics, 12(3), 309–316. https://doi.org/10.1080/07350015.1994.10524546
- Deelstra, G., Latouche, G., & Simon, M. (2020). On barrier option pricing by erlangization in a regime-switching model with jumps. Journal of Computational and Applied Mathematics, 371, 112606. https://doi.org/10.1016/j.cam.2019.112606
- Dieobold, F. X. (1986). Modeling the persistence of conditional variances: A comment. Econometric Reviews, 5(1), 51–56. https://doi.org/10.1080/07474938608800096
- Duan, J. C. (1995). The garch option pricing model. Mathematical Finance, 5(1), 13–32. https://doi.org/10.1111/j.1467-9965.1995.tb00099.x
- Duan, J. C., Popova, I., & Ritchken, P. (2002). Option pricing under regime switching. Quantitative Finance, 2(2), 116. https://doi.org/10.1088/1469-7688/2/2/303
- Duan, J. C., & Zhang, H. (2001). Pricing hang seng index options around the asian financial crisis–a garch approach. Journal of Banking & Finance, 25(11), 1989–2014. https://doi.org/10.1016/S0378-4266(00)00166-7
- Engle, R. F. (1982). Autoregressive conditional heteroscedasticity with estimates of the variance of united kingdom inflation. Econometrica: Journal of the Econometric Society, 50(4), 987–1007. https://doi.org/10.2307/1912773
- Godin, F., & Trottier, D. A. (2019). Option pricing under regime-switching models: Novel approaches removing path-dependence. Insurance, Mathematics & Economics, 87, 130–142.
- Godin, F., & Trottier, D. A. (2021). Option pricing in regime-switching frameworks with the extended girsanov principle. Mathematics and Economics.
- Gray, S. F. (1996). Modeling the conditional distribution of interest rates as a regime-switching process. Journal of Financial Economics, 42(1), 27–62. https://doi.org/10.1016/0304-405X(96)00875-6
- Hamilton, J. D. (1989). A new approach to the economic analysis of nonstationary time series and the business cycle. Econometrica: Journal of the Econometric Society, 57(2), 357–384. https://doi.org/10.2307/1912559
- Hamilton, J. D. (1990). Analysis of time series subject to changes in regime. Journal of Econometrics, 45(1–2), 39–70. https://doi.org/10.1016/0304-4076(90)90093-9
- Hamilton, J. D., & Susmel, R. (1994). Autoregressive conditional heteroskedasticity and changes in regime. Journal of Econometrics, 64(1–2), 307–333. https://doi.org/10.1016/0304-4076(94)90067-1
- Hardy, M. R. (2001). A regime-switching model of long-term stock returns. North American Actuarial Journal, 5(2), 41–53. https://doi.org/10.1080/10920277.2001.10595984
- Jarque, C. M., & Bera, A. K. (1980). Efficient tests for normality, homoscedasticity and serial independence of regression residuals. Economics Letters, 6(3), 255–259. https://doi.org/10.1016/0165-1765(80)90024-5
- Kirkby, J. L., & Nguyen, D. (2020). Efficient asian option pricing under regime switching jump diffusions and stochastic volatility models. Annals of Finance, 16(3), 307–351. https://doi.org/10.1007/s10436-020-00366-0
- Lamoureux, C. G., & Lastrapes, W. D. (1990). Persistence in variance, structural change, and the garch model. Journal of Business & Economic Statistics, 8(2), 225–234. https://doi.org/10.1080/07350015.1990.10509794
- Lin, S., & He, X. J. (2020). A regime switching fractional black–Scholes model and european option pricing. Communications in Nonlinear Science & Numerical Simulation, 85, 105222. https://doi.org/10.1016/j.cnsns.2020.105222
- Mitsui, H. (2000). An analysis of nikkei 225 option pricing by a garch model. Gendai-Finance, 7, 57–73.
- Mitsui, H., & Satoyoshi, K. (2010). Empirical study of nikkei 225 option with markov switching garch model. Asia-Pacific Financial Markets.
- Mitsui, H., & Watanabe, T. (2003). Bayesian analysis of garch option pricing models. Journal Japan Statistical Society(japanese Issue), 33, 307–324.
- Nasri, B. R., Remillard, B. N., & Thioub, M. Y. (2020). Goodness-of-fit for regime-switching copula models with application to option pricing. Canadian Journal of Statistics, 48(1), 79–96. https://doi.org/10.1002/cjs.11534
- Siu, T. K., Tong, H., Yang, H. (2004). On pricing derivatives under garch models: A dynamic gerber-shiu approach. North American Actuarial Journal, 8(3), 17–31. https://doi.org/10.1080/10920277.2004.12254408
- Yao, D. D., Zhang, Q., & Zhou, X. Y. (2006). A regime-switching model for european options. In Stochastic processes, optimization, and control theory: Applications in financial engineering, queueing networks, and manufacturing systems (pp. 281–300). Springer.