Next Article in Journal
Adaptive Fixed-Time Neural Network Tracking Control of Nonlinear Interconnected Systems
Next Article in Special Issue
Do Co-Worker Networks Increase or Decrease Productivity Differences?
Previous Article in Journal
On the Statistical Discrepancy and Affinity of Priority Vector Heuristics in Pairwise-Comparison-Based Methods
Previous Article in Special Issue
Research on Financial Systemic Risk in ASEAN Region
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

How to Promote the Performance of Parametric Volatility Forecasts in the Stock Market? A Neural Networks Approach

School of Finance, Qilu University of Technology, No. 3501, Daxue Road, Changqing District, Jinan 250353, China
Entropy 2021, 23(9), 1151; https://doi.org/10.3390/e23091151
Submission received: 28 July 2021 / Revised: 30 August 2021 / Accepted: 30 August 2021 / Published: 1 September 2021
(This article belongs to the Special Issue Structures and Dynamics of Economic Complex Networks)

Abstract

:
This study uses the fourteen stock indices as the sample and then utilizes eight parametric volatility forecasting models and eight composed volatility forecasting models to explore whether the neural network approach and the settings of leverage effect and non-normal return distribution can promote the performance of volatility forecasting, and which one of the sixteen models possesses the best volatility forecasting performance. The eight parametric volatility forecasts models are composed of the generalized autoregressive conditional heteroskedasticity (GARCH) or GJR-GARCH volatility specification combining with the normal, Student’s t, skewed Student’s t, and generalized skewed Student’s t distributions. Empirical results show that, the performance for the composed volatility forecasting approach is significantly superior to that for the parametric volatility forecasting approach. Furthermore, the GJR-GARCH volatility specification has better performance than the GARCH one. In addition, the non-normal distribution does not have better forecasting performance than the normal distribution. In addition, the GJR-GARCH model combined with both the normal distribution and a neural network approach has the best performance of volatility forecasting among sixteen models. Thus, a neural network approach significantly promotes the performance of volatility forecasting. On the other hand, the setting of leverage effect can encourage the performance of volatility forecasting whereas the setting of non-normal distribution cannot.

1. Introduction

Volatility is a statistical measure of the dispersion of returns for a given asset. A higher volatility means that an asset’s price can change dramatically over a short time period in either direction, and thus is expected to be less predictable. On the other hand, a lower volatility means that an asset’s price does not fluctuate dramatically, and then tends to be more steady (volatility is often measured by either the standard deviation or variance between returns from that same asset). Hence, volatility can be used to measure the amount of uncertainty or risk related to the size of changes in an asset’s price, and it obeys the criteria: ‘the higher the volatility and then the riskier the asset’. Because of the above property for volatility, volatility is usually used in asset allocation [1,2,3], option pricing [4,5], risk management [6,7,8,9] and hedge strategy [10,11]. Thus, how to accurately predict the volatility of an asset is a very important issue in the actual investment process in the financial field. As to the issue of volatility forecasting, most of literatures used the generalized autoregressive conditional heteroskedasticity (GARCH) family models, a parametric volatility forecasting approach, to predict the volatility of an asset [12,13,14,15,16,17,18,19]. Because this type of model can capture most common features of financial assets such as both the linear dependence and strong autoregressive conditional heteroskedasticity (ARCH) effect subsisting on the return series, and both the volatility clustering and leverage effect usually existing at the volatility of financial asset returns series [8,20,21,22,23] (Volatility clustering means that large changes tend to be followed by large changes, of either sign, and small changes tend to be followed by small changes [23]. Conversely, the leverage effect is the extra increase of volatility caused by the bad news [20]. Notably, the volatility clustering and leverage effect appear significantly for financial assets, especially for the stock market [8,22]).
Regarding past literatures about volatility forecasting, they tried to use a more flexible model or a complex approach to increase the performance of volatility forecasting [15,16,17,19]. For example, Aliyev et al. [12] used the GARCH, exponential GARCH (EGARCH) and the GJR-GARCH model of Glosten, Jagannathan and Runkle [24] models with normal distribution (i.e., the GARCH-N, EGARCH-N and GJR-GARCH-N) to estimate the volatility of the Nasdaq-100. Chun et al. [14] used the above three models, the GARCH-N, EGARCH-N and GJR-GARCH-N, to forecast the volatility in the KOSPI200 in the Korean market. They found that the GJR-GARCH model has the best performance of volatility forecasting among the above three GARCH-family models. Lee and Pai [15] used the GARCH with the normal, skewed Student’s t and skewed generalized error distributions (SGED) (i.e., the GARCH-N, GARCH-ST and GARCH-SGED) to forecast the volatility of REIT in the United States. They found that the GARCH-SGED model is superior to the GARCH-N and GARCH-ST models. Liu and Hung [16] applied the GARCH-N, GARCH-T, GARCH-HT, GARCH-SGT, GJR-GARCH-N and EGARCH-N models to forecast the volatility of the Standard & Poor’s 100 stock index. They found that the GJR-GARCH-N model achieves the most accurate volatility forecasts, followed by the EGARCH-N model (the GARCH-T, GARCH-HT and GARCH-SGT models are the GARCH model with the Student’s t, heavy-tailed and skewed generalized Student’s t distributions, respectively). Lv and Shan [17] used the RiskMetrics, GARCH, IGARCH, GJR-GARCH, EGARCH, FIGARCH and FIEGARCH with the normal and skewed Student’s t distributions to forecast the spot and futures price volatilities of natural gas. They found that the simple linear GARCH-class models overwhelmingly outperform nonlinear models in forecasting spot price volatility. However, nonlinear models are also superior to linear models in forecasting future price volatility (the simple linear GARCH-class models include the RiskMetrics, GARCH and IGARCH models whereas the nonlinear models contain the GJR-GARCH, EGARCH, FIGARCH and FIEGARCH models. Moreover, the IGARCH is the integrated GARCH model whereas the FIGARCH is the fractionally integrated GARCH model. In addition, the FIEGARCH is the fractionally integrated EGARCH model). Su [18] applied the GARCH and GJR-GARCH model with the SGED distribution (GJR-GARCH-SGED) to forecast the volatility of six stock indices such as the NYSE, Brussels, CAC40, DAX, SWISS and NIKKEI. He found that asymmetric volatility specification (GJR-GARCH model) forecasts out-of-sample volatility more accurately than symmetric volatility specification (GARCH model). From the above discussion, I discovered the following phenomena: First, most of literatures only used one sample asset such as the Nasdaq-100 [12], KOSPI200 [14], REIT [15], Standard & Poor’s 100 [16] and natural gas [17], indicating that the results may not be credible because the number of sample assets is not enough. Second, the empirical model may be a model with a flexible volatility specification [12,14], generalized return distribution [15] or with both flexible volatility specification and generalized return distribution [16,17,18]. Moreover, among the flexible volatility specifications the GJR–GARCH model has the best performance of volatility forecasting [13,14,16,18]. In addition, the GJR-GARCH-SGED model in Su [18] is the most flexible model among the models with both flexible volatility specification and generalized return distribution [16,17,18]. Third, the volatility forecasting models in the above literatures all focus on the GARCH-family with alternative distributions—a parametric approach.
Hence, in order to fill the gap of the above literatures, this study uses the fourteen stock indices as the sample assets and then utilizes a more flexible parametric volatility forecasting model combining with a neural network approach, an artificial intelligence technique, to promote the performance of volatility forecasting in the stock markets. The fourteen stock indices are the stock indices in the developed and emerging markets. This is the first contribution of this study because the number of sample assets is sufficient and the fourteen stock indices are distributed between countries with different degrees of economic development, indicating that the results obtained from this study are credible and representative in the stock market. The parametric volatility forecasting models in this study are mainly the asymmetric type of GARCH model, the GJR-GARCH of Glosten et al. [24], combined with the skewed generalized Student’s t (SGT) distribution of Theodossiou [25]. This is the second contribution of this study because the GJR-GARCH model with the SGT distribution (hereafter, the GJR-GARCH-SGT) is more flexible than the GJR-GARCH-SGED model in Su [18] and it can seize most of common features of financial assets such as, the volatility clustering and leverage effect existing on the volatility, and the skewness and kurtosis appearing on the return distribution (the GJR-GARCH model of Glosten et al. [24] can significantly seize the volatility clustering and leverage effect appearing at financial assets especially for the stock market. On the contrary, the skewed generalized Student’s t (SGT) distribution of Theodossiou [25] is able to capture the features of non-normal distribution such as the distribution of returns being left-skewed and having a larger and thicker tail than the normal distribution. In addition, the standardized SGT distribution degenerates into the standardized SGED distribution if n ). A neural network (NN) model is a multilayer perceptron, with one input layer, at least one hidden layer and one output layer, and one or more nodes within each layer. Moreover, each input node is connected to each hidden node, and in turn, each hidden node is connected to each output node (an artificial neural network is an interconnected group of nodes, inspired by a simplification of neurons in the brain. Here, each circular node represents an artificial neuron and an arrow represents a connection from the output of one artificial neuron to the input of another. Please see Figure 1 for more details). In other words, a neural network model is composed of an interconnected group of nodes. Hence, the structure of neural network models is similar to that of the network topology, and thus they both have the same concept (network topology is the topological structure of a network and may be depicted physically or logically. Physical topology is the placement of the various components of a network including several device locations and the cable installation between them. The device locations are regarded as the nodes of a network and the cable installation connects the device locations. Thus, the cable installation can be considered as the links or lines between the nodes in that network. On the other hand, logical topology illustrates how data flows within a network). Notably, because of the special training process of neural networks, a neural network approach is particularly useful for handling complex, non-linear univariate and multivariate relationships that would be difficult to fit using other techniques (in the training process of neural networks, these rules are initially determined by a set of initial weight values, and their weights are adjusted during the learning process to increase efficiency. Through continuous adjustment and learning, the true network output and target value are achieved. After reaching the same value, the weighted value in the network is fixed, and the training is considered complete). Hence, this study utilizes a neural network approach to promote the performance of volatility forecasts for alternative GARCH family models in the stock market. This is the third contribution of this study because the above literatures of volatility forecasting never use this approach.
Thus, this study uses the stock indices in the developed and emerging markets as the sample and then utilizes eight parametric volatility forecasting models and eight composed volatility forecasting models to explore whether the neural network approach can promote the performance of volatility forecasts, whether the settings of the leverage effect and return distribution can encourage the performance of volatility forecasts and which one of the 16 models possesses the best volatility forecasting performance. In addition, for each of fourteen stock indices, this study also investigates which model is the most suitable for it. The stock indices in the Group of Seven (G7) and Emerging Seven (E7) are used to represent the developed and emerging markets in this study, respectively (the Group of Seven (G7) is an intergovernmental organization consisting of the United States, Canada, the United Kingdom, France, Germany, Italy and Japan. Its members are the world’s largest IMF-advanced economies and wealthiest liberal democracies. On the other hand, the Emerging 7 (E7) is the seven countries Brazil, Russia, India, China, Mexico, Indonesia and Turkey. They have the highest economic performance in the class of emerging economies or they are the seven biggest emerging countries in terms of economic growth). The eight parametric volatility forecasting models are composed of the GJR-GARCH-SGT model and its seven degenerate models. The eight composed volatility forecasting models are obtained by eight parametric volatility forecasting models combined with a neural network approach. From the empirical results of both the preliminary analysis of data and the performance comparison of the 16 models, I discovered the following phenomena. The preliminary analysis of data includes the descriptive statistics of data and estimation results of the GJR-GARCH-SGT model. First, the stock indices in the E7 have higher return and higher risk than those in the G7. Second, all the stock indices in the G7 and E7 for the forecasting period have higher risk than those for the overall period because of COVID-19 spreading throughout the world in the last year. Third, the leverage effect is significant in the stock indices in the G7 and E7, especially for the G7. Fourth, the distribution of returns is left-skewed and has a larger and thicker tail than the normal distribution. Fifth, the performance for the composed volatility forecasting models is significantly superior to that for the parametric volatility forecasting models, and thus the neural network approach can significantly promote the performance of volatility forecasting. Sixth, the performance for the GJR-based models is significantly superior to that of the GARCH-based models, and thus the setting of the leverage effect can significantly encourage the performance of volatility forecasting. Seventh, the performance of the models with non-normal distribution is not superior to that of the models with normal distribution, and thus the setting of the non-normal return distribution cannot promote the performance of volatility forecasting. Eighth, among the 16 models in this study, the performance of the GJR-GARCH-N-NN models is the best followed by the GJR-GARCH-N, GJR-GARCH-T-NN and GJR-GARCH-SGT-NN, and thus the GJR-GARCH model combining both the normal distribution and a neural network approach has the better performance of volatility forecasting. Finally, for each of the 14 stock indices, the most suitable models are not necessarily the same but they possess the setting of leverage effect and further combine with a neural network approach, and thus these results are the same as those obtained from the analysis of previous issues.
The remainder of this paper is organized as follows. Section 2 describes the empirical models utilized in this study, eight parametric volatility forecasting models and eight composed volatility forecasting models, and two types of loss function to evaluate the above models. Section 3 states the basic statistical features of the return series for the stock indices in the E7 and G7 during the overall period and two sub-periods: estimation period and forecast period. Section 4 analyzes the results of the empirical model and further explores the issues addressed in this study. Finally, Section 5 concludes the findings in Section 3 and Section 4.

