Next Article in Journal
Improving DFT-Based Image Watermarking Using Particle Swarm Optimization Algorithm
Next Article in Special Issue
Ensemble Model of the Financial Distress Prediction in Visegrad Group Countries
Previous Article in Journal
Cross-Platform GPU-Based Implementation of Lattice Boltzmann Method Solver Using ArrayFire Library
Previous Article in Special Issue
Development of an Intelligent Decision Support System for Attaining Sustainable Growth within a Life Insurance Company
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Transformer: A New Neural Network-Based Architecture for Forecasting S&P Volatility

by
Eduardo Ramos-Pérez
1,
Pablo J. Alonso-González
2,* and
José Javier Núñez-Velázquez
2
1
Faculty of Economics, Universidad de Alcalá, Plaza de la Victoria 2, 28802 Alcalá de Henares, Madrid, Spain
2
Economics Department, Universidad de Alcalá, Plaza de la Victoria 2, 28802 Alcalá de Henares, Madrid, Spain
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(15), 1794; https://doi.org/10.3390/math9151794
Submission received: 26 June 2021 / Revised: 21 July 2021 / Accepted: 26 July 2021 / Published: 28 July 2021

Abstract

:
Events such as the Financial Crisis of 2007–2008 or the COVID-19 pandemic caused significant losses to banks and insurance entities. They also demonstrated the importance of using accurate equity risk models and having a risk management function able to implement effective hedging strategies. Stock volatility forecasts play a key role in the estimation of equity risk and, thus, in the management actions carried out by financial institutions. Therefore, this paper has the aim of proposing more accurate stock volatility models based on novel machine and deep learning techniques. This paper introduces a neural network-based architecture, called Multi-Transformer. Multi-Transformer is a variant of Transformer models, which have already been successfully applied in the field of natural language processing. Indeed, this paper also adapts traditional Transformer layers in order to be used in volatility forecasting models. The empirical results obtained in this paper suggest that the hybrid models based on Multi-Transformer and Transformer layers are more accurate and, hence, they lead to more appropriate risk measures than other autoregressive algorithms or hybrid models based on feed forward layers or long short term memory cells.

1. Introduction

Since the Financial Crisis of 2007–2008, financial institutions have enhanced their risk management framework in order to meet the new regulatory requirements set by Solvency II or Basel III. These regulations have the aim of measuring the risk profile of financial institutions and minimizing losses from unexpected events such as the European sovereign debt crisis or COVID-19 pandemic. Even though banks and insurance entities have reduced their losses thanks to the efforts made in the last years, unexpected events still cause remarkable losses to financial institutions. Thus, efforts are still required to further enhance market and equity risk models in which stock volatility forecasts play a fundamental role. Volatility, understood as a measure of an asset uncertainty [1,2], is not directly observed in stock markets. Thus, taking into consideration the stock market movements, a statistical model is applied in order to compute the volatility of a security.
GARCH-based models [3,4] are widely used for volatility forecasting purposes. This family of models is especially relevant because it takes into consideration the volatility clustering observed by [5]. Nevertheless, as the persistence of conditional variance tends to be close to zero, Refs. [6,7,8,9] developed more flexible variations of the traditional GARCH models. In addition, the models introduced by [10] (EGARCH) and [11] (GJR-GARCH) take into consideration that stocks volatility behaves differently depending on the market trend, bearish or bullish. Multivariate GARCH models were developed by [12,13]. Bollerslev et al. [14] applied the previous model to financial time series, while [15] introduced a time-varying multivariate GARCH. Dynamic conditional correlation GARCH, BEKK-GARCH and Factor-GARCH were other variants of this family that were developed by [16,17,18], respectively. Finally, it is worth mentioning that, in contrast to classical GARCH, the first-order zero-drift GARCH model (ZD-GARCH) proposed by [19] is non-stationary regardless of the sign of Lyapunov exponent and, thus, it can be used for studying heteroscedasticity and conditional heteroscedasticity together.
Another relevant family is composed by stochastic volatility models. As they assume that volatility follows its own stochastic process, these models are widely used in combination with Black–Scholes formula to assess derivatives price. The most popular process of this family is the [20] model which assumes that volatility follows an Cox-Ingersoll-Ross [21] process and stock returns a Brownian motion. The main challenge of the Heston model is the estimation of its parameters. Refs. [22,23] proposed a generalized method of moments to obtain the parameters of the stochastic process, while [24,25,26,27] used a simulation approach to estimate them. Other relevant stochastic volatility processes are Hull–White [28] and SABR [29] models.
The last relevant family is composed of those models based on machine and deep learning techniques. Even though GARCH models are considered part of the machine learning tool-kit, these models are considered another different family due to the significant importance that they have in the field of stock volatility. Thus, this family takes into consideration the models based on the rest of the machine and deep learning algorithms such as artificial neural networks [30], gradient boosting with regression trees [31], random forests [32] or support vector machines [33]. Refs. [34,35,36] applied machine learning techniques such as Support Vector Machines or hidded Markov models to forecast financial time series. Hamid and Iqbid [37] applied Artificial Neural Networks (ANNs) to demonstrate that the implied volatility forecasted by this algorithm is more accurate than Barone–Adesi and Whaley models.
ANNs have been also combined with other statistical models with the aim of improving the forecasting power of individual ANNs. The most common approach applied in the field of stocks volatility is merging GARCH-based models with ANNs. Refs. [38,39,40,41,42,43,44] developed different architectures based in the previous approach for stock volatility forecasting purposes. All these authors demonstrated that hybrid models overcome the performance of traditional GARCH models in the field of stock volatility forecasting. It is also worth mentioning the contribution of [45], who combined different GARCH models with ANNs in order to compare their predictive power. ANN-GARCH models have been also applied to forecast other financial time series such as metals [46,47] or oil [48,49] volatility. Apart from the combination with GARCH-based models, ANNs have been merged with other models for volatility forecasting purposes. Ramos-Pérez et al. [50] merged ANNs, random forests, support vector machines (SVM) and gradient boosting with regression trees in order to forecast S&P500 volatility. This model overcame the performance of a hybrid model based on feed forward layers and GARCH. Vidal and Kristjanpoller [51] proposed an architecture based on convolutional neural networks (CNNs) and long-short term memory (LSTM) units to forecast gold volatility. LSTMs were also used by [52] to forecast currency exchange rates volatility. It is also worth mentioning that GARCH models have not been only merged with ANNs, Peng et al. [53] combined SVM with GARCH-based models in order to predict cryptocurrencies volatility.
The aim of this paper is to introduce a more accurate stock volatility model based on an innovative machine and deep learning technique. For this purpose, hybrid models based on merging Transformer and Multi-Transformer layers with other approaches such as GARCH-based algorithms or LSTM units are introduced by this paper. Multi-Transformer layers, which are also introduced in this paper, are based on the Transformer architecture developed by [54]. Transformer layers have been successfully implemented in the field of natural language processing (NLP). Indeed, the models developed by [55,56] demonstrated that Transformer layers are able to overcome the performance of traditional NLP models. Thus, this recently developed architecture is currently considered the state-of-the-art in the field of NLP. In contrast to LSTM, Transformer layers do not incorporate recurrence in their structure. This novel structure relies on a multi-head attention mechanism and positional embeddings in order to forecast time series. As [54] developed Transformer for NLP purposes, positional embeddings are used in combination with word embeddings. The problem faced in this paper is the forecasting of stock volatility and, thus, the word embedding is not needed and the positional embedding has been modified as it is explained in Section 2.4.
In contrast to Transformer, Multi-Transformer randomly selects different subsets of training data and merges several multi-head attention mechanisms to produce the final output. Following the intuition of bagging, the aim of this architecture is to improve the stability and accurateness of the attention mechanism. It is worth mentioning that the GARCH-based algorithms used in combination with Transformer and Multi-Transformer layers are GARCH, EGARCH, GJR-GARCH, TrGARCH, FIGARCH and AVGARCH.
Therefore, three main contributions are provided by this study. First, Transformer layers are adapted in order to forecast stocks volatility. In addition, an extension of the previous structure is presented (Multi-Transformer). Second, this paper demonstrates that merging Transformer and Multi-Transformer layers with other models lead to more accurate volatility forecasting models. Third, the proposed stock volatility models generate appropriate risk measures in low and high volatility regimes. The Python implementation of the volatility models proposed in this paper is available in this repository.
As it is shown by the extensive literature included in this section, stock volatility forecasting has been a relevant topic not only for financial institutions and regulators but also for the academia. As financial markets can suffer drastic sudden drops, it is highly desirable to use models that can adequately forecast volatility. It is also useful to have indicators that can accurately measure risk.This paper makes use of recent deep and machine learning techniques to create more accurate stock volatility models and appropriate equity risk measures.
The rest of the paper is organized as follows: Section 2 describes the dataset, the measures used for validating the volatility forecasts and provides a look at the volatility models used as benchmark. Then, this section presents the volatility forecasting models proposed in this paper, which are based on Transformer and Multi-Transformer layers. As NLP Transformers need to be adapted in order to be used for volatility forecasting purposes and Multi-Transformer layers are introduced by this paper, explanations about the theoretical background of these structures are also given. The analysis of empirical results is presented in Section 3. Finally, the results are discussed in Section 4, followed by concluding remarks in Section 5.

2. Materials and Methods

This section is divided in five different subsections. The first one (Section 2.1) describes the data for fitting the models. The measures for validating the accuracy and value at risk (VaR) of each stock volatility model are explained in Section 2.2. Section 2.3 presents the stock volatility models and algorithms used for benchmarking purposes. Section 2.4 explains the adaptation of Transformer layers in order to be used for volatility forecasting purposes and, finally, the Multi-Transformer layers and the models based on them are presented in Section 2.5.

2.1. Data and Model Inputs

The proposed architectures and benchmark models are fitted using the rolling window approach (see Figure 1). This widely used methodology has been applied in finance, among others, by [57,58,59,60]. Rolling window uses a fixed sample length for fitting the model and, then, the following step is forecasted. As in this paper the window size is set to 650 and the forecast horizon to 1, the proposed and benchmark models are fitted using the last 650 S&P trading days and, then, the next day volatility is forecasted. This process is repeated until the whole period under analysis is forecasted. The periods used as training and testing set will be defined at the end of this subsection.
The input variables of the models proposed are the daily logarithmic returns ( r t i ) and the standard deviation of the last five daily logarithmic returns:
σ t 1 = i = 1 n ( r t i E [ r ] ) 2 n 1
As Multi-Transformer, Transformer and LSTM layers are able to manage time series, a lag of the last 10 observations of the previous variables are taken into consideration for fitting these layers. Thus, the input variables are:
X 1 = ( σ t 1 , σ t 2 , , σ t 10 )
X 2 = ( r t 1 , r t 2 , , r t 10 )
In accordance with other studies such as [38] or [50], the realized volatility is used as response variable for the models based on ANNs;
Y = σ ^ i , t = n = 0 i 1 ( r t + n E [ r f ] ) 2 i 1
where E [ r f ] = n = 0 i 1 r t + n / i and i = 5 . As shown in the previous formula, the realized volatility can be defined as the standard deviation of future logarithmic returns.
The dataset for fitting and evaluating the volatility forecasting models contains market data of S&P from 1 January 2008 to 31 December 2020. The optimum configuration of the models is obtained by applying the rolling window approach and selecting the configuration which minimizes the error (RMSE) in the period going from 1 January 2008 to 31 December 2015. The optimum configuration in combination with the rolling window methodology is applied in order to forecast the volatility contained in the testing set (from 1 January 2016 to 31 December 2020). The empirical results presented in Section 3.2 are based on the forecasts of the testing set.

2.2. Models Validation

This subsection presents the measures selected for validating and comparing the performance of the benchmark models with the algorithms proposed in this paper.
The mean absolute value ( M A E ) and the root mean squared error ( R M S E ) have been selected for validating the forecasting power of the different stock volatility models:
M A E = t = 1 N σ i , t σ ^ i , t N / R M S E = t = 1 N ( σ i , t σ ^ i , t ) 2 N
where N is the total number of observations.
The validation carried out by this study is not only interested on the accuracy, but also on the appropriateness of the risk measures generated by the different stock volatility forecasting models. In accordance with Solvency II Directive, 99.5 % VaR has been selected as risk measure. Although Solvency II has the aim of obtaining the yearly VaR, the calculations carried out in this paper will be based on a daily VaR in order to have more data points and, thus, more robust conclusions on the performance of the different models. The parametric approach developed by [61] is used for validating the different VaR estimations. The aim of this test is accepting (or rejecting) the hypothesis that the number of VaR exceedances are aligned with the confidence level selected for calculating the risk measure. In addition to the previous test, the approach suggested by [62] is also applied in order to validate the appropriateness of VaR.

2.3. Benchmark Models

This subsection introduces the benchmark models used in this paper: GARCH, EGARCH, AVGARCH, GJR-GARCH, TrARCH, FIGARCH and two architectures that combine GARCH-based algorithms with ANN and LSTM, respectively. The GARCH-based algorithms will be fitted assuming that innovations, ϵ t , follow a Student’s t-distribution. Thus, the returns generated by these models follow a conditional t-distribution [63].
The generalized autoregressive conditional heteroskedasticity (GARCH) model developed by [4] has been widely used for stock volatility forecasting purposes. GARCH(p,q) has the following expression:
σ ^ t 2 = ω + i = 1 q α i r t i 2 + i = 1 p β i σ t i 2 / r ^ t = σ ^ t ϵ t
where ω i , α i and β i are the parameters to be estimated, r t i the previous returns and σ t i 2 the last observed volatility. As previously stated, innovations ( ϵ t ) follow a Student’s t-distribution.
The absolute value GARCH [64], AVGARCH(p,q), is similar to the traditional GARCH model. In this case, the absolute value of previous return and volatility is taken into consideration to forecast volatility:
σ ^ t = ω + i = 1 q α i r t i + i = 1 p β i σ t i
As volatility behaves differently depending on the market tendency, models such as EGARCH, GJR-GARCH or TrGARCH were developed. EGARCH(p,q) [10] has the following expression for the logarithm of stocks volatility:
log σ ^ t 2 = ω + i = 1 p α i log σ ^ t i 2 + i = 1 q ( β i e t i + γ i ( e t i E ( e t i ) ) )
where ω i , α i , β i and γ i are the parameters to be estimated and e t = r t / σ t . The GJR-GARCH(p,o,q) developed by [11] has the following expression:
σ ^ t 2 = ω + i = 1 q α i r t i 2 + i = 1 o γ i r t i 2 I [ r t 1 < 0 ] + i = 1 p β i σ t i 2
As with EGARCH model, ω i , α i , β i and γ i are the parameters to be estimated. I [ r t 1 < 0 ] takes the value of 1 when the subscript condition is met. Otherwise I [ r t 1 < 0 ] = 0 . The volatility of the Threshold GARCH(p,o,q) (TrGARCH) model is obtained as follows:
σ ^ t = ω + i = 1 q α i r t i + i = 1 o γ i r t i I [ r t i < 0 ] + i = 1 p β i σ t i
As with the previous two architectures, ω i , α i , β i and γ i are the model parameters. The last GARCH-based algorithm used in this paper is the fractionally integrated GARCH (FIGARCH) model developed by [65]. The conditional variance dynamic is
σ ^ t = ω + 1 β L ϕ L ( 1 L ) d ϵ t 2 + σ h t 1
where L is the lag operator and d the fractional differencing parameter.
In addition to the previous approaches, two other hybrid models based on merging autoregressive algorithms with ANNs and LSTMs are also used as benchmark. Figure 2 shows the architecture of ANN-GARCH and LSTM-GARCH. The inputs of the algorithms are the following:
  • The last daily logarithmic return, r t 1 , for the ANN-GARCH and the last ten in the case of the LSTM-GARCH (as explained in Section 2.1).
  • The standard deviation of the last five daily logarithmic returns:
    σ t 1 = i = 1 n ( r t i E [ r ] ) 2 n 1
    where E [ r ] = i = 1 n r t i / n and n = 5 . As with the previous input variable, the last standard deviation is considered in the ANN-GARCH, whereas the last ten are taken into consideration by the LSTM-GARCH architecture.
The GARCH-based algorithms included within the ANN-GARCH and LSTM-GARCH models are the six algorithms previously presented in this same subsection (GARCH, EGARCH, AVGARCH, GJR-GARCH, TrARCH, FIGARCH).
As explained in Section 2.1, the true implied volatility, σ i , t , is used as response variable to train the models. This variable is the standard deviation of the future logarithmic returns:
σ ^ i , t = n = 0 i 1 ( r t + n E [ r f ] ) 2 i 1
where E [ r f ] = n = 0 i 1 r t + n / i . In this paper, i = 5 .
As it is shown in Figure 2, the input of the ANN-GARCH model is processed by two feed forward layers with dropout regularization. These layers have 16 and 8 neurons, respectively. The final output is produced by a feed forward layer with one neuron. In the case of the LSTM-GARCH, inputs are processed by a LSTM layer with 32 units and two feed forward layers with 8 and 1 neurons, respectively, in order to produce the final forecast.

2.4. Transformer-Based Models

Before explaining the volatility models based on Transformer layers (see Figure 3), all the modifications applied to their architecture are presented in this subsection. As previously stated, Transformer layers [54] were developed for NLP purposes. Thus, some modifications are needed in order to apply this layer for volatility forecasting purposes.
In contrast to LSTM, recurrence is not present in the architecture of Transformer layers. The two main components used by these layers in order to deal with time series are the following:
  • Positional encoder. As previously stated, Transformer layers have no recurrence structure. Thus, the information about the relative position of the observations within the time series needs to be included in the model. To do so, a positional encoding is added to the input data. In the context of NLP, Vaswani et al. [54] suggested the following wave functions as positional encoders:
    P E ( p o s , 2 i ) = sin ( p o s / 1000 2 i / d i m )
    P E ( p o s , 2 i + 1 ) = cos ( p o s / 1000 2 i / d i m )
    where d i m is the total number of explanatory variables (or word embedding dimension in NLP) used as input in the model, p o s is the position of the observation within the time series and i = ( 1 , 2 , , d i m 1 ) . This positional encoder modifies the input data depending on the lag of the time series and the embedding dimension used for the words.
    As volatility models do not use words as inputs, the positional encoder is modified in order to avoid any variation of the inputs depending on the number of time series used as input. Thus, the positional encoder suggested in this paper changes depending on the lag, but it remains the same across the different explanatory variables introduced in the model. As in the previous case, a wave function plays the role of positional encoder:
    P E p o s = cos π p o s N p o s 1 = sin π 2 + π p o s N p o s 1
    where p o s = ( 0 , 1 , , N p o s 1 ) is the position of the observation within the time series and N p o s maximum lag.
  • Multi-Head attention. It can be considered the key component of the Transformer layers proposed by [54]. As shown in Figure 3, Multi-Head attention is composed of several scaled dot-product attention units running in parallel. Scaled dot-product attention is computed as follows:
    A t t e n t i o n ( Q , K , V ) = s o f t m a x Q K T d k V
    where Q, K and V are input matrices and d k the number of input variables taken into consideration within the dot-product attention mechanism. Multi-Head attention splits the explicative variables in different groups or `heads’ in order to run the different scaled dot-product attention units in parallel. Once the different heads are calculated, the outputs are concatenated (Concat operator) and connected to a feed forward layer with linear activation. Thus, the Multi-Head attention mechanism has the following expression:
    M u l t i H e a d ( Q , K , V ) = C o n c a t h e a d 1 , , h e a d h W O
    h e a d i = A t t e n t i o n ( Q W i Q , K W i K , V W i V )
    where h is the number of heads. It is also worth mentioning that all the matrices of parameters ( W i Q , W i K , W i V and W O ) are trained using feed forward layers with linear activations.
In addition to the scaled dot-product and the Multi-Head attention mechanisms, Figure 3 shows the Transformer layers used in this paper. As suggested by [54], the Multi-Head attention is followed by a normalization, a feed forward layer with ReLU activation and, again, a normalization layer. Transformer layers also include two residual connections [66]. Thanks to these connections, the model will decide by itself if the training of some layers needs to be skipped during some phases of the fitting process.
The modified version of Transformer layers explained in the previous paragraphs are used in the volatility models presented in Figure 4. The T-GARCH architecture proposed in this paper merges the six GARCH algorithms presented in Section 2.3 with Transformer and feed forward layers in order to forecast σ ^ i , t . In addition to the previous algorithms and layers, TL-GARCH includes a LSTM with 32 units. In this last model, the temporal structure of the data is recognized and modelled by the LSTM layer and, thus, no positional encoder is needed in the Transformer layer. Both models have the following characteristics:
  • Adaptative Moment Estimator (ADAM) is the algorithm used for updating the weights of the feed forward, LSTM and Transformer layers. This algorithm takes into consideration current and previous gradients in order to implement a progressive adaptation of the initial learning rate. The values suggested by [67] for the ADAM parameters are used in this paper and the initial learning rate is set to δ = 0.01 .
  • The feed forward layers with dropout present in both models have 8 neurons, while the output layer has just one.
  • The level of dropout regularization θ [68] is optimized with the training set mentioned in Section 2.1.
  • The loss function used for weights optimization and back propagation purposes is the mean squared error.
  • Batch size is equal to 64 and the models are trained during 5000 epochs in order to obtain the final weights.

2.5. Multi-Transformer-Based Models

This subsection presents the Multi-Transformer layers and the volatility models based on them. The Multi-Transformer architecture proposed in this paper is a variant of the Transformer layers proposed by [54]. The main differences between both architectures are the following:
  • As shown in Figure 5, Multi-Transformer layers generate T different random samples of the input data. In the volatility models proposed in this paper, 90 % of the observations of the database are randomly selected in order to compute the different samples.
  • Multi-Transformer architecture is composed of T Multi-Head attention units (in this paper T = 5 ), one per each random sample of the input data. Then, the average of the different units is computed in order to obtain the final attention matrix. Thus, the Average Multi-Head (AMH) mechanism present in Multi-Transformer can be defined as follows:
    A M H ( Q , K , V ) = t = 1 T C o n c a t h e a d 1 , t , , h e a d h , t W t O T
    h e a d i , t = A t t e n t i o n ( Q t W i , t Q , K t W i , t K , V t W i , t V )
As with the Transformer architecture applied in this paper, the positional encoder used is P E p o s instead of P E ( p o s , 2 i ) and P E ( p o s , 2 i + 1 ) .
The aim of the Multi-Transformer layers introduced in the paper is to improve the stability and accuracy by applying bagging [69] to the attention mechanism. This technique is usually applied to algorithms such as linear regression, neural networks or decision trees. Instead of applying the procedure on all the data that are input into the model, the proposed methodology uses bagging only to the attention mechanism of the layer architecture.
The computational power required by bagging is one of the main limitations of this technique. As Multi-Transformer applies bagging to the attention mechanisms, their weights are trained several times in each epoch. Nevertheless, bagging is not applied to the rest of the layer weights and, thus, this offsets partially the previous limitation. It is also worth mentioning that bagging preserves the bias and this may result in underfitting.
On the other hand, this technique should bring two main advantages to the Multi-Transformer layer. First, bagging reduces significantly the error variance. Second, the aggregation of learners using this technique leads to a higher accuracy and reduces the risk of overfitting.
The structure of the volatility models based on Multi-Transformer layers (Figure 6) is similar to the architectures presented in Section 2.4. The MT-GARCH merges Multi-Transformer and feed forward layers with the six GARCH models presented in Section 2.3. In addition to the previous algorithms and layers, MTL-GARCH adds a LSTM with 32 units. The rest of the characteristics such as the optimizer, the number of neurons of the feed forward layers or the level of dropout regularization are the same than those presented in the previous section for T-GARCH and TL-GARCH.
The risk measures of ANN-GARCH, LSTM-GARCH and all the models introduced by this paper (Section 2.4 and Section 2.5) are calculated assuming that daily log-returns follow a non-standardize Student’s t-distribution with standard deviation equal to the forecasts made by the volatility models. It is worth mentioning that Student’s t-distribution generates more appropriate risk measures than normal distribution due to the shape of its tail [70,71]. In addition, this assumption is in line with the GARCH-based models used as benchmark and the inputs of the hybrid models presented in this paper.

3. Results

In this section, the forecasts and the risk measures of the volatility models presented in previous sections are compared with the ones obtained from the benchmark models. In addition, the following subsection shows the optimum hyperparameters of the benchmark and proposed hybrid volatility models.

3.1. Fitting of Models Based on Neural Networks

As explained in Section 2.1, rolling window approach ([57,58,59,60] among others) is applied for fitting the algorithms. The training set used for optimizing the level of dropout regularization contains S&P returns and observed volatilities from 1 January 2008 to 31 December 2015. Table 1 presents the error by model and level of θ .
The results of the optimization process reveals that θ = 0 generates higher error rates than the rest of the possible values regardless of the model. This means that models based on architectures such as Transformer, LSTM or feed forward layers need an appropriate level of regularization in order to avoid overfitting. According to the results, this is especially relevant for ANN-GARCH, where the error strongly depends on the level of regularization. The dropout level that minimizes the error of each model is selected.

3.2. Comparison against Benchmark Models

Once the optimum dropout level of each of the proposed volatility forecasting models based on Transformer and Multi-Transformer is selected, their performance is compared with the benchmark models (traditional GARCH processes, ANN-GARCH and LSTM-GARCH) presented in Section 2.3.
Table 2 and Table 3 present the validation error (RMSE and MAE) by year and model. The column `Total’ shows the error of the whole test period (from 1 January 2016 to 31 December 2020). The main conclusions drawn from the these tables are the following:
  • Traditional GARCH processes are outperformed by models based on merging artificial neural network architectures such as feed forward, LSTM or Transformer layers with the outcomes of autoregressive algorithms (also named hybrid models).
  • The comparison between ANN-GARCH and the rest of the volatility forecasting models based on artificial neural networks (LSTM-GARCH, T-GARCH, TL-GARCH, MT-GARCH and MTL-GARCH) reveals that feed forward layers lead to less accurate forecasts than other architectures. Multi-Transformer, Transformer and LSTM were specially created to forecast time series and, thus, the volatility models based on these layers are more accurate than ANN-GARCH.
  • Merging Multi-Transformer and Transformer layers with LSTMs leads to more accurate predictions than traditional LSTM-based architectures. Indeed, TL-GARCH achieves better results than LSTM-GARCH, even though the number of weights of TL-GARCH is significantly lower. Thus, the novel Transformer and Multi-Transformer layers introduced for NLPs purposes can be adapted as described in Section 2.4 and Section 2.5 in order to generate more accurate volatility forecasting models. It is also worth mentioning that Multi-Transformer layers, which were also introduced in this paper, lead to more accurate forecasts thanks to their ability to average several attention mechanisms. In fact, the model that achieves the lower MAE and RMSE is a mixture of Multi-Transformer and LSTM layers (MTL-GARCH).
To enhance the analysis of the results shown in Table 2 and Table 3, Figure 7 collects the RMSE and the observed volatility by year. Notice that only the most accurate GARCH-based model is shown in order to improve the visualization of the graph. The black dashed line shows that the observed volatility of 2020 was significantly higher than the rest of the years due to the turmoil caused by COVID-19 outbreak. As expected, the error of every model is also higher in 2020 because the market volatility was more unpredictable than the rest of the years. Nevertheless, it has to be mentioned that the 2020 forecasts of traditional autoregressive algorithms are significantly less accurate than hybrid models based on architectures such as LSTM, Transformer or Multi-Transformer layers.
Although the observed volatility is lower in years before 2020, autoregressive models are also outperformed by hybrid models. Nevertheless, the difference between both sets of models is remarkably lower.
The p-values of the Kupiec and Christoffersen tests by volatility model and year are shown in Table 4 and Table 5, respectively. In contrast to the approach suggested by Kupiec, Christoffersen test is not only focused on the total number of exceedances, but it also takes into consideration the number of consecutive VaR exceedances. As stated in Section 2.2, the risk measure and confidence level ( 99.5 % VaR) selected are in line with Solvency II Directive. This regulation sets the principles for calculating the capital requirements and assessing the risk profile of the insurance companies based in the European Union. This law covers not only the underwriting risks but also financial risks such as the potential losses due to variations on the interest rate curves or the equity prices.
The column ‘Total’ of Table 4 and Table 5 reveal that only TL-GARCH, MT-GARCH and MTL-GARCH produce appropriate risk measures (p-value higher than 0.05 in both tests) for the period 2016–2020. The rest of the models fail both tests and, thus, their risk measures can not be considered to be appropriate for that period.
As with any other statistical test, the higher the number of data points the more relevant are the outcomes obtained from the test. That is the reason why the previous paragraph focuses on the ‘Total’ column and not on the specific results obtained by year. The results by year show that most of the models fail the test in 2020 due to the high level of volatility produced by COVID-19 pandemic.
According to these results, the stock volatility models introduced in this paper (T-GARCH, TL-GARCH, MT-GARCH and MTL-GARCH) produce more accurate estimations and appropriate risk measures in most of the cases. Regarding the models accuracy, it is specially remarkable the difference observed in 2020, where COVID-19 caused a significant turmoil in the stock market. Concerning the appropriateness of equity risk measures, three out of four models based on Transformer and Multi-Transformer pass Kupiec and Christofferesen test for the period 2016–2020, while all the benchmark models fail at least one of them. Notice that the proposed models are compared with other approaches belonging to its own family (ANN-GARCH and LSTM-GARCH) and autoregressive models belonging to the GARCH family.

4. Discussion

This paper introduced a set of volatility forecasting models based on Transformer and Multi-Transformer layers. As Transformer layers were developed for NLP purposes [54], their architecture is adapted in order to generate stock volatility forecasting models. Multi-Transformer layers, which are introduced by this paper, have the aim of improving the stability and accuracy of Transformer layers by applying bagging to the attention mechanism. The predictive power and risk measures generated by the proposed volatility forecasting models (T-GARCH, TL-GARCH, MT-GARCH and MTL-GARCH) are compared with traditional GARCH processes and other hybrid models based on LSTM and feed forward layers.
Three main outcomes were drawn from the empirical results. First, hybrid models based on LSTM, Transformer or Multi-Transformer layers outperform traditional autoregressive algorithms and hybrid models based on feed forward layers. The validation error by year shows that this difference is more relevant in 2020, when the volatility of S&P500 was significantly higher than in the previous years due to COVID-19 pandemic. Volatility forecasting models are mainly used for pricing derivatives and assessing the risk profile of financial institutions. As the more relevant shocks on the solvency position of financial institutions and derivatives prices are observed in high volatility regimes, the accurateness of these models is particularly important in years such as 2020.
The higher performance of hybrid models have also been demonstrated by [38,39,40,41,42,43,44]. These papers merged traditional GARCH models with feed forward layers to predict stock market volatility. This type of models have shown also a superior performance in other financial fields such as oil market volatility [48,49] and metals price volatility [46,47]. Notice that this paper does not only present a comparison with traditional autoregressive models, but it also shows that Transformer and Multi-Transformer can lead to more accurate volatility estimations than other hybrid models.
Second, Multi-Transformer layers lead to more accurate volatility forecasting models than Transformer layers. As expected, applying bagging to the attention mechanism has a positive impact on the performance of the models presented in this paper. It is also remarkable that empirical results demonstrate that merging LSTM with Transformer or Multi-Transformer layers has also a positive impact on the models performance. On one hand, the volatility forecasting model based on Multi-Transformer and LSTM (named MTL-GARCH) achieves the best results in the period 2016–2020. On the other hand, the merging of Transfomer with LSTM (TL-GARCH) leads to a lower error rate than the hybrid model based only on LSTM layers (LSTM-GARCH) even though the number of weights of the first model is significantly lower. Thus, the use of Transfomer layers can lead to simpler and more accurate volatility forecasting models. Notice that Transformer layers are already considered the state of art thanks to BERT [55] and GPT-3 [56]. These models have been successfully used for sentence prediction, conversational response generation, sentiment classification, coding and writing fiction, among others.
Third, the results of Kupiec and Christoffersen tests revealed that only the risk estimations made by MTL-GARCH, TL-GARCH and MT-GARCH can be considered as appropriate for the period 2016–2020, whereas traditional autoregressive algorithms and hybrid models based on feed forward and LSTM layers failed, at least, one of the tests. As previously stated, volatility does not play only a key role in risk management but also in derivative valuation models. Thus, using a volatility model that generates appropriate risk measures can lead to more accurate derivatives valuation.

5. Conclusions

Transformer layers are the state of the art in natural language processing. Indeed, the performance of this layer have overcome the performance of any other previous model in this field [56]. As Transformer layers were specially created for natural language processing, they need to be modified in order to be used for other purposes. Probably, this is one of the main reasons why this layer have not been already extended to other fields. This paper provides the modifications needed to apply this layer for stock volatility forecasting purposes. The results shown in this paper demonstrates that Transformer layers can overcome also the performance of the main stock volatility models.
Following the intuition of bagging [69], this paper introduces Multi-Transformer layers. This novel architecture has the aim of improving the stability and accuracy of the attention mechanism, which is the core of Transformer layers. According to the results, it can be concluded that this procedure improves the accuracy of stock volatility models based on Transformer layers.
Leaving aside the comparisons between Transformer and Multi-Transformer layers, the hybrid models based on them have overcome the performance of autoregressive algorithms and other models based on feed forward layers and LSTMs. The architecture of these hybrid models (T-GARCH, TL-GARCH, MT-GARCH and MTL-GARCH) based on Transformer and Multi-Transformer layers is also provided in this paper.
According to the results, it is also worth noticing that the risk estimations based on the previous models are specially appropriate. The VaR of most of these models can be considered accurate even in years such as 2020, when the COVID-19 pandemic caused a remarkable turmoil in the stock market.
Consequently, the empirical results obtained with the hybrid models based on Transfomer and Multi-Transformer layers suggest that further investigation should be conducted about the possible application of them for derivative valuation purposes. Notice that volatility plays a key role in the financial derivatives valuation. In addition, the models can be extended by merging Transformer or Multi-Transformer layers with other algorithms (such as gradient boosting with trees or random forest) or modifying some key assumptions of the attention mechanism.

Author Contributions

Conceptualization, E.R.-P.; methodology, E.R.-P., P.J.A.-G. and J.J.N.-V.; software, E.R.-P.; validation, P.J.A.-G. and J.J.N.-V.; formal analysis, E.R.-P.; investigation, E.R.-P., P.J.A.-G. and J.J.N.-V.; writing—both original draft preparation, review and editing, E.R.-P., P.J.A.-G. and J.J.N.-V.; supervision, P.J.A.-G. and J.J.N.-V.; project administration, P.J.A.-G. and J.J.N.-V.; funding acquisition, P.J.A.-G. and J.J.N.-V. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by Economics Department of Universidad de Alcalá.

Data Availability Statement

The Python implementation of the volatility models proposed in this paper is available in https://github.com/EduardoRamosP/MultiTransformer (accessed on 26 June 2021).

Conflicts of Interest

The authors declare that they have no conflict of interest regarding the publication of the research article.

References

  1. Hull, J. Risk Management and Financial Institutions, 4th ed.; Wiley and Sons: London, UK, 2015. [Google Scholar]
  2. Rajashree, P.; Ranjeeeta, B. A differential harmony search based hybrid internal type2 fuzzy EGARCH model for stock market volatility prediction. Int. J. Approx. Reason. 2015, 59, 81–104. [Google Scholar]
  3. Engle, R. Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation. Econometrica 1982, 50, 987–1007. [Google Scholar] [CrossRef]
  4. Bollerslev, T. Generalized autoregressive conditional heteroskedasticity. J. Econom. 1986, 31, 307–327. [Google Scholar] [CrossRef] [Green Version]
  5. Mandelbrot, B. The variation of certain speculative prices. J. Bus. 1963, 36, 394–419. [Google Scholar] [CrossRef]
  6. Engle, R.; Lee, G. A permanent and transitory component model of stock return volatility. In Cointegration, Causality, and Forecasting: A Festschrift in Honor of Clive W. J. Granger; Engle, R., White, H., Eds.; Oxford University Press: Oxford, UK, 1999; pp. 475–497. [Google Scholar]
  7. Haas, M.; Mittnik, S.; Paolella, M. Mixed normal conditional heteroskedasticity. J. Financ. Econom. 2004, 2, 211–250. [Google Scholar] [CrossRef]
  8. Haas, M.; Mittnik, S.; Paolella, M. A new approach to Markov-switching GARCH models. J. Financ. Econom. 2004, 2, 493–530. [Google Scholar] [CrossRef]
  9. Haas, M.; Paolella, M. Mixture and regime-switching GARCH models. In Handbook of Volatility Models and Their Applications; Bauwens, L., Hafner, C., Laurent, S., Eds.; John Wiley and Sons: Hoboken, NJ, USA, 2012; pp. 71–102. [Google Scholar]
  10. Nelson, D.B. Conditional Heteroskedasticity in Asset Returns: A New Approach. Econometrica 1991, 59, 347–370. [Google Scholar] [CrossRef]
  11. Glosten, L.; Jagannathan, R.; Runkle, D. 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]
  12. Kraft, D.; Engle, R. Autoregressive Conditional Heteroskedasticity in Multiple Time Series; Department of Economics, UCSD: San Diego, CA, USA, 1982. [Google Scholar]
  13. Engle, R.; Granger, C.; Kraft, D. Combining competing forecasts of inflation with a bivariate ARCH model. J. Econ. Dyn. Control. 1984, 8, 151–165. [Google Scholar] [CrossRef]
  14. Bollerslev, T.; Engle, R.; Wooldridge, J. A Capital Asset Pricing Model with time-varying covariances. J. Political Econ. 1988, 96, 116–131. [Google Scholar] [CrossRef]
  15. Tse, Y.; Tsui, K. A multivariate GARCH model with time-varying correlations. J. Bus. Econ. Stat. 2002, 20, 351–362. [Google Scholar] [CrossRef]
  16. Engle, R. Dynamic conditional correlation: A simple class of multivariate generalized autoregressive conditional heteroskedasticity models. J. Bus. Econ. Stat. 2002, 20, 339–350. [Google Scholar] [CrossRef]
  17. Engle, R.; Kroner, F. Multivariate simultaneous generalized ARCH. Econom. Theory 1995, 11, 122–150. [Google Scholar] [CrossRef]
  18. Engle, R.; Ng, V.; Rotschild, M. Asset pricing with a factor-ARCH covariance structure: Empirical estimates for Treasury Bills. J. Econom. 1990, 45, 213–238. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, L.; Zhu, K.; Ling, S. The ZD-GARCH model: A new way to study heteroscedasticity. J. Econom. 2018, 202, 1–17. [Google Scholar] [CrossRef]
  20. Heston, S.L. A closed-form solution for options with stochastic volatility with applications to bond and currency options. Rev. Financ. Stud. 1993, 6, 327–343. [Google Scholar] [CrossRef] [Green Version]
  21. Cox, J.; Ingersoll, J.; Ross, S. A Theory of the Term Structure of Interest Rates. Econometrica 1985, 53, 385–407. [Google Scholar] [CrossRef]
  22. Melino, A.; Turnbull, S. Pricing foreign currency options with stochastic volatility. J. Econom. 1990, 45, 239–265. [Google Scholar] [CrossRef]
  23. Andersen, T.; Sorensen, B. GMM estimation of a stochastic volatility model: A Monte Carlo study. J. Bus. Econ. Stat. 1999, 14, 329–352. [Google Scholar]
  24. Durbin, J.; Koopman, S. Monte Carlo maximum likelihood estimation for non-Gaussian state space models. Biometrika 1997, 84, 669–684. [Google Scholar] [CrossRef] [Green Version]
  25. Broto, C.; Ruiz, E. Estimation methods for stochastic volatility models: A survey. J. Econ. Surv. 2004, 18, 613–649. [Google Scholar] [CrossRef] [Green Version]
  26. Danielsson, J. Stochastic volatility in asset prices: Estimation by simulated maximum likelihood. J. Econom. 2004, 64, 375–400. [Google Scholar] [CrossRef]
  27. Andersen, T. Encyclopedia of Complexity and System Sciences; Chapter Stochastic Volatility; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  28. Hull, J.C.; White, A. The Pricing of Options on Assets with Stochastic Volatilities. J. Financ. 1987, 42, 281–300. [Google Scholar] [CrossRef]
  29. Hagan, P.; Kumar, D.; Lesniewski, A.; Woodward, D. Managing Smile Risk. Wilmott Mag. 2002, 1, 84–108. [Google Scholar]
  30. Mcculloch, W.; Pitts, W. A Logical Calculus of Ideas Immanent in Nervous Activity. Bull. Math. Biophys. 1943, 5, 127–147. [Google Scholar] [CrossRef]
  31. Friedman, J.H. Greedy Function Approximation: A Gradient Boosting Machine. Ann. Stat. 2000, 29, 1189–1232. [Google Scholar]
  32. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  33. Cortes, C.; Vapnik, V. Support-Vector Networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  34. Gestel, T.; Suykens, J.; Baestens, D.; Lambrechts, A.; Laneknet, G. Financial time series prediction using least squares Support Vector Machines within the evidence framework. IEEE Trans. Neural Netw. 2001, 12, 809–821. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Gupta, A.; Dhinga, B. Stock markets prediction using hidden Markov models. In Proceedings of the 2012 Students Conference on Engineering and Systems, Allahabad, India, 16–18 March 2012; pp. 1–4. [Google Scholar]
  36. Dias, F.; Nogueira, R.; Peixoto, G.; Moreira, W. Decision-making for financial trading: A fusion approach of Machine Learning and Portfolio Selection. Expert Syst. Appl. 2019, 115, 635–655. [Google Scholar]
  37. Hamid, S.; Iqbid, Z. Using neural networks for forecasting volatility of S&P 500 Index futures prices. J. Bus. Res. 2002, 57, 1116–1125. [Google Scholar]
  38. Roh, T. Forecasting the Volatility of Stock Price Index. Expert Syst. Appl. 2006, 33, 916–922. [Google Scholar]
  39. Hajizadeh, E.; Seifi, A.; Zarandi, F.; Turksen, I. A hybrid modeling approach for forecasting the volatility of S&P 500 Index return. Expert Syst. Appl. 2012, 39, 531–536. [Google Scholar]
  40. Kristjanpoller, W.; Fadic, A.; Minutolo, M. Volatility forecast using hybrid neural network models. Expert Syst. Appl. 2014, 41, 2437–2442. [Google Scholar] [CrossRef]
  41. Monfared, S.A.; Enke, D. Volatility Forecasting Using a Hybrid GJR-GARCH Neural Network Model. Procedia Comput. Sci. 2014, 36, 246–253. [Google Scholar] [CrossRef] [Green Version]
  42. 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]
  43. Kim, H.; Won, C. Forecasting the volatility of stock price index: A hybrid model integrating LSTM with multiple GARCH-type models. Expert Syst. Appl. 2018, 103, 25–37. [Google Scholar] [CrossRef]
  44. Back, Y.; Kim, H. ModAugNet: A new forecasting framework for stock market index value with an overfitting prevention LSTM module and a prediction LSTM module. Expert Syst. Appl. 2018, 113, 457–480. [Google Scholar] [CrossRef]
  45. Bildirici, M.; Ersin, O. Improving forecasts of GARCH family models with the artificial neural networks: An applicaiton to the daily returns in Istanbul Stock Exchange. Expert Syst. Appl. 2009, 36, 7355–7362. [Google Scholar] [CrossRef]
  46. Kristjanpoller, W.; Minutolo, M. Gold price volatility: A Forecasting approach using the Artificial Neural Network-GARCH model. Expert Syst. Appl. 2015, 42, 7245–7251. [Google Scholar] [CrossRef]
  47. Kristjanpoller, W.; Hernández, E. Volatility of main metals forecasted by a hybrid ANN-GARCH model with regressors. Expert Syst. Appl. 2017, 84, 290–300. [Google Scholar] [CrossRef]
  48. Kristjanpoller, W.; Minutolo, M. Forecasting volatility of oil price using an Artificial Neural Network-GARCH model. Expert Syst. Appl. 2016, 65, 233–241. [Google Scholar] [CrossRef]
  49. Verma, S. Forecasting volatility of crude oil futures using a GARCH–RNN hybrid approach. Intell. Syst. Accounting, Financ. Manag. 2021. [Google Scholar] [CrossRef]
  50. Ramos-Pérez, E.; Alonso-González, P.; Núñez-Velázquez, J. Forecasting volatility with a stacked model based on a hybridized Artificial Neural Network. Expert Syst. Appl. 2019, 129, 1–9. [Google Scholar] [CrossRef]
  51. Vidal, A.; Kristjanpoller, W. Gold volatility prediction using a CNN-LSTM approach. Expert Syst. Appl. 2020, 157. [Google Scholar] [CrossRef]
  52. Jung, G.; Choi, S.Y. Forecasting Foreign Exchange Volatility Using Deep Learning Autoencoder-LSTM Techniques. Complexity 2021, 2021, 1–16. [Google Scholar] [CrossRef]
  53. Peng, Y.; Melo, P.; Camboim de Sá, J.; Akaishi, A.; Montenegro, M. The best of two worlds: Forecasting high frequency volatility for cryptocurrencies and traditional currencies with Support Vector Regression. Expert Syst. Appl. 2018, 97, 177–192. [Google Scholar] [CrossRef]
  54. Vaswani, A.; Shazeer, N.; Parmar, N.; Uszkoreit, J.; Jones, L.; Gomez, A.N.; Kaiser, L.; Polosukhin, I. Attention Is All You Need. Adv. Neural Inf. Process. Syst. 2017, 2017, 5998–6008. [Google Scholar]
  55. Devlin, J.; Chang, M.; Lee, K.; Toutanova, K. BERT: Pre-training of Deep Bidirectional Transformers for Language Understanding. arXiv 2018, arXiv:1810.04805. [Google Scholar]
  56. Brown, T.B.; Mann, B.; Ryder, N.; Subbiah, M.; Kaplan, J.; Dhariwal, P.; Neelakantan, A.; Shyam, P.; Sastry, G.; Askell, A.; et al. Language Models are Few-Shot Learners. In Advances in Neural Information Processing Systems; Larochelle, H., Ranzato, M., Hadsell, R., Balcan, M.F., Lin, H., Eds.; Curran Associates, Inc.: New York, NY, USA, 2020; Volume 33, pp. 1877–1901. [Google Scholar]
  57. Swanson, N.R. Money and output viewed through a rolling window. J. Monet. Econ. 1998, 41, 455–474. [Google Scholar] [CrossRef]
  58. Goyal, A.; Welch, I. Predicting the Equity Premium With Dividend Ratios; NBER Working Papers 8788; National Bureau of Economic Research, Inc.: Cambridge, MA, USA, 2002. [Google Scholar]
  59. Zivot, E.; Wang, J. Modeling Financial Time Series with S-PLUS®; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  60. Molodtsova, T.; Papell, D. Taylor Rule Exchange Rate Forecasting during the Financial Crisis. NBER Int. Semin. Macroecon. 2012, 9, 55–97. [Google Scholar] [CrossRef] [Green Version]
  61. Kupiec, P.H. Techniques for Verifying the Accuracy of Risk Measurement Models. J. Deriv. 1995, 3, 73–84. [Google Scholar] [CrossRef]
  62. Christoffersen, P.F.; Bera, A.; Berkowitz, J.; Bollerslev, T.; Diebold, F.; Giorgianni, L.; Hahn, J.; Lopez, J.; Mariano, R. Evaluating Interval Forecasts. Int. Econ. Rev. 1997, 39, 841–862. [Google Scholar] [CrossRef]
  63. Bauwens, L.; Hafner, C.; Laurent, S. Handbook of Volatility Models and Their Applications; Wiley Handbooks in Financial E; Wiley: Hoboken, NJ, USA, 2012. [Google Scholar]
  64. Taylor, S.J. Modelling Financial Time Series; Wiley: Hoboken, NJ, USA, 1986. [Google Scholar]
  65. Baillie, R.T.; Bollerslev, T.; Mikkelsen, H.O. Fractionally integrated generalized autoregressive conditional heteroskedasticity. J. Econom. 1996, 74, 3–30. [Google Scholar] [CrossRef]
  66. He, K.; Zhang, X.; Ren, S.; Sun, J. Deep Residual Learning for Image Recognition. In Proceedings of the 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Las Vegas, NV, USA, 27–30 June 2016; pp. 770–778. [Google Scholar] [CrossRef] [Green Version]
  67. Kingma, D.P.; Ba, J. Adam: A Method for Stochastic Optimization. arXiv 2014, arXiv:1412.6980. [Google Scholar]
  68. Srivastava, N.; Hinton, G.; Krizhevsky, A.; Sutskever, I.; Salakhutdinov, R. Dropout: A Simple Way to Prevent Neural Networks from Overfitting. J. Mach. Learn. Res. 2014, 15, 1929–1958. [Google Scholar]
  69. Breiman, L. Bagging Predictors. Mach. Learn. 1996, 24, 123–140. [Google Scholar] [CrossRef] [Green Version]
  70. Jeremic, Z.; Terzić, I. Empirical estimation and comparison of Normal and Student T linear VaR on the Belgrade Stock Exchange. In Proceedings of the Sinteza 2014—Impact of the Internet on Business Activities in Serbia and Worldwide, Belgrade, Serbia, 25–26 April 2014; pp. 298–302. [Google Scholar] [CrossRef] [Green Version]
  71. McNeil, A.J.; Frey, R.; Embrechts, P. Quantitative Risk Management: Concepts, Techniques and Tools; Princeton University Press: Princeton, NJ, USA, 2015. [Google Scholar]
Figure 1. Rolling window methodology.
Figure 1. Rolling window methodology.
Mathematics 09 01794 g001
Figure 2. ANN-GARCH and LSTM-GARCH architectures.
Figure 2. ANN-GARCH and LSTM-GARCH architectures.
Mathematics 09 01794 g002
Figure 3. Transformer and Multi-Head attention mechanism.
Figure 3. Transformer and Multi-Head attention mechanism.
Mathematics 09 01794 g003
Figure 4. T-GARCH and TL-GARCH volatility models.
Figure 4. T-GARCH and TL-GARCH volatility models.
Mathematics 09 01794 g004
Figure 5. Multi-Transformer architecture.
Figure 5. Multi-Transformer architecture.
Mathematics 09 01794 g005
Figure 6. MT-GARCH and MTL-GARCH volatility models.
Figure 6. MT-GARCH and MTL-GARCH volatility models.
Mathematics 09 01794 g006
Figure 7. Observed volatility and RMSE by year.
Figure 7. Observed volatility and RMSE by year.
Mathematics 09 01794 g007
Table 1. RMSE by level of θ .
Table 1. RMSE by level of θ .
Model θ = 0 θ = 0.05 θ = 0.10 θ = 0.15
ANN-GARCH0.03510.00920.00850.0082
LSTM-GARCH0.00650.00570.00560.0054
T-GARCH0.00890.00760.00720.0074
TL-GARCH0.00500.00450.00440.0045
MT-GARCH0.00680.00620.00640.0064
MTL-GARCH0.00470.00450.00420.0044
Source: own elaboration.
Table 2. RMSE by volatility model and year.
Table 2. RMSE by volatility model and year.
Model20162017201820192020Total
GARCH(1,1)0.00580.00260.00950.00730.10260.0464
AVGARCH(1,1)0.00530.00270.00760.00560.08470.0383
EGARCH(1,1)0.00560.00280.00930.00780.08800.0399
GJR-GARCH(1,1,1)0.00900.00280.01260.00680.12480.0565
TrGARCH(1,1,1)0.00740.00270.01150.00580.11530.0521
FIGARCH(1,1)0.00620.00290.00950.00660.10110.0457
ANN-GARCH0.00420.00230.00600.00440.01710.0086
LSTM-GARCH0.00320.00210.00430.00300.01010.0054
T-GARCH0.00480.00290.00580.00440.01170.0067
TL-GARCH0.00300.00190.00330.00260.00700.0040
MT-GARCH0.00360.00210.00460.00330.00960.0054
MTL-GARCH0.00300.00160.00330.00260.00660.0038
Source: own elaboration.
Table 3. MAE by volatility model and year.
Table 3. MAE by volatility model and year.
Model20162017201820192020Total
GARCH(1,1)0.00370.00190.00580.00440.03630.0105
AVGARCH(1,1)0.00340.00190.00490.00370.02960.0087
EGARCH(1,1)0.00350.00200.00600.00480.03330.0100
GJR-GARCH(1,1,1)0.00480.00200.00740.00420.04040.0118
TrGARCH(1,1,1)0.00420.00200.00690.00380.03650.0107
FIGARCH(1,1)0.00380.00210.00550.00410.03610.0104
ANN-GARCH0.00290.00190.00380.00290.00950.0042
LSTM-GARCH0.00220.00150.00270.00210.00600.0029
T-GARCH0.00350.00210.00410.00310.00700.0040
TL-GARCH0.00200.00140.00210.00180.00440.0023
MT-GARCH0.00240.00160.00310.00230.00570.0030
MTL-GARCH0.00190.00120.00210.00180.00410.0022
Source: own elaboration.
Table 4. Kupiec test (p-values) by volatility model and year.
Table 4. Kupiec test (p-values) by volatility model and year.
Model20162017201820192020Total
GARCH(1,1)0.5430.5400.0510.5430.0520.008
AVGARCH(1,1)0.5430.5400.0510.5430.0520.008
EGARCH(1,1)0.5430.5400.0510.5430.0520.008
GJR-GARCH(1,1,1)0.5430.5400.0110.5430.1900.008
TrGARCH(1,1,1)0.5430.5400.0510.8100.1900.042
FIGARCH(1,1)0.5430.5400.0510.5430.0520.008
ANN-GARCH0.5430.5400.0010.0020.0120.001
LSTM-GARCH0.8100.1860.5400.1880.1900.042
T-GARCH0.1880.5400.0020.5430.0520.001
TL-GARCH0.5430.5400.8130.8100.8100.782
MT-GARCH0.1120.5400.5400.1880.0520.089
MTL-GARCH0.5430.1130.1130.8100.1900.910
Source: own elaboration.
Table 5. Christoffersen test (p-values) by volatility model and year.
Table 5. Christoffersen test (p-values) by volatility model and year.
Model20162017201820192020Total
GARCH(1,1)0.5220.5200.0040.5230.0480.002
AVGARCH(1,1)0.5220.5200.0040.5230.0480.002
EGARCH(1,1)0.5220.5200.0040.5230.0480.002
GJR-GARCH(1,1,1)0.5220.5200.0020.5230.1790.002
TrGARCH(1,1,1)0.5220.5200.0040.8000.1790.009
FIGARCH(1,1)0.5220.5200.0040.5230.0480.002
ANN-GARCH0.5220.5200.0010.0020.0020.001
LSTM-GARCH0.8000.1800.5200.1770.1790.037
T-GARCH0.1760.5200.0010.5230.0480.001
TL-GARCH0.5220.5200.8030.8000.7970.693
MT-GARCH0.1130.5200.5200.1770.0480.079
MTL-GARCH0.5220.1130.1130.8000.1790.790
Source: own elaboration.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ramos-Pérez, E.; Alonso-González, P.J.; Núñez-Velázquez, J.J. Multi-Transformer: A New Neural Network-Based Architecture for Forecasting S&P Volatility. Mathematics 2021, 9, 1794. https://doi.org/10.3390/math9151794

AMA Style

Ramos-Pérez E, Alonso-González PJ, Núñez-Velázquez JJ. Multi-Transformer: A New Neural Network-Based Architecture for Forecasting S&P Volatility. Mathematics. 2021; 9(15):1794. https://doi.org/10.3390/math9151794

Chicago/Turabian Style

Ramos-Pérez, Eduardo, Pablo J. Alonso-González, and José Javier Núñez-Velázquez. 2021. "Multi-Transformer: A New Neural Network-Based Architecture for Forecasting S&P Volatility" Mathematics 9, no. 15: 1794. https://doi.org/10.3390/math9151794

APA Style

Ramos-Pérez, E., Alonso-González, P. J., & Núñez-Velázquez, J. J. (2021). Multi-Transformer: A New Neural Network-Based Architecture for Forecasting S&P Volatility. Mathematics, 9(15), 1794. https://doi.org/10.3390/math9151794

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