# Model Calibration and Validation for the Fuzzy-EGARCH-ANN Model

Authors: Geleta T. Mohammed; Jane A. Aduda; Ananda O. Kube
Journal: Applied Computational Intelligence and Soft Computing
 (2021)

---

## Abstract

This work shown as the fuzzy-EGARCH-ANN (fuzzy-exponential generalized autoregressive conditional heteroscedastic-artificial neural network) model does not require continuous model calibration if the corresponding DE algorithm is used appropriately, but other models such as GARCH, EGARCH, and EGARCH-ANN need continuous model calibration and validation so they fit the data and reality very well up to the desired accuracy. Also, a robust analysis of volatility forecasting of the daily S&P 500 data collected from Yahoo Finance for the daily spanning period 1/3/2006 to 20/2/2020. To our knowledge, this is the first study that focuses on the daily S&P 500 data using high-frequency data and the fuzzy-EGARCH-ANN econometric model. Finally, the research finds that the best performing model in terms of one-step-ahead forecasts based on realized volatility computed from the underlying daily data series is the fuzzy-EGARCH-ANN (1,1,2,1) model with Student’st-distribution.


---

## Body

## 1. Introduction

Mathematical models can provide a cost-effective, objective, and flexible approach for assessing management decisions, mainly when these decisions are strategic alternatives. In some instances, a mathematical model is the only means available for evaluating and testing alternatives. However, for this potential to be realized, models must be valid for the application and must provide credible and reliable results. The process of ensuring validity, credibility, and reliability typically consists of model validation and calibration.Model calibration is the process of determining to what extent the model user can or is required to modify the default input parameter values that describe the underlying mechanics to obtain the better values of the problem. Model users can often not determine the impact that input parameter values have on the selected model. This may arise from several causes, including a lack of understanding of principles, a lack of understanding of the model, poor model documentation, or a combination of these [1]. Moreover, given observed data, adjust model parameters until the model reproduces observation “closely enough.” In fact, perfect accuracy is not possible since some errors cannot be reduced in the model, e.g., measurement errors inherent in well test results (see Hill and Tiedeman for more on this topic [2]).Model validation is the process of determining the degree to which a theory, an approach, or a model is a “good enough” representation of the reality from the viewpoint of the intended uses of the theory, the approach, or the model [3–5]. It is concerned with quantifying the accuracy of the model by comparing model data to the observed data [6]. It is a complex process because the concept of “good enough” includes subjective judgments of what constitutes a reasonable degree of “good enough,” and it differs from the point of view of different individuals. Subjective judgments prevent to make a general validation approach of theories, approaches, or models. Therefore, absolute validation is philosophically not possible because it requires an infinite number of tests [3].A theory, a methodology, an approach, or a mathematical model can never be proven to be valid; instead, we can say that there is not enough evidence to reject them. Therefore, as long as a methodology or an approach does not have sufficient evidence to reject it, we can accept it. Validation relative to a specific series of tests may be legitimate as a basis for making decisions [7]. In these situations, relative validation can be possible when validation involves a comparison of observed events with those predicted by models [8, 9].Finally, the main objective of this paper is to deal with model calibration and validation of the fuzzy-EGARCH-ANN model. A lot of scholars tried to develop a hybrid model which fits the given data in reality very well. Even if most of them are better than the base models, they still need continuous model calibration and validation. For these cases, one can look at the references listed from [10–13]. Therefore, we will see whether the fuzzy-EGARCH-ANN model needs continuous model calibration and validation or not.This paper is organized as follows. Introduction is given in Section1. Research methodology is presented in Section 2. In Section 3, details of results and discussion are given. In Section 4, a conclusion part of this study is presented.


## 2. Research Methodology

### 2.1. Data Source

To achieve the objectives of the study, daily S&P 500 data were collected from Yahoo Finance. As explained in the literature, the study split the given dataset into an in-sample and out-sample period. The in-sample dataset was used for initial parameter estimation and model selection. The out-sample dataset was used to evaluate forecast performance. In this study, the in-sample period ranges from 1/3/2006 to 3/5/2019 (800 observations), whereas the out-of-sample period ranges from 4/5/2019 to 20/2/2020 (201 observations).


### 2.2. Model Specification and Econometric Tests

#### 2.2.1. Unit Root Test of the Variables

For determining whether a time series is stationary or nonstationary, there are many tests. This study used the augmented Dickey–Fuller (ADF) test [14] and the Phillips–Perron (PP) test [15]. Ljung–Box statistics Qm and Lagrange multiplier (LM) tests are the appropriate ARCH effect tests and were used in this study. The Ljung–Box test statistic Q was used to assess the independency among the residuals [16]. When dealing with ARCH/GARCH−type models, we first examine the characteristics of the unconditional distribution of the stock price. This will enable us to explore and explain some stylized facts that exist in the financial time series. In statistics, the Bera and Jarque test is a test of departure from normality based on the sample skewness and kurtosis [17].


#### 2.2.2. Fuzzy-EGARCH-ANN Model

As proposed by the authors, the fuzzy-EGARCH-ANN model is described by a collection of fuzzy rules in the form of If-Then statements in order to describe the stock fluctuations with volatility clustering overlooked by EGARCH-ANN model via fuzzy rules and the asymmetric responses of volatility to positive and negative shocks via an EGARCH-ANN model [18].Thekth rule of the fuzzy system for EGARCH-ANN p,q,m,s model is described by the following.Rulek: if x1t is Fk1 and …xdt is Fkd, then(1)ykt=σktεkt,lnσkt2=Ek+∑j=1qβkjlnσt−j2+∑i=1pγkiyt−iσt−i−Eyt−iσt−i+∑i=1pαkiyt−iσt−i+∑h=1sξkhGZtWh+wh,0,where(2)GZtWh+wh,0=tanhZtWh+wh,0,Wh=wh,1,…wh,m′,zt−d=yt−d−EyEy2,d=1:m,xt=yt−1,yt−2,…yt−p,σt−12,σt−22,…σt−q2T=x1t,x2t,…,xdt′.As given by the authors,xt is the input vector with d=p+q at instance t. Fkl is the fuzzy set to describe the stock market return and volatility for l=1,2,3,…,d, R is the number of rules, and xlt is the premise variable. The authors assumed that the residual series distribution follows either the Gaussian normal or Student’s t-distribution. If it is not following the Gaussian normal distribution, it follows Student’s t-distribution directly.The authors had described the output of this fuzzy-EGARCH-ANN model as the weighted average of each individual rule, and it is obtained by using FIS fundamental steps as in [18, 19] as follows.Step 1 (fuzzification layer).
using the Gaussian membership function, find the grade of membership of the inputxlt in Fkl as follows:(3)Fklxlt=exp−12xlt−cklakl2,where ckl is the center and akl is the spread of kth rule membership function corresponding to the lth premise variable.Step 2 (firing strength layer).
find the firing strength ofkth rule by assuming the product T-norm of the antecedent fuzzy sets as(4)ukxt=∏l=1dFklxlt=∏l=1d−12xlt−cklakl2.Step 3 (normalization layer).
find the normalized firing strengths, i.e., the ratio of thekth rule’s firing strength to the sum of all rule’s firing strengths:(5)w¯k=ukxt∑k=1Rukxt.Step 4 (consequent and defuzzification layer).
combine the normalized firing strengths and the consequent corresponding rule to produce the model output as the weighted average of each individual rule:(6)lnσt2=∑k=1Rw¯kfk,where fk is the output of the kth rule. That is,(7)fk=Ek+∑j=1qβkjlnσt−j2+∑i=1pγkiyt−iσt−i−Eyt−iσt−i+∑i=1pαkiyt−iσt−i+∑h=1sξkhGZtWh+wh,0.
This implies that(8)lnσt2=∑k=1Rw¯kfk=∑k=1Rukxt∑k=1Rukxtfk.
The collection of theR rules assembles a model as a combination of local models. As a result, the exponential of this weighted average value gives us the predicted stock market return volatility using the fuzzy-EGARCH-ANN model [18].


### 2.3. Choosing the Optimal Lag Length and Model Selection Criterion

The optimal lag length for ARIMA and EGARCH-ANN models was chosen using the ACF, PACF plot, and information criteria. Some of the most popular information criteria are AIC and BIC, which are used to select the appropriate models. The PACF of a time series is a function of its ACF and is a useful tool for determining the orderp of an AR model because PACF cuts off at lag p for an AR (P) process. For MA models, ACF is useful in specifying the order because ACF cuts off at lag q for an MA (q) [20].


### 2.4. Parameter Estimation

As briefly explained by Jingqiao Zhang and Sanderson, the ordinary least squares (OLS) method works great (assuming we meet some initial conditions), but one assumption that must be made for OLS to work is that the disturbance termϵt is homoscedastic [21]. However, this is not always a realistic assumption in real life since variance is not always constant. Under the presence of ARCH effects, OLS estimation is not efficient. Therefore, the study employed maximum likelihood estimation (MLE) for estimating unknown parameters in EGARCH-ANN models.For the estimation of the fuzzy-EGARCH-ANN model parameters without suffering from local optimum, the differential evolution (DE) algorithm with Archive is employed [21]. Furthermore, the step-by-step self-explanatory summary of the algorithm employed in this study is given in Table 1. For easy understanding, the software we have used throughout this study is MATLAB R2018a.Table 1
Procedure of JADE with Archive and its description.
Procedure of JADE with Archive as in [18, 21]StartInitialize:μCR=0.5; μF=0.5; A=∅; create a random initial population Xi,0:i=1:NPForg=1:G; generation loopSF=∅; SCR=∅;Fori=1:NP; population size loopGenerate:CRi= randni μCR,0.1, Fi= randciμF,0.1.Randomly choose:Xbest,gp as one of the 100 p% best vectors, Xr1,g≠Xi,g from current population P,andX˜r2,g≠Xr1,g≠Xi,g from P∪A.ComputeVi,g=Xi,g+Fi∗Xbest,gp−Xi,g+Fi∗Xr1,g−X˜r2,gGenerate:jrand= randint 1,DForj=1:D dimension loopIfj=jrand or rand 0,1<CRi then uj,i,g=vj,i,g Else uj,i,g=xj,i,g. End ifEnd for dimension loopIfEXi,g≤EUi,g then Xi,g+1=Xi,g Else Xi,g+1=Ui,g; Xi,g⟶A; CRi⟶SCR, Fi⟶SF.End ifEnd for population size loopRandomly remove solutions fromA so that the cardinal of A is less or equal to NPμCR=1−c∗μCR+c∗meanASCRμF=1−c∗μF+c∗meanLSFEnd for Generation loopEndc is a positive constant between 0 and 1, meanA(.) is the usual arithmetic mean, and meanL (·) is the Lehmer mean.The Lehmer mean is given by(9)meanLSF=ΣFεSFF2ΣFεSFF.


## 2.1. Data Source

To achieve the objectives of the study, daily S&P 500 data were collected from Yahoo Finance. As explained in the literature, the study split the given dataset into an in-sample and out-sample period. The in-sample dataset was used for initial parameter estimation and model selection. The out-sample dataset was used to evaluate forecast performance. In this study, the in-sample period ranges from 1/3/2006 to 3/5/2019 (800 observations), whereas the out-of-sample period ranges from 4/5/2019 to 20/2/2020 (201 observations).


## 2.2. Model Specification and Econometric Tests

### 2.2.1. Unit Root Test of the Variables

For determining whether a time series is stationary or nonstationary, there are many tests. This study used the augmented Dickey–Fuller (ADF) test [14] and the Phillips–Perron (PP) test [15]. Ljung–Box statistics Qm and Lagrange multiplier (LM) tests are the appropriate ARCH effect tests and were used in this study. The Ljung–Box test statistic Q was used to assess the independency among the residuals [16]. When dealing with ARCH/GARCH−type models, we first examine the characteristics of the unconditional distribution of the stock price. This will enable us to explore and explain some stylized facts that exist in the financial time series. In statistics, the Bera and Jarque test is a test of departure from normality based on the sample skewness and kurtosis [17].


### 2.2.2. Fuzzy-EGARCH-ANN Model

As proposed by the authors, the fuzzy-EGARCH-ANN model is described by a collection of fuzzy rules in the form of If-Then statements in order to describe the stock fluctuations with volatility clustering overlooked by EGARCH-ANN model via fuzzy rules and the asymmetric responses of volatility to positive and negative shocks via an EGARCH-ANN model [18].Thekth rule of the fuzzy system for EGARCH-ANN p,q,m,s model is described by the following.Rulek: if x1t is Fk1 and …xdt is Fkd, then(1)ykt=σktεkt,lnσkt2=Ek+∑j=1qβkjlnσt−j2+∑i=1pγkiyt−iσt−i−Eyt−iσt−i+∑i=1pαkiyt−iσt−i+∑h=1sξkhGZtWh+wh,0,where(2)GZtWh+wh,0=tanhZtWh+wh,0,Wh=wh,1,…wh,m′,zt−d=yt−d−EyEy2,d=1:m,xt=yt−1,yt−2,…yt−p,σt−12,σt−22,…σt−q2T=x1t,x2t,…,xdt′.As given by the authors,xt is the input vector with d=p+q at instance t. Fkl is the fuzzy set to describe the stock market return and volatility for l=1,2,3,…,d, R is the number of rules, and xlt is the premise variable. The authors assumed that the residual series distribution follows either the Gaussian normal or Student’s t-distribution. If it is not following the Gaussian normal distribution, it follows Student’s t-distribution directly.The authors had described the output of this fuzzy-EGARCH-ANN model as the weighted average of each individual rule, and it is obtained by using FIS fundamental steps as in [18, 19] as follows.Step 1 (fuzzification layer).
using the Gaussian membership function, find the grade of membership of the inputxlt in Fkl as follows:(3)Fklxlt=exp−12xlt−cklakl2,where ckl is the center and akl is the spread of kth rule membership function corresponding to the lth premise variable.Step 2 (firing strength layer).
find the firing strength ofkth rule by assuming the product T-norm of the antecedent fuzzy sets as(4)ukxt=∏l=1dFklxlt=∏l=1d−12xlt−cklakl2.Step 3 (normalization layer).
find the normalized firing strengths, i.e., the ratio of thekth rule’s firing strength to the sum of all rule’s firing strengths:(5)w¯k=ukxt∑k=1Rukxt.Step 4 (consequent and defuzzification layer).
combine the normalized firing strengths and the consequent corresponding rule to produce the model output as the weighted average of each individual rule:(6)lnσt2=∑k=1Rw¯kfk,where fk is the output of the kth rule. That is,(7)fk=Ek+∑j=1qβkjlnσt−j2+∑i=1pγkiyt−iσt−i−Eyt−iσt−i+∑i=1pαkiyt−iσt−i+∑h=1sξkhGZtWh+wh,0.
This implies that(8)lnσt2=∑k=1Rw¯kfk=∑k=1Rukxt∑k=1Rukxtfk.
The collection of theR rules assembles a model as a combination of local models. As a result, the exponential of this weighted average value gives us the predicted stock market return volatility using the fuzzy-EGARCH-ANN model [18].


## 2.2.1. Unit Root Test of the Variables

For determining whether a time series is stationary or nonstationary, there are many tests. This study used the augmented Dickey–Fuller (ADF) test [14] and the Phillips–Perron (PP) test [15]. Ljung–Box statistics Qm and Lagrange multiplier (LM) tests are the appropriate ARCH effect tests and were used in this study. The Ljung–Box test statistic Q was used to assess the independency among the residuals [16]. When dealing with ARCH/GARCH−type models, we first examine the characteristics of the unconditional distribution of the stock price. This will enable us to explore and explain some stylized facts that exist in the financial time series. In statistics, the Bera and Jarque test is a test of departure from normality based on the sample skewness and kurtosis [17].


## 2.2.2. Fuzzy-EGARCH-ANN Model

As proposed by the authors, the fuzzy-EGARCH-ANN model is described by a collection of fuzzy rules in the form of If-Then statements in order to describe the stock fluctuations with volatility clustering overlooked by EGARCH-ANN model via fuzzy rules and the asymmetric responses of volatility to positive and negative shocks via an EGARCH-ANN model [18].Thekth rule of the fuzzy system for EGARCH-ANN p,q,m,s model is described by the following.Rulek: if x1t is Fk1 and …xdt is Fkd, then(1)ykt=σktεkt,lnσkt2=Ek+∑j=1qβkjlnσt−j2+∑i=1pγkiyt−iσt−i−Eyt−iσt−i+∑i=1pαkiyt−iσt−i+∑h=1sξkhGZtWh+wh,0,where(2)GZtWh+wh,0=tanhZtWh+wh,0,Wh=wh,1,…wh,m′,zt−d=yt−d−EyEy2,d=1:m,xt=yt−1,yt−2,…yt−p,σt−12,σt−22,…σt−q2T=x1t,x2t,…,xdt′.As given by the authors,xt is the input vector with d=p+q at instance t. Fkl is the fuzzy set to describe the stock market return and volatility for l=1,2,3,…,d, R is the number of rules, and xlt is the premise variable. The authors assumed that the residual series distribution follows either the Gaussian normal or Student’s t-distribution. If it is not following the Gaussian normal distribution, it follows Student’s t-distribution directly.The authors had described the output of this fuzzy-EGARCH-ANN model as the weighted average of each individual rule, and it is obtained by using FIS fundamental steps as in [18, 19] as follows.Step 1 (fuzzification layer).
using the Gaussian membership function, find the grade of membership of the inputxlt in Fkl as follows:(3)Fklxlt=exp−12xlt−cklakl2,where ckl is the center and akl is the spread of kth rule membership function corresponding to the lth premise variable.Step 2 (firing strength layer).
find the firing strength ofkth rule by assuming the product T-norm of the antecedent fuzzy sets as(4)ukxt=∏l=1dFklxlt=∏l=1d−12xlt−cklakl2.Step 3 (normalization layer).
find the normalized firing strengths, i.e., the ratio of thekth rule’s firing strength to the sum of all rule’s firing strengths:(5)w¯k=ukxt∑k=1Rukxt.Step 4 (consequent and defuzzification layer).
combine the normalized firing strengths and the consequent corresponding rule to produce the model output as the weighted average of each individual rule:(6)lnσt2=∑k=1Rw¯kfk,where fk is the output of the kth rule. That is,(7)fk=Ek+∑j=1qβkjlnσt−j2+∑i=1pγkiyt−iσt−i−Eyt−iσt−i+∑i=1pαkiyt−iσt−i+∑h=1sξkhGZtWh+wh,0.
This implies that(8)lnσt2=∑k=1Rw¯kfk=∑k=1Rukxt∑k=1Rukxtfk.
The collection of theR rules assembles a model as a combination of local models. As a result, the exponential of this weighted average value gives us the predicted stock market return volatility using the fuzzy-EGARCH-ANN model [18].


## 2.3. Choosing the Optimal Lag Length and Model Selection Criterion

The optimal lag length for ARIMA and EGARCH-ANN models was chosen using the ACF, PACF plot, and information criteria. Some of the most popular information criteria are AIC and BIC, which are used to select the appropriate models. The PACF of a time series is a function of its ACF and is a useful tool for determining the orderp of an AR model because PACF cuts off at lag p for an AR (P) process. For MA models, ACF is useful in specifying the order because ACF cuts off at lag q for an MA (q) [20].


## 2.4. Parameter Estimation

As briefly explained by Jingqiao Zhang and Sanderson, the ordinary least squares (OLS) method works great (assuming we meet some initial conditions), but one assumption that must be made for OLS to work is that the disturbance termϵt is homoscedastic [21]. However, this is not always a realistic assumption in real life since variance is not always constant. Under the presence of ARCH effects, OLS estimation is not efficient. Therefore, the study employed maximum likelihood estimation (MLE) for estimating unknown parameters in EGARCH-ANN models.For the estimation of the fuzzy-EGARCH-ANN model parameters without suffering from local optimum, the differential evolution (DE) algorithm with Archive is employed [21]. Furthermore, the step-by-step self-explanatory summary of the algorithm employed in this study is given in Table 1. For easy understanding, the software we have used throughout this study is MATLAB R2018a.Table 1
Procedure of JADE with Archive and its description.
Procedure of JADE with Archive as in [18, 21]StartInitialize:μCR=0.5; μF=0.5; A=∅; create a random initial population Xi,0:i=1:NPForg=1:G; generation loopSF=∅; SCR=∅;Fori=1:NP; population size loopGenerate:CRi= randni μCR,0.1, Fi= randciμF,0.1.Randomly choose:Xbest,gp as one of the 100 p% best vectors, Xr1,g≠Xi,g from current population P,andX˜r2,g≠Xr1,g≠Xi,g from P∪A.ComputeVi,g=Xi,g+Fi∗Xbest,gp−Xi,g+Fi∗Xr1,g−X˜r2,gGenerate:jrand= randint 1,DForj=1:D dimension loopIfj=jrand or rand 0,1<CRi then uj,i,g=vj,i,g Else uj,i,g=xj,i,g. End ifEnd for dimension loopIfEXi,g≤EUi,g then Xi,g+1=Xi,g Else Xi,g+1=Ui,g; Xi,g⟶A; CRi⟶SCR, Fi⟶SF.End ifEnd for population size loopRandomly remove solutions fromA so that the cardinal of A is less or equal to NPμCR=1−c∗μCR+c∗meanASCRμF=1−c∗μF+c∗meanLSFEnd for Generation loopEndc is a positive constant between 0 and 1, meanA(.) is the usual arithmetic mean, and meanL (·) is the Lehmer mean.The Lehmer mean is given by(9)meanLSF=ΣFεSFF2ΣFεSFF.


## 3. Results and Discussion

### 3.1. Calibration and Validation of the Fuzzy-EGARCH-ANN Model

This work investigates a global optimization algorithm for the calibration of stochastic volatility models. It is shown that commonly used gradient-based optimization procedures may not lead to a good solution and often converge to a local optimum. Many evolutionary algorithms have been introduced to explore the whole parameter domain thoroughly. The number of iterations, population size, and other user-defined inputs influence solution quality and computation time. Practical information from the literature shows that the evolutionary algorithm approach outperforms the standard gradient-based method.When using optimization techniques to find model parameters that minimize the error between the stock’s model prices and market prices, the calibration’s success can substantially depend on the optimization algorithm. The calibration performed with an unsuitable algorithm can lead to parameter instability. As will be demonstrated, the optimization method used to calibrate a model can become as crucial as the model itself. A loss function describes the difference between empirical and model prices. For the same reason mentioned by the authors [18], this study selected the most commonly used loss function, the so-called mean squared forecast error (MSFE), and used it for the fuzzy-EGARCH-ANN model calibration and validation.Generally, numerical models and field data are imperfect representations of the real world. The goal is to obtain a reasonably accurate representation that is at least internally consistent. This consistency is achieved by calibrating or adjusting the model until errors are minimized, as measured by qualitative or quantitative means. Figure1 describes the general model calibration and validation of the GARCH-type models.Figure 1
Flowchart for calibration and validation of EGARCH and EGARCH-ANN models.


#### 3.1.1. Simulation Model Setup

This step consists of the definition of study scope and purpose, determination of the performance measure, field data collection and reduction, and network coding [1, 9].


#### 3.1.2. Initial Evaluation

This step tests whether the default calibration parameters are acceptable by comparing the distribution of a selected performance measure [1, 9].Since the parameter estimation technique for both EGARCH and EGARCH-ANN models is MLE which is a gradient-based optimization method [18], it might not give us a globally optimal solution to these volatility models’ parameter estimation problems. That is why the study is required to perform model calibration, that is, model calibration using evolutionary algorithms.From the estimation of EGARCH-ANN model parameters, as explained by the authors, the study is only interested in the best orders of the EGARCH-ANN model but not to the parameters’ values model [18]. As a result, at the stage of the EGARCH-ANN model parameter estimation, we are not required to perform model calibration for EGARCH-ANN. This is because we are required to perform model calibration to adjust parameters to the one that can take us to reality after we have identified the parameter’s order. Hence, performing model calibration for the EGARCH-ANN model cannot change the fitted EGARCH-ANN model’s order. Therefore, for our case, there is no need to perform model calibration under the EGARCH-ANN model fitting.As explained by the authors, the parameter estimation method for the fuzzy-EGARCH-ANN model was one of the evolutionary algorithms, the so-called differential evolution (DE) algorithm, which has advantages over the other evolutionary algorithms [18, 21]. By itself, it gives a globally optimum solution to the optimization problem of the fuzzy-EGARCH-ANN model parameter estimation [18]. Its accuracy depends upon some user-provided inputs like the number of iterations, population size, and so on. The recommended amount of iteration for JADE is 500, and the population size is between 5D and 10D, where D is the dimension of the given parameter space of the fuzzy-EGARCH-ANN model [18, 21, 22]. Regarding the recommended number of iterations and population size employed, the JADE with Archive algorithm can provide us globally optimum parameters of the fuzzy-EGARCH-ANN model. Therefore, besides using the recommended amount of iteration and population size, there is no need to perform continuous model calibration for the fuzzy-EGARCH-ANN volatility model. Hence, for justification of this, we decided to apply this model to one of the well-known financial time series data.


### 3.2. The Fuzzy-EGARCH-ANN Model Applications to the S&P 500 Data

We applied the fuzzy-EGARCH-ANN model onto the S&P 500 time series data to show that the fuzzy-EGARCH-ANN model with its parameter estimation method does not need model calibration and validation.The plot of daily S&P 500 data for the period of 1/3/2006 to 20/2/2020 is shown in Figure2. From this figure, it is evident that unconditional mean and variance are changing over time, and the series has an increasing trend over time. The changing mean and variance of daily S&P 500 data over time is an indication of the nonstationarity of the level series. This implies that it is challenging to model ARIMA and ARCH/GARCH-type models for nonstationary series. Therefore, to achieve stationarity, logarithmic transformation was applied.Figure 2
The daily S&P 500 stock price series.Table2 presents the findings of the ADF test and PP test and formally confirms that the transformation of natural logarithms is stationary.Table 2
Statistical properties of daily S&P 500 log return data.
Statistical properties of daily S&P 500 log return dataMin.−0.043 945Max.0.045 966Mean−0.000 355Median−0.000 228Variance0.000 060Skewness−0.802 855Kurtosis8.347 911Arch testp value = 2.2 204e − 16Jarque–Bera test (for alpha = 0.05)p value = 0ADF testp value = 1.0 000e-03PP testp value = 1.0 000e-03Table2 portrays the ARCH effect tests on the residuals of ARIMA (1, 0, 0) model. These results show the presence of ARCH effects in the daily S&P 500 data. This suggests that modeling ARCH/GARCH types is appropriate for the dataset. As a result, the obtained basic statistical properties of these data are given in Table 2.Note thatMean< median usually flags left skewness.Negative skewness means left skewed data.Kurtosis>3 means larger peakedness than Gaussian.Arch test result with very smallp value implies existence of strong heteroscedasticity.Jarque–Bera test result with very smallp value implies that the residual series strongly follows Student’s t-distribution.Furthermore, Figures2–5 are the corresponding figures of the daily S&P 500 data series, log return series, residual series, and histogram of the residual series, respectively.Figure 3
The daily S&P 500 log return series.Figure 4
The S&P 500 residual series.Figure 5
The histogram of the S&P 500 residual series.Now by using the first 800 (that is, forT = 800) simulated data points, the study has obtained the following results. Since the study has identified that the residual series follows Student’s t-distribution, then by using this distribution type, EGARCH−ANN1,1,2,1 model is obtained as the best fitted model to these selected data. Moreover, its fitness summary is presented in Table 3.Table 3
The EGARCH-ANN (1,1,2,1) model fitted to the given S&P 500 data.
ParametersCoefficientsStd. errorst-StatsK−0.308 20.084 6−3.642 0ARCH10.217 60.050 84.282 3GARCH10.983 50.009 998.999Leverage1−0.094 50.040 7−2.323 4leveragebyANN1−0.521 00.241 3−2.159 0DoF4.226 10.673 76.273 2With addition information like log likelihood: 2956, Akaike Information Criterion: -5900 Bayesian Information Criterion: -5872 Hannan-Quinn Information Criterion: -5889.Now, to estimate the fuzzy-EGARCH-ANN model parameters, the number of rules (R) and the realized volatility are the ones to be determined. For our case, as mentioned earlier,R is going to be determined by the subtractive clustering algorithm (SCA), and the realized volatility at time t is the square root of the square of residual series at time t. R = 2 is obtained, and the graph of the realized volatility is given in Figure 6.Figure 6
The realized volatility of daily S&P 500 series.Therefore, forp=1, q=1, m=2, and s=1, by JADE with Achieve, the estimated parameter vector is given by Table 4. In addition, the corresponding mean squared forecast error (MSFE) value is 2.4907e−05, and MAE is 0.0034.Table 4
The estimated parameters of fuzzy-EGARCH-ANN model on the first 800 given S&P 500 data points.
ParametersCoefficientsParametersCoefficientsE1−2.724 7E2−0.833 9α11−3.168 7α12−2.789 7γ11−0.811 8γ121.360 3β110.670 8β120.526 2ξ11−0.052 8ξ12−0.077 3Figure7 presents a summary and analysis of the MSFE obtained by JADE with Archive during the parameter estimation of the proposed model.Figure 7
The MSFE obtained by JADE with Archive for S&P 500 data.Finally, Figures8 and 9 present a summary of estimated volatility by the proposed model for T = 800 and validation for T = 1000, respectively. Here, the validation dataset refers to the last 200 data points.Figure 8
The estimated volatility of daily S&P 500 data forT = 800.Figure 9
The validation of volatility of daily S&P 500 data forT = 1000.Note. The mean squared forecast error (MSFE) for this validation is 2.6 045e − 05, which is very near to that of the mean squared forecast error (MSFE) obtained during the parameter estimation. In addition, its corresponding MAE value is 0.0036. Here, the model has fitted the data very well. As a result, the following forecasted values were obtained for different periods, as shown in Table 5.Table 5
One-step-ahead volatility prediction analysis of the fuzzy-EGARCH-ANN(1,1,2,1) model.
The realized volatility valueT = 900T = 990T = 995T = 999T = 10000.006 380.005400.003 020.007140.002 190.0052070.004 480.004420.00403


## 3.1. Calibration and Validation of the Fuzzy-EGARCH-ANN Model

This work investigates a global optimization algorithm for the calibration of stochastic volatility models. It is shown that commonly used gradient-based optimization procedures may not lead to a good solution and often converge to a local optimum. Many evolutionary algorithms have been introduced to explore the whole parameter domain thoroughly. The number of iterations, population size, and other user-defined inputs influence solution quality and computation time. Practical information from the literature shows that the evolutionary algorithm approach outperforms the standard gradient-based method.When using optimization techniques to find model parameters that minimize the error between the stock’s model prices and market prices, the calibration’s success can substantially depend on the optimization algorithm. The calibration performed with an unsuitable algorithm can lead to parameter instability. As will be demonstrated, the optimization method used to calibrate a model can become as crucial as the model itself. A loss function describes the difference between empirical and model prices. For the same reason mentioned by the authors [18], this study selected the most commonly used loss function, the so-called mean squared forecast error (MSFE), and used it for the fuzzy-EGARCH-ANN model calibration and validation.Generally, numerical models and field data are imperfect representations of the real world. The goal is to obtain a reasonably accurate representation that is at least internally consistent. This consistency is achieved by calibrating or adjusting the model until errors are minimized, as measured by qualitative or quantitative means. Figure1 describes the general model calibration and validation of the GARCH-type models.Figure 1
Flowchart for calibration and validation of EGARCH and EGARCH-ANN models.


### 3.1.1. Simulation Model Setup

This step consists of the definition of study scope and purpose, determination of the performance measure, field data collection and reduction, and network coding [1, 9].


### 3.1.2. Initial Evaluation

This step tests whether the default calibration parameters are acceptable by comparing the distribution of a selected performance measure [1, 9].Since the parameter estimation technique for both EGARCH and EGARCH-ANN models is MLE which is a gradient-based optimization method [18], it might not give us a globally optimal solution to these volatility models’ parameter estimation problems. That is why the study is required to perform model calibration, that is, model calibration using evolutionary algorithms.From the estimation of EGARCH-ANN model parameters, as explained by the authors, the study is only interested in the best orders of the EGARCH-ANN model but not to the parameters’ values model [18]. As a result, at the stage of the EGARCH-ANN model parameter estimation, we are not required to perform model calibration for EGARCH-ANN. This is because we are required to perform model calibration to adjust parameters to the one that can take us to reality after we have identified the parameter’s order. Hence, performing model calibration for the EGARCH-ANN model cannot change the fitted EGARCH-ANN model’s order. Therefore, for our case, there is no need to perform model calibration under the EGARCH-ANN model fitting.As explained by the authors, the parameter estimation method for the fuzzy-EGARCH-ANN model was one of the evolutionary algorithms, the so-called differential evolution (DE) algorithm, which has advantages over the other evolutionary algorithms [18, 21]. By itself, it gives a globally optimum solution to the optimization problem of the fuzzy-EGARCH-ANN model parameter estimation [18]. Its accuracy depends upon some user-provided inputs like the number of iterations, population size, and so on. The recommended amount of iteration for JADE is 500, and the population size is between 5D and 10D, where D is the dimension of the given parameter space of the fuzzy-EGARCH-ANN model [18, 21, 22]. Regarding the recommended number of iterations and population size employed, the JADE with Archive algorithm can provide us globally optimum parameters of the fuzzy-EGARCH-ANN model. Therefore, besides using the recommended amount of iteration and population size, there is no need to perform continuous model calibration for the fuzzy-EGARCH-ANN volatility model. Hence, for justification of this, we decided to apply this model to one of the well-known financial time series data.


## 3.1.1. Simulation Model Setup

This step consists of the definition of study scope and purpose, determination of the performance measure, field data collection and reduction, and network coding [1, 9].


## 3.1.2. Initial Evaluation

This step tests whether the default calibration parameters are acceptable by comparing the distribution of a selected performance measure [1, 9].Since the parameter estimation technique for both EGARCH and EGARCH-ANN models is MLE which is a gradient-based optimization method [18], it might not give us a globally optimal solution to these volatility models’ parameter estimation problems. That is why the study is required to perform model calibration, that is, model calibration using evolutionary algorithms.From the estimation of EGARCH-ANN model parameters, as explained by the authors, the study is only interested in the best orders of the EGARCH-ANN model but not to the parameters’ values model [18]. As a result, at the stage of the EGARCH-ANN model parameter estimation, we are not required to perform model calibration for EGARCH-ANN. This is because we are required to perform model calibration to adjust parameters to the one that can take us to reality after we have identified the parameter’s order. Hence, performing model calibration for the EGARCH-ANN model cannot change the fitted EGARCH-ANN model’s order. Therefore, for our case, there is no need to perform model calibration under the EGARCH-ANN model fitting.As explained by the authors, the parameter estimation method for the fuzzy-EGARCH-ANN model was one of the evolutionary algorithms, the so-called differential evolution (DE) algorithm, which has advantages over the other evolutionary algorithms [18, 21]. By itself, it gives a globally optimum solution to the optimization problem of the fuzzy-EGARCH-ANN model parameter estimation [18]. Its accuracy depends upon some user-provided inputs like the number of iterations, population size, and so on. The recommended amount of iteration for JADE is 500, and the population size is between 5D and 10D, where D is the dimension of the given parameter space of the fuzzy-EGARCH-ANN model [18, 21, 22]. Regarding the recommended number of iterations and population size employed, the JADE with Archive algorithm can provide us globally optimum parameters of the fuzzy-EGARCH-ANN model. Therefore, besides using the recommended amount of iteration and population size, there is no need to perform continuous model calibration for the fuzzy-EGARCH-ANN volatility model. Hence, for justification of this, we decided to apply this model to one of the well-known financial time series data.


## 3.2. The Fuzzy-EGARCH-ANN Model Applications to the S&P 500 Data

We applied the fuzzy-EGARCH-ANN model onto the S&P 500 time series data to show that the fuzzy-EGARCH-ANN model with its parameter estimation method does not need model calibration and validation.The plot of daily S&P 500 data for the period of 1/3/2006 to 20/2/2020 is shown in Figure2. From this figure, it is evident that unconditional mean and variance are changing over time, and the series has an increasing trend over time. The changing mean and variance of daily S&P 500 data over time is an indication of the nonstationarity of the level series. This implies that it is challenging to model ARIMA and ARCH/GARCH-type models for nonstationary series. Therefore, to achieve stationarity, logarithmic transformation was applied.Figure 2
The daily S&P 500 stock price series.Table2 presents the findings of the ADF test and PP test and formally confirms that the transformation of natural logarithms is stationary.Table 2
Statistical properties of daily S&P 500 log return data.
Statistical properties of daily S&P 500 log return dataMin.−0.043 945Max.0.045 966Mean−0.000 355Median−0.000 228Variance0.000 060Skewness−0.802 855Kurtosis8.347 911Arch testp value = 2.2 204e − 16Jarque–Bera test (for alpha = 0.05)p value = 0ADF testp value = 1.0 000e-03PP testp value = 1.0 000e-03Table2 portrays the ARCH effect tests on the residuals of ARIMA (1, 0, 0) model. These results show the presence of ARCH effects in the daily S&P 500 data. This suggests that modeling ARCH/GARCH types is appropriate for the dataset. As a result, the obtained basic statistical properties of these data are given in Table 2.Note thatMean< median usually flags left skewness.Negative skewness means left skewed data.Kurtosis>3 means larger peakedness than Gaussian.Arch test result with very smallp value implies existence of strong heteroscedasticity.Jarque–Bera test result with very smallp value implies that the residual series strongly follows Student’s t-distribution.Furthermore, Figures2–5 are the corresponding figures of the daily S&P 500 data series, log return series, residual series, and histogram of the residual series, respectively.Figure 3
The daily S&P 500 log return series.Figure 4
The S&P 500 residual series.Figure 5
The histogram of the S&P 500 residual series.Now by using the first 800 (that is, forT = 800) simulated data points, the study has obtained the following results. Since the study has identified that the residual series follows Student’s t-distribution, then by using this distribution type, EGARCH−ANN1,1,2,1 model is obtained as the best fitted model to these selected data. Moreover, its fitness summary is presented in Table 3.Table 3
The EGARCH-ANN (1,1,2,1) model fitted to the given S&P 500 data.
ParametersCoefficientsStd. errorst-StatsK−0.308 20.084 6−3.642 0ARCH10.217 60.050 84.282 3GARCH10.983 50.009 998.999Leverage1−0.094 50.040 7−2.323 4leveragebyANN1−0.521 00.241 3−2.159 0DoF4.226 10.673 76.273 2With addition information like log likelihood: 2956, Akaike Information Criterion: -5900 Bayesian Information Criterion: -5872 Hannan-Quinn Information Criterion: -5889.Now, to estimate the fuzzy-EGARCH-ANN model parameters, the number of rules (R) and the realized volatility are the ones to be determined. For our case, as mentioned earlier,R is going to be determined by the subtractive clustering algorithm (SCA), and the realized volatility at time t is the square root of the square of residual series at time t. R = 2 is obtained, and the graph of the realized volatility is given in Figure 6.Figure 6
The realized volatility of daily S&P 500 series.Therefore, forp=1, q=1, m=2, and s=1, by JADE with Achieve, the estimated parameter vector is given by Table 4. In addition, the corresponding mean squared forecast error (MSFE) value is 2.4907e−05, and MAE is 0.0034.Table 4
The estimated parameters of fuzzy-EGARCH-ANN model on the first 800 given S&P 500 data points.
ParametersCoefficientsParametersCoefficientsE1−2.724 7E2−0.833 9α11−3.168 7α12−2.789 7γ11−0.811 8γ121.360 3β110.670 8β120.526 2ξ11−0.052 8ξ12−0.077 3Figure7 presents a summary and analysis of the MSFE obtained by JADE with Archive during the parameter estimation of the proposed model.Figure 7
The MSFE obtained by JADE with Archive for S&P 500 data.Finally, Figures8 and 9 present a summary of estimated volatility by the proposed model for T = 800 and validation for T = 1000, respectively. Here, the validation dataset refers to the last 200 data points.Figure 8
The estimated volatility of daily S&P 500 data forT = 800.Figure 9
The validation of volatility of daily S&P 500 data forT = 1000.Note. The mean squared forecast error (MSFE) for this validation is 2.6 045e − 05, which is very near to that of the mean squared forecast error (MSFE) obtained during the parameter estimation. In addition, its corresponding MAE value is 0.0036. Here, the model has fitted the data very well. As a result, the following forecasted values were obtained for different periods, as shown in Table 5.Table 5
One-step-ahead volatility prediction analysis of the fuzzy-EGARCH-ANN(1,1,2,1) model.
The realized volatility valueT = 900T = 990T = 995T = 999T = 10000.006 380.005400.003 020.007140.002 190.0052070.004 480.004420.00403


## 4. Conclusions

This study provides a robust analysis of volatility forecasting of daily S&P 500 data using weekly data spanning from 1/3/2006 to 19/2/2020. The forecasting performance of the semiparametric nonlinear fuzzy-EGARCH-ANN model is investigated based on forecasting performance criteria such as MSE and MAE-based tests and alternative measures of realized volatility. It is evident that fuzzy-EGARCH-ANN(1,1,2,1) with Student’st-distribution is found to perform best in terms of one-step-ahead forecasting. To our knowledge, this is the first study that focuses on daily S&P 500 data using high-frequency data and the fuzzy-EGARCH-ANN econometric model. From the results and discussion section of this study and the concepts of calibration that have been stated in this study, it is concluded that there is no need for model calibration for the proposed model as far as the practitioner uses the corresponding DE algorithm suitably. Finally, except the computational time, the model is novel.



---

*Source: 6637091-2021-12-24.xml*