2. Methodology

In order to accurately forecast volatility, the selected empirical model should capture the common features of financial assets. For example, the distribution of returns is skewed to the right or left and has a larger and thicker tail than the normal distribution. In other words, the return series is not normally distributed. Moreover, the return series exhibits linear dependence and strong ARCH effect. In addition, the volatility clustering and leverage effect usually exists in the volatility of financial asset return series [8,20,23]. Hence, the empirical models include a symmetric type of volatility specification, the GARCH, and an asymmetric one, the GJR-GARCH, combined with the normal (N), Student’s t (T), skewed Student’s t (ST) and skewed generalized Student’s t (SGT) distributions, totaling eight different models (the GJR-GARCH model of Glosten et al. [24] is an asymmetric type of GARCH-based model, and it can seize the financial features: the volatility clustering and leverage effect. On the contrary, the asymmetric type of distributions, skewed Student’s t (ST) and skewed generalized Student’s t (SGT), can capture the skewness and fat-tails on the distribution of return). The eight models above can be divided into two categories. The first category of model includes the GARCH-N, GARCH-T, GARCH-ST and GARCH-SGT models, named the GARCH-based models. The second category of model consists of the GJR-GARCH-N, GJR-GARCH-T, GJR-GARCH-ST and GJR-GARCH-SGT models, called the GJR-based models. The eight models above are called parametric volatility forecasting models. Notably, this study combines the above parametric volatility forecasting models with a neural network (NN) approach to promote the volatility forecasting performance in the stock market. Hence, there are an additional eight models, the GARCH-N-NN, GARCH-T-NN, GARCH-ST-NN, GARCH-SGT-NN, GJR-GARCH-N-NN, GJR-GARCH-T-NN, GJR-GARCH-ST-NN and GJR-GARCH-SGT-NN models, respectively, representing the GARCH-N, GARCH-T, GARCH-ST, GARCH-SGT, GJR-GARCH-N, GJR-GARCH-T, GJR-GARCH-ST and GJR-GARCH-SGT models combined with a neural network approach. The above eight models are named as composed volatility forecasting models.

2.1. Parametric Volatility Forecasting Models

Among the eight parametric volatility forecasting models in this study, the GJR-GARCH-SGT model can degenerate into the other seven models. Hence, in this subsection, I mainly illustrate the mean and variance equations of the GJR-GARCH-SGT model and then describe how the GJR-GARCH-SGT model be degenerated into the other seven models. The mean and variance equations of the GJR-GARCH-SGT model are expressed as follows:
r t = μ t + e t ,   μ t = ϕ 0 + ϕ 1 r t 1 ,   e t = z t h t ,   z t ~ IID   SGT 0 , 1 ; κ , λ , n
h t = ω + α + η I t 1 e t 1 2 + β h t 1
where r t represents the return of the stock indices in the emerging and developed markets where r t = lnP t lnP t 1 × 100 . P t is the close price of the stock index at time t and e t is the current error. μ t and h t represent the conditional mean and variance of return, respectively. Moreover, I t 1 is an indicator dummy that takes the value of 1 if e t 1 < 0 and zero otherwise, and thus parameter η is used to capture the leverage effect of volatility. Furthermore, ω , α and β are the parameters of variance equation and they obey the constraints ω , α ,   β > 0 and β + α + 0.5 η < 1 (according to Example 2.1 of Ling and McAleer [26], the necessary and sufficient condition for the existence of the second moment condition for all GJR-based models is β + α + 0.5 η < 1 ). Notably, if η = 0 in Equation (2), then the GJR-GARCH volatility specification degenerates into the GARCH volatility specification. IID denotes that the standardized errors z t are independent and identically distributed. Because z t is drawn from the standardized SGT distribution which allows returns innovation to follow a flexible treatment of both skewness and excess kurtosis in the conditional distribution of returns. The probability density function for the standardized SGT distribution is derived in Lee and Su [7] and can be represented as follows (the standardized SGT distribution, which has zero mean and unit variance, was checked by Mathematica software and another analogous standardized SGT distribution was proposed by Bali and Theodossiou [27]):
f z t = C 1 + z t + δ κ 1 + sign z t + δ λ κ θ κ n + 1 κ
where θ = 1 S λ B 1 κ , n κ 1 2 B 3 κ , n 2 κ 1 2 , S λ = 1 + 3 λ 2 4 A 2 λ 2 , A = B 2 κ , n 1 κ B 1 κ , n κ 1 2 B 3 κ , n 2 κ 1 2 , δ = 2 λ A S λ , C = κ 2 θ B 1 κ , n κ 1 .
where κ , n and λ are the scaling parameters and C and θ are the normalizing constants ensuring that f · is a proper probability density function. The parameters κ and n control the height and tails of density with the constraints κ > 0 and n > 2 , respectively. The skewness parameter λ controls the rate of descent of the density around the mode of z t with 1 < λ < 1 . In the case of positive and negative skewness, the density function skews toward to the right and left, respectively. B · is the beta function whereas ‘sign’ denotes a sign function. The parameter n has the degrees of freedom interpretation in case λ = 0 and κ = 2 . The log-likelihood function of the GJR-GARCH-SGT model thus can be expressed as:
L Ψ = ln f r t Ω t 1 ; Ψ = ln C ln h t n + 1 κ ln 1 + r t μ t h t + δ κ 1 + sign r t μ t h t + δ λ κ θ κ
where Ψ = ϕ 0 , ϕ 1 ,   ω , α , β , η , κ , λ , n is the vector of parameters to be estimated, and Ω t 1 denotes the information set of all observed returns up to time t−1. Notably, if κ = 2 in Equation (3), then the standardized SGT distribution degenerates into the standardized ST distribution. Using the same inference process, the standardized SGT distribution degenerates into the standardized Student’s t distribution if κ = 2 and λ = 0 and it also degenerates into the standardized normal distribution if κ = 2 , λ = 0 and n (regarding the process of the SGT distribution degenerating into the normal, Student’s t, skewed Student’s t (ST) distributions, please see Lee and Su [7] for more details).

2.2. Composed Volatility Forecasting Models

The eight composed volatility forecasting models are the GARCH-N-NN, GARCH-T-NN, GARCH-ST-NN, GARCH-SGT-NN, GJR-GARCH-N-NN, GJR-GARCH-T-NN, GJR-GARCH-ST-NN and GJR-GARCH-SGT-NN models. They are obtained by eight parametric volatility forecasting models combined with a neural network approach. A neural network (NN) approach is particularly useful for handling the complex, non-linear univariate and multivariate relationships that would be difficult to fit by using other techniques. A neural network model is composed of a multilayer perceptron with an interconnected group of nodes. For example, Figure 1a shows one backpropagation neural network with three layers, an input layer including two input nodes ( x 1   and   x 2 ), a hidden layer of three hidden nodes ( h ˜ 1 , h ˜ 2   and   h ˜ 3 ) and an output layer containing an output node ( y ). On the other hand, Figure 1b displays the other backpropagation neural network with an input layer including one input node ( x ), a hidden layer of five hidden nodes ( h ˜ 1 , h ˜ 2 , h ˜ 3 , h ˜ 4 , and   h ˜ 5 ) and an output layer containing one output node ( y ). Notably, the input nodes and the output nodes are analogous to the explanatory variables and the dependent variables in a regression model, respectively. The theory of neural network (NN) models is illustrated as follows (regarding the theory of neural network (NN) models, please see Lu et al. [28] and chapter 12 in the user’s guide of RATS version 6 or Doan [29]): Subsequently, I used a vector to represent the nodes of each layer and the total number of elements in a vector denote that of node in a layer. For example, the vectors X = x 1 , x 2 , , x d and Y = y 1 , y 2 , , y c represent the nodes in the input layer and output layer, respectively. On the other hand, the vector H = h ˜ 1 , h ˜ 2 , , h ˜ m denotes the nodes in the hidden layer. The three layers’ perceptron model is obtained by a weighted linear combination of the d input values from the d input nodes, X = x 1 , x 2 , , x d , and is expressed as follows:
a j = i = 1 d w ji 1 x i ·
The activation of hidden unit j can be achieved by transforming the linear sum via using a logistic activation function g a j = 1 / 1 + e a j :
h j ˜ = g a j = g i = 1 d w ji 1 x i
Thus, the node of output layer is defined as:
y k = g ˜ j = 1 m w jk 2 g i = 1 d w ji 1 x i
If the output function is taken linear, g ˜ a = a , the output model reduces to:
y k = j = 1 m w jk 2 g i = 1 d w ji 1 x i ·
where w ji 1 and w jk 2 are the weights of the hidden node and output node, respectively. Notably, fitting a neural networks model involves a training process of this model. The training process is executed by supplying a set of known input and output values, and then allowing the neural networks algorithm to adjust the hidden node and output node weights until the output produced by the networks matches the actual output in the training sample to the desired degree of accuracy (in this study, the input values are a volatility forecasting series h ^ t obtained by the variance equation of parametric volatility forecasting models as shown in Equation (2). On the contrary, the output values are the true values of variance and are replaced by the squared intraday return series ( r t 2 ) as the proxies. In addition, the neural networks algorithm in the training process is the instruction of ‘NNLEARN’ in RATS version 6, manufacturer, city and country). Once the training process is completed, the model can be used to generate new output data from the other sets of inputs (the new output data may be the fitted values for the in-sample volatility forecasts or the forecast values for the out-of-sample volatility forecasts in this study. Notably, the new output data is obtained by the instruction of ‘NNTEST’ in RATS version 6). Assuming that the fit is good, the relationships represented by the sample input and output data can generalize to other samples, thus the model can produce good predictions. Thus, the volatility forecasts of the GJR-GARCH-SGT-NN model are obtained by assigning a volatility forecasting series h ^ t obtained from the GJR-GARCH-SGT model as the input values and using a backpropagation neural network with one input node, five hidden nodes and an output node as shown in Figure 1b. Thus, I took an example of the GJR-GARCH-SGT-NN model to explain the process of the volatility forecasting for a neural network approach. Figure 2 lists the schematic diagram for the process of the volatility forecasting of the GJR-GARCH-SGT-NN model, and is illustrated as follows:
Step 1
 Fit a return series ( r t ) into the GJR-GARCH-SGT model (see, Equations (1)–(4)). Then, get a forecasted variance series ( h ^ t ) (see, Equation (2))
Step 2
 Input a forecasted variance series ( h ^ t ) obtained from step 1 and the true variance series ( r t 2 ) into the ‘NNLEARN’ function. Then, obtain the weights of the hidden nodes ( w ji 1 ) and the weights of the output node ( w jk 2 ). The ‘NNLEARN’ function respectively regards h ^ t and r t 2 as the input value and output value of Figure 1b, and then substitutes them ( h ^ t and r t 2 ) into Equation (8) to execute a training process of the neural network in order to obtain the weights w ji 1 and w jk 2 (in the training process of neural networks, the weights of the hidden nodes and the output node ( w ji 1 and w jk 2 ) are initially determined by a set of initial weights values, and their weights are adjusted until the output produced by the networks ( h ^ t ) matches the actual output ( r t 2 ) in the training sample to the desired degree of accuracy. Te training is then considered complete, the weights w ji 1 and w jk 2 are obtained). Figure 1b is a structure of a backpropagation neural network with one node in the input layer, five nodes in the hidden layer and one node in the output layer.
Step 3
 Input a forecasted variance series ( h ^ t ) obtained from step 1 and the weights of the hidden nodes ( w ji 1 ) and the weights of the output node ( w jk 2 ) obtained from step 2 into the ‘NNTEST’ function. Furthermore, obtain the forecasted variance series for a neural network approach ( h ^ t NN ). That is, the ‘NNTEST’ function substitutes the input value h ^ t and weights of the hidden nodes and output node ( w ji 1 and w jk 2 ) into Equation (8), and then the output value h ^ t NN is obtained.
Step 4
 Given a forecasted variance series for a neural network approach ( h ^ t NN ) and the true variance series ( r t 2 ), then calculate the values of loss functions, MAE and RMSE (see Equations (9) and (10)).
In addition, I used two types of loss function to perform the performance comparison of volatility forecasting for the sixteen models or four categories of model. The first category of model is composed of one parametric volatility forecasting model and one corresponding composed volatility forecasting model (for example, the GARCH-N and GARCH-N-NN, the GARCH-T and GARCH-T-NN, the GARCH-ST and GARCH-ST-NN, the GARCH-SGT and GARCH-SGT-NN, the GJR-GARCH-N and GJR-GARCH-N-NN, the GJR-GARCH-T and GJR-GARCH-T-NN, the GJR-GARCH-ST and GJR-GARCH-ST-NN and the GJR-GARCH-SGT and GJR-GARCH-SGT-NN, totaling 8 paired models in the first category of model). The second category of model is composed of one GARCH-based model and one corresponding GJR-based model (for instance, the GARCH-N and GJR-GARCH-N, the GARCH-T and GJR-GARCH-T, the GARCH-ST and GJR-GARCH-ST, the GARCH-SGT and GJR-GARCH-SGT, the GARCH-N-NN and GJR-GARCH-N-NN, the GARCH-T-NN and GJR-GARCH-T-NN, the GARCH-ST-NN and GJR-GARCH-ST-NN and the GARCH-SGT-NN and GJR-GARCH-SGT-NN, totaling 8 paired models in the second category of model). The third category of model is a group of models with the same volatility specification and volatility forecasting approach but different distribution (for example, the GARCH-N, GARCH-T, GARCH-ST and GARCH-SGT; the GJR-GARCH-N, GJR-GARCH-T, GJR-GARCH-ST and GJR-GARCH-SGT; the GARCH-N-NN, GARCH-T-NN, GARCH-ST-NN and GARCH-SGT-NN; and the GJR-GARCH-N-NN, GJR-GARCH-T-NN, GJR-GARCH-ST-NN and GJR-GARCH-SGT-NN, totaling four groups of model in the third category of model). The fourth category of model is composed of all sixteen models in this study. Hence, the performance comparison of volatility forecasting for the first category of model can explore whether the neural network approach can promote the performance of volatility forecasting. On the contrary, the performance comparison of volatility forecasting for the second and third categories of model can investigate whether the settings of the leverage effect and return distribution can encourage the performance of volatility forecasts, respectively. In addition, the performance comparison of volatility forecasting for the fourth category of model can explore which one of the 16 models possesses the best volatility forecasting performance. Two types of loss function are the mean absolute error (MAE) and root mean squared error (RMSE). The MAE measures the average magnitude of the errors in a set of forecasts, without considering their direction. The MAE for the out-of-sample volatility forecasting can be evaluated by the following equation.
MAE = 1 T i = 1 T ε i = 1 T i = 1 T h ^ t + i t + i 1 h t + i
where ε i denotes the forecast error; h ^ t + i t + i 1 is the one-step-ahead forecast of the variance of the returns dependent on all information upon the time t + i -1 and can be estimated by one of the sixteen models in this study; h t + i is the true value of variance and is replaced by the squared intraday returns as the proxies; T is the number of computing 1-day-ahead variance and is equal to 250 in this study. Moreover, the RMSE for the out-of-sample volatility forecasting can be obtained by the following equation (the MAE and RMSE for in-sample volatility forecasting can be evaluated by the equation: MAE = 1 T i = 1 T ε i = 1 T i = 1 T h ^ i h i and RMSE = 1 T i = 1 T ε i 2 = 1 T i = 1 T h ^ i h i 2 where T is the total number of observations for overall period and is equal to 3501 in this study).
RMSE = 1 T i = 1 T ε i 2 = 1 T i = 1 T h ^ t + i t + i 1 h t + i 2
where ε i , h ^ t + i t + i 1 , h t + i and T are defined the same as for Equation (9). Hence, the MAE and RMSE can measure the differences between the values predicted by a model and the values actually observed from the thing being modeled or estimated.

3. Data and Descriptive Statistics

This study uses the stock indices in the developed and emerging markets as the sample to explore the issue of ‘how to promote the performance of parametric volatility forecasting models? A neural networks approach’. I used the stock indices in the G7 and E7 to represent the stock markets in the developed and emerging markets, respectively. The stock indices in the G7 include the Dow Jones (DJ), TSX, FTSE, CAC40, DAX, MIB and N225, respectively corresponding to the United States, Canada, the United Kingdom, France, Germany, Italy and Japan. On the contrary, the stock indices in the E7 are the BVSP, RTSI, BSE, SSE, MXX, JKSE and XU100, respectively corresponding to Brazil, Russia, India, China, Mexico, Indonesia and Turkey.
Table 1 reports basic descriptive statistics of daily return of stock indices in the G7 and E7 during the overall period and two sub-periods. The overall period starts from 10 January 2001 and ends 27 July 2020 and is used to perform the in-sample volatility forecasting. The overall period is divided into two sub-periods, the estimation period and forecasting period, to execute the out-of-sample volatility forecasting. The estimation period is from 10 January 2001 to 8 April 2019 whereas the forecasting period is from 9 April 2019 to 27 July 2020. Notably, all the daily close price data of the 14 stock indices were obtained from the Yahoo finance website. As shown in panel A of Table 1, the mean values of G7 are between −0.0215 (MIB) and 0.0263 (DJ) whereas those of E7 range from 0.0254 (SSE) to 0.0717 (JKSE). Conversely, the values of standard deviation for the G7 are between 1.2631 (TSX) and 1.8306 (MIB) whereas those for the E7 range from 1.4361 (MXX) to 2.3155 (RTSI). These results indicate that during the overall period, the stock indices in the E7 have the higher return and higher risk than those in the G7 (the reason is that the maximum value of mean for the E7 (0.0717) is greater than that for the G7 (0.0263) and the minimum value of mean for the E7 (0.0254) is also greater than that for the G7 (−0.0215). On the other hand, the maximum value of standard deviation for the E7 (2.3155) is greater than that for the G7 (1.8306) and the minimum value of standard deviation for the E7 (1.4361) is also greater than that for the G7 (1.2631)). The above finding is consistent with that found in Su [9]. As illustrated in panel B and panel C of Table 1, I discovered that the stock indices in the G7 and E7 for the estimation and forecasting periods possess the same phenomena found from the overall period. That is, the stock indices in the E7 have higher return and higher risk than those in the G7. Notably, regarding the forecasting period, the values of standard deviation for the G7 are between 1.7269 (N225) and 2.3058 (MIB) whereas those for the E7 range from 1.6095 (MXX) to 3.1898 (SSE). These results indicate that the stock indices in the G7 and E7 for the forecasting period have higher risk than those for the overall period (the reason is, regarding the G7, the maximum value of standard deviation for the forecasting period (2.3058) is greater than that for the overall period (1.8306) and the minimum value of standard deviation for the forecasting period (1.7269) is greater than that for the overall period (1.2631). Conversely, regarding the E7, the maximum value of standard deviation for the forecasting period (3.1898) is greater than that for the overall period (2.3155) and the minimum value of standard deviation for the forecasting period (1.6095) is greater than that for the overall period (1.4361)). This phenomenon is attributed to COVID-19 spreading throughout the world in the last year. Figure 3 illustrates the trend of price levels, and the variation of return for the 14 stock indices during the overall period. From Figure 3, I also discovered that in the last year, the price of stock indices underwent a severe decline and its return experienced a serious variation, indicating that the high risk appears at the forecasting period. This phenomenon is the same as that found from the above analysis. In addition, the volatility clustering occurs significantly in the overall period.
Finally, regarding the other descriptive statistics, they have the same features as those for most of the financial return series. For example, the distribution of returns is left-skewed and has a larger and thicker tail than the normal distribution, indicating that the return series is not normally distributed. The above results are found by coefficient of skewness, excess kurtosis and the J-B normality test statistics [30]. In addition, the return series exhibit linear dependence and strong ARCH effect as shown by the Ljung–Box Q 2 20 statistics for the squared returns. From the above findings, a GARCH family model is very suitable to seize the fat tails and time-varying volatility found in these asset return series.

4. Empirical Results

As described in Section 2, the empirical models in this study can be divided into two categories: the parametric volatility forecasting models and composed volatility forecasting models. The composed volatility forecasting models are the parametric volatility forecasting models combined with a neural network approach. The parametric volatility forecasting models include the GARCH-N, GARCH-T, GARCH-ST, GARCH-SGT, GJR-GARCH-N, GJR-GARCH-T, GJR-GARCH-ST and GJR-GARCH-SGT models. Among the above eight parametric volatility forecasting models, the GJR-GARCH-SGT model is the most flexible because this model can capture most of the common features of financial assets and this model can degenerate into the other seven models under the setting of some restrictions. Hence, I will report some financial features for the 14 stock indices via using the empirical results of the GJR-GARCH-SGT model.

4.1. Estimation Results for the GJR-GARCH-SGT Model

Table 2 illustrates the empirical results of the GJR-GARCH-SGT model for the stock indices in the G7 and E7. As shown in Table 2, parameters ω , α and β are significantly positive for most stock indices and all stock indices obey the constraint β + α + 0.5 η < 1 as reported by the numbers listed in row ‘C’ of Table 2. For example, the constraint β + α + 0.5 η for BVSP is equal to 0.9756, and is less than 1 because the values of parameters β , α and η are equal to 0.9197, 0.0143 and 0.0832, respectively. Moreover, parameter η is significantly positive at 1% for all stock indices. The values of parameter η for the stock indices in the G7 are greater than 0.1421 (TSX) whereas those for the stock indices in the E7 are smaller than 0.0873 (MXX) except for BSE. These results indicate that the leverage effect exists significantly in the stock indices in the G7 and E7, especially for the G7 because the values of parameter η for the case of G7 are greater than those for the E7. In addition, the shape parameters κ , n and λ are significant for most of stock indices and obey the constraints κ > 0 , n > 2 and 1 < λ < 1 . Notably, the values of parameter λ are significantly negative for most of stock indices. These results indicate that the distribution of returns is left-skewed and has a larger and thicker tail than the normal distribution. Finally, regarding the Ljung–Box Q2 (20) statistics for the squared returns, they are not significant at the 10% level for most of cases and they are all far smaller than the same statistics appearing in Table 1. These results indicate that the serial correlation does not exist in standard residuals and the GJR-GARCH-SGT model is sufficient to correct the serial correlation of these returns series in the conditional variance equation for the stock indices in the G7 and E7.

4.2. The Performance Assessment of Volatility Forecasts

Via the analysis of empirical results in Table 2, the selected empirical model in this study can capture the common features of financial assets well. I then executed the in-sample and out-of-sample volatility forecasts of the 16 models for the the 14 stock indices in the G7 and E7. Table 3 and Table 4 report the results of in-sample volatility forecasts, respectively, based on the MAE and RMSE loss functions for the overall period. On the contrary, Table 5 and Table 6 list the results of out-of-sample volatility forecasts, respectively, based on the MAE and RMSE loss functions for the forecasting period via using a rolling window approach (for each data series, the eight parametric volatility forecasting models and eight composed volatility forecasting models, totaling sixteen models, were first estimated using a sample of 3250 daily returns, and a volatility forecast for the next period was obtained. Subsequently, the estimation period was rolled forward by adding one new day and omitting the most distant day. By repeating this procedure, the out-of-sample volatility forecasts were calculated for the next 250 days). Figure 4 shows the trend of actual variance and its two out-of-sample variance forecasts obtained by the GJR-GARCH-SGT and GJR-GARCH-SGT-NN models. From Figure 4, I observed that there is a significantly sharp value of variance on March of last year owing to COVID-19 spreading througout the world.
Subsequently, regarding Table 3, Table 4, Table 5 and Table 6, I performed the volatility forecasting performance comparison for four categories of model to explore whether the neural network approach can promote the performance of volatility forecasting, whether the settings of leverage effect and non-normal return distribution can encourage the performance of volatility forecasting and which one of the 16 models possesses the best volatility forecasting performance, and then record the results of performance comparison in columns S1, S2, S3 and S4 of each table (regarding four categories of model, please see Section 2.2 for more details). The results in columns S1, S2, S3 and S4 of Table 3, Table 4, Table 5 and Table 6 are also summarized in Table 7 in order to easily explore the four main issues of this study. I used the data in Table 3 to illustrate the performance comparison of volatility forecasting for the four categories of the model. Regarding the performance comparison of volatility forecasting for the first category of model, I took the following two examples in Table 3 to illustrate it. First, regarding the paired models, ‘the GARCH-N and GARCH-N-NN’, the GARCH-N model has the lower value of MAE for the cases of DJ, TSX, FTSE and BVSP whereas the GARCH-N-NN model possesses the lower value of MAE for the cases of the other ten stock indices. For example, regarding the DJ, the value of MAE for the GARCH-N model (1.94301) is lower than that for the GARCH-N-NN model (1.94312). Furthermore, in Table 3 the results ‘4’ for the GARCH-N model and ‘10’ for the GARCH-N-NN model are recorded in column ‘S1’, respectively corresponding to the rows ‘GARCH-N’ and ‘GARCH-N-NN’. In Table 7, the above results, ‘4’ and ‘10’, are also recorded in column ‘MAE’ below ‘In-sample’ of S1, respectively corresponding to the rows ‘GARCH-N’ and ‘GARCH-N-NN’. Second, regarding the paired models, ‘the GJR-GARCH-SGT and GJR-GARCH-SGT-NN’, the GJR-GARCH-SGT-NN model possesses the lower value of MAE for all 14 stock indices but the GJR-GARCH-SGT model does not obtain the lower value of MAE. In Table 3 the results ‘0’ and ‘14’ are then recorded in column ‘S1’, respectively corresponding to the rows ‘GJR-GARCH-SGT’ and ‘GJR-GARCH-SGT-NN’. In Table 7, the above results, ‘0’ and ‘14’, are also recorded in column ‘MAE’ below ‘In-sample’ of S1, respectively corresponding to the rows ‘GJR-GARCH-SGT’ and ‘GJR-GARCH-SGT-NN’. Regarding the performance comparison of volatility forecasting for the second category of model, I took the following example in Table 3 to explain it. Regarding the paired models, ‘the GARCH-N and GJR-GARCH-N models’, the GARCH-N model has the lower value of MAE for the cases of MIB, RTSI, BSE and XU100 whereas the GJR-GARCH-N model possesses the lower value of MAE for the other ten stock indices. For example, regarding the MIB, the value of MAE for the GARCH-N model (3.79126) is lower than that for the GJR-GARCH-N model (3.79163). In Table 3 the results ‘4’ and ‘10’ are then recorded in column ‘S2’, respectively corresponding to the rows ‘GARCH-N’ and ‘GJR-GARCH-N’. In Table 7, the above results, ‘4’ and ‘10’, are also recorded in column ‘MAE’ below ‘In-sample’ of S2, respectively corresponding to the rows ‘GARCH-N’ and ‘GJR-GARCH-N’. Regarding the performance comparison of volatility forecasting for the third category of model, I took the following example in Table 3 to explain it. Among a group of models, ‘the GARCH-N, GARCH-T, GARCH-ST and GARCH-SGT’, the GARCH-N model possesses the lowest value of MAE for the cases of DJ, TSX, FTSE, CAC40, DAX, BVSP and RTSI but the GARCH-T model never obtains the lowest value of MAE. For instance, regarding the DJ, the value of MAE for the GARCH-N model (1.94301) is lower than that for the other three GARCH-based models such as the GARCH-T (1.99808), GARCH-ST (1.98032) and GARCH-SGT (1.95943). In other words, regarding the DJ, the GARCH-N model possesses the lowest value of MAE among four GARCH-based models. In addition, the GARCH-ST model has the lowest value of MAE for the cases of BSE, MXX and XU100 whereas the GARCH-SGT model obtains the lowest value of MAE for the cases of MIB, N225, SSE and JKSE. In Table 3 the results ‘7’, ‘0’, ‘3’, and ‘4’ are then recorded in column ‘S3’, respectively corresponding to the rows ‘GARCH-N’, ‘GARCH-T’, ‘GARCH-ST’ and ‘GARCH-SGT’. In Table 7, the above results ‘7’, ‘0’, ‘3’, and ‘4’ are also recorded in column ‘MAE’ below ‘In-sample’ of S3, respectively corresponding to the rows ‘GARCH-N’, ‘GARCH-T’, ‘GARCH-ST’ and ‘GARCH-SGT’. Regarding the performance comparison of volatility forecasting for the fourth category of model, I took the following example in Table 3 to explain it. Among all 16 models, the GJR-GARCH-N model has the lowest value of MAE for the cases of TSX and BVSP. For example, regarding the TSX, the value of MAE for the GJR-GARCH-N model (1.64781) is lower than that for the other fifteen models such as the GARCH-N (1.67592), GARCH-T (1.72006), GARCH-ST (1.70677), GARCH-SGT (1.70920), GJR-GARCH-T (1.69383), GJR-GARCH-ST (1.68906), GJR-GARCH-SGT (1.69515), GARCH-N-NN (1.71390), GARCH-T-NN (1.71509), GARCH-ST-NN (1.71786), GARCH-SGT-NN (1.71812), GJR-GARCH-N-NN (1.68169), GJR-GARCH-T-NN (1.69363), GJR-GARCH-ST-NN (1.69340) and GJR-GARCH-SGT-NN (1.69350). In other words, regarding the TSX, the GJR-GARCH-N model has the lowest value of MAE among all sixteen models. In Table 3 the result ‘2’ is then recorded in column ‘S4’, corresponding to the row ‘GJR-GARCH-N’. In Table 7, the above result ‘2’ is also recorded in column ‘MAE’ below ‘In-sample’ of S4, corresponding to the row ‘GJR-GARCH-N’.
Table 7 lists the summary results of performance comparison for the in-sample and out-of-sample volatility forecasts based on the MAE and RMSE loss functions. In other words, the numbers in column ‘MAE’ below ‘In-sample’ of S1, S2, S3 and S4 are respectively summarized from those in columns ‘S1’, ‘S2’, ‘S3’ and ‘S4’ of Table 3. On the contrary, the numbers in column ‘RMSE’ below ‘In-sample’ of S1, S2, S3 and S4 are respectively summarized from those in columns ‘S1’, ‘S2’, ‘S3’ and ‘S4’ of Table 4. The numbers in column ‘MAE’ below ‘Out-of-sample’ of S1, S2, S3 and S4 are respectively summarized from those in columns ‘S1’, ‘S2’, ‘S3’ and ‘S4’ of Table 5. Conversely, the numbers in column ‘RMSE’ below ‘Out-of-sample’ of S1, S2, S3 and S4 are respectively summarized from those in columns ‘S1’, ‘S2’, ‘S3’ and ‘S4’ of Table 6. In order to easily explore the four main issues of this study, I performed calculations for the summation of all four numbers in column ‘S1’ for each model, as well as columns ‘S2’, ‘S3’ and ‘S4’. For example, regarding the GARCH-N model, the numbers ‘4′ and ‘0’ are respectively in columns ‘MAE’ and ‘RMSE’ below ‘In-sample’ of S1. Moreover, the numbers ‘3′ and ‘5′ are respectively in columns ‘MAE’ and ‘RMSE’ below ‘Out-of-sample’ of S1. Hence, in Table 7, the summation of all four numbers in column ‘S1’ is equal to 12, and is recorded in the column ‘Sum’ below ‘S1’ and the row ‘GARCH-N’. Regarding the other 15 models, the summation of all four numbers in columns ‘S1’ must be done with the same inference process. With regard to the 16 models, the summation of all four numbers in columns ‘S2’, ‘S3’ or ‘S4’ must also be evaluated with the same inference process. The above summation results in columns ‘S2’, ‘S3’ or ‘S4’ are recorded in the column ‘Sum’ below ‘S1’, ‘S3’ or ‘S4’, respectively. Subsequently, I used all 16 numbers in column ‘Sum’ below ‘S1’, ‘S2’, ‘S3’, and ‘S4’ of Table 7 to execute the performance comparison of volatility forecasting for four categories of model. As shown by the numbers at column ‘Sum’ below S1 of Table 7, I found that the numbers for all eight composed volatility forecasting models are far greater than those for all eight corresponding parametric volatility forecasting models. For example, the number for the GARCH-N-NN model (44) is far greater than that for the GARCH-N model (12). These results indicate that the performance for the composed volatility forecasting models is significantly superior to that for the parametric volatility forecasting models. In other words, the neural network approach can significantly improve the performance of volatility forecasting. As reported by the numbers in column ‘Sum’ below S2 of Table 7, I found that, regarding the parametric volatility forecasting approach, the numbers for all four GJR-based models are far greater than those for all four corresponding GARCH-based models, as shown in panel A of this table. For instance, the number for the GJR-GARCH-N model (38) is far greater than that for the GARCH-N model (18). I also found that, regarding the composed volatility forecasting approach, the numbers for all four GJR-based models are far greater than those for all four corresponding GARCH-based models, as shown in panel B of this table. These results imply that irrespective of the parametric forecasting approach or composed forecasting approach, the performance for the GJR-based models is significantly superior to that of the GARCH-based models. That is to say, the setting of the leverage effect can significantly encourage the performance of volatility forecasting (as shown in Section 2, the GJR-based model can seize the leverage effect appearing at the financial assets whereas the GARCH-based model cannot.) As illustrated by the numbers in column ‘Sum’ below S3 of Table 7, I found that the numbers for the models with non-normal distribution are not greater than those for the models with normal distribution based on the same volatility forecasting approach and volatility specification. For example, the number for the GARCH-N model (35) is far greater than those for the GARCH-T (1), GARCH-ST (5) and GARCH-SGT (15) models. Moreover, the number for the GJR-GARCH-N model (45) is far greater than those for the GJR-GARCH-T (2), GJR-GARCH-ST (3) and GJR-GARCH-SGT (6) models. Furthermore, the number for the GARCH-N-NN model (37) is far greater than those for the GARCH-T-NN (13), GARCH-ST-NN (4) and GARCH-SGT-NN (2) models. In addition, the number for the GJR-GARCH-N-NN model (32) is far greater than those for the GJR-GARCH-T-NN (12), GJR-GARCH-ST-NN (3) and GJR-GARCH-SGT-NN (10) models. The above results indicate that irrespective of volatility forecasting approach or volatility specification, the performance of the models with the non-normal distribution is not superior to that of the models with the normal distribution. In other words, the setting of the non-normal return distribution cannot promote the performance of volatility forecasting. As listed by the 16 numbers in column ‘Sum’ below S4 of Table 7, I found the number for the GJR-GARCH-N-NN model (15) is the greatest. On the contrary, the numbers for the GJR-GARCH-N, GJR-GARCH-T-NN and GJR-GARCH-SGT-NN are all equal to 9, the second greatest among the 16 numbers. The above result indicates that, among the 16 models in this study, the performance of the GJR-GARCH-N-NN models is the best followed by GJR-GARCH-N, GJR-GARCH-T-NN and GJR-GARCH-SGT-NN. In other words, the GJR-GARCH model combined with both the normal distribution and a neural networks approach has the best performance of volatility forecasting among the sixteen models in this study.
In addition, this study also investigates which model is the most suitable for each of the fourteen stock indices. That is, regarding each stock index, which model has the best performance of volatility forecasting in order to find the most suitable model for each stock index. In order to easily explore this issue, I summarized the most superior model for each stock index based on two types of volatility forecasts (in-sample and out-of-sample) and two types of loss function (MAE and RMSE). Taking an example of ‘DJ’ stock index, among the 16 models, the GJR-GARCH-SGT-NN model has the best performance for in-sample volatility forecast based on MAE (respectively, RMSE) as shown in the column ‘DJ’ of Table 3 (respectively, Table 4). These results are recorded in column ‘DJ’ and rows ‘MAE’ and ‘RMSE’ of ‘In-sample’ in Table 8. Conversely, among the 16 models, the GJR-GARCH-N-NN model has the best performance for out-of-sample volatility forecast based on MAE (respectively, RMSE) as shown in the column ‘DJ’ of Table 5 (respectively, Table 6). These results are recorded in column ‘DJ’ and rows ‘MAE’ and ‘RMSE’ of ‘Out-of-sample’ in Table 8. Hence, Table 8 summarizes the results of the most suitable mode for alternative stock indices. In other words, the results listed in row ‘MAE’ (respectively, ‘RMSE’) of ‘In-Sample’ in Table 8 are summarized from the results of the performance comparison for the fourth category of model in Table 3 (respectively, Table 4). On the other hand, the results listed in row ‘MAE’ (respectively, ‘RMSE’) of ‘Out-of-Sample’ in Table 8 are summarized from the results of the performance comparison for the fourth category of model in Table 5 (respectively, Table 6). From Table 8, I found that both GJR-GARCH-N-NN and GJR-GARCH-SGT-NN are the most suitable models for the DJ stock index because both the GJR-GARCH-N-NN and GJR-GARCH-SGT-NN appear twice among four cases in column ‘DJ’ in Table 8 (the four cases in Table 8 are composed of two types of volatility forecasts (in-sample and out-of-sample) and two types of loss function (MAE and RMSE) when the volatility forecasting of a specific stock index is executed). These results are recorded in row ‘Best model’, corresponding to column ‘DJ’ in Table 8. On the contrary, GJR-GARCH-N is the most suitable model for the TSX stock index because GJR-GARCH-N is the most relevant for the four cases. These results are recorded in row ‘Best model’, corresponding to column ‘TSX’ in Table 8. In the same inference process, I found the most suitable models for the others stock indices, and I recorded them in row ‘Best model’ and the column corresponding to the specific stock index in Table 8. From the results listed in the row ‘Best model’ of Table 8, I obtained the following conclusion. First, GJR-GARCH-N is the most suitable model for the TSX, BVSP and MXX. Second, GARCH-N-NN is the most suitable model only for the RTSI. Third, GJR-GARCH-N-NN is the most suitable model for the DJ, FTSE, MIB and SSE. Fourth, GJR-GARCH-T-NN is the most suitable model for the DAX, JKSE and XU100. Fifth, the GJR-GARCH-SGT-NN is the most suitable model for the DJ, N225 and BSE. To sum up, the most suitable models for the 14 stock indices are distributed at the GJR-GARCH-N, GARCH-N-NN, GJR-GARCH-N-NN, GJR-GARCH-T-NN and GJR-GARCH-SGT-NN models. These results indicate that the most suitable models are not necessarily the same for each of the 14 stock indices. Regarding the most suitable models above, they possess the setting of leverage effect and further combine with a neural networks approach. As to the setting of distribution, they are randomly distributed at the normal, Student’s t and SGT. Hence, the above conclusions are the same as those obtained from the analysis of previous issues. That is, a neural network approach and the setting of leverage effect can significantly promote the performance of volatility forecasting but the setting of non-normal distribution cannot.

5. Conclusions

This study uses the stock indices in the developed and emerging markets as the sample and then utilizes eight parametric volatility forecasting models and eight composed volatility forecasting models to explore whether the neural networks approach can promote the performance of volatility forecasting, whether the settings of leverage effect and non-normal return distribution can encourage the performance of volatility forecasting and which one of the 16 models posseses the best volatility forecasting performance. In addition, this study also investigates which model is the most suitable for each of 14 stock indices. The stock indices in the G7 and E7 are used to represent the stock markets in the developed and emerging markets, respectively. The eight parametric volatility forecasts models are composed of the GJR-GARCH or GARCH models with the normal, Student’s t, skewed Student’s t and generalized skewed Student’s t distributions. The eight composed volatility forecasting models are the eight parametric volatility forecasting models combined with a neural network approach. Notably, the neural network model has the same concepts as the network topology.
The empirical findings can be summarized as follows. From the descriptive statistics of data and estimation results of the GJR-GARCH-SGT model, I obtained the following conclusions. First, the stock indices in the E7 have higher return and higher risk than those in the G7. Second, all the stock indices in the G7 and E7 for the forecasting period have higher risk than those for the overall period because of COVID-19 spreading throughout the world in the last year. Third, the leverage effect exists significantly in the stock indices in the G7 and E7, especially for the G7. Fourth, the distribution of returns is left-skewed and has a larger and thicker tail than the normal distribution. From the performance comparison for the four categories of model, I drew the following conclusions, irrespective of in-sample or out-of-sample volatility forecasting. First, the performance for the composed volatility forecasting models is significantly superior to that of the parametric volatility forecasting models, indicating that the neural network approach can significantly improve the performance of volatility forecasting. Second, irrespective of parametric forecasting approach or composed forecasting approach, the performance for the GJR-based models is significantly superior to that of the GARCH-based models, implying that the setting of the leverage effect can significantly encourage the performance of volatility forecasting. Third, irrespective of volatility forecasting approach or volatility specification, the performance of the models with non-normal distribution are not superior to that of the models with the normal distribution, indicating that the setting of the non-normal return distribution cannot promote the performance of volatility forecasting. Fourth, among the 16 models in this study, the performance of the GJR-GARCH-N-NN models is the best followed by GJR-GARCH-N, GJR-GARCH-T-NN and GJR-GARCH-SGT-NN, inferring that the GJR-GARCH model combined with both the normal distribution and a neural network approach has the better performance of volatility forecasting. Finally, for each of the 14 stock indices, the most suitable models are not necessarily the same but they possess the setting of leverage effect and further combine with a neural network approach, and thus these results are the same a those obtained from the analysis of previous issues.
Based on the above empirical results, I propose the following important policy implications for investors and fund managers. First, the investors should use the asymmetric GARCH model to precisely forecast the volatility of stock indices because the asymmetric GARCH model such as the GJR-based model can seize the leverage effect existing in the stock indices and the setting of the leverage effect can significantly encourage the performance of volatility forecasting. Second, the investors should utilize the composed volatility forecasting models to accurately forecast the volatility of stock indices because the neural networks approach can handle the complex, non-linear univariate and multivariate relationships that would be difficult to fit using other techniques, and therefore it can significantly promote the performance of volatility forecasting. Hence, the two policy implications above can enable investors and fund managers to precisely execute subsequent investment processes such as the asset allocation, option pricing, risk management and hedge strategy in the stock market.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the daily close price data of 14 stock indices was obtained from the Yahoo finance website. http://finance.yahoo.com (accessed on 28 July 2021).

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Aziz, N.S.A.; Vrontos, S.; Hasim, H.M. Evaluation of multivariate GARCH models in an optimal asset allocation framework. N. Am. J. Econ. Financ. 2019, 47, 568–596. [Google Scholar] [CrossRef] [Green Version]
  2. Harris, R.D.F.; Nguyen, A. Long memory conditional volatility and asset allocation. Int. J. Forecast. 2013, 29, 258–273. [Google Scholar] [CrossRef] [Green Version]
  3. Su, J.B. The implementation of asset allocation approaches: Theory and evidence. Sustainability 2020, 12, 7162. [Google Scholar] [CrossRef]
  4. Escobar-Anel, M.; Rastegari, J.; Stentoft, L. Option pricing with conditional GARCH models. Eur. J. Oper. Res. 2021, 289, 350–363. [Google Scholar] [CrossRef]
  5. Yoo, E.G.; Yoon, S.J. CBOE VIX and Jump-GARCH option pricing models. Int. Rev. Econ. Financ. 2020, 69, 839–859. [Google Scholar] [CrossRef]
  6. Ardia, D.; Bluteau, K.; Boudt, K.; Catania, L. Forecasting risk with Markov-switching GARCH models: A large-scale performance study. Int. J. Forecast. 2018, 34, 733–747. [Google Scholar] [CrossRef]
  7. Lee, C.F.; Su, J.B. Alternative statistical distributions for estimating Value-at-Risk: Theory and evidence. Rev. Quant. Financ. Account. 2012, 39, 309–331. [Google Scholar] [CrossRef]
  8. Su, J.B. Empirical analysis of long memory, leverage, and distribution effects for stock market risk estimates. N. Am. J. Econ. Financ. 2014, 30, 1–39. [Google Scholar] [CrossRef]
  9. Su, J.B. Value-at-risk estimates of the stock indices in developed and emerging markets including the spillover effects of currency market. Econ. Model. 2015, 204–224. [Google Scholar] [CrossRef]
  10. Chang, C.L.; McAleer, M.; Tansuchat, R. Crude oil hedging strategies using dynamic multivariate GARCH. Energy Econ. 2011, 33, 912–923. [Google Scholar] [CrossRef] [Green Version]
  11. Chang, C.L.; González-Serrano, L.; Jimenez-Martin, J.A. Currency hedging strategies using dynamic multivariate GARCH. Math. Comput. Simul. 2013, 94, 164–182. [Google Scholar] [CrossRef] [Green Version]
  12. Aliyev, F.; Ajayi, R.; Gasim, N. Modelling asymmetric market volatility with univariate GARCH models: Evidence from Nasdaq-100. J. Econ. Asymmetries 2020, 22, e00167. [Google Scholar] [CrossRef]
  13. Bhowmik, R.; Wang, S. Stock market volatility and return analysis: A systematic literature review. Entropy 2020, 22, 522. [Google Scholar] [CrossRef]
  14. Chun, D.; Cho, H.; Ryu, D. Forecasting the KOSPI200 spot volatility using various volatility measures. Phys. A 2019, 514, 156–166. [Google Scholar] [CrossRef]
  15. Lee, Y.H.; Pai, T.Y. REIT volatility prediction for skew-GED distribution of the GARCH model. Expert Syst. Appl. 2010, 37, 4737–4741. [Google Scholar] [CrossRef]
  16. Liu, H.C.; Hung, J.C. Forecasting S&P-100 stock index volatility: The role of volatility asymmetry and distributional assumption in GARCH models. Expert Syst. Appl. 2010, 37, 4928–4934. [Google Scholar] [CrossRef]
  17. Lv, X.; Shan, X. Modeling natural gas market volatility using GARCH with different distributions. Phys. A 2013, 392, 5685–5699. [Google Scholar] [CrossRef]
  18. Su, J.B. How candlestick features affect the performance of volatility forecasts: Evidence from the stock market. Eur. J. Financ. 2015, 21, 486–506. [Google Scholar] [CrossRef]
  19. Zhang, Y.J.; Yao, T.; He, L.Y.; Ripple, R. Volatility forecasting of crude oil market: Can the regime switching GARCH model beat the single-regime GARCH models? Int. Rev. Econ. Financ. 2019, 59, 302–317. [Google Scholar] [CrossRef]
  20. Christie, A. The stochastic behavior of common stock variances: Value, leverage, and interest rate effects. J. Financ. Econ. 1982, 10, 407–432. [Google Scholar] [CrossRef]
  21. Chu, G.; Li, X.; Shen, D.; Zhang, Y. Unexpected information demand and volatility clustering of Chinese stock returns: Evidence from Baidu index. Entropy 2020, 22, 44. [Google Scholar] [CrossRef] [Green Version]
  22. Inglada-Perez, L. A comprehensive framework for uncovering non-linearity and Chaos in financial markets: Empirical evidence for four major stock market indices. Entropy 2020, 22, 1435. [Google Scholar] [CrossRef]
  23. Mandelbrot, B. The variation of certain speculative prices. J. Bus. 1963, 36, 394–419. [Google Scholar] [CrossRef]
  24. Glosten, L.R.; Jagannathan, R.; Runkle, D.E. On the relation between the expected value and the volatility of the nominal excess return on stocks. J. Financ. 1993, 48, 1779–1801. [Google Scholar] [CrossRef]
  25. Theodossiou, P. Financial data and the skewed generalized t distribution. Manag. Sci. 1998, 44, 1650–1661. [Google Scholar] [CrossRef]
  26. Ling, S.; McAleer, M. Stationarity and the existence of moments of a family of GARCH process. J. Econom. 2002, 106, 109–117. [Google Scholar] [CrossRef] [Green Version]
  27. Bali, T.G.; Theodossiou, P. A conditional-SGT-VaR approach with alternative GARCH models. Ann. Oper. Res. 2007, 151, 241–267. [Google Scholar] [CrossRef]
  28. Lu, X.; Que, D.; Cao, G. Volatility forecast based on the hybrid artificial neural network and GARCH-type models. Procedia Comput. Sci. 2016, 91, 1044–1049. [Google Scholar] [CrossRef] [Green Version]
  29. Doan, T.A. RATS User’s Guide, 6th ed.; Estima: Fullerton, CA, USA, 2004; pp. 448–455. [Google Scholar]
  30. Jarque, C.M.; Bera, A.K. A test for normality of observations and regression residuals. Int. Stat. Rev. 1987, 55, 163–172. [Google Scholar] [CrossRef]
Figure 1. The structure of backpropagation neural network. (a) One backpropagation neural networks with an input layer including two input nodes ( x 1   and   x 2 ), a hidden layer of three hidden nodes ( h ˜ 1 , h ˜ 2   and   h ˜ 3 ), and an output layer containing an output node ( y ). (b) One backpropagation neural networks with an input layer including one input node ( x ), a hidden layer of five hidden nodes ( h ˜ 1 , h ˜ 2 , h ˜ 3 , h ˜ 4 , and   h ˜ 5 ), and an output layer containing one output node ( y ).
Figure 1. The structure of backpropagation neural network. (a) One backpropagation neural networks with an input layer including two input nodes ( x 1   and   x 2 ), a hidden layer of three hidden nodes ( h ˜ 1 , h ˜ 2   and   h ˜ 3 ), and an output layer containing an output node ( y ). (b) One backpropagation neural networks with an input layer including one input node ( x ), a hidden layer of five hidden nodes ( h ˜ 1 , h ˜ 2 , h ˜ 3 , h ˜ 4 , and   h ˜ 5 ), and an output layer containing one output node ( y ).
Entropy 23 01151 g001
Figure 2. The schematic diagram for the volatility forecasting process of a neural network approach.
Figure 2. The schematic diagram for the volatility forecasting process of a neural network approach.
Entropy 23 01151 g002
Figure 3. The price levels and returns of spot prices for overall period.
Figure 3. The price levels and returns of spot prices for overall period.
Entropy 23 01151 g003aEntropy 23 01151 g003b
Figure 4. The actual variance and its two out-of-sample forecasts (the GJR-GARCH-SGT and GJR-GARCH-SGT-NN models).
Figure 4. The actual variance and its two out-of-sample forecasts (the GJR-GARCH-SGT and GJR-GARCH-SGT-NN models).
Entropy 23 01151 g004
Table 1. Descriptive statistics of daily return for the G7 and E7 stock indices.
Table 1. Descriptive statistics of daily return for the G7 and E7 stock indices.
MeanSDMax.Min.SKKURJ-B Q 2 20 Obs.
Panel A. Overall period
DJ0.02631.356310.764−12.577−0.6470 c11.298 c18,866.79 c4031.63 c3501
TSX0.01811.263111.294−16.998−1.3781 c24.159 c86,251.62 c2302.80 c3501
FTSE7 × 10−51.397710.870−11.512−0.3820 c8.844 c11,495.47 c2956.18 c3501
CAC40−0.00391.680911.285−13.098−0.3211 c6.579 c6374.96 c2364.63 c3501
DAX0.01981.738211.588−13.054−0.3698 c6.076 c5466.49 c2320.15 c3501
MIB−0.02151.830613.761−18.541−0.8287 c10.120 c15,341.27 c1130.00 c3501
N2250.01461.701511.805−12.715−0.3951 c6.136 c5584.40 c1438.09 c3501
BVSP0.05192.088513.022−18.749−0.5818 c7.231 c7826.8 c2207.54 c3501
RTSI0.06652.315515.642−20.779−0.5733 c8.425 c10,546.1 c803.26 c3501
BSE0.06331.725914.412−17.183−0.5704 c11.296 c18,804.7 c777.52 c3501
SSE0.02541.999645.160−12.7633.1136 c77.546 c882,866.9 c8.473501
MXX0.05291.436110.153−16.277−0.5429 c10.026 c14,837.2 c1372.47 c3501
JKSE0.07171.615013.624−12.925−0.3663 c11.109 c18,081.8 c922.40 c3501
XU1000.06852.276415.642−20.330−0.4234 c8.447 c10,514.0 c582.22 c3501
Panel B. Estimation period
DJ0.02801.25547.877−11.269−0.5146 c6.8384 c6478.18 c2837.12 c3251
TSX0.01991.18008.709−16.998−1.1476 c20.079 c55,328.6 c1291.52 c3251
FTSE0.00621.350610.870−10.327−0.1799 c7.4255 c7486.55 c3469.63 c3251
CAC40−0.00101.644111.285−11.476−0.1231 c5.4008 c3959.36 c2637.30 c3251
DAX0.01921.705811.588−11.829−0.2596 c4.8478 c3220.02 c2745.92 c3251
MIB−0.02061.788813.761−13.331−0.4568 c6.8566 c6481.45 c1414.19 c3251
N2250.01441.699511.805−12.715−0.3750 c6.2068 c5294.73 c1379.55 c3251
BVSP0.05372.018210.967−18.749−0.4689 c5.598 c4365.59 c1502.78 c3251
RTSI0.06782.338315.642−20.779−0.5240 c8.347 c9587.39 c805.08 c3251
BSE0.06881.688614.412−17.183−0.5087 c11.009 c16,558.95 c632.53 c3251
SSE0.01331.876512.950−12.763−0.1915 c5.443 c4033.75 c914.03 c3251
MXX0.06271.422710.153−16.277−0.5190 c10.645 c15,495.60 c1239.05 c3251
JKSE0.08421.607913.624−12.925−0.3487 c11.231 c17,152.38 c810.85 c3251
XU1000.06732.306715.642−20.330−0.3862 c8.431 c9710.74 c537.72 c3251
Panel C. Forecasting period
DJ0.00302.294810.764−12.577−0.7590 c10.203 c1113.02 c355.686 c251
TSX−0.00362.057711.294−13.176−1.6405 c18.926 c3859.00 c376.123 c251
FTSE−0.07671.90638.666−11.512−1.2277 c10.858 c1296.15 c198.535 c251
CAC40−0.04012.10098.056−13.098−1.5277 c11.175 c1403.71 c171.667 c251
DAX0.02732.115610.414−13.054−1.1140 c12.113 c1586.48 c144.527 c251
MIB−0.03222.30588.549−18.541−3.0502 c24.071 c6449.01 c96.842 c251
N2250.01771.72696.889−7.142−0.6490 c5.439 c327.02 c109.667 c251
BVSP0.03242.849213.022−15.993−1.0412 c10.454 c1188.4 c343.631 c251
RTSI0.05871.99868.280−12.049−1.6007 c9.016 c957.3 c15.876251
BSE−0.00772.152111.573−11.3274−0.8935 c10.961 c1290.0 c113.595 c251
SSE0.18123.189845.160−6.712411.2323 c159.326 c270,762 c0.325251
MXX−0.06051.60955.618−7.3361−0.6836 c4.408 c222.7 c206.215 c251
JKSE−0.09501.69799.704−9.6359−0.5312 c10.033 c1064.5 c193.267 c251
XU1000.07781.84115.810−9.6871−1.3201 c5.618 c403.0 c80.058 c251
Notes: 1. The superscripts a, b and c denote significance at the 10%, 5% and 1% levels, respectively. 2. SK and KUR denote the skewness and excess kurtosis, respectively. 3. J-B statistics are based on Jarque and Bera (1987) and are asymptotically chi-squared-distributed with 2 degrees of freedom. 4. Q 2   20 statistics are asymptotically chi-squared-distributed with 20 degrees of freedom. 5. Obs. denotes the number of observations. 6. The overall period starts from 10 January 2001 to 27 July 2020. On the other hand, the estimation period is from 10 January 2001 to 8 April 2019 whereas the forecasting period is from 9 April 2019 to 27 July 2020. 7. Bold font in column ‘Mean’ (resp. ‘SD’) of each panel denotes the largest value when all seven numbers in column ‘Mean’ (resp. ‘SD’) corresponding to the G7 or E7 are compared with each other. 8. Underline font in column ‘Mean’ (resp. ‘SD’) of each panel denotes the smallest value when all seven numbers in column ‘Mean’ (resp. ‘SD’) corresponding to the G7 or E7 are compared with each other.
Table 2. The empirical results of the GJR-GARCH-SGT model for the overall period.
Table 2. The empirical results of the GJR-GARCH-SGT model for the overall period.
DJTSXFTSECAC40DAXMIBN225
ϕ 0 0.0391 (0.013) c0.0367 (0.012) c0.0018 (0.015)0.0187 (0.018)0.0377 (0.019) a0.0177 (0.020)0.0380 (0.022) a
ϕ 1 −0.0472 (0.016) c0.0037 (0.016)−0.0306 (0.016) a−0.0303 (0.016) a0.0012 (0.016)−0.0217 (0.016)−0.0252 (0.016)
ω 0.0290 (0.004) c0.0227 (0.004) c0.0309 (0.005) c0.0398 (0.004) c0.0387 (0.006) c0.0326 (0.006) c0.0861 (0.016) c
α −0.0218 (0.008) c0.0241 (0.009) b−0.0168 (0.007) b−0.0222 (0.006) c−0.0162 (0.005) c−0.0007 (0.008)0.0232 (0.009) b
β 0.8929 (0.010) c0.8837 (0.011) c0.9008 (0.012) c0.9026 (0.003) c0.9121 (0.008) c0.9150 (0.010) c0.8732 (0.014) c
η 0.2176 (0.021) c0.1421 (0.020) c0.1943 (0.022) c0.2046 (0.010) c0.1759 (0.017) c0.1433 (0.015) c0.1488 (0.023) c
C0.97990.978850.981150.98270.983850.985950.9708
n 7.0798 (1.423) c5.1560 (0.701) c5.7090 (0.905) c5.9134 (0.917) c9.5401 (2.326) c5.7027 (0.886) c7.1754 (1.542) c
λ −0.1229 (0.022) c−0.1548 (0.023) c−0.1143 (0.023) c−0.1380 (0.023) c−0.1090 (0.021) c−0.1445 (0.023) c−0.0554 (0.021) b
κ 1.7311 (0.126) c2.1707 (0.162) c2.1041 (0.163) c2.1338 (0.159) c1.6614 (0.119) c2.0680 (0.157) c1.7058 (0.139) c
Q2 (20)3.38039.941 c12.51010.8139.60548.903 c7.029
LL−4974.32−4631.04−5207.87−5922.27−6058.64−6215.33−6306.45
BVSPRTSIBSESSEMXXJKSEXU100
ϕ 0 0.0565 (0.026) b0.0991 (0.026) c0.0778 (0.017) c0.0116 (0.018)0.0426 (0.018) b0.0850 (0.016) c0.0958 (0.030) c
ϕ 1 −0.0144 (0.016)0.0102 (0.015)0.0463 (0.016) c0.0142 (0.011)0.0460 (0.016) c−0.0024 (0.015)−0.0004 (0.016)
ω 0.0946 (0.020) c0.1205 (0.011) c0.0671 (0.012) c0.0351 (0.011)c0.0360 (0.005) c0.0613 (0.015) c0.0960 (0.012) c
α 0.0143 (0.007) a0.0449 (0.004) c0.0363 (0.009) c0.0460 (0.009)c0.0296 (0.008) c0.0698 (0.013) c0.0425 (0.006) c
β 0.9197 (0.010) c0.9002 (0.002) c0.8701 (0.013) c0.9235 (0.012)c0.9089 (0.003) c0.8706 (0.016) c0.9161 (0.003) c
η 0.0832 (0.015) c0.0672 (0.007) c0.1534 (0.026) c0.0520 (0.016)c0.0873 (0.012) c0.0824 (0.021) c0.0455 (0.013) c
C0.97560.97870.98310.99550.982150.98160.98135
n 4.9911 (0.709) c4.2549 (0.547) c3.8139 (0.386) c6.9225 (1.712) c4.8618 (0.651) c4.5756 (0.669) c4.1622 (0.503) c
λ −0.0808 (0.022) c−0.0262 (0.020)−0.0428 (0.020) b−0.0070 (0.007)−0.0784 (0.023) c−0.0603 (0.018) c−0.0593 (0.022) c
κ 2.1926 (0.180) c1.9903 (0.165) c2.4594 (0.220) c1.3106 (0.104) c2.0735 (0.158) c1.7512 (0.146) c2.1166 (0.169) c
Q2 (20)27.4384.7582.46531.257 a21.58113.7416.592
LL−7048.77−7280.81−6022.33−6451.92−5582.73−5838.77−7254.09
Note: 1. The superscripts a, b and c denote significance at the 10%, 5% and 1% levels, respectively. 2. Numbers in parentheses at the rows of parameters are standard errors. 3. LL indicates the log-likelihood value. 4. Q2 (20) statistics are asymptotically chi-squared-distributed with 20 degrees of freedom. 5. The numbers in row ‘C’ denote the value of the expression ‘ β + α + 0.5 η ‘ and are used to check the constraint ‘ β + α + 0.5 η < 1 , the necessary and sufficient condition for the existence of the second moment condition for all GJR-based models. Please see Example 2.1 of Ling and McAleer [26] for more details.
Table 3. The results of in-sample volatility forecasts based on the MAE loss function for the overall period.
Table 3. The results of in-sample volatility forecasts based on the MAE loss function for the overall period.
DJTSXFTSECAC40DAXMIBN225BVSPRTSIBSESSEMXXJKSEXU100S1S2S3S4
Panel A. Parametric volatility forecasting model
G-n1.943011.675922.100583.024613.229763.791263.207124.611726.044523.364025.507372.241333.126755.817864470
G-t1.998081.720062.117733.067663.269713.736403.239954.654636.081473.366805.352932.242153.176115.811971400
G-st1.980321.706772.110963.046093.249803.706463.227534.649096.070693.361795.357352.240803.170685.799203530
G-sgt1.959431.709202.109963.044633.237863.701503.202704.662656.064733.388825.074302.243763.115945.820813540
GJR-n1.909321.647812.046732.970393.199643.791633.188554.569136.103943.43455.351302.212153.099895.82530210102
GJR-t1.932221.693832.065063.017713.236443.750083.237764.616846.152223.40735.212952.223283.144035.8246201000
GJR-st1.929621.689062.068173.014583.230193.736373.227874.616366.142803.40615.218592.226403.144685.816431930
GJR-sgt1.917561.695152.072053.019553.224193.740663.208024.633356.140093.45854.940372.230543.100015.843870910
Panel B. Composed volatility forecasting model
G-n-NN1.943121.713902.113613.016333.218023.647543.172104.658496.014013.334014.874502.238293.036185.7818410050
G-t-NN1.947261.715092.112153.010603.217023.651683.171564.667796.005473.338304.830832.237063.028265.7787913281
G-st-NN1.949741.717862.113023.013243.219943.655783.172754.672516.006913.338984.830202.237953.029365.7798011210
G-sgt-NN1.949581.718122.113483.013213.219883.655673.172884.679576.006913.338894.833812.237913.028965.7798911200
GJR-n-NN1.881921.681692.033212.933903.143093.593973.141764.597186.012553.32464.828762.209813.012935.77303121453
GJR-t-NN1.855251.693632.034632.932733.142973.598303.113034.602196.008003.30044.861492.209713.004185.76911141265
GJR-st-NN1.855101.693492.034042.932653.143573.600063.112004.603966.008373.30074.858112.209803.005265.76948131211
GJR-sgt-NN1.854891.693502.034002.932923.143333.600153.111654.604266.008333.30044.855522.209833.005515.76948141233
Note: 1. G-n, G-t, G-st, G-sgt, GJR-n, GJR-t, GJR-st and GJR-sgt respectively denote the GARCH-N, GARCH-T, GARCH-ST, GARCH-SGT, GJR-GARCH-N, GJR-GARCH-T, GJR-GARCH-ST and GJR-GARCH-SGT models. On the other hand, G-n-NN, G-t-NN, G-st-NN, G-sgt-NN, GJR-n-NN, GJR-t-NN, GJR-st-NN and GJR-sgt-NN respectively denote the GARCH-N, GARCH-T, GARCH-ST, GARCH-SGT, GJR-GARCH-N, GJR-GARCH-T, GJR-GARCH-ST and GJR-GARCH-SGT models with a neural network approach. 2. Bold font denotes the lower value of loss function when the predictive accuracies of one parametric volatility forecasting model and one corresponding composed volatility forecasting model are compared with each other based on the same volatility specification and return distribution. 3. Italic font denotes the lower value of loss function when the predictive accuracies of one GARCH-based model and one corresponding GJR-based model are compared with each other based on the same volatility forecasting approach and return distribution. 4. Underline font denotes the lowest value of loss function when the predictive accuracies of a group of models are compared with each other based on the same volatility specification and volatility forecasting approach but the different distribution. 5. Shade font denotes the lowest value of loss function when the predictive accuracies of all sixteen models are compared with each other. 6. The number in column S1 denotes the total number of indices that a specific model possesses, the lower value of loss function when the predictive accuracies of two models are compared with each other as shown by note 2. This result can illustrate the performance comparison of two volatility forecasting models with different volatility forecasting approach, a parametric or composed volatility forecasting approach. 7. The number in column S2 denotes the total number of indices that a specific model possesses, the lower value of loss function when the predictive accuracies of two models are compared with each other as shown by note 3. This result can illustrate the performance comparison of two volatility forecasting models with different volatility specifications, the GJR-GARCH or GARCH. 8. The number in column S3 denotes the total number of indices that a specific model possesses, the lowest value of loss function when the predictive accuracies of four models are compared with each other as shown by note 4. This result can illustrate the performance comparison of a group of volatility forecasting models with different return distributions, the normal, Student’s t, ST, and SGT. 9. The number in column S4 denotes the total number of indices that a specific model possesses, the lowest value of loss function when the predictive accuracies of all sixteen models are compared with each other as shown by note 5.
Table 4. The results of in-sample volatility forecasts based on the RMSE loss function for the overall period.
Table 4. The results of in-sample volatility forecasts based on the RMSE loss function for the overall period.
DJTSXFTSECAC40DAXMIBN225BVSPRTSIBSESSEMXXJKSEXU100S1S2S3S4
Panel A. Parametric volatility forecasting model
G-n5.889877.460735.895147.656917.9684611.21877.7820512.247616.640510.422236.83346.780589.0755616.215003100
G-t5.929257.491195.903157.675667.9751811.19617.7909212.275716.659410.445836.15816.780759.0713216.21140400
G-st5.925797.490695.903797.672987.9729611.18897.7905912.287016.658710.445136.16026.782409.0713816.21130410
G-sgt5.917197.491415.903557.672577.9712011.18717.7867512.286316.658110.450236.07886.782699.0573616.21300430
GJR-n5.704247.382305.779537.571927.8695611.25217.6850112.062216.689210.433636.27586.747638.9801616.2136011110
GJR-t5.743837.438025.802087.615117.8825011.21017.7097612.092516.708410.489935.81146.751028.9693716.218301010
GJR-st5.747907.438475.807587.623077.8826511.21707.7058312.100916.705810.490735.81656.753448.9734916.218601000
GJR-sgt5.738787.443765.809797.625847.8825111.21957.6986012.102216.705310.519835.78926.754238.9618316.222301020
Panel B. Composed volatility forecasting model
G-n-NN5.889737.444495.815967.616567.9217511.06077.7694512.205816.627910.371435.54756.771049.0215416.1876141101
G-t-NN5.908287.455535.824817.623017.9246811.09697.7713112.255016.635310.381135.53336.770549.0113516.183814030
G-st-NN5.911837.462335.827017.627667.9284711.09897.7732712.265616.636010.382135.53296.772009.0131316.184414000
G-sgt-NN5.909987.461885.826497.627487.9274111.09787.7728712.261816.636010.381835.53276.771979.0123216.184514010
GJR-n-NN5.668387.329585.734527.471307.7946610.96587.6484712.041216.630010.370435.49886.730878.9436416.1678141376
GJR-t-NN5.656907.344515.730457.473387.7933410.99647.5469412.092016.634710.334635.50856.731468.9166416.1633141433
GJR-st-NN5.655757.342915.729907.473117.7939010.99817.5450012.100716.635110.335235.50836.731698.9194416.1634141400
GJR-sgt-NN5.654947.341845.729157.472967.7934410.99857.5446812.102116.635110.334335.50446.731768.9206416.1634141444
Note: 1. G-n, G-t, G-st, G-sgt, GJR-n, GJR-t, GJR-st and GJR-sgt respectively denote the GARCH-N, GARCH-T, GARCH-ST, GARCH-SGT, GJR-GARCH-N, GJR-GARCH-T, GJR-GARCH-ST and GJR-GARCH-SGT models. On the other hand, G-n-NN, G-t-NN, G-st-NN, G-sgt-NN, GJR-n-NN, GJR-t-NN, GJR-st-NN and GJR-sgt-NN respectively denote the GARCH-N, GARCH-T, GARCH-ST, GARCH-SGT, GJR-GARCH-N, GJR-GARCH-T, GJR-GARCH-ST and GJR-GARCH-SGT models with a neural network approach. 2. Bold font denotes the lower value of loss function when the predictive accuracies of one parametric volatility forecasting model and one corresponding composed volatility forecasting model are compared with each other based on the same volatility specification and return distribution. 3. Italic font denotes the lower value of loss function when the predictive accuracies of one GARCH-based model and one corresponding GJR-based model are compared with each other based on the same volatility forecasting approach and return distribution. 4. Underline font denotes the lowest value of loss function when the predictive accuracies of a group of models are compared with each other based on the same volatility specification and volatility forecasting approach but the different distribution. 5. Shade font denotes the lowest value of loss function when the predictive accuracies of all sixteen models are compared with each other. 6. The number in column S1 denotes the total number of indices that a specific model possesses, the lower value of loss function when the predictive accuracies of two models are compared with each other as shown by note 2. This result can illustrate the performance comparison of two volatility forecasting models with different volatility forecasting approach, a parametric or composed volatility forecasting approach. 7. The number in column S2 denotes the total number of indices that a specific model possesses, the lower value of loss function when the predictive accuracies of two models are compared with each other as shown by note 3. This result can illustrate the performance comparison of two volatility forecasting models with different volatility specifications, the GJR-GARCH or GARCH. 8. The number in column S3 denotes the total number of indices that a specific model possesses, the lowest value of loss function when the predictive accuracies of four models are compared with each other as shown by note 4. This result can illustrate the performance comparison of a group of volatility forecasting models with different return distributions, the normal, Student’s t, ST, and SGT. 9. The number in column S4 denotes the total number of indices that a specific model possesses, the lowest value of loss function when the predictive accuracies of all sixteen models are compared with each other as shown by note 5.
Table 5. The results of out-of-sample volatility forecasts based on the MAE loss function for the forecasting period.
Table 5. The results of out-of-sample volatility forecasts based on the MAE loss function for the forecasting period.
DJTSXFTSECAC40DAXMIBN225BVSPRTSIBSESSEMXXJKSEXU100S1S2S3S4
Panel A. Parametric volatility forecasting model
G-n5.238354.086144.007234.887465.255296.873913.580238.454085.210865.5440919.71192.538843.218634.251123880
G-t5.487814.345794.107545.058975.415906.803223.619238.705835.201985.5735517.80972.535973.284064.281564900
G-st5.450564.336984.099355.040405.394626.770973.606258.745395.189255.5775617.76332.533553.276784.274884910
G-sgt5.393034.332204.094075.029315.364386.759013.574348.751805.184625.6196717.17152.537163.206564.282063950
GJR-n4.813453.771814.369835.203515.440527.414043.533138.145525.457335.6659517.86802.523633.222404.3166346133
GJR-t4.974633.995254.406015.388835.602227.476803.568408.405495.573115.7450015.82102.546203.286134.373883500
GJR-st4.978683.994444.435725.420695.620177.518993.560418.460825.565375.7602415.80782.550583.288664.373283500
GJR-sgt4.941633.997924.440685.427105.603097.524523.542978.477175.565925.8556515.48292.559563.232124.390854510
Panel B. Composed volatility forecasting model
G-n-NN5.241354.265283.921624.689674.941446.132393.543648.781855.148615.1494411.61752.538733.037694.23808118126
G-t-NN5.344394.394024.005364.724235.031076.204663.532079.006705.182405.3023911.62622.545943.037224.3038210610
G-st-NN5.352384.432733.995404.781525.009096.239983.579059.175925.169465.3035611.37412.545483.065664.2930110811
G-sgt-NN5.369604.453644.024414.781245.045966.205043.541739.074225.166165.3094511.59822.547083.113894.2802811700
GJR-n-NN4.401753.945534.009384.729844.941266.244503.606948.362625.390255.2354911.89652.528822.996074.29655106113
GJR-t-NN4.553684.082224.044344.746605.029596.396463.495998.681415.475665.3331211.51892.545443.021024.3976311821
GJR-st-NN4.615194.079344.114414.742465.049896.489823.500908.739675.474155.3398312.54742.541553.165454.3912211600
GJR-sgt-NN4.604644.071514.116784.701485.004596.588093.590518.811775.478065.3347711.67272.540633.015034.3973410710
Note: 1. G-n, G-t, G-st, G-sgt, GJR-n, GJR-t, GJR-st and GJR-sgt respectively denote the GARCH-N, GARCH-T, GARCH-ST, GARCH-SGT, GJR-GARCH-N, GJR-GARCH-T, GJR-GARCH-ST and GJR-GARCH-SGT models. On the other hand, G-n-NN, G-t-NN, G-st-NN, G-sgt-NN, GJR-n-NN, GJR-t-NN, GJR-st-NN and GJR-sgt-NN respectively denote the GARCH-N, GARCH-T, GARCH-ST, GARCH-SGT, GJR-GARCH-N, GJR-GARCH-T, GJR-GARCH-ST and GJR-GARCH-SGT models with a neural network approach. 2. Bold font denotes the lower value of loss function when the predictive accuracies of one parametric volatility forecasting model and one corresponding composed volatility forecasting model are compared with each other based on the same volatility specification and return distribution. 3. Italic font denotes the lower value of loss function when the predictive accuracies of one GARCH-based model and one corresponding GJR-based model are compared with each other based on the same volatility forecasting approach and return distribution. 4. Underline font denotes the lowest value of loss function when the predictive accuracies of a group of models are compared with each other based on the same volatility specification and volatility forecasting approach but the different distribution. 5. Shade font denotes the lowest value of loss function when the predictive accuracies of all sixteen models are compared with each other. 6. The number in column S1 denotes the total number of indices that a specific model possesses, the lower value of loss function when the predictive accuracies of two models are compared with each other as shown by note 2. This result can illustrate the performance comparison of two volatility forecasting models with different volatility forecasting approach, a parametric or composed volatility forecasting approach. 7. The number in column S2 denotes the total number of indices that a specific model possesses, the lower value of loss function when the predictive accuracies of two models are compared with each other as shown by note 3. This result can illustrate the performance comparison of two volatility forecasting models with different volatility specifications, the GJR-GARCH or GARCH. 8. The number in column S3 denotes the total number of indices that a specific model possesses, the lowest value of loss function when the predictive accuracies of four models are compared with each other as shown by note 4. This result can illustrate the performance comparison of a group of volatility forecasting models with different return distributions, the normal, Student’s t, ST, and SGT. 9. The number in column S4 denotes the total number of indices that a specific model possesses, the lowest value of loss function when the predictive accuracies of all sixteen models are compared with each other as shown by note 5.
Table 6. The results of out-of-sample volatility forecasts based on the RMSE loss function for the forecasting period.
Table 6. The results of out-of-sample volatility forecasts based on the RMSE loss function for the forecasting period.
DJTSXFTSECAC40DAXMIBN225BVSPRTSIBSESSEMXXJKSEXU100S1S2S3S4
Panel A. Parametric volatility forecasting model
G-n15.761616.435712.202814.953115.845326.28277.6526025.207313.047515.5771133.4826.003909.186198.7904653101
G-t15.832616.544412.255815.049115.900126.32767.6575625.320813.090115.6034130.9535.996779.175358.843556410
G-st15.834516.579312.256315.055615.905426.31987.6568925.399413.087415.6054130.9205.997529.178648.846926400
G-sgt15.817216.579512.255015.051415.898526.31377.6522625.397613.086715.6098130.7125.997719.164108.848336430
GJR-n14.708215.463712.195715.146015.664626.73127.5435124.108013.219815.3434131.8525.947918.995128.73118611114
GJR-t14.823415.632812.264815.419315.753526.76657.5834624.172313.320115.2647129.8955.958278.954438.7751751011
GJR-st14.821815.623412.299915.473715.769826.83197.5757424.256013.314215.2693129.8885.961448.961298.7760651000
GJR-sgt14.805115.623612.304815.481515.766026.83667.5685324.252213.314415.2992129.8575.965658.947538.7810351021
Panel B. Composed volatility forecasting model
G-n-NN15.715816.795411.935714.865815.720325.91997.6486625.235513.073415.5384128.7806.010519.169538.8495192100
G-t-NN15.809617.029012.176714.858615.765126.02827.6502525.724013.093215.6256128.7756.013319.161248.897378110
G-st-NN15.820517.090912.034214.923815.638526.00647.6616425.757013.094715.5870128.7426.014629.151288.902858220
G-sgt-NN15.887217.090612.186314.910915.769626.02387.6504625.795913.087715.6027128.7366.014869.154798.903208211
GJR-n-NN14.588016.096111.605214.509515.307425.70607.5809324.299013.151915.3102128.8135.949589.085638.7672381293
GJR-t-NN14.719716.260111.727214.503415.350225.83607.5258524.846413.194715.3500128.7415.952299.053188.8224491310
GJR-st-NN14.689616.282611.839714.459715.356325.81057.5258024.949913.199515.3505128.9305.949968.990938.8216791221
GJR-sgt-NN14.679116.271711.854014.258315.176525.80097.5506624.954713.201415.3489128.7595.950119.052888.8230491222
Note: 1. G-n, G-t, G-st, G-sgt, GJR-n, GJR-t, GJR-st and GJR-sgt respectively denote the GARCH-N, GARCH-T, GARCH-ST, GARCH-SGT, GJR-GARCH-N, GJR-GARCH-T, GJR-GARCH-ST and GJR-GARCH-SGT models. On the other hand, G-n-NN, G-t-NN, G-st-NN, G-sgt-NN, GJR-n-NN, GJR-t-NN, GJR-st-NN and GJR-sgt-NN respectively denote the GARCH-N, GARCH-T, GARCH-ST, GARCH-SGT, GJR-GARCH-N, GJR-GARCH-T, GJR-GARCH-ST and GJR-GARCH-SGT models with a neural network approach. 2. Bold font denotes the lower value of loss function when the predictive accuracies of one parametric volatility forecasting model and one corresponding composed volatility forecasting model are compared with each other based on the same volatility specification and return distribution. 3. Italic font denotes the lower value of loss function when the predictive accuracies of one GARCH-based model and one corresponding GJR-based model are compared with each other based on the same volatility forecasting approach and return distribution. 4. Underline font denotes the lowest value of loss function when the predictive accuracies of a group of models are compared with each other based on the same volatility specification and volatility forecasting approach but the different distribution. 5. Shade font denotes the lowest value of loss function when the predictive accuracies of all sixteen models are compared with each other. 6. The number in column S1 denotes the total number of indices that a specific model possesses, the lower value of loss function when the predictive accuracies of two models are compared with each other as shown by note 2. This result can illustrate the performance comparison of two volatility forecasting models with different volatility forecasting approach, a parametric or composed volatility forecasting approach. 7. The number in column S2 denotes the total number of indices that a specific model possesses, the lower value of loss function when the predictive accuracies of two models are compared with each other as shown by note 3. This result can illustrate the performance comparison of two volatility forecasting models with different volatility specifications, the GJR-GARCH or GARCH. 8. The number in column S3 denotes the total number of indices that a specific model possesses, the lowest value of loss function when the predictive accuracies of four models are compared with each other as shown by note 4. This result can illustrate the performance comparison of a group of volatility forecasting models with different return distributions, the normal, Student’s t, ST, and SGT. 9. The number in column S4 denotes the total number of indices that a specific model possesses, the lowest value of loss function when the predictive accuracies of all sixteen models are compared with each other as shown by note 5.
Table 7. The summary results of performance comparison for the in-sample and out-of-sample volatility forecasts.
Table 7. The summary results of performance comparison for the in-sample and out-of-sample volatility forecasts.
S1S2S3S4
In-SampleOut-of-SampleSumIn-SampleOut-of-SampleSumIn-SampleOut-of-SampleSumIn-SampleOut-of-SampleSum
MAERMSEMAERMSEMAERMSEMAERMSEMAERMSEMAERMSEMAERMSEMAERMSE
Panel A. Parametric volatility forecasting model
G-n4035124383187108103500011
G-t1046114494210001100000
G-st3046135494223110500000
G-sgt30361254942243531500000
GJR-n204612101161138101113114520349
GJR-t003581010510350101200011
GJR-st10359910510343000300000
GJR-sgt00459910510341212600011
Panel B. Composed volatility forecasting model
G-n-NN10141194401821151012103701607
G-t-NN1314108452061983111310001
G-st-NN1114108432082121012400101
G-sgt-NN1114118442072110101200011
GJR-n-NN1214108441413612455711932363315
GJR-t-NN14141194812148134763211253109
GJR-st-NN1314119471214612441002310012
GJR-sgt-NN14141094712147124534121034029
Note: 1. The numbers in column ‘MAE’ below ‘In-sample’ of S1, S2, S3 and S4, respectively are summarized from those in columns ‘S1’, ‘S2’, ‘S3’ and ‘S4’ of Table 3. 2. The numbers in column ‘RMSE’ below ‘In-sample’ of S1, S2, S3 and S4, respectively are summarized from those in columns ‘S1’, ‘S2’, ‘S3’ and ‘S4’ of Table 4. 3. The numbers in column ‘MAE’ below ‘Out-of-sample’ of S1, S2, S3 and S4, respectively are summarized from those in columns ‘S1’, ‘S2’, ‘S3’ and ‘S4’ of Table 5. 4. The numbers in column ‘RMSE’ below ‘Out-of-sample’ of S1, S2, S3 and S4, respectively are summarized from those in columns ‘S1’, ‘S2’, ‘S3’ and ‘S4’ of Table 6. 5. The numbers in column ‘Sum’ below ‘S1’ denote the summation of four numbers in columns ‘MAE’ and ‘RMSE’ below ‘In-sample’ and ‘Out-of-sample’ of S1. 6. The numbers in column ‘Sum’ below ‘S2’ denote the summation of four numbers in columns ‘MAE’ and ‘RMSE’ below ‘In-sample’ and ‘Out-of-sample’ of S2. 7. The numbers in column ‘Sum’ below ‘S3’ denote the summation of four numbers in columns ‘MAE’ and ‘RMSE’ below ‘In-sample’ and ‘Out-of-sample’ of S3. 8. The numbers in column ‘Sum’ below ‘S4’ denote the summation of four numbers in columns ‘MAE’ and ‘RMSE’ below ‘In-sample’ and ‘Out-of-sample’ of S4. 9. Bold font in column ‘S1’ denotes the larger number when two numbers at the same column are compared with each other where two numbers respectively correspond to a parametric volatility forecasting model and its corresponding composed volatility forecasting model, such as the GARCH-N and GARCH-N-N models. 10. Bold font in column ‘S2’ denotes the larger number when two numbers at the same column are compared with each other where two numbers respectively correspond to a GARCH-based model and its corresponding GJR-based model, such as the GARCH-N and GJR-GARCH-N models. 11. Bold font in column ‘S3’ denotes the largest number when four numbers at the same column are compared with each other where four numbers respectively correspond to four models with the same volatility specification and volatility forecasting approach but different return distributions. 12. Bold font in column ‘S4’ denotes the largest number when all sixteen numbers in the same column are compared with each other where sixteen numbers respectively correspond to sixteen models.
Table 8. The summary results of the best volatility forecasting model for alternative stock indices.
Table 8. The summary results of the best volatility forecasting model for alternative stock indices.
DJTSXFTSECAC40DAXMIBN225
In-SampleMAEGJR-sgt-NNGJR-nGJR-n-NNGJR-st-NNGJR-t-NNGJR-n-NNGJR-sgt-NN
RMSEGJR-sgt-NNGJR-n-NNGJR-sgt-NNGJR-n-NNGJR-t-NNGJR-n-NNGJR-sgt-NN
Out-of-SampleMAEGJR-n-NNGJR-nG-n-NNG-n-NNGJR-n-NNG-n-NNGJR-t-NN
RMSEGJR-n-NNGJR-nGJR-n-NNGJR-sgt-NNGJR-sgt-NNGJR-n-NNGJR-st-NN
Best model GJR-n-NN; GJR-sgt-NNGJR-nGJR-n-NN?GJR-t-NNGJR-n-NNGJR-sgt-NN
BVSPRTSIBSESSEMXXJKSEXU100
In-SampleMAEGJR-nG-t-NNGJR-t-NN; GJR-sgt-NNGJR-n-NNGJR-t-NNGJR-t-NNGJR-t-NN
RMSEGJR-n-NNG-n-NNGJR-sgt-NNGJR-n-NNGJR-n-NNGJR-t-NNGJR-t-NN
Out-of-SampleMAEGJR-nG-n-NNG-n-NNG-st-NNGJR-nGJR-n-NNG-n-NN
RMSEGJR-nG-nGJR-tG-sgt-NNGJR-nGJR-sgtGJR-n
Best model GJR-nG-n-NNGJR-sgt-NNGJR-n-NNGJR-nGJR-t-NNGJR-t-NN
Note: 1. This table summarizes the results of the best volatility forecasting model for alternative stock indices in Table 3, Table 4, Table 5 and Table 6. In other words, the results listed in row ‘MAE’ (respectively, ‘RMSE’) of ‘In-Sample’ in Table 8 are summarized from the results of the performance comparison for the fourth category of model in Table 3 (respectively, Table 4). On the other hand, the results listed in row ‘MAE’ (respectively, ‘RMSE’) of ‘Out-of-Sample’ in Table 8 are summarized from the results of the performance comparison for the fourth category of model in Table 5 (respectively, Table 6). 2. GJR-n represents the GJR-GARCH-N model. On the contrary, G-n-NN, GJR-n-NN, GJR-t-NN and GJR-sgt-NN denote the GARCH-N, GJR-GARCH-N, GJR-GARCH-T and GJR-GARCH-SGT models with a neural networks approach. 3. The models listed in row ‘Best model’ and the column corresponding to a specific stock index denote the suitable models for this specific stock index. 4. The symbol ‘?’ in the row ‘Best model’ denotes I could not obtain a suitable model because there is no model that can simultaneously appear twice or more. 5. The bold font in rows ‘MAE’ and ‘RMSE’ of ‘In-sample’ (or ‘Out-of-sample’) denote the suitable models for a specific stock index because for the above four cases this model can simultaneously appear twice or more.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Su, J.-B. How to Promote the Performance of Parametric Volatility Forecasts in the Stock Market? A Neural Networks Approach. Entropy 2021, 23, 1151. https://doi.org/10.3390/e23091151

AMA Style

Su J-B. How to Promote the Performance of Parametric Volatility Forecasts in the Stock Market? A Neural Networks Approach. Entropy. 2021; 23(9):1151. https://doi.org/10.3390/e23091151

Chicago/Turabian Style

Su, Jung-Bin. 2021. "How to Promote the Performance of Parametric Volatility Forecasts in the Stock Market? A Neural Networks Approach" Entropy 23, no. 9: 1151. https://doi.org/10.3390/e23091151

APA Style

Su, J. -B. (2021). How to Promote the Performance of Parametric Volatility Forecasts in the Stock Market? A Neural Networks Approach. Entropy, 23(9), 1151. https://doi.org/10.3390/e23091151

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop