![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
Abstract
An option is a financial contract that can be used to reduce risks in an investment. It is widely known that a fair price of this contract depends significantly on the volatility of an underlying asset price, which may be affected differently by positive and negative information. Therefore, the fair price of option has been studied through various methods. In this research, an analytical formula for European option pricing via the TGARCH model is derived based on an Edgeworth expansion of the density of cumulative asset return. Furthermore, the accuracy of the proposed method is investigated by comparing numerical results with the GARCH model, TGARCH model, analytical approximation via the GARCH model and Monte Carlo technique. The results reveal that in the case of in the money (ITM) with the proposed method performed better than the others. The behaviour of the proposed method is also discussed.
1. Introduction
An option is a financial contract that is widely used by both investors and business companies to reduce risks in an underlying asset investment. The risks of option holders are transferred to the option seller after paying a premium. Therefore, many researchers are interested in the fair premium, such as Black and Scholes (Citation1973); Duan (Citation1995); Duan et al. (Citation1999); Harding (Citation2013).
Black and Scholes (Citation1973) proposed a model, based on a partial differential equation, for option pricing. To investigate the pricing formula, several strong assumptions are assumed. For example, the volatility of the asset return is assumed to be a constant and the asset price follows a lognormal distribution. Motivated by these limitations, Finucane (Citation1989); Hull and White (Citation1987) improved the Black-Scholes model by providing the formula without assuming a constant volatility. For a volatility depending on time, the GARCH model is often used to estimate the volatility in various asset markets. Alternatively, this model is unable to explain some feature of asset return (Cai & Li, Citation2019) such as a correlation between stock returns and changes in return volatility (Black, Citation1976). Therefore, several versions of modified GARCH models have been proposed such as the EGARCH, GJR-GARCH and TGARCH models. The literature review of GARCH and modified GARCH models are provided in the next section. For the lognormal distribution assumption of asset price, Jarrow and Rudd (Citation1982) noticed that the distribution of asset prices may not be lognormal distribution. They relaxed this condition by using an approximating distribution approach and formulated an analytical approximation of option prices with constant volatility. In 1999, Duan et al. improved the analytical formula to a stochastic volatility case via the GARCH model.
In this research, an analytical approximation for the option prices via the TGARCH model is investigated. This analytical approximation can be formulated through a combination of the first four moments of cumulative asset return. The terms of such moments are approximated by using the TGARCH model. Apart from this step, the estimated option prices can be computed by using parameters from the TGARCH model fitting. Also, the numerical performance of the proposed method is studied by comparing the performance to that of Monte Carlo simulation and other models. This study contributes to both academic knowledge and actual practice. For an academic knowledge, although the Black-Scholes formula has been widely used in options markets, it requires a number of unreasonable assumptions. This study proposes a more general formula to relax some conditions. For an empirical practice, a new efficient method for pricing options in the case of ITM is provided. This formula is applied to a real data from the stock exchange of Thailand (SET). The study results may suggest a new alternative approach for option investors.
The contents of this paper are organized into six sections. In section 2, a literatures review and information about the TGARCH model are provided. An analytical approximation of option prices is examined in section 3 and a computational part of the analytical approximation is explained in section 4. This was followed by a numerical performance and discussion section. The paper concludes with section 6.
2. Literature review
The GARCH model, introduced by Bollerslev (Citation1986), is widely used to estimate the non-constant volatility, depending on time. The GARCH model provides a good approximation for smooth and persistent change volatility (Chen et al., Citation2013). Therefore, the application of this model to estimate the volatility in time series data is widely studied in various sectors, see Agnolucci (Citation2009); Angabini and Wasiuzzama (Citation2011); Mcmillan and Speight (Citation2004). However, the GARCH model is unable to describe some features of asset return (Cai & Li, Citation2019). To overcome the limitations of GARCH models, many researchers have developed modified GARCH models. For example, Nelson (Citation1991) and Glosten et al. (Citation1993) introduced the EGARCH model and the GJR-GARCH model, respectively. The EGARCH and GJR-GARCH models aim at studying the impact of negative and positive returns on conditional volatility. Zakoian (Citation1994) also proposed the TGARCH model for purposes similar to the GJR-GARCH model, but it focuses on conditional standard deviation (volatility) instead of conditional variance, and used the TGARCH model to approximate the volatility of a stock index in France (CAC index).
Following Zakoian (Citation1994), some researchers have studied and used the TGARCH model to estimate the volatility of other underlying asset prices, such as carbon, crude oil, ethanol, natural gas, coal-three, corn and sugar; see Alberola et al. (Citation2008); Hasan et al. (Citation2013); Trujillo-Barrera et al. (Citation2012) and Godana et al. (Citation2014) Moreover, the TGARCH model is used to describe mortgage risk factors of housing price and capture the house price dynamic on the logarithmic return and to estimate the housing price volatility for the pricing reverse mortgage derivatives (Lee et al., Citation2015). For financial market research, Sabiruzzaman et al. (Citation2010) investigated the pattern of volatility in the daily trading volume index of the Hong Kong stock exchange by using two approaches; GARCH and TGARCH models. They found that the TGARCH specification is superior to the GARCH specification. Moreover, the TGARCH model is used to estimate forex market volatility. In 2014, Ghosh found evidence of significant volatility spill-over effect between financial markets and the forex market for the Indian economy. His results showed asymmetric reactions in the forex market volatility. Furthermore, the TGARCH model is also used to study arbitrage trading. In 2015, Cui et al. proposed a new statistical arbitrage model called TGARCH-WNN. In this model, the TGARCH model is applied to capture and predict the standard deviation of the price spread and wavelet neural networks (WNNs) are used to predict the trading threshold.
Let p and q be positive integers. The TGARCH model with parameter p and q is given by
(1)
(1)
where
is a sequence of independent identically distributed (i.i.d.) random variables with zero mean and unit variance,
and
are the positive and negative parts of
is the information set (σ-field) of all information up to including time t, ht is the conditional variance which is
α0 is a positive real number,
and βj are non-negative real numbers for
and
The model in (1) can be rewritten in the form
(2)
(2)
where αi and γi are non-negative real numbers for all and
Under a locally risk neutral valuation relationship (LRNVR) situation, Duan (Citation1995) used the GARCH model to estimate the volatility of the S&P100 index and then used his result to approximate the option price. Over the past few decades, many researchers, such as Harding (Citation2013); Hsieh and Ritchken (Citation2006), have applied the idea of Duan to different types of GARCH models, such as NGARCH and HN-GARCH models, respectively, to estimate the option price. In 2016, Hongwiengjan andThongtha applied the idea of Duan (Citation1995) to the TGARCH model and used it to approximate the option price of the SET50 index. The TGARCH model under LRNVR can be formulated by replacing in (2) by
The model becomes
(3)
(3)
Apart from the constant volatility, Jarrow and Rudd (Citation1982) studied the way to relax the assumption of the Black-Scholes model with regards to the lognormal distribution of stock prices. They relaxed this condition by approximating a non-lognormal distribution using a series expansion technique, the Edgeworth expansion. They also constructed an analytical formula, with constant volatility, for option prices. Their formula is similar to a Black-Scholes formula adjusted by skewness and kurtosis of the cumulative asset return. In 1999, Duan et al. used Jarrow & Rudd’s idea to developed an analytical approximation with stochastic volatility via the GARCH model. The performance of this approximation was studied by comparing the obtained results with the prices from the Monte Carlo simulation. The pricing error of their analytical approximation was small, with respect to those from Monte Carlo technique, for all cases depending on time to maturity date. Following this research, Duan et al. (Citation2006) provided analytical formulae via the EGARCH and GJR-GARCH models. Furthermore, Huang et al. (Citation2017) derived the option pricing formula via the Realized GARCH model with more terms in the Edgeworth expansion.
3. An analytical approximation of option prices
Duan et al. (Citation1999) developed an analytical approximation for option pricing by computing an option price as a Black-Scholes formula adjusted by skewness and kurtosis of cumulative asset return with GARCH model under LRNVR. Let ST be an underlying asset price at expiration date and St be the underlying asset price at time t. Define
(4)
(4)
and let
where
Q is a pricing measure satisfying the LRNVR condition,
and
are the conditional expectation and variance on
under probability measure Q, respectively. Dual et al. established an analytical approximation of the option price (
) using the Edgeworth expansion around the normal density function. The formula is
(5)
(5)
K is the strike price, r is the risk free interest rate, and and
are the density function and cumulative distribution function of a standard normal random variable, respectively. The analytical approximation of the option price can be computed by estimating
The next section describes the calculation of these values.
4. Computing the analytical approximation of option prices
In computing the analytical approximation of the option price by using (5), the values of κ3 and κ4 are needed. These values can be computed directly by using the first four moments of the asset return under LRNVR. Duan et al. (Citation1999) provided the formula of the first four moments of ρT under LRNVR as
(6)
(6)
where for
Expanding the right side of (6), each term is in the form of
where γ is a constant,
The expression terms of the first four moments are shown in (Duan et al., Citation1999) and are evaluated into various expression terms. The above summation can be rewritten as
(7)
(7)
In this equation, some expression terms can be computed directly using the TGARCH model under LRNVR (shown in Appendix A). However, for the other expression terms, it is hard to formulate to their explicit forms. In this case, the Taylor expansion is used to approximate such expression terms. In the first case, we compute the terms by transforming ht into TGARCH form under LRNVR with The model can be rewritten as
(8)
(8)
(9)
(9)
where
(10)
(10)
and
is a standard normal random variable. After that, a general form of some expressions in (7) is derived by using (9). However, some random variables, such as Yt and
are dependent. Then, the expectations of the product of Yt and
are formulated in the following proposition.
Proposition 4.1.
Let and
Then
Proof.
See Appendix B. □
Some expressions in (7) cannot be formulated into explicit forms, such as the expression terms ending with 2, power of ht. These expressions can be approximated by using the Taylor expansion. Expanding the Taylor expansion of
around
we have
Therefore,
(11)
(11)
From the above two cases, the expression in (7) can be approximated. Therefore, the approximate option prices in (5) can be computed. The results of this estimation are provided in the next section, and a numerical performance is also studied.
5. Numerical performance and discussion
From the previous section, we note that the analytical approximation of option prices depends on the value of the first four moments of asset return. Therefore, in this section, our study focuses on two main points. In section 5.1, the behaviours of the first four moments of cumulative asset return and the analytical approximation of option prices are investigated. The results of these values are computed by using set data. For the second point, the option prices estimated by using the analytical formula and those obtained from other methods, such as Monte Carlo simulation and the GARCH model, are compared. The accuracy of such techniques with the real data is investigated in section 5.2. The SET data is used as observational data. For the Monte Carlo simulation, the 500,000 samples are generated with 500,000 paths simulation.
5.1. Simulation results
In this section, the behaviour of the approximated formula is studied in different directions. First, mean, variance, skewness and kurtosis of the cumulative return are computed with the TGARCH model by using (6). Second, the approximated option prices with skewness and kurtosis are investigated in various cases. In this study, some parameters and data are set. From (5), the values of and St are needed. These following values are specified based on the value in real market. The values of strike price (K) and time to maturity date (T – t) are set in different cases, which are 900, 1000, 1100 and 10, 30, 90, 150, respectively. The risk-free interest rate (r) is set to be 2% as the government bond yield rate, and the stock price at the current period (St) is 1000. The unconditional variance
is also required and then set, as in Eriksson (Eriksson, Citation2013):
(12)
(12)
whenever
Considering the results in section 5.2, the 100 sets of parameters of the TGARCH model in (8) are randomly selected and associated to the condition in (12) and the following conditions:
Then, the mean and variance of the cumulative asset return can be computed. The standardized version of mean and variance is defined as and
where
with
In , the average of the standardized mean and variance of cumulative asset return from the simulation results are shown when the unconditional variance is set to be 20% below and above h1. The result shows that the values of mean and variance of cumulative asset return, obtained by using (6), are not significantly different from those obtained by Monte Carlo simulation.
Table 1. The average of standardized mean and variance of the cumulative asset return from the simulation results.
For the skewness and kurtosis, their standardized version are defined as and
respectively. Note that we need to compute the third and fourth moments. Some expression terms in (7) can be computed directly, as shown in Appendix A. The others are approximated by using the Taylor expansion given in (11). As can be seen in , the results obtained by using analytical technique are significantly different from those received by Monte Carlo simulation, especially when the time to maturity increases.
Table 2. The average of standardized skewness and kurtosis of the cumulative asset return from the simulation results.
Next, we want to study the behaviour of an analytical approximation of call option prices with skewness and kurtosis. In this study, the moneyness ratio () is set to be 0.90, 1.00 and 1.10. In , the results from an analytical approximation prices and Monte Carlo simulation have the same direction, especially when the time to maturity increases.
Figure 1. Average of simulation results of an analytical approximation for call option prices. Source: Created by the authors.
![Figure 1. Average of simulation results of an analytical approximation for call option prices. Source: Created by the authors.](/cms/asset/89e6b648-a6aa-425c-a7ec-e7fb9df71ff5/rero_a_1805636_f0001_c.jpg)
From , the results of ATM and OTM cases from the proposed method appear unreasonable. In this simulation, the strike price K is set to be 900, 1000 and 1100 and the current stock price is 1000. However, the results of the ATM case is negative while those of the OTM case seem to be explosive when the time to maturity increases. This shows that the proposed method cannot estimate the option prices for both ATM and OTM cases. For the ITM case, the results from both methods show a similar pattern. Therefore, to deepen this study, more ITM cases, with and 0.95, are studied. From , the results show that the behaviour of two option prices have the same trend when
=0.80 and 0.85. However, the analytical approximation prices are too high when
=0.95. Therefore, the results from these two methods seem to be uncorrelated when
increases to 1.
Figure 2. Average of simulation results of an analytical approximation for call option prices: ‘ITM Cases’. Source: Created by the authors.
![Figure 2. Average of simulation results of an analytical approximation for call option prices: ‘ITM Cases’. Source: Created by the authors.](/cms/asset/8f5814e3-49db-4e44-9d30-099e970147c2/rero_a_1805636_f0002_c.jpg)
After getting the analytical option prices and analysing the behaviour of the approximated formula in comparison to the Monte Carlo approach, the impact of skewness and kurtosis is considered. In this point, the formulas are set into three cases as
where
and A4 are given in (5). Note that the Pricenc does not focus on the impact of both skewness and kurtosis while the Pricenc adjusted only the skewness is defined as the Pricesc. For the Priceskc, it is the formula in (5) that is the Pricenc adjusted by both skewness and kurtosis. One case of simulation result is shown in . The results for other cases show the same trend.
Figure 3. Impact of skewness and kurtosis correspondence on call option prices. Source: Created by the authors.
![Figure 3. Impact of skewness and kurtosis correspondence on call option prices. Source: Created by the authors.](/cms/asset/f542af56-535c-461e-a37b-24632b0afae3/rero_a_1805636_f0003_c.jpg)
It can be seen in that the kurtosis affects the analytical approximation of option prices extremely because the in (5) is so small. For the skewness, it does not provide a significant impact on the approximation.
5.2. Approximated option prices in real market
In this section, the analytical formula of option prices is applied to the real SET data. The accuracy of the analytical approximation is investigated by using four different measures in the case of ITM. In this study, the risk-free interest rate is set to be 1.99% per annum as the Thai government bond rate. The time to expiration (T – t) is considered to three cases: 30, 60 and 90 days. For other information, such as the asset prices at time t (St) and strike prices (K), the SET data are used as observation data. To approximate the option prices, we first compute the parameter values and λ in (8). In this step, the maximum likelihood is used to estimate these appropriate parameters. The results of the new analytical approximation of call option prices are compared with the observation prices and prices obtained from the other methods such as the GARCH model, TGARCH model, GJR-GARCH model, analytical approximation via GARCH model and Monte Carlo simulation. The tools used to measure the accuracy of the results are correlation coefficient (CORR), mean absolute percentage Error (MAPE), root mean square error (RMSE) and maximum of absolute percentage error (MaxAPE). The formulas are given as follows:
where the yi is the observed price, xi is the estimated price,
is the mean of observed price,
is the mean of estimated price and n is the number of observations. The results in reveal that, in the case that
the proposed method provides an excellent approximated price. From the CORR measure, the results from these six methods have a high positive relation to the observation prices. From MAPE, RMSE and MaxAPE, the proposed method gives more accuracy than the other approaches. This conclusion satisfies the result from the previous section. The results are shown in and .
Figure 4. CORR, MAPE, RMSE and MaxAPE of call option prices. Source: Created by the authors based on SET50 data.
![Figure 4. CORR, MAPE, RMSE and MaxAPE of call option prices. Source: Created by the authors based on SET50 data.](/cms/asset/8b5aea3e-98d4-4deb-84c8-76195c9284f3/rero_a_1805636_f0004_c.jpg)
Table 3. The performance of option prices with
For the 11 ITM cases when approaches to 1, the results of the proposed method do not show a good performance, as can be seen in .
Figure 5. Impact of correspondence on call option prices. Source: Created by the authors based on SET50 data.
![Figure 5. Impact of K/St correspondence on call option prices. Source: Created by the authors based on SET50 data.](/cms/asset/395fa0bf-8132-4e5b-80e8-acbeec9e6bb3/rero_a_1805636_f0005_c.jpg)
However, after computing the 27 OTM cases with real-data option prices, the results of the proposed method reveal a poor performance. This method cannot capture the option prices. In most cases, the estimated prices are negative and then are set as 0.
6. Conclusion
An analytical approximation formula via the TGARCH model has been developed in this study. The formula is similar to the Black-Scholes formula, adjusted by the skewness and kurtosis of the cumulative asset return. The first four moments of the cumulative asset return are derived and computed by using the TGARCH model. The new proposed method generalises the Black-Scholes by reducing two assumptions: constant volatility and the lognormal distribution asset prices. After formulating the proposed method, both simulation and empirical results are studied. With these results, it is shown that in the case of ITM with the behaviour of the analytical approximation via the TGARCH model does not significantly differ from the compared methods, such as GARCH, TGARCH, GJR-GARCH, analytical approximation via the GARCH model and Monte Carlo simulation. However, with regards to accuracy, the empirical results demonstrate that the analytical approximation via the TGARCH model is more accurate than the other methods. Furthermore, the analytical approximation via the TGARCH model required fewer assumptions than the Black-Scholes formula. Future research can be devoted to improving a new method with more general assumptions in the real market.
Acknowledgment
The authors would like to thank the referees for their careful reading and helpful comments.
Disclosure statement
No potential conflict of interest was reported by the authors.
Additional information
Funding
References
- Agnolucci, P. (2009). Volatility in crude oil futures: A comparison of the predictive ability of GARCH and implied volatility models. Energy Economics, 31(2), 316–321. https://doi.org/https://doi.org/10.1016/j.eneco.2008.11.001
- Alberola, E., Chevallier, J., & Chèze, B. (2008). The EU emissions trading scheme: The effect of industrial production and CO2 emissions on carbon prices. International Economics, 116, 93–126.
- Angabini, A., & Wasiuzzama, S. (2011). GARCH models and the financial crisis-A study of the Malaysian stock market. The International Journal of Applied Economics and Finance, 5(3), 226–236. https://doi.org/https://doi.org/10.3923/ijaef.2011.226.236
- Black, F. (1976). Studies of stock market volatility changes. Proceedings of the American Statistical Association, Business and Economic Statistics Section, 177–181.
- Black, F., & Scholes, M. (1973). The pricing of options and corporate liabilities. Journal of Political Economy, 81(3), 637–654. https://doi.org/https://doi.org/10.1086/260062
- Bollerslev, T. (1986). Generalized autoregressive conditional heteroskedasticity. Journal of Econometrics, 31(3), 307–327. https://doi.org/https://doi.org/10.1016/0304-4076(86)90063-1
- Cai, Y., & Li, G. (2019). A quantile function approach to the distribution of financial returns following TGARCH models. Statistical Modelling. https://doi.org/https://doi.org/10.1177/1471082X19876371
- Chen, C. W. S., Lin, E. M. H., & Lin, Y. R. (2013). Bayesian perspective on mixed GARCH models with jumps. Uncertainty Analysis in Econometrics with Applications, 31, 141–154.
- Cui, L., Huang, K., & Cai, K. J. (2015). Application of a TGARCH-wavelet neural network to arbitrage trading in the metal futures market in China. Quantitative Finance, 15(2), 371–384. https://doi.org/https://doi.org/10.1080/14697688.2013.819987
- Duan, J. C. (1995). The GARCH option pricing model. Mathematical Finance, 5(1), 13–32. https://doi.org/https://doi.org/10.1111/j.1467-9965.1995.tb00099.x
- Duan, J. C., Gauthier, G., & Simonato, J. G. (1999). An analytical approximation for the GARCH option pricing model. The Journal of Computational Finance, 2(4), 75–116. https://doi.org/https://doi.org/10.21314/JCF.1999.033
- Duan, J. C., Gauthier, G., Simonato, J. G., & Sasseville, C. (2006). Approximating the GJR-GARCH and EGARCH option pricing models analytically. The Journal of Computational Finance, 9(3), 41–69. https://doi.org/https://doi.org/10.21314/JCF.2006.156
- Eriksson, K. (2013). On return innovation distribution in GARCH volatility modelling: Empirical evidence from the Stockholm Stock Exchange [Unpublished bachelor’s thesis]. Umeȧ University.
- Finucane, T. J. (1989). Black-Scholes approximations of call option prices with stochastic volatilities: A note. The Journal of Financial and Quantitative Analysis, 24(4), 527–532. https://doi.org/https://doi.org/10.2307/2330984
- Glosten, L. R., Jagannathan, R., & Runkle, D. E. (1993). On the relation between the expected value and the volatility of the nominal excess return on stocks. The Journal of Finance, 48(5), 1779–1801. https://doi.org/https://doi.org/10.1111/j.1540-6261.1993.tb05128.x
- Godana, A. A., Ashebir, Y. A., & Yirtaw, T. G. (2014). Statistical analysis of domestic price volatility of sugar in Ethiopia. American Journal of Theoretical and Applied Statistics, 3(6), 177–183. https://doi.org/https://doi.org/10.11648/j.ajtas.20140306.12
- Harding, D. (2013). The performance of GARCH option pricing models - An empirical study on Swedish OMXS30 Call Options [Unpublished master’s thesis]. Jönköping International Business School.
- Hasan, M. Z., Akhter, S., & Rabbi, F. (2013). Asymmetry and persistence of energy price volatility. International Journal of Finance and Accounting, 7(2), 373–378.
- Hongwiengjan, W., & Thongtha, D. (Eds.). (2016, May). Call Option Pricing via TGARCH Approach: The Case Study in SET50. AMM&APAM 2016. Proceedings of the 21st Annual Meeting in Mathematics (AMM 2016) and Annual Pure and Applied Mathematics Conference 2016 (APAM 2016) (pp. 217–226). Bangkok, Thailand: Chulalongkorn University.
- Huang, Z., Wang, T., & Hansen, P. R. (2017). An analytical approximation approach: Option pricing with the realized GARCH model. Journal of Futures Markets, 37(4), 328–358. https://doi.org/https://doi.org/10.1002/fut.21821
- Hull, J., & White, A. (1987). The pricing of options on assets with stochastic volatilities. The Journal of Finance, 42(2), 281–300. https://doi.org/https://doi.org/10.1111/j.1540-6261.1987.tb02568.x
- Hsieh, K. C., & Ritchken, P. (2006). An empirical comparison of GARCH option pricing models. Review of Derivatives Research, 8(3), 129–150. https://doi.org/https://doi.org/10.1007/s11147-006-9001-3
- Jarrow, R., & Rudd, A. (1982). Approximate option valuation for arbitrary stochastic processes. Journal of Financial Economics, 10(3), 347–369. https://doi.org/https://doi.org/10.1016/0304-405X(82)90007-1
- Lee, C. C., Chen, K., & Shyu, D. S. (2015). A study on the main determination of mortgage risk: Evidence from reverse mortgage markets. International Journal of Financial Research, 6(2), 84–94. https://doi.org/https://doi.org/10.5430/ijfr.v6n2p84
- Mcmillan, D. G., & Speight, A. E. H. (2004). Daily volatility forecasts: Reassessing the performance of GARCH models. Journal of Forecasting, 23(6), 449–460. https://doi.org/https://doi.org/10.1002/for.926
- Nelson, D. B. (1991). Conditional heteroskedasticity in asset returns: A new approach. Econometrica, 59(2), 347–370. https://doi.org/https://doi.org/10.2307/2938260
- Sabiruzzaman, M., Monimul Huq, M., Beg, R. A., & Anwar, S. (2010). Modelling and forecasting trading volume index: GARCH versus TGARCH approach. The Quarterly Review of Economics and Finance, 50(2), 141–145. https://doi.org/https://doi.org/10.1016/j.qref.2009.11.006
- Trujillo-Barrera, A., Mallory, M., & Garcia, P. (2012). Volatility spillovers in U.S. crude oil, ethanol, and corn futures markets. International Journal of Financial Research, 37(2), 247–262.
- Zakoian, J. M. (1994). Threshold heteroskedastic models. Journal of Economic Dynamics and Control, 18(5), 931–955. https://doi.org/https://doi.org/10.1016/0165-1889(94)90039-6
Appendix A:
Some general form of terms in the first four moments of cumulative asset return
This appendix is reserved for providing general forms of some expression terms in (7). We can derive these terms by
A.1 Some terms ending with ![](//:0)
![](//:0)
If q = 0, we get
If q = 1, we get
If q = 2, we get
A.2 Some terms ending with ![](//:0)
![](//:0)
If q = 0, we get
If q = 1, we get
If q = 2, we get
If q = 3, we get
A.3 Some terms ending with ![](//:0)
![](//:0)
If s = 0, we get
If s = 1, we get
A.4 Some terms ending with ![](//:0)
![](//:0)
If s = 0, we get
If s = 1, we get
If s = 2, we get
A.5 Some terms ending with ![](//:0)
![](//:0)
If u = 0, we get
If u = 1, we get
A.6 Some terms ending with ![](//:0)
![](//:0)
If q = 0, we get
If q = 1, we get
If q = 2, we get
A.7 Some terms ending with ![](//:0)
![](//:0)
If s = 0, we get
If s = 1, we get
Appendix B:
Proof of Proposition 4.1
Proof.
It can be derived by using an integration by parts that
With this facts and fact that
where I is an indicator function, we have
The proof can be completed by applying the above equations, the fact that when n is odd,
when k is even to the equations: