An Empirical Evaluation in GARCH Volatility Modeling: Evidence from the Stockholm Stock Exchange

Show more

1. Introduction

The development of econometrics led to the invention of adjusted methodologies for the modeling of mean value and variance. Models of generalized conditional autoregressive heteroscedasticity (GARCH) are based on the assumption that random components in models present changes on volatility. These models were developed by Engle [2] , in a simple form, and they were generalized later by Bollerslev [3] .

The models of autoregressive conditional heteroscedasticity (GARCH) have a long and noteworthy history but they are not free of limitations. For example, Black [4] on his paper claims that stock market returns are negatively correlated with changes on volatility returns implying that volatility tends to rise in response to bad news and fall in response to good news. On the other hand, on GARCH models we assume that only the size of return of the conditional variance is defined and not the positivity or negativity of volatility’s return, which are unpredicted. Another crucial limitation of GARCH models is the non-negativity of parameters in order to ensure the positivity of the conditional variance. All these limitations cause difficulties in the estimation of GARCH models.

GARCH models were applied with great success on the modeling of changing variability or the variance volatility on time series for measuring financial investments. After the determination of an asymmetric relationship between conditional volatility and conditional mean value, econometricians focused their efforts on planning methodologies for modeling this phenomenon.

Nelson [5] suggested an exponential GARCH model (EGARCH) expressed in logarithms of the conditional variance volatility. The EGARCH model has become popular as it presents asymmetric volatility on positive and negative returns. A number of modifications on this model were made over the years. Glosten, Jagannathan and Runkle [6] suggest another asymmetric model known as GJR-GARCH which deals with the limitations of the symmetric GARCH models.

The purpose of this paper is to quantify two asymmetric models using prices from the Stockholm stock market for the period 30 September 1986 until 11 May 2016, representing 7,434 observations. The first 7,000 values in the model were used for quantification and statistical verification and the last 434 values for the forecast demonstration in retrospect.

The remainder of the paper is organized as follows: Section 2 provides a brief literature review. Section 3 discusses the symmetric and asymmetric GARCH models. Section 4 summarizes the data. The results are discussed in Section 5 and Section 6 proposes the forecasting methodology. Finally, the last section offers the concluding remarks.

2. Literature Review

The ability of GARCH models that study the relationship between risk and return has been validated in many studies. For example, Donaldson and Kamstra [7] made a nonlinear GARCH model based on neural networks. They evaluated the model’s ability in forecasting the volatility of returns on the stock markets of London, New York, Tokyo and Toronto. The results of their paper showed that neural network models captures volatility effects and forecasting better than GARCH, EGARCH and GJR models.

Nam, Pyun and Arize [8] used the asymmetric non-linear GARCH-M model for US market indices for the period 1926:01-1997:12. The results of their paper showed that negative returns on average reverted more quickly in the long term rather than positive returns.

Tudor [9] uses GARCH and GARCH-M models to examine the volatility of US and Romanian stock markets for the period January, 03 2001 until February 09, 2008 or a total of 1,853 daily returns. The results showed that the GARCH-M model performs better and confirms the results between volatility and expected returns on both markets.

Panait and Slavescu [10] use daily, weekly and monthly data for seven Romanian listed companies on the Bucharest stock market for the period 1997-2012. Using the GARCH-in-mean model they compare the volatility of companies in three phases. The results of their paper showed that persistency is more evident in the daily returns rather than in the weekly and monthly series. Furthermore, the GARCH-in-mean model failed to confirm that an increase in volatility leads to a rise in future returns.

Gao, Zhang and Zhang [11] use the Markov chain Monte Carlo (MCMC) method instead of the Maximum Likelihood Ratio method for the estimation of the coefficients of GARCH models. Using daily data from the stock market of China for the period 1 January 2000 until 29 April 2011, they compare the results of the volatility of the data with three different models and two different distributions. The results showed that the GED-GARCH model is better than the t-GARCH, and that the t-GARCH is better than N-GARCH.

Dutta [12] examined exchange rate parities of USA and Japan for the period 1 January 2000 until 31 January 2012. The data are estimated not only with symmetric but also with asymmetric GARCH models. Findings indicate that positive shocks are more common than the negative shocks in this return series. Also, asymmetric tests for volatility show a size effect on news, which is stronger for good news than for bad news.

Given the different backgrounds for each market it is expected that risk and return differ from country to country. This paper attempts to examine the volatility and return of the Stockholm stock market using symmetric and asymmetric models.

3. Methodology

One of the fundamental hypotheses for a stationary time series is the stable variance. But there are time series, mostly financial, that have intervals of large volatility. These series are characterized with periods of sharp increases and downturns during which their variance is varying. Thus, researchers are not interested in examining the variance of such series throughout the sample period but only interested in the varying or conditional variance. Based on this division between conditional and unconditional variance, we can characterize time series models with conditional variance as conditional heteroskedastic models. The notion of “conditional” heteroscedasticity was first introduced by Engle [2] . Engle suggested that varying variance can be explained through an autoregressive scheme as a function of previous values. For this reason this model is called Autoregressive Conditional Heteroskedastic Model, known as ARCH.

The varying GARCH models consist of two equations. The first one (the equation of mean) describes the data as a function of other variables adding an error term. The second equation (the equation of variance) determines the evolution of the conditional variance of the error from the mean equation as a function of past conditional variances and lagged errors. The first equation (equation of the mean) on GARCH varying models is not of great interest, contrary to the second equation (equation of variance) which is the one that we pay more attention to in order to compare different variances on the same equation of mean.

3.1. Symmetric GARCH Models

The traditional measuring methods of volatility (variance or standard deviation) are absolute and cannot conceive the characteristics of financial data (time series) like volatility clustering, asymmetries, leverage effect and long memory. The basic model suggested by Engel [2] is the following:

${\epsilon}_{t}={z}_{t}{\sigma}_{t}$

where ${z}_{t}$ is an independent identically distributed (i.i.d.) process with mean zero and variance 1. ${\sigma}_{t}$ is the volatility that evolves over time. The volatility ${\sigma}_{t}^{2}$ in the basic ARCH (q) model is defined as:

${\sigma}_{t}^{2}=\omega +{\displaystyle \underset{i=1}{\overset{q}{\sum}}{\alpha}_{i}{\epsilon}_{t-i}^{2}}$ (1)

where ${\sigma}_{t}^{2}$ is the conditional variance, $\omega >0$ and ${\alpha}_{i}\ge 0$ for ${\sigma}_{t}^{2}$ to be positive.

This model shows that after a large (small) shock, it is likely that a large (small) shock will follow. In other words, a large (small) ${\epsilon}_{t-1}^{2}$ implies a large (small) ${\epsilon}_{t}^{2}$ on the current period according to Equation (1) thus a large (small) variance (volatility).

Finally, we can point out that large values on time lags on ARCH models presuppose large periods of volatility contrary to small values of time lags that foresee smooth periods. This may not occur in reality. In order to overcome this problem, Bollerslev, Chou and Kroner [13] suggested a new model where the conditional variance does not depend only on previous square error values but also on previous values of the same variance. The model suggested is known as Generalized Autoregressive Conditional Heteroskedasticity (GARCH) model.

GARCH (p, q) Model

The generalized form of GARCH(p, q) is given as:

${R}_{t}=\mu +{\epsilon}_{t}$ (mean equation)

${\sigma}_{t}^{2}=\omega +{\displaystyle \underset{i=1}{\overset{q}{\sum}}{\alpha}_{i}{\epsilon}_{t-i}^{2}}+{\displaystyle \underset{j=1}{\overset{p}{\sum}}{\beta}_{j}{\sigma}_{t-j}^{2}}$ (variance equation) (2)

where ${R}_{t}$ are returns of time series at time t, $\mu $ is the mean value of the returns, ${\epsilon}_{t}$ is the error term at time t, which is assumed to be normally distributed with zero mean and conditional variance ${\sigma}_{t}^{2}$ , p is the order of GARCH and q is the order of ARCH process, $\mu $ , $\omega $ , ${\alpha}_{i}$ and ${\beta}_{j}$ are parameters for estimation. All parameters in variance equation must be positive ( $\mu >0$ , $\omega >0$ , ${\alpha}_{i}\ge 0$ , and ${\beta}_{j}\ge 0$ for ${\sigma}_{t}^{2}$ to be positive). Also, we expect the value of parameter ω to be small. Parameter ${\alpha}_{i}$ measures the response of volatility on market variances and parameter ${\beta}_{j}$ expresses the difference which was caused from outliers on conditional variance. Finally, we expect the sum ${\alpha}_{i}+{\beta}_{j}<1$

The model GARCH (1,1) has the following form:

${R}_{t}=\mu +{\epsilon}_{t}$ (mean equation)

${\sigma}_{t}^{2}=\omega +{\alpha}_{1}{\epsilon}_{t-1}^{2}+{\beta}_{1}{\sigma}_{t-1}^{2}$ (variance equation) (3)

As we expect a positive variance, we can argue that regression coefficients are always positive $\omega \ge 0$ , ${\alpha}_{1}\ge 0$ and ${\beta}_{1}\ge 0$ . Also, we should point out that in order to achieve stationarity on the variance, regression coefficients ${\alpha}_{1}$ and ${\beta}_{1}$ should be less than one $\left({\alpha}_{1}<1\right)$ and $\left({\beta}_{1}<1\right)$ . Thus, on the previous model, the following relations are valid: $\omega \ge 0$ ${\alpha}_{1}\ge 0$ and ${\beta}_{1}\ge 0$ for a positive value of ${\sigma}_{t}^{2}$ and ${\alpha}_{1}<1$ and ${\beta}_{1}<1$ .

The conditional variance of the returns of Equation (3) is defined from three outcomes:

・ The constant given by $\omega $ coefficient.

・ The variance part expressed from the relationship ${\alpha}_{1}{\epsilon}_{t-1}^{2}$ defined as ARCH component.

・ The part of predicted variance from past period expressed by ${\beta}_{1}{\sigma}_{t-1}^{2}$ and is called GARCH.

The sum of regression coefficients ${\alpha}_{1}+{\beta}_{1}$ expresses the impact of variables’ variance of the previous period regarding the current value of volatility. This value is usually near to one and is regarded as a sign of increasing inactivity of shocks of the volatility of returns on the financial assets.

3.2. Asymmetric GARCH Models

The main disadvantage of GARCH models is their inappropriateness in the cases where an asymmetric effect is usually observed and is registered from a different instability in the case of good and bad news. In the asymmetric models, upward and downward trends of returns are interpreted as bad and good news. If the decline of a return is accompanied with an increase of instability larger than the instability caused by the increase then it is said to have a leverage effect.

Given that all terms in a GARCH model are squared, there will always be an asymmetric response in positive and negative periods. However, due to natural leverage in most companies, a negative shock is more damaging than a positive shock because it produces larger volatility.

Among the most widely known asymmetric models are the Exponential GARCH model (EGARCH) and the asymmetric GJR model.

3.2.1. Asymmetric GARCH Models

One of the most popular asymmetric ARCH models is the EGARCH model proposed by Nelson [5] . The EGARCH (p, q) model is given by

$\mathrm{log}{\sigma}_{t}^{2}=\omega +{\displaystyle \underset{i=1}{\overset{p}{\sum}}{\alpha}_{i}\left|\frac{{\epsilon}_{t-i}}{{\sigma}_{t-i}}\right|}+{\displaystyle \underset{j=1}{\overset{q}{\sum}}{\beta}_{j}\mathrm{log}{\sigma}_{t-j}^{2}+{\displaystyle \underset{k=1}{\overset{r}{\sum}}{\gamma}_{k}}}\frac{{\epsilon}_{t-k}}{{\sigma}_{t-k}}$ (4)

where $\omega $ , ${\alpha}_{i}$ , ${\beta}_{j}$ , and ${\gamma}_{k}$ are parameters which can be estimated using the maximum likelihood method. We should also point out that $\left|{\beta}_{j}\right|<1$ and ${\gamma}_{k}$ parameter is the one that gives the result of leverage effect. In other words, we consider that ${\epsilon}_{t-k}$ term is the one that establishes the asymmetry of EGARCH (p, q) when parameter ${\gamma}_{k}\ne 0$ . Also, when parameter ${\gamma}_{k}<0$ , then positive shocks cause short volatility in relation to negative shocks. Furthermore, we expect that parameters ${\gamma}_{k}+{\alpha}_{i}>0$ , given that parameter ${\gamma}_{k}<0$ .

The conditional variance of the above model is expressed in logarithmic form which ensures the non-negativity without imposing more constraints of non-

negativity. The term $\frac{{\epsilon}_{t-k}}{{\sigma}_{t-k}}$ on the above equation represents the asymmetric

effect of shocks. According to Poon and Granger [14] , a negative shock that leads to the largest conditional variance will not be the same on a positive shock in the next period.

The EGARCH (1,1) model is often used for the estimation of variance ${\sigma}^{2}$ and has the following form:

$\mathrm{log}{\sigma}_{t}^{2}=\omega +{\alpha}_{1}\left|\frac{{\epsilon}_{t-1}}{{\sigma}_{t-1}}\right|+{\beta}_{1}\mathrm{log}{\sigma}_{t-1}^{2}+{\gamma}_{1}\frac{{\epsilon}_{t-1}}{{\sigma}_{t-1}}$ (5)

For a positive shock $\frac{{\epsilon}_{t-1}}{{\sigma}_{t-1}}>0$ the above equation becomes

$\mathrm{log}{\sigma}_{t}^{2}=\omega +{\beta}_{1}\mathrm{log}{\sigma}_{t-1}^{2}+\left({\alpha}_{1}+{\gamma}_{1}\right)\frac{{\epsilon}_{t-1}}{{\sigma}_{t-1}}$ (6)

whereas for a negative shock $\frac{{\epsilon}_{t-1}}{{\sigma}_{t-1}}<0$ the above equation becomes

$\mathrm{log}{\sigma}_{t}^{2}=\omega +{\beta}_{1}\mathrm{log}{\sigma}_{t-1}^{2}+\left({\alpha}_{1}-{\gamma}_{1}\right)\frac{{\epsilon}_{t-1}}{{\sigma}_{t-1}}$ (7)

The EGARCH model has many advantages when compared to the GARCH (p, q) model.

・ The first is the logarithmic form which does not allow the positive constraint among parameters.

・ Another advantage of EGARCH model is that it incorporates asymmetries in the change of volatility of returns.

・ Parameters $\alpha $ and $\gamma $ define two important asymmetries in the conditional variance. If ${\gamma}_{1}<0$ then negative changes increase volatility (instability) more than positive changes of the same size.

・ EGARCH model can successfully define the change of volatility.

The form of the EGARCH model denotes that the conditional variance is an exponential function of the examined variables which ensures a positive character. In other words, conditional variance ensures the exponential nature of the EGARCH model where external changes will have a stronger influence on the predicted volatility than TGARCH. An asymmetric influence is indicated by the no null value of γ_{1} coefficient whereas the presence of leverage is indicated by the negative value of the same coefficient.

3.2.2. The GJR-GARCH Model

The GJR-GARCH(p, q) model is another asymmetric GARCH model proposed by Glosten, Jagannatahan and Runkle [6] . The generalized form of the GJR- GARCH (p, q) model is given in the following form:

${\sigma}_{t}^{\delta}=\omega +{\displaystyle \underset{i=1}{\overset{p}{\sum}}{\alpha}_{i}{\epsilon}_{t-i}^{2}+{\displaystyle \underset{j=1}{\overset{q}{\sum}}{\beta}_{j}{\sigma}_{t-j}^{2}+{\gamma}_{i}{I}_{t-i}{\epsilon}_{t-i}^{2}}}$ (8)

${I}_{t-i}=\{\begin{array}{l}1\text{when}{\epsilon}_{t-i}<0\\ 0\text{when}{\epsilon}_{t-i}\ge 0\end{array}$

when $\omega $ , ${\alpha}_{i}$ , ${\beta}_{j}$ and ${\gamma}_{i}$ are parameters under estimation

${I}_{t-i}$ is a dummy variable, meaning that ${I}_{t-i}$ is a functional index which takes zero value when ${\epsilon}_{t-i}$ is positive and value one when ${\epsilon}_{t-i}$ is negative. If parameter ${\gamma}_{i}>0$ then negative errors are leveraged meaning that negative innovations or bad news have larger impact than good news. Finally, we assume that on the GJR-GARCH model parameters are positive and the relationship

${\alpha}_{i}+{\beta}_{j}+\frac{{\gamma}_{i}}{2}<1$ is valid.

The GJR-GARCH (1,1) model is the one that is more often used for the estimation of ${\sigma}^{2}$ variance and has the following form:

${\sigma}_{t}^{2}=\omega +{\alpha}_{1}{\epsilon}_{t-1}^{2}+{\beta}_{1}{\sigma}_{t}^{2}+{\gamma}_{1}{I}_{t-1}{\epsilon}_{t-1}^{2}$ (9)

${I}_{t-1}=\{\begin{array}{l}1\text{when}{\epsilon}_{t-1}<0\\ 0\text{when}{\epsilon}_{t-1}\ge 0\end{array}$

Engle and Ng [15] compare the response of conditional variance to shocks implied by various econometric models and find evidence that the GJR model fits stock return data the best.

3.3. Estimation of the GARCH Model

The estimation of GARCH models can be done with the Ordinary Least Squares method. Due to the fact that error terms are not independently and identically distributed iid(0,1), it is better to avoid using the OLS method mainly on small samples. In this case, it is better to use the maximum likelihood method(see Greene, [16] ). The parameters of GARCH models maximize the log likelihood function. The estimation of parameters on the log likelihood function derives through nonlinear least squares using Marquardt’s algorithm [1] . The log likelihood function is given below:

$\mathrm{ln}L\left[\left({y}_{t}\right),\theta \right]={\displaystyle \underset{t=1}{\overset{T}{\sum}}\left[\mathrm{ln}\left[D\left({z}_{t}\left(\theta \right),\upsilon \right)\right]-\frac{1}{2}\mathrm{ln}\left[{\sigma}_{t}^{2}\left(\theta \right)\right]\right]}$ (10)

where $\theta $ is the vector of the parameters that have to be estimated for the conditional mean, conditional variance and density function, ${z}_{t}$ denoting their density function, $D\left({z}_{t}\left(\theta \right),\upsilon \right)$ , is the log-likelihood function of $\left[{y}_{t}\left(\theta \right)\right]$ , for a sample of T observation. The maximum likelihood estimator $\stackrel{^}{\theta}$ for the true parameter vector is found by maximizing (10).

The models GARCH assumed Gaussian innovations, but nonetheless imply non-Gaussian unconditional distributions. However, time-varying volatility models with Gaussian innovations generally do not generate sufficient unconditional non-Gaussianity to match certain financial asset return data (see, Poon and Granger [14] ).

3.3.1. Conditional Distributions

In this section we describe the log-likelihood functions used for the estimation of parameters on volatility models for all theoretical distributions.

1) Normal Distribution

In the case of a standard normal distribution for the i.i.d. random variables $\left\{{z}_{t}\right\}$ , the following log-likelihood function needs to be maximized.

$\mathrm{ln}L\left[\left({y}_{t}\right),\theta \right]=-\frac{1}{2}\left[T\mathrm{ln}\left(2\text{\pi}\right)+{\displaystyle \underset{t=1}{\overset{T}{\sum}}{z}_{t}^{2}+{\displaystyle \underset{t=1}{\overset{T}{\sum}}\mathrm{ln}\left({\sigma}_{t}^{2}\right)}}\right]$ (11)

where $\theta $ is the vector of the parameters that have to be estimated for the conditional mean, conditional variance and density function, T is observations.

2) Student-t Distribution

The Student-t distribution can handle more severe leptokurtosis. The log- likelihood function is defined as

$\begin{array}{c}\mathrm{ln}L\left[\left({y}_{t}\right),\theta \right]=T\left[\mathrm{ln}\Gamma \left(\frac{\upsilon +1}{2}\right)-\mathrm{ln}\Gamma \left(\frac{\upsilon}{2}\right)-\frac{1}{2}\mathrm{ln}\left[\pi \left(\upsilon -2\right)\right]\right]\\ \text{}-\frac{1}{2}{\displaystyle \underset{t=1}{\overset{T}{\sum}}\left[\mathrm{ln}\left({\sigma}_{t}^{2}\right)+\left(1+\upsilon \right)\mathrm{ln}\left(1+\frac{{z}_{t}^{2}}{\upsilon -2}\right)\right]}\end{array}$ (12)

where $\Gamma \left(\upsilon \right)={\displaystyle {\int}_{0}^{\infty}{\text{e}}^{-x}{x}^{\upsilon -1}\text{d}x}$ is the gamma function and $\upsilon $ is the degree of freedom.

The t-Student is symmetric around zero. The Student-t distribution incorporates the standard normal distribution as a special case when $\upsilon =\infty $ and the Cauchy distribution when $\upsilon =1$ . Hence, a lower value, $\upsilon $ yields a distribution with “fatter tails”.

3.3.2. Generalized Error Distribution

Nelson [5] proposed the use of GED when estimating EGARCH since it is more appealing in terms of fulfilling stationarity compared tο the Student-t distribution. In the case of a t-Student distribution the unconditional means and variances may not be finite in the EGARCH. The log-likelihood function for the standard GED is defined as

$\mathrm{ln}L\left[\left({y}_{t}\right),\theta \right]={\displaystyle \underset{t=1}{\overset{T}{\sum}}\left[\mathrm{ln}\left(\frac{\upsilon}{\lambda}\right)-\frac{1}{2}{\left|\frac{{z}_{t}}{\lambda}\right|}^{\upsilon}-\left(1+{\upsilon}^{-1}\right)\mathrm{ln}\left(2\right)-\mathrm{ln}\Gamma \left(\frac{1}{\upsilon}\right)-\frac{1}{2}\mathrm{ln}\left({\sigma}_{t}^{2}\right)\right]}$ (13)

where $\lambda ={\left[{2}^{-2/\upsilon}\frac{\Gamma \left(\frac{1}{\upsilon}\right)}{\Gamma \left(\frac{3}{\upsilon}\right)}\right]}^{1/2}$

The distribution of generalized error (GED) incorporates both normal distribution when (
$\upsilon =2$ ), Laplace distribution when (
$\upsilon =1$ ), and the unique distribution for
$\upsilon =\infty $ . Specifically, we would say when
$\upsilon =2$ the distribution of the random variable z_{t} would be the standard normal distribution. When
$\upsilon <2$ , the distribution of the random variable z_{t}, will have thicker tails than that of normal distribution. For
$\upsilon =1$ the distribution of the random variable z_{t} will have a double exponential distribution. For
$\upsilon >2$ the distribution of the random variable z_{t} will have thinner tails than normal distribution, and for
$\upsilon =\infty $ the distribution of the random variable z_{t} will be a uniform distribution.

In our paper we estimate the conditional volatility using normal distribution, t-student and generalized error distribution. Engle [2] , who introduced ARCH models on the estimations of models used normal distribution. However, in the literature it is stated that the returns of assets do not follow a normal distribution thus Engle’s estimations could have been biased on these models ignoring volatility. Many authors such as Brooks, Clare and Persand [17] and Vilasuso [18] proved that estimations on GARCH models following normal distribution had the lowest forecasting performance on models that reflected skewness and kurtosis in innovations. Bollerslev [19] , in order to measure the kyrtosis on assets returns, introduced a standardized Student’s t distribution with $\upsilon >2$ degrees of freedom on GARCH models.

4. Data and Descriptive Statistics

The data in our study are collected from the official website www.nasdaqomxnordic.com. The data span is the period from 30 of September 1986 to 11 May 2016 and comprise 7434 observations. The daily stock return is calculated as:

${R}_{t}=\mathrm{ln}\left(\frac{{I}_{t}}{{I}_{t-1}}\right)\times 100=\left(\mathrm{ln}{I}_{t}-\mathrm{ln}{I}_{t-1}\right)\times 100$ (14)

where ${I}_{t}$ is the daily closing value of the stock market on day t, and ${R}_{t}$ is the daily stock return.

The daily closing values of OMX Stockholm 30 Index and its returns are displayed in Figure 1 and Figure 2, respectively.

As it can be seen in Figure 1, the closing values of OMX Stockholm 30 Index show a random walk.

As it can be seen in Figure 2, the daily returns of OMX Stockholm 30 Index are stationary. The return data is tested for autocorrelation both in returns as well as in squared returns and are displayed in Figure 3 and Figure 4, respectively.

The Ljung and Box Q-statistics [20] [21] on the 1st, 10th, 20th and 36th lags of the sample autocorrelations functions of the return series indicate significant serial correlations. When autocorrelation has been detected on data that we examine (Figure 3), we should find the autocorrelation form (ARMA (p, q)).

In addition, in Figure 4 (Squared Daily Stock Returns), the Ljung and Box Q-statistics [20] [21] on the 1st, 10th, 20th and 36th lags of the sample autocorrelations functions are statistically significant, thus indicating that there is an

Figure 1. Daily closing values of OMX stockholm 30 Index in the period from 30 September 1986 to 11 May 2016.

Figure 2. Daily stock returns of OMX stockholm 30 Index in the period from 30 September 1986 to 11 May 2016.

Figure 3. Correlogram of daily stock returns of the OMX Stockholm 30 Index.

ARCH effect. Since we have detected that there is an ARCH effect (Figure 4), we have to find the most suitable GARCH model which can adjust the data on the autocorrelation form.

The summary of the descriptive statistics for the daily logarithmic stock index returns of the OMX Stockholm 30 Index is presented in Figure 5.

The results in Figure 5 show that the daily return rates do not follow normal distribution. In other words, the returns of the Stockholm stock market present positive asymmetry and kurtosis (leptokurtic), suggesting that the return distribution is a fat-tailed one. In Figure 6 the Q-Q plot is presented, which displays the quantiles of return data series against the quantiles of the normal distribution.

Figure 4. Correlogram of squared daily stock returns of the OMX Stockholm 30 Index.

The Q-Q plot, which displays the quantiles of return data series against the quantiles of the normal distribution, shows that there is a low degree of fit of the empirical distribution to the normal distribution.

The leptokurtic behavior of the data is confirmed by the normal quintile and empirical density graph presented in Figure 7.

The summary of the descriptive statistics, normality tests, ARCH tests and unit root tests for the daily stock index returns of the OMX Stockholm 30 Index is presented in Table 1.

After the detection of series stationarity, we define the form of the ARMA (p, q) model from the correlogram of Figure 3. Parameters p and q can be determined from partial autocorrelation coefficients and autocorrelation coefficients

Figure 5. Summary descriptive statistics for the daily returns of the OMX Stockholm 30 Index.

Figure 6. Q-Q plot of daily stock returns of the OMX Stockholm 30 Index.

Figure 7. Normal density graphs of daily stock returns of the OMX Stockholm 30 Index.

Table 1. Summary descriptive statistics, normality tests, ARCH tests and unit root tests for the daily returns of the OMX Stockholm 30 Index.

respectively, compared to the critical value $\pm \frac{2}{\sqrt{n}}=\pm \frac{2}{\sqrt{7433}}=\pm 0.023$ . There-,

fore the value of p will be between $0\le p\le 3$ , and respectively, the value of q will be between $0\le q\le 3$ . Thereafter, we create Table 2 with the values of p and q as follows:

The results from Table 2 indicate that according to the criteria of Akaike (AIC), Schwartz (SIC) and Hannan-Quinn (HQ), the most suitable model is the ARIMA (0,0,1) model.

After the estimation of the above model in Figure 8, we test for the existence

Table 2. Comparison of models within the range of exploration using AIC, SIC and HQ.

Figure 8. Estimation of the ARIMA (0, 0, 1) model.

of conditional heteroscedasticity (ARCH(q) test) from the squared residuals of the above model. Figure 9 gives these results.

From the results of Figure 9 we can see that autocorrelation coefficients and partial autocorrelation coefficients are statistically significant. Consequently, the null hypothesis for the absence of ARCH or GARCH procedure is rejected.

5. Empirical Results

Since there are ARCH effects in the Stockholm stock return data, we can proceed with the estimation of the ARIMA(0,0,1)-GARCH models.

Figure 9. Q-statistics for the standardized squared residuals.

First of all, we estimate the symmetric ARIMA(0,0,1)-GARCH(1,1) model with normal distribution, t-student distribution as well as the Generalized error distribution (GED). The estimation of parameters is done with the maximum likelihood method using the Broyden-Fletcher-Goldfarb-Shanno (BFGS) (see Press et al. [22] ) algorithm which is a repeating method for solving non-linear optimization problems without constraint. The parameters (coefficients) of estimated models and the residuals’ test of normality, autocorrelation and conditional heteroskedasticity are provided in Table 3. A higher log-likelihood value yields a better fit.

Table 3 gives both the estimation of parameters along with the value of the log-likelihood function as well as the residual tests of normality, autocorrelation and conditional heteroskedasticity. From the above table we point out that coefficients are statistically significant with all distributions. Also, there is no autocorrelation and conditional heteroskedasticity problem. Furthermore, the ARIMA(0,0,1)-GARCH(1,1) model has the largest value in logarithmic likelihood (LL) with t-student distribution. Thus, we can use this model for forecasting.

We proceed with the following asymmetric (non linear) GARCH models such as the ARIMA(0,0,1)-EGARCH (1,1) model as well as the ARIMA(0,0,1)-GJR- GARCH(1,1) model with the normal distribution, t-student distribution and Generalized error distribution (GED). The estimation of parameters is done with the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm Marquardt [1] . The parameters (coefficients) of estimated models and the residuals test of normality, autocorrelation and conditional heteroskedasticity are provided in Table 4.

From the results of Table 4 we can see that all the coefficients of non-linear GARCH models are statistically significant. These results show that asymmetry exists. Furthermore, diagnostics tests of non-linear GARCH models seem to be

Table 3. Estimated symmetric GARCH models for the daily returns of the OMX Stockholm 30 Index.

Notes: 1. The persistence is calculated as α_{1} + β_{1} for the ARMA(0,0,1)-GARCH(1,1) model. 2. Values in parentheses denote the p-values. 3. LL is the value of the log-likelihood.

Table 4. Estimated asymmetric GARCH models for the Daily Returns of the OMX Stockholm 30 Index.

Notes: 1. The persistence is calculated as ${\beta}_{\text{1}}$ for ARIMA(0,0,1)-EGARCH(1,1) model, and ${\alpha}_{1}+{\gamma}_{1}/2+{\beta}_{1}$ for ARIMA(0,0,1)-GJR-GARCH(1,1) model. 2. Values in parentheses denote the p-values. 3. LL is the value of the log-likelihood.

satisfactory. Also, the results from the models show that Q-statistics for the standardized squared residuals and the ARCH-LM test are insignificant with high p values. From the above table we can see that ARIMA(0,0,1)-GJR- GARCH(1,1) model has the largest logarithmic likelihood (LL) value with t-student distribution. Thus, this model can be used for forecasting.

5.1. Asymmetric and Leverage Effects

Asymmetry and leverage effects results are examined for non-linear variances of ARIMA(0,0,1)-EGARCH(1,1) and ARIMA(0,0,1)-GJR-GARCH(1,1) models from three different distributions. Since the coefficients are statistically significant in all cases, asymmetry exists. Positive signs of the coefficients on the ARIMA(0,0,1)-GJR-GARCH(1,1) models as well as negative signs on the ARIMA(0,0,1)-EGARCH(1,1) models indicate that there are leverage effects. In addition, bad news has more impact on volatility than good news in all distributions that we used. In the following Table 5, we present the models with the three distributions indicating that bad news have more impact on volatility. For example, on the ARIMA(0,0,1)-GJR-GARCH(1,1) model and t-student the effect of bad news on conditional volatility is 6.03 times higher than good news.

5.2. Test of Asymmetries

In order to examine if an asymmetric model is suitable for forecasting, Engle and Ng [15] created a test known as the sign and size bias test defining if an asymmetric model is suitable for the examined series or to what extent the symmetric GARCH model is considered adequate. The Engel-Ng [15] test is usually applied to the residuals of a GARCH fit to the returns data. The sign and size bias test is based on the significance of b_{1} coefficient of the following regression:

${\stackrel{^}{\epsilon}}_{t}^{2}={b}_{0}+{b}_{1}{D}_{i,t-1}^{-}+{v}_{t}$ (15)

where

${\stackrel{^}{\epsilon}}_{t}^{2}$ are the squared residuals from the symmetric GARCH model.

${D}_{i,t-1}^{-}$ is a dummy variable which takes the value 1 if ${\epsilon}_{t-i}$ is negative and 0 otherwise, and gives the slope dummy value.

${v}_{t}$ is an i.i.d. error term. (see Dutta [14] ).

If b_{1} coefficient is statistically significant in positive and negative changes relatively to conditional variance then there is asymmetry on the GARCH model.

A test for sign bias can also be conducted using the following regression:

${\stackrel{^}{\epsilon}}_{t}^{2}={b}_{0}+{b}_{1}{D}_{t-1}^{-}{\epsilon}_{t-1}+{v}_{t}$ (16)

Like regression (15), the statistical significance of b_{1} coefficienton regression (16) indicates that the size of a shock will have an asymmetric impact on volatility. Regression (16) tests the negative bias size. For the positive bias size we use the following regression: (see Dutta [14] )

${\stackrel{^}{\epsilon}}_{t}^{2}={b}_{0}+{b}_{1}\left(1-{D}_{t-1}^{-}\right){\epsilon}_{t-1}+{v}_{t}$ (17)

A joint test can be conducted through defining ${D}_{t-1}^{+}$ as $1-{D}_{t-1}^{-}$ which indicates a positive size bias. The joint test for positive sign bias and positive or negative size bias is presented on the following regression:

Table 5. The magnitude of news impact on volatility.

Notes: The asymmetry is calculated as
$\text{1}-{\gamma}_{\text{1}}$ , and
$\text{1}+{\gamma}_{\text{1}}$ for the ARIMA(0,0,1)-EGARCH(1,1) model,
${\alpha}_{\text{1}}+{\gamma}_{\text{1}}$ _{ , }and
${\alpha}_{\text{1}}$ for the ARIMA(0,0,1)-GJR-GARCH(1,1) model.

${\stackrel{^}{\epsilon}}_{t}^{2}={b}_{0}+{b}_{1}{D}_{t-1}^{-}+{b}_{2}{D}_{t-1}^{-}{\epsilon}_{t-1}+{b}_{3}{D}_{t-1}^{+}{\epsilon}_{t-1}+{v}_{t}$ (18)

The significance of b_{1} coefficient of regression (17) shows the existence of sign bias where positive and negative changes have different consequences in volatility compared to the symmetric GARCH model. On the other hand, the significance of b_{2} and b_{3} coefficients of regression (18) indicates not only size bias but also if the size of change is significant. The test follows a χ^{2} distribution with degrees of freedom equal to 3. The joint test statistic is given from the formula TR^{2}. The null hypothesis for the joint test is that there is no asymmetric result. (see Brooks, [23] , pp. 474-475). Table 6 presents the results of the asymmetry and volatility tests.

The results of Table 6 show that the sign bias test is statistically significant on both models. Thus, there is asymmetry. This result is also confirmed from two size bias tests having large statistical significance. Also, from the results of the above table we can see that the size effect of bad news is stronger than that of good news.

5.3. Likelihood Ratio Tests

The Likelihood ratio (LR) tests consist of estimations on two models, (an unrestricted model and a restricted one). The null hypothesis that is examined is ${H}_{0}:{\gamma}_{1}=0$ . The maximized values of log-likelihood function (LLF) are used on this test according to the following:

$LR=-2\left(LL{F}_{r}-LL{F}_{u}\right)\to {\chi}^{2}\left(m\right)$ (19)

where

$LL{F}_{r}$ is the value of maximum likelihood function from the constrained model

$LL{F}_{u}$ is the value of maximum likelihood function from the unconstrained model

$m$ the number of constraints.

The Likelihood ratio (LR) tests follow asymptotically χ^{2} distribution with m degrees of freedom. Table 7 presents the maximized values from all estimated models with the corresponding distributions.

Table 6. Tests of asymmetries.

Notes: 1. (*) denotes significance at the (1%); 2. The t-statistics for the sign bias, negative size bias and positive size bias tests are those of coefficient b_{1} in regression (15), (16) and (17), respectively. 3. The F-statistic is basedon regression (18). 4. Values in parentheses denote the p-values.

Table 7. Likelihood-ratio test results.

Note: (*) denotes significance at the (1%) level, LLF_{r} is the value from maximum likelihood function from the constrained model, LLF_{u} is the value from maximum likelihood function from the unconstrained model and LR is Likelihood ratiotest.

The results from Table 7 show that null hypothesis is rejected so asymmetry exists. Thus, we can use both models of forecasting due to existence of asymmetry (Engle and Ng, and LR test). Furthermore, the negative value of ${\gamma}_{\text{1}}$ coefficient on the ARIMA(0,0,1)-EGARCH(1,1) model and the positive value of ${\gamma}_{\text{1}}$ coefficient on the ARIMA(0,0,1)-GJR-GARCH(1,1) model denote the presence of leverage.

6. Forecasting the Volatility of the Stockholm Stock Exchange Index Returns ex Post

In this section we present the forecasting results from the two asymmetric models. In our paper we forecast for future values on rate of return and volatility on the Stockholm stock market using the static1-step ahead method based on estimated parameters of the two asymmetric models. The last 434 series observations were used for an ex-post forecast, with the main focus on the forecast of volatility. In the literature, a variety of statistics has been used which evaluates and compares the forecasts of returns. The optimal forecasting value is evaluated through Mean Squared Error. Other statistical indices usually used for the return of forecasting are the Mean Absolute Error (MAE), the Root Mean Square Error (RMSE), the Mean Absolute Percentage Error (MAPE) and the (Theil U-Theil index [24] ). The lower the values of the RMSE, MAE and MAPE indices, the better the forecast of models. Figure 10 and Figure 11 of the two models for the period we analyze are given below:

The above figures show that estimation intervals are stable on both models. However, there are some indices that help us to test the forecast of a model. The first one is referred to as the Mean Square Error. According to Patton [25] this criterion is the most powerful for the evaluation of models. However, the mean square error has some constraints on the forecasting of variance. According to Vilhelmsson [26] , this index is sensitive to outliers. So, Vilhelmsson suggests the

Figure 10. Ex-post forecast of the volatility of the Stockholm Stock Exchange index returns (ARIMA(0,0,1)-EGARCH(1,1).

Figure 11. Ex-post forecast of the volatility of the Stockholm Stock Exchange index returns (ARIMA(0,0,1)-GJR-GARCH(1,1).

mean absolute error index which is more robust to outliers. According to the Mean Absolute Error (MAE), the Root Mean Square Error(RMSE) and Mean Absolute Percentage Error (MAPE), the ARIMA(0,0,1)-EGARCH(1,1) model provides more exact forecasts on the returns of the Stockholm stock market.

7. Conclusions

The modeling and forecasting of volatility in financial markets used to be a fundamental issue for many researchers. The importance of this problem increased over the last years as there is upheaval in the financial world. The aim of this paper is to compare various volatility models and forecast for future values regarding the rate of return and volatility of the Stockholm stock market using 1-step ahead. Measuring the period from 30 September 1986 until 11 May 2016 and using as a sample 7,434 daily observations for different models we concluded that the asymmetric models give better results on the returns and volatility of the Stockholm stock market. More specifically, we estimated the symmetric ARIMA(0,0,1)-GARCH(1,1) model, as well as the asymmetric models ARIMA(0,0,1)-EGARCH(1,1) and ARIMA(0,0,1)-GJR-GARCH(1,1) models with different residual distributions. The analysis of estimations indicated that t-student distribution is considered the most suitable on the estimation of parameters for all models. These results are in accordance with the empirical works by (Hamilton and Susmel [27] ), (Poon and Granger [28] ) and many others documented conditional non-Gaussianities in financial data. Moreover, an effort was made to test the asymmetric response of volatility on the positive and negative changes of models. The results of our paper, as far as the asymmetric tests are concerned, showed that negative changes are more frequent and more powerful (robust) on the returns of the Stockholm stock market.

To sum up, the results of our paper confirm previous findings that GARCH models with normal errors do not seem to fully capture the leptokurtosis in empirical time series (see, e.g. Kim and White [29] ). Contrary to t-student distribution and GED, which provide a better frame on conditional volatility, we can better test the time-varying heteroskedasticity, skewness and kurtosis of the series. This is accomplished by allowing the parameters of the two distributions to vary through time. Finally, the asymmetric models appear to better formulate the different responses to different past shocks and to explain conditional volatility.

References

[1] Marquardt, D.W. (1963) An Algorithm for Least Squares Estimation of Nonlinear Parameters. Journal of the Society for Industrial and Applied Mathematics, 11, 431-441.

https://doi.org/10.1137/0111030

[2] Engle, R.F. (1982) Autoregressive Conditional Heteroscedasticity with Estimates of the Variance of UK Inflation. Econometrica, 50, 987-1008.

https://doi.org/10.2307/1912773

[3] Bollerslev, T. (1986) Generalized Autoregressive Conditional Heteroskedasticity. Journal of Econometrics, 31, 307-327.

https://doi.org/10.1016/0304-4076(86)90063-1

[4] Black, F. (1976) Studies in Stock Price Volatility Changes of the Nominal Excess Return on Stocks. Proceedings of the American Statistical Association, Business and Economics Statistics Section, 177-181.

[5] Nelson, D.B. (1991) Conditional Heteroskedasticity in Asset Returns: A New Approach. Econometrica, 59, 347-370.

https://doi.org/10.2307/2938260

[6] Glosten, L.R., Jagannathan, R. and Runkle, D.E. (1993) On the Relation between the Expected Value and the Volatility of the Nominal Excess Return on Stocks. The Journal of Finance, 48, 1779-1801.

https://doi.org/10.1111/j.1540-6261.1993.tb05128.x

[7] Donaldson, R.G. and Kamstra, M. (1997) An Artificial Neural Network—GARCH Model for International Stock Return Volatility. Journal of Empirical Finance, 4, 17-46.

https://doi.org/10.1016/s0927-5398(96)00011-4

[8] Nam, K., Pyun, C.S. and Arize, C.A. (2002) Asymmetric Mean-Reversion and Contrarian Profits: ANSTGARCH Approach. Journal of Empirical Finance, 9, 563-588.

https://doi.org/10.1016/S0927-5398(02)00011-7

[9] Tudor, C. (2008) Modeling Time Series Volatilities Using Symmetrical GARCH Models. The Romanian Economic Journal, 30, 183-208.

[10] Panait, I. and Slavescu, E.O. (2012) Using Garch-in-Mean Model to Investigate Volatility and Persistence at Different Frequencies for Bucharest Stock Exchange during 1997-2012. Theoretical and Applied Economics, 19, 55-76.

[11] Gao, Y., Zhang, C. and Zhang, L. (2012) Comparison of GARCH Models Based on Different Distributions. Journal of Computers, 7, 1967-1973.

https://doi.org/10.4304/jcp.7.8.1967-1973

[12] Dutta, A. (2014) Modelling Volatility: Symmetric or Asymmetric GARCH Models? Journal of Statistics: Advances in Theory and Applications, 12, 99-108.

[13] Bollerslev, T., Chou, R.Y. and Kroner, K.F. (1992) ARCH Modeling in Finance: A Review of the Theory and Empirical Evidence. Journal of Econometrics, 52, 5-59.

https://doi.org/10.1016/0304-4076(92)90064-X

[14] Poon, S.H. and Granger C.W.J. (2003) Forecasting Volatility in Financial Markets: A Review. Journal of Economic Literature, 41, 478-539.

https://doi.org/10.1257/.41.2.478

[15] Engle, R.F. and Ng, V.K. (1993) Measuring and Testing the Impact of News on Volatility. Journal of Finance, 48, 1749-1778.

https://doi.org/10.1111/j.1540-6261.1993.tb05127.x

[16] Greene, W.H. (2012) Econometric Analysis. 7th Edition, Prentice Hall, Upper Saddle River.

[17] Brooks, C., Clare, A.D. and Persand G. (2000) A Word of Caution on Calculating Market Based Minimum Capital Risk Requirements. Journal of Banking and Finance, 24, 1557-1574.

https://doi.org/10.1016/S0378-4266(99)00092-8

[18] Vilasuso, J. (2002) Forecasting Exchange Rate Volatility, Economics Letters, 76, 59-64.

https://doi.org/10.1016/S0165-1765(02)00036-8

[19] Bollerslev, T. (1987) Conditionally Heteroskedastic Time Series Model for Speculative Prices and Rates of Return. Review of Economics and Statistics, 69, 542-547.

https://doi.org/10.2307/1925546

[20] Ljung, G. and Box, G. (1978) On a Measure of Lack of Fit in Time Series Models. Biometrika, 65, 297-303.

https://doi.org/10.1093/biomet/65.2.297

[21] Ljung, G. and Box, G. (1979) The Likelihood Function of Stationary Autoregressive-Moving Average Models. Biometrika, 66, 265-270.

https://doi.org/10.1093/biomet/66.2.265

[22] Press, W., Flannery, B., Teukolsky, S. and Vettering, W. (1988) Numerical Recipes in C. Cambridge University Press, New York.

[23] Brooks, C. (2008) Introductory Econometrics for Finance. 2nd Edition, Cambridge University Press, New York.

https://doi.org/10.1017/CBO9780511841644

[24] Theil, H. (1961) Economic Forecasts and Policy. North-Holland Publishing Company, Amsterdam.

[25] Patton, A. (2006) Volatility Forecast Comparison Using Imperfect Volatility Proxies. University of Technology, Sydney.

[26] Vilhelmsson, A. (2006) Garch Forecasting Performance under Different Distribution Assumptions. Journal of Forecasting, 25, 561-578.

https://doi.org/10.1002/for.1009

[27] Hamilton, J.D. and Susmelb, R. (1994) Autoregressive Conditional Heteroskedasticity and Changes in Regime. Journal of Econometrics, 64, 307-333.

https://doi.org/10.1016/0304-4076(94)90067-1

[28] Poon, S.H. and Granger, C. (2003) Forecasting Volatility in Financial Markets: A Review. Journal of Economic Literature, 41, 478-539.

https://doi.org/10.1257/.41.2.478

[29] Kim, T.H. and White, H. (2004) On More Robust Estimation of Skewness and Kurtosis. Finance Research Letters, 1, 56-73.

https://doi.org/10.1016/S1544-6123(03)00003-5