Next Article in Journal
Heteroskedasticity in One-Way Error Component Probit Models
Next Article in Special Issue
Simultaneous Indirect Inference, Impulse Responses and ARMA Models
Previous Article in Journal
A Comparison of Some Bayesian and Classical Procedures for Simultaneous Equation Models with Weak Instruments
Previous Article in Special Issue
Monte Carlo Inference on Two-Sided Matching Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Multi-Step-Ahead Prediction of ARCH/GARCH Models and NoVaS Transformation

Department of Mathematics, University of California, San Diego, CA 92093, USA
*
Author to whom correspondence should be addressed.
Submission received: 31 October 2018 / Revised: 9 July 2019 / Accepted: 21 July 2019 / Published: 8 August 2019
(This article belongs to the Special Issue Resampling Methods in Econometrics)

Abstract

:
This paper gives a computer-intensive approach to multi-step-ahead prediction of volatility in financial returns series under an ARCH/GARCH model and also under a model-free setting, namely employing the NoVaS transformation. Our model-based approach only assumes i . i . d innovations without requiring knowledge/assumption of the error distribution and is computationally straightforward. The model-free approach is formally quite similar, albeit a GARCH model is not assumed. We conducted a number of simulations to show that the proposed approach works well for both point prediction (under L 1 and/or L 2 measures) and prediction intervals that were constructed using bootstrapping. The performance of GARCH models and the model-free approach for multi-step ahead prediction was also compared under different data generating processes.

1. Introduction

Multi-step-ahead prediction in a time series amounts to predicting a sequence of future values using only the values observed over a finite time interval. Examples of time series for which multi-step-ahead prediction is useful include crop yields, stock prices, traffic volume, and electrical power consumption. In the paper at hand, we focus on multi-step-ahead prediction of squared financial returns, which is related to the so-called volatility, i.e., the conditional expectation of the squared returns. A typical approach to solve this problem, known as multi-stage or iterated prediction, is to construct a single model from the past observed time series data and then apply the model step by step to predict its future values. The iterated method uses the predicted value of the current time step to determine its value in the next time step. However, empirical evidence points to the fact that multi-stage prediction is susceptible to the error accumulation problem, i.e., errors committed in the past are propagated into future predictions.
The benchmark model for financial returns has been the GARCH(1,1); see Bollerslev et al. (1992) and the references therein. Notably, the work in Andersen et al. (2006), page 811, stated that, beyond the one-step-ahead case, there is no analytical form of the predictive density for multi-step-ahead predictions of volatility in GARCH models. However, the analytical form of the multi-step ahead predictive probability density function of a GARCH(1,1) process under Gaussian or Student t innovations has been recently derived in the working paper of Abadir et al. (2018).
The work of Abadir et al. (2018) was an important breakthrough, but it hinged on knowing the error distribution that drove the GARCH(1,1) model. By contrast, the paper at hand proposes an alternative computer-intensive approach to multi-step-ahead prediction of the squared returns in ARCH/GARCH and related models that does not require knowledge of the error distribution. Furthermore, an analogous methodology can be employed in a model-free setting using the normalizing and variance-stabilizing transformation (NoVaS) approach; for more details on NoVaS, see Politis (2003, 2007, 2015). Notably, our method is not of the multi-stage/iterated type, and therefore eliminates the errors accumulation issue.
The remainder of the paper is organized as follows. Section 2 presents the proposed method for the optimal multi-step-ahead point predictions for ARCH/GARCH processes, as well as in a model-free setting using the NoVaS transformation. Section 3 addresses the associated methods for the construction of multi-step-ahead prediction intervals using bootstrapping. Section 4 illustrates the numerical performance of the proposed methods by means of simulated examples; some concluding remarks are provided in Section 5.

2. Optimal Multi-Step-Ahead Point Prediction

Consider data X 1 , , X n from a zero mean and (strictly) stationary financial returns time series { X t } . Our goal is to predict the future squared returns X n + h 2 for any h 2 ; the case h = 1 was treated in Politis (2007, 2015).
Let F n be a short-hand for the observed information set, i.e., F n = { X t , 1 t n } . In the L 2 sense, the optimal predictor of X n + h 2 based on F n is the conditional mean and given by:
X n + h 2 ^ = E ( X n + h 2 | F n ) .
Similarly, the optimal L 1 predictor is the conditional median:
X n + h 2 ^ = M e d i a n ( X n + h 2 | F n ) .
In the following parts of this section, we study the multi-step-ahead prediction in the nonlinear financial models ARCH/GARCH and in the NoVaS setting, which is an application of the model-free approach to financial returns data.

2.1. L 2 Optimal Prediction for ARCH(p) and GARCH(1,1) Models

Suppose the data follow the ARCH(p) process of Engle (1982) defined by the recursion:
X t = σ t ϵ t , a n d σ t 2 = α + a 1 X t 1 2 + + a p X t p 2
where α 0 , a j 0 for all j = 1 , , p , and { ϵ t } i . i . d . N ( 0 , 1 ) .
First, consider the simplest case h = 2 . Based on ARCH(P) Model in (3), we can express X n + 1 and X n + 2 in the following way:
X n + 1 = ϵ n + 1 σ n + 1 2 , a n d σ n + 1 2 = α + a 1 X n 2 + + a p X n p + 1 2 ,
X n + 2 = ϵ n + 2 σ n + 2 2 , a n d σ n + 2 2 = α + a 1 X n + 1 2 + + a p X n p + 2 2 .
Obviously, X n + 1 can be easily written as a function of the past observations X n , , X n + 1 p and the unknown future error ϵ n + 1 . Furthermore, we can also rewrite X n + 2 to be a function of X n , , X n + 1 p and the unknown future errors ϵ n + 1 and ϵ n + 2 . The notations are as follows:
X n + 1 = ϵ n + 1 α + a 1 X n + a 2 X n 1 2 + + a p X n p + 1 2 = f 1 ( X n , , X n p + 1 ; ϵ n + 1 )
and:
X n + 2 = ϵ n + 2 α + a 1 X n + 1 2 + a 2 X n 2 + + a p X n p + 2 2 = ϵ n + 2 α + a 1 ϵ n + 1 2 ( α + a 1 X n + a 2 X n 1 2 + + a p X n p + 1 2 ) + a 2 X n 2 + + a p X n p + 2 2 = f 2 ( X n , , X n p + 1 ; ϵ n + 1 , ϵ n + 2 ) .
Recursively, we can express X n + h for any h 1 as a function of past observations { X 1 , , X n } and the unknown future innovations { ϵ n + 1 , , ϵ n + h } in the form:
X n + h = f h ( X 1 , , X n ; ϵ n + 1 , , ϵ n + h ) .
Since { X 1 , , X n } are given and known, we can write (6) simply as:
X n + h = f h ( ϵ n + 1 , , ϵ n + h ) , for any h 1 .
The squared financial returns can be rewritten as f h 2 . Based on the assumption that ϵ t is i . i . d N ( 0 , 1 ) , the conditional distribution function F f h 2 of the future squared returns f h 2 ( · ) can be derived. Hence, the optimal predictor (conditional median for L 1 or conditional mean for L 2 ) of x n + h 2 is easy to calculate by F f h 2 .
Take h = 1 and h = 2 as examples. By (4) and (5), the L 2 optimal predictors of X n + 1 2 and X n + 2 2 are:
X n + 1 2 ^ = E { ϵ n + 1 2 ( α + a 1 X n + a 2 X n 1 2 + + a p X n p + 1 2 ) | F n } = α + a 1 X n 2 + a 2 X n 1 2 + + a p X n p + 1 2 ,
and:
X n + 2 2 ^ = E [ ϵ n + 2 2 ( α + a 1 σ n + 1 2 ϵ n + 1 2 + a 2 X n 2 + + a p X n p + 2 2 ) | X 1 , , X n ] = α + a 1 σ n + 1 2 + a 2 X n 2 + + a p X n p + 2 2
since E ( ϵ n + 1 2 | F n ) = 1 and E ( ϵ n + 2 2 | F n ) = 1 by assumption. First, we can note that X n + 1 2 ^ = f 1 2 ( ϵ n + 1 2 = 1 ) and X n + 2 2 ^ = f 2 2 ( ϵ n + 1 2 = 1 , ϵ n + 2 2 = 1 ) . Actually, we can easily verify that the L 2 optimal predictor for any h 1 is given by:
X n + h 2 ^ = f h 2 ( ϵ n + 1 2 = 1 , , ϵ n + h 2 = 1 ) .
Note that because all ϵ t ’s are independent of each other, as well as of the past values of the X series, the h-step-ahead predictor in (10) is equivalent to the method of multi-stage/iterated prediction that uses the predicted values of the current time step to determine its value in the next time step. However, for the L 1 case, because the median function is not a linear operator, this equivalence breaks down; see Section 2.3 in what follows.
Remark 1.
As already mentioned, the benchmark for fitting financial returns is the GARCH(1,1). Nevertheless, a GARCH(1,1) model is tantamount to an ARCH(p) model with p = and an exponentially-decreasing coefficient; see, e.g., Francq and Zakoian (2011). Because of the exponentially decrease of the ARCH coefficients, it is customary to approximate the GARCH(1,1) models with an ARCH(p) where p is finite, albeit large. In this sense, all the above results apply verbatim to a GARCH(1,1) process as well. Of course, in fitting a GARCH(1,1) model, the GARCH equation is used to fit four parameters, which are then expanded to the p coefficients of the approximating ARCH(p) model.

2.2. L 2 Optimal Prediction for NoVaS

Given a sequence of observations { X 1 , , X n } , we can fit the data by a special application of the model-free methodology, NoVaS, which was introduced by Politis (2003, 2007) for stationary data in prediction of squared financial returns. Let us continue considering a zero mean and (strictly) stationary financial return time series { X t } . The NoVaS methodology is trying to map the dataset X 1 , , X n to an i . i . d Gaussian dataset { W t , t n } .
The starting point is the ARCH model defined by:
X t = Z t a + i = 1 p a i X t i 2
under which, the residual:
X t a + i = 1 p a i X t i 2
is thought of as perfectly normalized and variance-stabilized, as it is assumed to be i . i . d . N ( 0 , 1 ) , which is actually not true here. This ratio can be interpreted as an attempt to “Studentize” the return X t by dividing with a time-localized measure of the standard deviation of X t . However, there seems to be no reason to exclude the value of X t from an empirical, causal estimate of the standard deviation of X t ; recall that a causal estimate is one involving present and past data only, i.e., the data { X s , s t } .
Hence, the work in Politis (2003) defined a new “Studentized” quantity as follows:
W t : = X t α s t 1 2 + a 0 X t 2 + i = 1 p a i X t i 2 for t = p + 1 , p + 2 , , n .
In the above, s t 1 2 is an estimator of σ X 2 = V a r ( X 1 ) based on the data up to (but not including1) time t; under the zero mean assumption for X 1 , the natural estimator is s t 1 2 = ( t 1 ) 1 k = 1 t 1 X k 2 .
The definition in Equation (13) describes the proposed normalizing and variance-stabilizing transformation under which the data series { X t } is mapped to the new series { W t } . The order p ( 0 ) and the vector of nonnegative parameters ( α , a 0 , , a p ) are chosen by the practitioner with the twin goals of normalization and variance stabilization.
Furthermore, the NoVaS transformation Equation (13) can be re-arranged to yield:
X t = W t α s t 1 2 + a 0 X t 2 + i = 1 p a i X t i 2 .
Formally, the only real difference between the NoVaS of Equation (14) and the ARCH of Equation (11) is the presence of the term X t 2 paired with the coefficient a 0 . Replacing the term a in Equation (11) by the term α s t 1 2 in Equation (14) is only natural since the former has, by necessity, units of variance; in other words, the term a in Equation (11) is not scale invariant, whereas the term α in Equation (14) is.
Given the assumed structure of the return series, the target of variance stabilization, which amounts to constructing a local estimator of scale for Studentization purposes, requires:
α 0 , a i 0 for all i 0 , and α + i = 0 p a i = 1 .
Equation (15) has the interesting implication that the { W t } series can be assumed to have a (unconditional) variance that is (approximately) unity. Nevertheless, note that p and α , a 0 , , a p must be carefully chosen to achieve a degree of conditional homoscedasticity as well; to do this, one must necessarily take p small enough, as well as α small enough or even equal to zero, so that a local (as opposed to global) estimator of scale is obtained. The work in Politis (2003) provided two structures for the a i coefficients satisfying Equation (15). One is to let α = 0 and a i = 1 / ( p + 1 ) for all 0 i p ; this specification is called the simple NoVaS transformation and involves only one parameter, namely the order p, to be chosen by the practitioner. The other one is given by the exponential (decay) NoVaS, where α = 0 and a i = c e c i for all 0 i p . The exponential scheme involves choosing two parameters: p and c > 0 , since c is determined by Equation (15). For more details of how to select the optimal parameters here, see Politis (2015).
The above Equation (14) can be used for one-step-ahead prediction in an analogous way to the ARCH/GARCH models already discussed. In fact, Equation (14) is formally analogous to an ARCH(p) model with i . i . d . errors given by ϵ t = W t 1 a 0 W t 2 . Hence, the construction of the previous subsection can be repeated to write X n + h as some function of { X 1 , , X n } and { W t , t = 1 , , h } for any h 1 , i.e.,
X n + h = f h ( X 1 , , X n ; W n + 1 , , W n + h ) .
Since the data { X 1 , , X n } are given, we can simplify (16) as:
X n + h = f h ( W n + 1 , , W n + h ) .
In the L 2 sense, the optimal predictor of X n + h 2 based on F n is given by:
X n + h 2 ^ = E ( X n + h 2 | F n ) = E { f h 2 ( W n + 1 , , W n + h ) | F n ) }
Since the W t are i . i . d , we can get analogous results with those concerning the ARCH/GARCH models, i.e.,
X n + h 2 ^ = f h 2 ( W n + 1 2 = 1 , , W n + h 2 = 1 ) )
Therefore, for any h 1 , we can use similar ideas to that in the ARCH/GARCH cases to conduct multi-step-ahead prediction in NoVaS by approximating the conditional mean or median from their conditional distribution functions.

2.3. L 1 Optimal Prediction and Generalizations

We can generalize the above prediction method to an interesting class of prediction functions g ( · ) , namely the power family where g ( x ) = x k for some fixed k and the power-absolute value family where g ( x ) = | x | k . So far, we have worked on the prediction of X n + h 2 , that is g ( x ) = x 2 . More generally, we can derive the best L 2 or L 1 predictor of g ( X n + h ) given F n .
Regarding L 1 optimal prediction, it was already mentioned that it is n o t equivalent to the multi-stage/iterated approach. For the h = 2 case, we can easily get an analytic formula of the conditional distribution of g ( X n + h ) assuming an ARCH(p) model with normal errors. An approximate analytic form of the conditional distribution function of X n + 2 2 is given as follows:
F f 2 2 ( x | { X n , , X 1 } ) = γ ( 1 2 , x 2 A ) B π , x > 0
where:
A = α + a 2 X n 2 + a 3 X n 1 2 + + a p X n p + 2 2 ,
B = a 1 σ n + 1 2 = a 1 ( α + a 1 X n 2 + + a p X n p + 1 2 ) ,
γ ( s , x ) = 0 x t s 1 e t d t , s > 0 .
Solving F f h 2 ( x | { X n , , X 1 } ) = 1 2 , we obtain X n + 2 2 ^ π 8 A B 2 + π 4 B 4 for the L 1 optimal predictor.
The general case h 1 has been recently worked out by Abadir et al. (2018). However, the analytical method crucially depends on the assumption for the error distribution. A more robust way to approximate the conditional distribution of g ( X n + h ) can be derived using bootstrapping2.
Before going into bootstrapping, note that a simple Monte Carlo simulation can re-produce the analytical calculations in a straightforward manner. For example, if the ARCH/GARCH errors are assumed independent with a particular distribution function F ϵ , we can easily generate pseudo-replicates of X n + h by using Equation (7), and simulating many sets of { ϵ n + 1 , , ϵ n + h } where each ϵ t is drawn i . i . d from F ϵ .
If F ϵ has a known structural form with some unknown parameters, e.g., Student-t with unknown degrees of freedom, one can use a data-based estimate of the unknown parameter in order to estimate F ϵ and then proceed with the simulation; this is then equivalent to a parametric bootstrap procedure. In the more realistic case where F ϵ is treated as unknown, it can be estimated by the empirical distribution of the ARCH/GARCH residuals and then used in the simulation; this is equivalent to the standard (nonparametric) bootstrap. Because of the formal analogy of NoVaS to ARCH models, a similar bootstrap method works for NoVaS as well; this is the so-called model-free bootstrap. Detailed Algorithms 1–5 are given in the following two subsections.

2.4. Bootstrap Algorithms for ARCH/GARCH Point Prediction

Assume an ARCH/GARCH model with all parameters known and errors { ϵ t } i . i . d . with known distribution F ϵ . Under the independence of { ϵ t } for all t 1 , we can generate many ϵ n + 1 * , , ϵ n + h * i . i . d . from F ϵ by Monte Carlo simulation and compute many pseudo-values of the quantity of interest g ( X n + h * ) .
Algorithm 1:h-step ahead prediction with parameters known: Monte Carlo.
Step 1.
Generate { ϵ n + 1 * , , ϵ n + h * } i . i . d . from F ϵ , and plug them into the function (7) to obtain the pseudo-value X n + h * . Repeat the above procedure M times, and denote the M pseudo-values by  { X n + h ( 1 ) , . , X n + h ( M ) } .
Step 2.
Calculate the optimal predictor g ( X n + h ) ^ of g ( X n + h ) by taking the sample median (under L 1 risk) or sample mean (under L 2 risk) of the set { g ( X n + h ( 1 ) ) , , g ( X n + h ( M ) ) } .
If the parameters and F ϵ in the ARCH model are unknown, estimates must be used, in which case the Monte Carlo simulation becomes bootstrapping. Let F ^ ϵ denote the estimator of F ϵ . The two cases, parametric and nonparametric bootstrap, depend on whether the parametric form of F ϵ is known or not. In the former case, F ^ ϵ uses the parametric form of F ϵ with parameters estimated and plugged in. In the latter, F ^ ϵ is typically taken to be the empirical distribution of the ARCH residuals normalized to unit variance.
Algorithm 2:h-step ahead prediction with parameters unknown: bootstrap.
Step 1.
Fit the data with an ARCH(p) or GARCH(1,1) model, and obtain the estimators { a ^ 0 , a ^ 1 , , a ^ p }  of  { a 0 , a 1 , , a p } . Furthermore, record the residuals { ϵ ^ 1 , , ϵ ^ n } with the distribution function F ^ ϵ . We will use F ^ ϵ to estimate F ϵ in the following steps.
Step 2.
Perform Step 1 and Step 2 of Algorithm 1 using F ^ ϵ instead F ϵ and { a ^ 0 , a ^ 1 , , a ^ p } instead of { a 0 , a 1 , , a p } .
Remark 2.
The work in Bose and Mukherjee (2009) proposed a weighted linear estimator (WLE) to estimate the ARCH parameters. This method does not involve nonlinear optimization and gives a closed-form expression, so it is computationally easier to obtain the estimator compared to maximum likelihood. In our numerical work, we used the WLE to obtain the estimators { a ^ 0 , a ^ 1 , , a ^ p } of { a 0 , a 1 , , a p } in Algorithm 2.

2.5. Bootstrap Algorithms for NoVaS-Based Point Prediction

We now go back to the model-free setting of Section 2.2. In order to estimate the conditional mean or conditional median in the NoVaS setting, we should first use one of the NoVaS methods (simple vs. exponential, generalized or not, etc.) to obtain the coefficients α , a 0 , a 1 , , a p . Based on the independence of the W t , we can use Monte Carlo and/or bootstrap to generate different W n + k * for k = 1 , , h , and consequently approximate the distribution of f h ( W n + 1 , , W n + h ) . Denote by F ^ W the empirical distribution of the transformed data W p + 1 , , W n . Similar to Algorithm 2, we can use either a (truncated)3 standard normal distribution or F ^ W to generate the pseudo-values W n + k * .
Algorithm 3:h-step ahead prediction for NoVaS: bootstrap.
Step 1.
Use one of the NoVaS methods (simple vs. exponential, generalized or not, etc.) to obtain the transformed data { W t for t = p + 1 , , n } and the coefficients α , p, and a 0 , a 1 , , a p .
Step 2.
Compute the analytic form of Equation (16), i.e., express X n + h as a function of { X 1 , , X n } and  { W n + 1 , , W n + h } using the values { a 0 , a 1 , , a p } obtained in Step 1.
Step 3.
Generate { W n + 1 * , , W n + h * } as i . i . d . either from a (truncated) standard normal distribution or from  F ^ W , and plug them into the function (16) to obtain the pseudo-value X n + h * . Repeat the above procedure M times and denote the M pseudo-values by { X n + h ( 1 ) , . , X n + h ( M ) } .
Step 4.
Calculate the optimal predictor g ( X n + h ) ^ of g ( X n + h ) by taking the sample median (under L 1 risk) or sample mean (under L 2 risk) of the set { g ( X n + h ( 1 ) ) , , g ( X n + h ( M ) ) } .

3. Optimal Multi-Step-Ahead Prediction Intervals

Going beyond point prediction, it may be desirable to construct prediction intervals for g ( X n + h ) with a target coverage level ( 1 β ) 100 % . One-step ahead prediction intervals have been discussed in detail in Pan and Politis (2016); see also Chapter 10 of Politis (2015). In this section, we will propose a construction of multi-step ahead prediction intervals in the given setting of financial returns data { X 1 , , X n } .
As explained in Politis (2015), the bootstrap is a sine qua non for the construction of prediction intervals as it allows us to incorporate the variability of estimated quantities in our estimate of the conditional distribution of g ( X n + h ) given { X 1 , , X n } . The variability of estimated quantities is not so important in point prediction when only the center of the conditional distribution is of interest. However, it is crucial in order to obtain an estimate of the conditional distribution that is not too narrow, yielding prediction intervals with accurate coverage.
In what follows, we give the algorithms for bootstrap prediction intervals for g ( X n + h ) in the two settings: model-based for ARCH/GARCH models and model-free based on NoVaS. The algorithms follow the “forward bootstrap” paradigm introduced in Pan and Politis (2016).
Algorithm 4: Bootstrap prediction intervals for g ( X n + h ) under ARCH/GARCH models.
Step 1.
Fit the ARCH(p) model to the data { X 1 , , X n } , i.e., obtain the estimators { a ^ 0 , a ^ 1 , , a ^ p } and the residuals { ϵ ^ p + 1 , , ϵ ^ n } .
Step 2.
Use Algorithm 1 or Algorithm 2 to compute g ( X n + h ) ^ , the point predictor of g ( X n + h ) of choice.
Step 3.
(a) Re-sample (with replacement) the residuals { ϵ ^ p + 1 , , ϵ ^ n } to create the pseudo-errors ϵ p + 1 * , , ϵ n * and ϵ n + 1 * , , ϵ n + h * .
(b) Let ( X 1 * , , X p * ) = ( X 1 + I , , X p + I ) where I is generated as a discrete random variable uniform on the values 0 , 1 , , n p . Now, use the fitted ARCH model of Step 1 to generate bootstrap pseudo-data X t * for t = p + 1 , , n in a recursive manner.
(c) Based on the bootstrap data X 1 * , , X n * , re-estimate the parameters obtaining { a ^ 0 * , a ^ 1 * , , a ^ p * } .
(d) Re-define the last p values of the bootstrap data to match the original, i.e., re-define X t * = X t  for t = n p + 1 , , n ; this is the “forward bootstrap”4 construction.
(e) Use the fitted ARCH model of Step 1, the bootstrap data X 1 * , , X n * , and the pseudo-errors  ϵ n + 1 * , , ϵ n + h * to generate recursively the future bootstrap data X t * for t = n + 1 , , n + h .
(f). Based on the bootstrap data X n p + 1 * , , X n * and the re-estimated parameters { a ^ 0 * , a ^ 1 * , , a ^ p * } , use Algorithm 1 or 2 (according to which one was used in Step 2 to calculate the bootstrap predictor denoted by g ( X n + h * ) ^ .
(g) Calculate the bootstrap root5: g ( X n + h * ) g ( X n + h * ) ^ .
Step 4.
Repeat Step 3 above B times; the B bootstrap root replicates are collected in an empirical distribution whose α -quantile is denoted q ( α ) . The ( 1 β ) 100 % equal-tailed prediction interval for g ( X n + h ) is then given by [ g ( X n + h ) ^ + q ( β / 2 ) , g ( X n + h ) ^ + q ( 1 β / 2 ) ] .
Algorithm 5: Model-free (MF) bootstrap prediction intervals for g ( X n + h ) based on NoVaS.
Step 1.
Use one of the NoVaS algorithms (simple vs. exponential, generalized or not, etc.) to obtain the transformed data { W t for t = p + 1 , , n } that are assumed to be approximately i . i . d . Let p, α , and a i denote the fitted NoVaS parameters.
Step 2.
Use Algorithm 3 to calculate g ( X n + h ) ^ , the chosen point predictor of g ( X n + h ) .
Step 3.
(a) Re-sample randomly (with replacement)6 the transformed variables { W t for t = p + 1 , , n }  to create the pseudo-data W p + 1 * , , W n 1 * , W n * and W n + 1 * , , W n + h * .
(b) Let ( X 1 * , , X p * ) = ( X 1 + I , , X p + I ) where I is generated as a discrete random variable uniform on the values 0 , 1 , , n p . Generate the bootstrap pseudo-data X t * for t = p + 1 , , n  using:
X t * = W t * 1 a 0 W t * 2 α s t 1 * 2 + i = 1 p a i X t i * 2 for t = p + 1 , , n
where s t 1 * 2 = ( t 1 ) 1 k = 1 t 1 X k * 2 .
(c) Based on the bootstrap data X 1 * , , X n * , re-estimate the NoVaS parameters, obtaining p * , α * , and a i * ; for simplicity, we can keep the same value for p, i.e., let p * equal p.
(d) Re-define the last p values of the bootstrap data to match the original, i.e., re-define X t * = X t for t = n p + 1 , , n ; this is the “forward bootstrap” construction.
(e) Calculate the bootstrap future value X n + h * by iteration as:
X n + 1 * = W n + 1 * 1 a 0 W n + 1 * 2 α s n 2 + i = 1 p a i X n i + 1 2
where s n 2 = n 1 i = 1 n X i 2 .
If h < p , for j = 2 , , h :
X n + j * = W n + j * 1 a 0 W n + j * 2 α s n + 1 j 2 + k = 1 j 1 a k X n k + j * 2 + i = j p a i X n i + j 2
If h p , for j = 2 , , h :
X n + j * = W n + j * 1 a 0 W n + j * 2 α s n + 1 j 2 + i = 1 p a i X n i + j * 2
where s n + 1 j 2 = ( n + j 1 ) 1 ( i = 1 n X i 2 + k = 1 j 1 X n + k * 2 ) .
(f) Based on the bootstrap data X n p + 1 * , , X n * and the parameters p * , α * , a 0 * , a 1 * , , a p * , use Algorithm 3 to calculate the bootstrap predictor g ( X n + h * ) ^ .
(g) Calculate the bootstrap root: g ( X n + h * ) g ( X n + h * ) ^ .
Step 4.
Repeat Step 3 above B times; the B bootstrap root replicates are collected in an empirical distribution whose α -quantile is denoted q ( α ) . The ( 1 β ) 100 % equal-tailed prediction interval for g ( X n + h ) is given by [ g ( X n + h ) ^ + q ( β / 2 ) , g ( X n + h ) ^ + q ( 1 β / 2 ) ] .

4. Simulations and Finite Sample Performance

In this section, we conduct simulations to examine the finite sample performance of our algorithms.

4.1. Settings

In the simulation, 200 datasets X n = ( X 1 , , X n ) , each of size n = 100 , are generated separately by the following seven different GARCH(1,1) models.
  • Model 1. Standard GARCH with Gaussian errors and finite fourth moment:
    X t = σ t ϵ t , σ t 2 = 0.00001 + 0.73 σ t 1 2 + 0.10 X t 1 2 , { ϵ t } i . i . d . N ( 0 , 1 ) .
  • Model 2. Standard GARCH with Gaussian errors and infinite fourth moment:
    X t = σ t ϵ t , σ t 2 = 0.00001 + 0.8895 σ t 1 2 + 0.10 X t 1 2 , { ϵ t } i . i . d . N ( 0 , 1 ) .
  • Model 3. Standard GARCH with Student-t errors:
    X t = σ t ϵ t , σ t 2 = 0.00001 + 0.73 σ t 1 2 + 0.10 X t 1 2 , { ϵ t } i . i . d . t distributed with five degrees of freedom.
  • Model 4. GARCH with time-varying parameters (TV-GARCH):
    The value of β decreases as a linear function of t, starting at β 1 = 0.10 for t = 1 , and ending at β = 0.05 for t = n . At the same time, the value of α increases as a linear function of t, starting at α = 0.73 for t = 1 , and ending at α = 0.93 for t = n . ω = 0.00001 and { ϵ t } i . i . d . N ( 0 , 1 ) .
  • Model 5. Two-state Markov Switching GARCH(1,1) (MS-GARCH):
    X t = σ t ϵ t , σ t 2 = s = 1 2 1 { P ( S t = s ) } [ ω s + α t σ t 1 2 + β s X t 1 2 ]
    In the first regime, we set α 1 = 0.9 , β 1 = 0.07 , ω 1 = 2.4 e 5 . In the second regime, we set α 2 = 0.7 , β 2 = 0.22 , ω 2 = 1.2 e 4 . The transition probabilities for the first regime are p 11 = 0.9 and p 12 = 0.1 , while for the second regime, we use p 21 = 0.3 and p 22 = 0.7 . { ϵ t } i . i . d . N ( 0 , 1 ) .
  • Model 6. Smooth transition GARCH (ST-GARCH):
    X t = [ a b ( t / T ) ] σ t ϵ t , σ t 2 = ω + α σ t 1 2 + β X t 1 2
    where { ϵ t } i . i . d . N ( 0 , 1 ) . ω = 1.2 e 5 , α = 0.9 , β = 0.07 , a = α + β = 0.97 , and b = β / α 0.078 .
  • Model 7. Stochastic volatility model (SV-GARCH):
    X t | h t N ( 0 , exp ( h t ) ) ,
    h t | h t 1 N ( μ + ϕ ( h t 1 μ ) , η 2 ) , h 0 N ( μ , η 2 / ( 1 ϕ 2 ) ) ,
    where μ = 10 , ϕ = 0.95 , η = 0.2 .
We performed up to five-step-ahead point predictions and interval predictions for each dataset. We used M = 5000 simulations to compute the point predictions. For the bootstrap prediction intervals, we used B = 300 replications; we would have liked to use a higher number of bootstrap replications, but it was practically infeasible having 200 simulated datasets in each of which the point predictors were computed in a computer-intensive way. However, a practitioner having a single dataset at hand could (and should) use a higher B, e.g., B = 1000 or more.
Five models or transformations were used to fit the data in both point predictions and interval predictions as follows: fitting a GARCH(1,1) model, simple-NoVaS, exponential NoVaS (Exp-NoVaS), generalized simple NoVaS (GS-NoVaS), and generalized exponential NoVaS (GE-NoVaS).
In point predictions, the mean absolute deviations (MAD) and mean squared errors (MSE) for five-step-ahead point predictions in both the L 1 and L 2 sense (the absolute value or the square of the prediction error at the updated time point averaged over the 200 replications) were recorded. Furthermore, the bootstrap prediction interval ( L i , U i ) with a nominal coverage 95 % was constructed for the future values X n + h with h = 1 , , 5 .
The corresponding empirical average coverage level (CVR) and the average length (LEN) of the constructed intervals and the standard error (St.err) associated with each length of the constructed intervals are calculated as:
C V R = 1 N i = 1 N 1 [ L i , U i ] X ( n + h , i )
L E N = 1 N i = 1 N L E N i and S t . e r r = 1 N i = 1 N ( L E N i L E N ) 2
where L E N i = U i L i .

4.2. Results and Discussions

The simulation results for point predictions are shown in Table 1, Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9, Table 10, Table 11, Table 12, Table 13, Table 14, Table 15, Table 16, Table 17, Table 18, Table 19, Table 20, Table 21, Table 22, Table 23, Table 24, Table 25, Table 26, Table 27 and Table 28. The following conclusions can be obtained from the results:
  • When comparing the MADs between the L 1 and L 2 predictions by fitting the same models, we can find that the MADs of L 1 predictions were always smaller than those of L 2 predictions. Furthermore, we can find that the MSEs of L 1 predictions were always bigger than those of L 2 predictions, when comparing the MSEs between the L 1 and L 2 predictions with the same model settings. This was expected since in the L 1 sense, we tried to minimize the mean absolute deviations, while for L 2 , the loss function to be minimized was the mean squared error.
  • Furthermore, for each model’s fitting results, there were no obvious error accumulation problems in the multi-step-ahead prediction for both L 1 and L 2 measures.
  • NoVaS methods consistently performed better than GARCH(1,1) for all data generating processes. When the prediction step h was higher, the difference of their respective performances became smaller.
As regards prediction intervals, the simulation results are summarized in Table 29, Table 30, Table 31, Table 32, Table 33, Table 34 and Table 35. Generally, the conclusions were similar to those of the point predictions. Furthermore, NoVaS methods gave more accurate coverage than GARCH(1,1) in the L 1 sense prediction for all seven data generating processes. In the L 2 sense, when the data were generated from a standard GARCH(1,1) with normal errors, GARCH(1,1) gave as good coverage as NoVaS. When we used other models to generate data, for example, GARCH(1,1) with t distributed errors or MS-GARCH(1,1) or ST-GARCH(1,1) etc., GARCH(1,1) performed very poorly, while NoVaS methods were still performing well.
Our results accentuate the drawbacks of GARCH(1,1) associated with one-step ahead predictions from previous work; see Politis (2015) and the references therein. In particular, NoVaS methods were invariably more robust than GARCH(1,1) fitting when the data have a stochastic structure that deviates from a stationary GARCH(1,1) model, e.g., a time-varying GARCH, a GARCH with structure breaks, etc.

5. Conclusions

In this paper, we derived a new way of multi-step-ahead predictions for ARCH/GARCH and NoVaS methods only based on the basic assumptions of models or transformation. This method has good properties based on our theoretical methodology and simulated results. To sum up:
  • The ARCH/GARCH version of our algorithms worked well for data that are generated by a stationary GARCH(1,1) model.
  • The NoVaS version of our algorithms worked well for time series data that are either GARCH or have a stochastic structure that deviates from a stationary GARCH(1,1) model.
  • There was no apparent error accumulation issue in the multi-step-ahead prediction.
  • The methods were theoretically and computationally straightforward.
  • Combined with the one-step ahead prediction results in Politis (2015), NoVaS was shown to outperform GARCH model fitting most of the time, for h-step-ahead prediction for all h 1 .

Author Contributions

Conceptualization, D.N.P.; methodology, J.C. and D.N.P.; software, J.C.; validation, J.C.; formal analysis, J.C.; investigation, J.C.; resources, J.C.; data curation, J.C.; writing original draft preparation, J.C.; writing review and editing, D.N.P.; visualization, J.C.; supervision, D.N.P.; project administration, D.N.P.; funding acquisition, D.N.P.

Funding

This research was partially supported by the National Science Foundation Grant DMS 16-13026.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; nor in the decision to publish the results.

References

  1. Abadir, Karim, Alessandra Luati, and Paolo Paruolo. 2018. The Forecast Density of a Garch(1,1). Working Paper. Bologna: University of Bologna. [Google Scholar]
  2. Andersen, Torben G., Tim Bollerslev, Peter F. Christoffersen, and Francis X. Diebold. 2006. Volatility and Correlation Forecasting. In Handbook of Economic Forecasting. Amsterdam: Elsevier, vol. 1, chp. 15. pp. 777–878. [Google Scholar]
  3. Bollerslev, Tim, Ray Y. Chou, and Kenneth F. Kroner. 1992. ARCH modeling in finance: A review of the theory and empirical evidence. Journal of Econometrics 52: 5–59. [Google Scholar] [CrossRef]
  4. Bose, Arup, and Kanchan Mukherjee. 2009. Bootstrapping a weighted linear estimator of the arch parameters. Journal of Time Series Analysis 30: 315–31. [Google Scholar] [CrossRef]
  5. Engle, Robert F. 1982. Autoregressive conditional heteroscedasticity with estimates of the variance of united kingdom inflation. Econometrica 50: 987–1007. [Google Scholar] [CrossRef]
  6. Francq, Christian, and Jean-Michel Zakoian. 2011. GARCH Models: Structure, Statistical Inference and Financial Applications. Chichester: John Wiley & Sons. [Google Scholar]
  7. Pan, Li, and Dimitris N. Politis. 2014. Bootstrap prediction intervals for markov processes. Computational Statistics & Data Analysis 100: 467–94. [Google Scholar]
  8. Pan, Li, and Dimitris N. Politis. 2016. Bootstrap prediction intervals for linear, nonlinear and nonparametric autoregressions. Journal of Statistical Planning and Inference 177: 1–27. [Google Scholar] [CrossRef] [Green Version]
  9. Politis, Dimitris N. 2003. A normalizing and variance–stabilizing transformation for financial time series. In Recent Advances and Trends in Nonparametric Statistics. Edited by Michael G. Akritas and Dimitris N. Politis. Amsterdam: JAI, pp. 335–47. [Google Scholar] [CrossRef]
  10. Politis, Dimitris N. 2007. Model-free versus Model-based Volatility Prediction. Journal of Financial Econometrics 5: 358–59. [Google Scholar] [CrossRef]
  11. Politis, Dimitris N. 2015. Model-Free Prediction and Regression: A Transformation-Based Approach to Inference. New York: Springer. [Google Scholar]
1
The reason for not including time t in the variance estimator is for purposes of notational clarity, as well as the easy identifiability of the effect of the coefficient a 0 associated with X t 2 in the denominator of Equation (13).
2
The bootstrap validity is not shown in this paper, because it is beyond the scope of the paper.
3
From (13), it follows that:
1 W t 2 = α s t 1 2 + a 0 X t 2 + i = 1 p a i X t i 2 X t 2 a 0
since all the parameters are nonnegative; thus, | W t | 1 / a 0 , i.e., the range of the W t is finite. Typically, the { W t } variables have a large enough range such that the boundedness is not seen as spoiling the normality from a practical perspective, but in any theoretical works and/or simulations, it is necessary to use the standard normal distribution truncated to the range ± 1 / a 0 .
4
Most of the bootstrap methods will not have their last p values X t *  where t = n , , n p + 1 , exactly equal to the original values as needed for the prediction process. Herein lies the problem, since the behavior of the predictors for future values needs to be captured conditionally on the original values. In this forward bootstrap step, we redefine the last p values of X t *  and make them equal to the original ones; see Pan and Politis (2014).
5
Since we do not have much information of the distribution of g ( X n + h ) g ( X n + h ) ^ , we could use the distribution of g ( X n + h * ) g ( X n + h * ) ^ to approximate it. Therefore, we can employ the quantile of the approximated distribution to calculate the prediction intervals in the following steps. See more theoretical inference in Politis (2015).
6
It is also possible to use other bootstrap methods, for example wild bootstrap or block bootstrap, to create pseudo-data. It depends on the information you have for W t .
Table 1. MADs of L 1 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.8895 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) . Exp, exponential; GS, generalized simple; GE, generalized exponential.
Table 1. MADs of L 1 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.8895 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) . Exp, exponential; GS, generalized simple; GE, generalized exponential.
Prediction Step12345
Fitting a GARCH8.23 × 10−57.35 × 10−56.93 × 10−58.86 × 10−51.49 × 10−4
Simple-NoVaS6.99 × 10−57.88 × 10−58.29 × 10−51.06 × 10−41.63 × 10−4
Exp-NoVaS7.21 × 10−58.28 × 10−58.72 × 10−51.14 × 10−41.70 × 10−4
GS-NoVaS6.30 × 10−57.31 × 10−58.07 × 10−58.71 × 10−59.94 × 10−5
GE-NoVaS7.02 × 10−58.44 × 10−58.76 × 10−51.16 × 10−41.71 × 10−4
Table 2. MADs of L 2 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.8895 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) .
Table 2. MADs of L 2 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.8895 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) .
Prediction Step12345
Fitting a GARCH1.48 × 10−42.07 × 10−42.67 × 10−43.20 × 10−44.60 × 10−4
Simple-NoVaS5.22 × 10−56.44 × 10−56.61 × 10−58.82 × 10−51.49 × 10−4
Exp-NoVaS5.13 × 10−56.41 × 10−56.58 × 10−58.84 × 10−51.50 × 10−4
GS-NoVaS4.71 × 10−55.97 × 10−56.06 × 10−58.48 × 10−51.46 × 10−4
GE-NoVaS4.84 × 10−56.26 × 10−56.38 × 10−58.73 × 10−51.48 × 10−4
Table 3. MSEs of L 1 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.8895 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) .
Table 3. MSEs of L 1 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.8895 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) .
Prediction Step12345
Fitting a GARCH2.48 × 10−84.42 × 10−84.36 × 10−81.94 × 10−71.12 × 10−6
Simple-NoVaS1.23 × 10−84.24 × 10−84.17 × 10−81.94 × 10−71.12 × 10−6
Exp-NoVaS1.18 × 10−84.21 × 10−84.14 × 10−81.94 × 10−71.12 × 10−6
GS-NoVaS1.08 × 10−84.09 × 10−83.94 × 10−81.92 × 10−71.13 × 10−6
GE-NoVaS1.10 × 10−84.14 × 10−84.06 × 10−81.92 × 10−71.12 × 10−6
Table 4. MSEs of L 2 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.8895 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) .
Table 4. MSEs of L 2 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.8895 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) .
Prediction Step12345
Fitting a GARCH1.03 × 10−72.11 × 10−73.92 × 10−76.26 × 10−71.86 × 10−6
Simple-NoVaS1.36 × 10−84.15 × 10−84.10 × 10−81.93 × 10−71.11 × 10−6
Exp-NoVaS1.31 × 10−84.13 × 10−84.07 × 10−81.93 × 10−71.12 × 10−6
GS-NoVaS1.03 × 10−83.88 × 10−83.73 × 10−81.91 × 10−71.13 × 10−6
GE-NoVaS1.13 × 10−83.98 × 10−83.93 × 10−81.88 × 10−71.11 × 10−6
Table 5. MADs of L 1 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) .
Table 5. MADs of L 1 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) .
Prediction Step12345
Fitting a GARCH4.84 × 10−54.69 × 10−55.18 × 10−55.45 × 10−55.99 × 10−5
Simple-NoVaS4.94 × 10−54.73 × 10−55.26 × 10−55.59 × 10−56.00 × 10−5
Exp-NoVaS4.87 × 10−54.69 × 10−55.24 × 10−55.55 × 10−56.00 × 10−5
GS-NoVaS4.83 × 10−54.67 × 10−55.17 × 10−55.44 × 10−55.94 × 10−5
GE-NoVaS4.84 × 10−54.69 × 10−55.25 × 10−55.43 × 10−56.00 × 10−5
Table 6. MADs of L 2 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) .
Table 6. MADs of L 2 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) .
Prediction Step12345
Fitting a GARCH5.68 × 10−55.45 × 10−55.54 × 10−56.25 × 10−5.26 × 10−5
Simple-NoVaS6.20 × 10−56.05 × 10−56.06 × 10−56.96 × 10−56.67 × 10−5
Exp-NoVaS6.20 × 10−56.09 × 10−56.15 × 10−57.05 × 10−56.87 × 10−5
GS-NoVaS5.98 × 10−55.81 × 10−55.84 × 10−56.55 × 10−56.37 × 10−5
GE-NoVaS5.50 × 10−55.30 × 10−55.75 × 10−56.22 × 10−56.23 × 10−5
Table 7. MSEs of L 1 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) .
Table 7. MSEs of L 1 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) .
Prediction step12345
Fitting a GARCH7.94 × 10−98.48 × 10−97.80 × 10−98.67 × 10−91.03 × 10−8
Simple-NoVaS8.00 × 10−98.35 × 10−97.85 × 10−98.88 × 10−91.03 × 10−8
Exp-NoVaS7.95 × 10−98.33 × 10−97.83 × 10−98.84 × 10−91.03 × 10−8
GS-NoVaS7.78 × 10−98.35 × 10−97.63 × 10−98.59 × 10−91.02 × 10−8
GE-NoVaS8.02 × 10−98.64 × 10−97.99 × 10−98.92 × 10−91.02 × 10−8
Table 8. MSEs of L 2 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) .
Table 8. MSEs of L 2 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . N ( 0 , 1 ) .
Prediction Step12345
Fitting a GARCH7.28 × 10−97.47 × 10−96.50 × 10−97.64 × 10−98.65 × 10−9
Simple-NoVaS7.79 × 10−97.50 × 10−96.84 × 10−98.67 × 10−98.85 × 10−9
Exp-NoVaS7.48 × 10−97.34 × 10−96.64 × 10−98.39 × 10−98.70 × 10−9
GS-NoVaS7.12 × 10−97.40 × 10−96.41 × 10−97.90 × 10−98.43 × 10−9
GE-NoVaS6.99 × 10−97.44 × 10−96.43 × 10−97.68 × 10−98.53 × 10−9
Table 9. MADs of L 1 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . t distributed with degrees of freedom of 5.
Table 9. MADs of L 1 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . t distributed with degrees of freedom of 5.
Prediction Step12345
Fitting a GARCH1.46 × 10−41.26 × 10−41.29 × 10−41.84 × 10−41.77 × 10−4
Simple-NoVaS1.44 × 10−41.24 × 10−41.27 × 10−41.82 × 10−41.76 × 10−4
Exp-NoVaS1.43 × 10−41.22 × 10−41.26 × 10−41.80 × 10−41.75 × 10−4
GS-NoVaS1.43 × 10−41.23 × 10−41.27 × 10−41.80 × 10−41.76 × 10−4
GE-NoVaS1.43 × 10−41.23 × 10−41.29 × 10−41.81 × 10−41.77 × 10−4
Table 10. MADs of L 2 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . t distributed with degrees of freedom of 5.
Table 10. MADs of L 2 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . t distributed with degrees of freedom of 5.
Prediction Step12345
Fitting a GARCH1.72 × 10−41.53 × 10−41.63 × 10−42.17 × 10−42.16 × 10−4
Simple-NoVaS1.58 × 10−41.35 × 10−41.44 × 10−41.97 × 10−41.92 × 10−4
Exp-NoVaS1.57 × 10−41.35 × 10−41.46 × 10−41.98 × 10−41.93 × 10−4
GS-NoVaS1.56 × 10−41.35 × 10−41.45 × 10−41.99 × 10−41.92 × 10−4
GE-NoVaS1.59 × 10−41.34 × 10−41.48 × 10−42.01 × 10−41.94 × 10−4
Table 11. MSEs of L 1 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . t distributed with degrees of freedom of 5.
Table 11. MSEs of L 1 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . t distributed with degrees of freedom of 5.
Prediction Step12345
Fitting a GARCH1.06 × 10−77.11 × 10−88.97 × 10−82.95 × 10−73.81 × 10−7
Simple-NoVaS1.03 × 10−76.86 × 10−89.00 × 10−82.93 × 10−73.83 × 10−7
Exp-NoVaS1.03 × 10−76.83 × 10−89.00 × 10−82.93 × 10−73.82 × 10−7
GS-NoVaS1.03 × 10−76.90 × 10−89.05 × 10−82.94 × 10−73.84 × 10−7
GE-NoVaS1.05 × 10−76.88 × 10−89.23 × 10−82.96 × 10−73.85 × 10−7
Table 12. MSEs of L 2 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . t distributed with degrees of freedom of 5.
Table 12. MSEs of L 2 predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and { ϵ t } i . i . d . t distributed with degrees of freedom of 5.
Prediction Step12345
Fitting a GARCH1.05 × 10−76.94 × 10−88.91 × 10−82.84 × 10−73.79 × 10−7
Simple-NoVaS9.24 × 10−86.03 × 10−88.32 × 10−82.77 × 10−73.66 × 10−7
Exp-NoVaS9.16 × 10−85.96 × 10−88.28 × 10−82.76 × 10−73.64 × 10−7
GS-NoVaS9.17 × 10−86.07 × 10−88.35 × 10−82.77 × 10−73.67 × 10−7
GE-NoVaS9.48 × 10−86.12 × 10−88.53 × 10−82.80 × 10−73.67 × 10−7
Table 13. MADs of L 1 predictions for data generated from GARCH(1,1) with slowing-varying parameters (TV-GARCH).
Table 13. MADs of L 1 predictions for data generated from GARCH(1,1) with slowing-varying parameters (TV-GARCH).
Prediction Step12345
Fitting a GARCH1.94 × 10−42.17 × 10−42.01 × 10−41.76 × 10−42.09 × 10−4
Simple-NoVaS1.91 × 10−42.12 × 10−42.03 × 10−41.72 × 10−42.07 × 10−4
Exp-NoVaS1.91 × 10−42.12 × 10−42.02 × 10−41.72 × 10−42.06 × 10−4
GS-NoVaS1.91 × 10−42.13 × 10−42.02 × 10−41.73 × 10−42.06 × 10−4
GE-NoVaS1.97 × 10−42.17 × 10−42.02 × 10−41.79 × 10−42.13 × 10−4
Table 14. MADs of L 2 predictions for data generated from GARCH(1,1) with slowing-varying parameters (TV-GARCH).
Table 14. MADs of L 2 predictions for data generated from GARCH(1,1) with slowing-varying parameters (TV-GARCH).
Prediction Step12345
Fitting a GARCH1.90 × 10−42.13 × 10−42.03 × 10−41.73 × 10−42.02 × 10−4
Simple-NoVaS1.94 × 10−42.11 × 10−42.07 × 10−41.74 × 10−42.10 × 10−4
Exp-NoVaS1.94 × 10−42.11 × 10−42.05 × 10−41.74 × 10−42.09 × 10−4
GS-NoVaS1.91 × 10−42.10 × 10−41.99 × 10−41.71 × 10−42.03 × 10−4
GE-NoVaS1.89 × 10−42.09 × 10−41.99 × 10−41.72 × 10−42.07 × 10−4
Table 15. MSEs of L 1 predictions for data generated from GARCH(1,1) with slowing-varying parameters (TV-GARCH).
Table 15. MSEs of L 1 predictions for data generated from GARCH(1,1) with slowing-varying parameters (TV-GARCH).
Prediction Step12345
Fitting a GARCH1.21 × 10−71.70 × 10−71.28 × 10−71.06 × 10−71.58 × 10−7
Simple-NoVaS1.15 × 10−71.61 × 10−71.23 × 10−79.86 × 10−81.53 × 10−7
Exp-NoVaS1.16 × 10−71.61 × 10−71.23 × 10−79.95 × 10−81.52 × 10−7
GS-NoVaS1.16 × 10−71.62 × 10−71.24 × 10−79.99 × 10−81.53 × 10−7
GE-NoVaS1.25 × 10−71.70 × 10−71.31 × 10−71.08 × 10−71.61 × 10−7
Table 16. MSEs of L 2 predictions for data generated from GARCH(1,1) with slowing-varying parameters (TV-GARCH).
Table 16. MSEs of L 2 predictions for data generated from GARCH(1,1) with slowing-varying parameters (TV-GARCH).
Prediction Step12345
Fitting a GARCH1.07 × 10−71.55 × 10−71.12 × 10−79.44 × 10−81.41 × 10−7
Simple-NoVaS9.60 × 10−81.41 × 10−71.06 × 10−78.36 × 10−81.32 × 10−7
Exp-NoVaS9.60 × 10−81.40 × 10−71.05 × 10−78.39 × 10−81.30 × 10−7
GS-NoVaS9.70 × 10−81.41 × 10−71.07 × 10−78.51 × 10−81.32 × 10−7
GE-NoVaS1.04 × 10−71.47 × 10−71.12 × 10−79.02 × 10−81.38 × 10−7
Table 17. MADs of L 1 predictions for data generated from two-state Markov switching GARCH(1,1) (MS-GARCH).
Table 17. MADs of L 1 predictions for data generated from two-state Markov switching GARCH(1,1) (MS-GARCH).
Prediction Step12345
Fitting a GARCH6.76 × 10−47.67 × 10−48.44 × 10−47.80 × 10−47.14 × 10−4
Simple-NoVaS7.00 × 10−47.74 × 10−48.85 × 10−47.90 × 10−47.26 × 10−4
Exp-NoVaS7.02 × 10−47.75 × 10−48.87 × 10−47.92 × 10−47.26 × 10−4
GS-NoVaS6.97 × 10−47.70 × 10−48.80 × 10−47.91 × 10−47.23 × 10−4
GE-NoVaS7.06 × 10−47.75 × 10−48.85 × 10−47.98 × 10−47.27 × 10−4
Table 18. MADs of L 2 predictions for data generated from two-state Markov switching GARCH(1,1) (MS-GARCH).
Table 18. MADs of L 2 predictions for data generated from two-state Markov switching GARCH(1,1) (MS-GARCH).
Prediction Step12345
Fitting a GARCH7.80 × 10−49.01 × 10−49.21 × 10−49.17 × 10−48.80 × 10−4
Simple-NoVaS6.77 × 10−47.50 × 10−48.59 × 10−47.73 × 10−47.10 × 10−4
Exp-NoVaS6.78 × 10−47.50 × 10−48.61 × 10−47.75 × 10−47.09 × 10−4
GS-NoVaS6.76 × 10−47.50 × 10−48.59 × 10−47.76 × 10−47.08 × 10−4
GE-NoVaS6.81 × 10−47.50 × 10−48.62 × 10−47.79 × 10−47.09 × 10−4
Table 19. MSEs of L 1 predictions for data generated from two-state Markov switching GARCH(1,1) (MS-GARCH).
Table 19. MSEs of L 1 predictions for data generated from two-state Markov switching GARCH(1,1) (MS-GARCH).
Prediction Step12345
Fitting a GARCH1.27 × 10−62.30 × 10−62.35 × 10−62.78 × 10−61.49 × 10−6
Simple-NoVaS1.45 × 10−62.46 × 10−62.73 × 10−63.09 × 10−61.76 × 10−6
Exp-NoVaS1.45 × 10−62.46 × 10−62.73 × 10−63.10 × 10−61.76 × 10−6
GS-NoVaS1.43 × 10−62.45 × 10−62.71 × 10−63.10 × 10−61.74 × 10−6
GE-NoVaS1.46 × 10−62.47 × 10−62.74 × 10−63.12 × 10−61.76 × 10−6
Table 20. MSEs of L 2 predictions for data generated from two-state Markov switching GARCH(1,1) (MS-GARCH).
Table 20. MSEs of L 2 predictions for data generated from two-state Markov switching GARCH(1,1) (MS-GARCH).
Prediction Step12345
Fitting a GARCH1.29 × 10−62.28 × 10−62.05 × 10−62.55 × 10−61.51 × 10−6
Simple-NoVaS1.34 × 10−62.33 × 10−62.58 × 10−62.97 × 10−61.64 × 10−6
Exp-NoVaS1.34 × 10−62.33 × 10−62.59 × 10−62.98 × 10−61.64 × 10−6
GS-NoVaS1.33 × 10−62.32 × 10−62.57 × 10−62.97 × 10−61.62 × 10−6
GE-NoVaS1.35 × 10−62.33 × 10−62.59 × 10−62.99 × 10−61.64 × 10−6
Table 21. MADs of L 1 predictions for data generated from smooth transition GARCH(1,1) (ST-GARCH).
Table 21. MADs of L 1 predictions for data generated from smooth transition GARCH(1,1) (ST-GARCH).
Prediction Step12345
Fitting a GARCH1.83 × 10−41.78 × 10−42.01 × 10−42.02 × 10−42.22 × 10−4
Simple-NoVaS1.83 × 10−41.79 × 10−42.03 × 10−42.02 × 10−42.24 × 10−4
Exp-NoVaS1.82 × 10−41.78 × 10−42.03 × 10−42.02 × 10−42.24 × 10−4
GS-NoVaS1.81 × 10−41.78 × 10−42.02 × 10−41.98 × 10−42.21 × 10−4
GE-NoVaS1.82 × 10−41.79 × 10−42.05 × 10−41.99 × 10−42.24 × 10−4
Table 22. MADs of L 2 predictions for data generated from smooth transition GARCH(1,1) (ST-GARCH).
Table 22. MADs of L 2 predictions for data generated from smooth transition GARCH(1,1) (ST-GARCH).
Prediction Step12345
Fitting a GARCH2.23 × 10−42.14 × 10−42.18 × 10−42.39 × 10−42.42 × 10−4
Simple-NoVaS1.91 × 10−41.88 × 10−42.02 × 10−42.09 × 10−42.19 × 10−4
Exp-NoVaS1.89 × 10−41.86 × 10−42.03 × 10−42.09 × 10−42.19 × 10−4
GS-NoVaS1.90 × 10−41.89 × 10−42.03 × 10−42.09 × 10−42.17 × 10−4
GE-NoVaS1.86 × 10−41.83 × 10−42.06 × 10−42.03 × 10−42.20 × 10−4
Table 23. MSEs of L 1 predictions for data generated from smooth transition GARCH(1,1) (ST-GARCH)
Table 23. MSEs of L 1 predictions for data generated from smooth transition GARCH(1,1) (ST-GARCH)
Prediction Step12345
Fitting a GARCH1.12 × 10−71.23 × 10−71.14 × 10−71.13 × 10−71.43 × 10−7
Simple-NoVaS1.16 × 10−71.24 × 10−71.19 × 10−71.21 × 10−71.48 × 10−7
Exp-NoVaS1.16 × 10−71.24 × 10−71.20 × 10−71.21 × 10−71.49 × 10−7
GS-NoVaS1.12 × 10−71.22 × 10−71.16 × 10−71.16 × 10−71.45 × 10−7
GE-NoVaS1.18 × 10−71.27 × 10−71.22 × 10−71.21 × 10−71.50 × 10−7
Table 24. MSEs of L 2 predictions for data generated from smooth transition GARCH(1,1) (ST-GARCH).
Table 24. MSEs of L 2 predictions for data generated from smooth transition GARCH(1,1) (ST-GARCH).
Prediction Step12345
Fitting a GARCH1.05 × 10−71.13 × 10−79.61 × 10−81.04 × 10−71.25 × 10−7
Simple-NoVaS1.01 × 10−71.09 × 10−79.96 × 10−81.04 × 10−71.23 × 10−7
Exp-NoVaS9.93 × 10−81.08 × 10−71.00 × 10−71.03 × 10−71.23 × 10−7
GS-NoVaS9.84 × 10−81.08 × 10−79.85 × 10−81.00 × 10−71.21 × 10−7
GE-NoVaS1.00 × 10−71.11 × 10−71.02 × 10−71.02 × 10−71.24 × 10−7
Table 25. MADs of L 1 predictions for data generated from the stochastic volatility model (SV-GARCH).
Table 25. MADs of L 1 predictions for data generated from the stochastic volatility model (SV-GARCH).
Prediction Step12345
Fitting a GARCH5.35 × 10−55.77 × 10−54.72 × 10−54.28 × 10−53.86 × 10−5
Simple-NoVaS5.34 × 10−55.69 × 10−54.79 × 10−54.45 × 10−53.78 × 10−5
Exp-NoVaS5.35 × 10−55.66 × 10−54.74 × 10−54.35 × 10−53.75 × 10−5
GS-NoVaS5.23 × 10−55.68 × 10−54.74 × 10−54.23 × 10−53.79 × 10−5
GE-NoVaS5.23 × 10−55.71 × 10−54.79 × 10−54.31 × 10−53.85 × 10−5
Table 26. MADs of L 2 predictions for data generated from the stochastic volatility model (SV-GARCH).
Table 26. MADs of L 2 predictions for data generated from the stochastic volatility model (SV-GARCH).
Prediction Step12345
Fitting a GARCH5.77 × 10−56.23 × 10−55.56 × 10−55.29 × 10−55.12 × 10−5
Simple-NoVaS6.02 × 10−56.17 × 10−55.85 × 10−55.99 × 10−55.41 × 10−5
Exp-NoVaS6.10 × 10−56.19 × 10−55.99 × 10−56.04 × 10−55.68 × 10−5
GS-NoVaS5.74 × 10−56.29 × 10−56.10 × 10−55.85 × 10−55.69 × 10−5
GE-NoVaS5.66 × 10−55.96 × 10−55.45 × 10−55.28 × 10−54.77 × 10−5
Table 27. MSEs of L 1 predictions for data generated from the stochastic volatility model (SV-GARCH).
Table 27. MSEs of L 1 predictions for data generated from the stochastic volatility model (SV-GARCH).
Prediction Step12345
Fitting a GARCH1.04 × 10−81.69 × 10−81.08 × 10−87.03 × 10−95.06 × 10−9
Simple-NoVaS1.02 × 10−81.65 × 10−81.05 × 10−87.06 × 10−94.82 × 10−9
Exp-NoVaS1.04 × 10−81.65 × 10−81.06 × 10−86.99 × 10−94.78 × 10−9
GS-NoVaS1.00 × 10−81.66 × 10−81.04 × 10−86.92 × 10−94.80 × 10−9
GE-NoVaS1.03 × 10−81.70 × 10−81.08 × 10−86.94 × 10−95.04 × 10−9
Table 28. MSEs of L 2 predictions for data generated from the stochastic volatility model (SV-GARCH).
Table 28. MSEs of L 2 predictions for data generated from the stochastic volatility model (SV-GARCH).
Prediction Step12345
Fitting a GARCH8.99 × 10−91.51 × 10−81.02 × 10−86.76 × 10−95.29 × 10−9
Simple-NoVaS8.70 × 10−91.42 × 10−89.99 × 10−97.62 × 10−95.49 × 10−9
Exp-NoVaS8.82 × 10−91.39 × 10−81.00 × 10−87.41 × 10−95.51 × 10−9
GS-NoVaS8.10 × 10−91.42 × 10−81.02 × 10−86.92 × 10−95.32 × 10−9
GE-NoVaS8.37 × 10−91.46 × 10−81.00 × 10−86.62 × 10−95.27 × 10−9
Table 29. Interval predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.8895 , θ = 0.10 , and ϵ i . i . d N ( 0 , 1 ) . CVR, average coverage level; LEN, average length.
Table 29. Interval predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.8895 , θ = 0.10 , and ϵ i . i . d N ( 0 , 1 ) . CVR, average coverage level; LEN, average length.
L2L1
GARCH(1,1)GARCH(1,1)
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.7141.28 × 10−21.77 × 10−210.7441.01 × 10−21.55 × 10−2
20.7461.28 × 10−21.79 × 10−220.7341.13 × 10−21.83 × 10−2
30.7461.37 × 10−21.87 × 10−230.7681.19 × 10−21.95 × 10−2
40.7661.18 × 10−21.59 × 10−240.7341.20 × 10−21.99 × 10−2
50.7861.28 × 10−21.75 × 10−250.7441.28 × 10−22.04 × 10−2
EXP-NoVaSEXP-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9582.00 × 10−21.35 × 10−210.9361.99 × 10−21.39 × 10−2
20.9522.07 × 10−21.44 × 10−220.9362.19 × 10−21.52 × 10−2
30.9521.98 × 10−21.27 × 10−230.9441.94 × 10−21.27 × 10−2
40.9462.17 × 10−21.43 × 10−240.9442.08 × 10−21.28 × 10−2
50.9502.19 × 10−21.32 × 10−250.9362.19 × 10−21.51 × 10−2
Simple-NoVaSSimple-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9461.98 × 10−21.35 × 10−210.9341.84 × 10−21.32 × 10−2
20.9421.78 × 10−21.41 × 10−220.9461.91 × 10−21.49 × 10−2
30.9461.92 × 10−21.35 × 10−230.9362.04 × 10−21.47 × 10−2
40.9462.07 × 10−21.34 × 10−240.9641.93 × 10−21.40 × 10−2
50.9562.21 × 10−21.43 × 10−250.9542.07 × 10−21.45 × 10−2
GS-NoVaSGS-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9482.05 × 10−21.39 × 10−210.951.61 × 10−21.21 × 10−2
20.9421.78 × 10−21.41 × 10−220.942.07 × 10−21.47 × 10−2
30.9521.93 × 10−21.43 × 10−230.9361.62 × 10−21.25 × 10−2
40.9482.14 × 10−21.30 × 10−240.9361.70 × 10−21.42 × 10−2
50.9542.26 × 10−21.40 × 10−250.941.78 × 10−21.32 × 10−2
GE-NoVaSGE-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9481.73 × 10−21.14 × 10−210.9581.98 × 10−21.13 × 10−2
20.9481.66 × 10−21.05 × 10−220.9421.68 × 10−21.16 × 10−2
30.9521.98 × 10−21.27 × 10−230.9441.94 × 10−21.27 × 10−2
40.9462.17 × 10−21.43 × 10−240.9442.08 × 10−21.28 × 10−2
50.9502.19 × 10−21.32 × 10−250.9422.32 × 10−21.13 × 10−2
Table 30. Interval predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and ϵ i . i . d N ( 0 , 1 ) .
Table 30. Interval predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and ϵ i . i . d N ( 0 , 1 ) .
L2L1
GARCH(1,1)GARCH(1,1)
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9484.39 × 10−32.11 × 10−310.923.76 × 10−33.18 × 10−3
20.9404.53 × 10−32.23 × 10−320.9365.49 × 10−35.21 × 10−3
30.9504.47 × 10−32.74 × 10−330.9385.99 × 10−35.52 × 10−3
40.9524.02 × 10−31.89 × 10−340.9227.16 × 10−36.31 × 10−3
50.9343.77 × 10−32.74 × 10−350.924.57 × 10−34.21 × 10−3
EXP-NoVaSEXP-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9544.46 × 10−32.69 × 10−310.9544.74 × 10−32.29 × 10−3
20.9724.40 × 10−32.55 × 10−320.954.62 × 10−32.19 × 10−3
30.9384.17 × 10−32.59 × 10−330.954.64 × 10−32.12 × 10−3
40.9584.62 × 10−32.59 × 10−340.9484.58 × 10−32.12 × 10−3
50.9504.58 × 10−32.49 × 10−350.9424.45 × 10−31.88 × 10−3
Simple-NoVaSSimple-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9604.54 × 10−32.50 × 10−310.9464.28 × 10−32.45 × 10−3
20.9584.27 × 10−32.97 × 10−320.954.26 × 10−32.33 × 10−3
30.9684.63 × 10−32.87 × 10−330.9524.21 × 10−32.62 × 10−3
40.9604.73 × 10−32.85 × 10−340.9544.25 × 10−32.55 × 10−3
50.9484.15 × 10−32.93 × 10−350.9484.19 × 10−32.32 × 10−3
GS-NoVaSGS-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9494.37 × 10−32.53 × 10−310.9464.26 × 10−32.37 × 10−3
20.9484.63 × 10−22.78 × 10−320.954.26 × 10−32.33 × 10−3
30.9384.17 × 10−32.59 × 10−330.954.22 × 10−32.42 × 10−3
40.9453.76 × 10−32.71 × 10−340.9484.20 × 10−31.91 × 10−3
50.9504.58 × 10−32.49 × 10−350.9484.19 × 10−32.32 × 10−3
GE-NoVaSGE-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9464.92 × 10−32.57 × 10−310.965.37 × 10−32.00 × 10−3
20.9464.68 × 10−32.38 × 10−320.9485.13 × 10−33.34 × 10−3
30.9584.39 × 10−32.35 × 10−330.9524.03 × 10−32.05 × 10−3
40.9544.30 × 10−32.03 × 10−340.954.80 × 10−32.03 × 10−3
50.9484.15 × 10−32.93 × 10−350.9444.42 × 10−32.78 × 10−3
Table 31. Results of interval predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and ϵ i . i . d t 5 .
Table 31. Results of interval predictions for data generated from GARCH(1,1) with ω = 0.00001 , α = 0.73 , θ = 0.10 , and ϵ i . i . d t 5 .
L2L1
GARCH(1,1)GARCH(1,1)
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9244.17 × 10−33.18 × 10−310.9362.35 × 10−38.60 × 10−3
20.9313.75 × 10−32.57 × 10−320.9282.15 × 10−37.50 × 10−3
30.9224.47 × 10−32.24 × 10−330.922.37 × 10−38.52 × 10−3
40.9254.02 × 10−32.63 × 10−340.9382.92 × 10−36.95 × 10−3
50.9224.56 × 10−32.79 × 10−350.922.79 × 10−38.20 × 10−3
EXP-NoVaSEXP-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9454.46 × 10−32.66 × 10−310.953.38 × 10−32.95 × 10−3
20.9424.25 × 10−32.72 × 10−320.9583.80 × 10−32.72 × 10−3
30.9434.54 × 10−32.77 × 10−330.9463.75 × 10−32.40 × 10−3
40.9494.94 × 10−32.72 × 10−340.9523.76 × 10−32.43 × 10−3
50.9544.72 × 10−33.08 × 10−350.9463.40 × 10−32.84 × 10−3
Simple-NoVaSSimple-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9504.82 × 10−32.35 × 10−310.9534.01 × 10−33.07 × 10−3
20.9634.86 × 10−33.26 × 10−320.9423.58 × 10−32.73 × 10−3
30.9664.85 × 10−32.82 × 10−330.9523.28 × 10−32.48 × 10−3
40.9545.04 × 10−33.05 × 10−340.9523.51 × 10−32.54 × 10−3
50.9444.36 × 10−32.51 × 10−350.9444.20 × 10−33.03 × 10−3
GS-NoVaSGS-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9504.82 × 10−32.35 × 10−310.9483.64 × 10−34.49 × 10−3
20.9544.47 × 10−32.93 × 10−320.9563.22 × 10−35.92 × 10−3
30.9524.69 × 10−33.03 × 10−330.9463.31 × 10−34.08 × 10−3
40.9504.57 × 10−32.99 × 10−340.9483.23 × 10−34.52 × 10−3
50.9544.50 × 10−33.06 × 10−350.953.62 × 10−34.54 × 10−3
GE-NoVaSGE-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9504.90 × 10−32.57 × 10−310.9543.65 × 10−32.78 × 10−3
20.9464.36 × 10−32.93 × 10−320.9543.86 × 10−32.84 × 10−3
30.9484.48 × 10−32.82 × 10−330.9463.64 × 10−32.81 × 10−3
40.9524.58 × 10−32.78 × 10−340.953.53 × 10−32.89 × 10−3
50.9524.53 × 10−32.94 × 10−350.9666.38 × 10−33.92 × 10−3
Table 32. Results of interval predictions for data generated from TV-GARCH(1,1).
Table 32. Results of interval predictions for data generated from TV-GARCH(1,1).
L2L1
GARCH(1,1)GARCH(1,1)
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.8501.57 × 10−31.07 × 10−210.7261.12 × 10−34.97 × 10−3
20.8481.20 × 10−31.60 × 10−320.7161.32 × 10−35.88 × 10−3
30.8581.84 × 10−32.77 × 10−330.7188.87 × 10−33.90 × 10−3
40.8442.32 × 10−32.02 × 10−340.7169.78 × 10−34.39 × 10−3
50.8562.29 × 10−31.68 × 10−350.721.24 × 10−35.77 × 10−3
EXP-NoVaSEXP-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9482.40 × 10−32.20 × 10−310.9522.76 × 10−32.45 × 10−3
20.9502.23 × 10−32.17 × 10−320.9522.74 × 10−32.50 × 10−3
30.9522.93 × 10−32.15 × 10−330.952.69 × 10−32.62 × 10−3
40.9543.02 × 10−32.18 × 10−340.9422.78 × 10−32.65 × 10−3
50.9502.86 × 10−32.12 × 10−350.9482.82 × 10−32.61 × 10−3
Simple-NoVaSSimple-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9502.32 × 10−32.37 × 10−310.9422.20 × 10−32.30 × 10−3
20.9562.20 × 10−32.36 × 10−320.9562.82 × 10−32.45 × 10−3
30.9602.88 × 10−32.15 × 10−330.9482.60 × 10−32.50 × 10−3
40.9522.50 × 10−32.27 × 10−340.9462.79 × 10−32.36 × 10−3
50.9542.80 × 10−32.01 × 10−350.9462.47 × 10−32.62 × 10−3
GS-NoVaSGS-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9502.32 × 10−32.37 × 10−310.9422.20 × 10−32.30 × 10−3
20.9502.52 × 10−32.57 × 10−320.9482.57 × 10−32.33 × 10−3
30.9522.59 × 10−32.53 × 10−330.9522.52 × 10−32.18 × 10−3
40.9502.37 × 10−32.05 × 10−340.9462.97 × 10−32.22 × 10−3
50.9502.62 × 10−32.12 × 10−350.952.64 × 10−32.21 × 10−3
GE-NoVaSGE-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9482.40 × 10−32.20 × 10−310.9482.77 × 10−32.66 × 10−3
20.9502.23 × 10−32.17 × 10−320.9522.71 × 10−32.83 × 10−3
30.9422.47 × 10−32.25 × 10−330.952.53 × 10−32.50 × 10−3
40.9482.44 × 10−32.17 × 10−340.952.75 × 10−32.49 × 10−3
50.9492.29 × 10−32.12 × 10−350.9542.64 × 10−32.47 × 10−3
Table 33. Results of interval predictions for data generated from MS-GARCH(1,1).
Table 33. Results of interval predictions for data generated from MS-GARCH(1,1).
L2L1
GARCH(1,1)GARCH(1,1)
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.8683.32 × 10−21.78 × 10−210.8563.27 × 10−21.25 × 10−2
20.8723.42 × 10−21.50 × 10−220.893.04 × 10−21.10 × 10−2
30.8683.58 × 10−21.62 × 10−230.8823.10 × 10−21.07 × 10−2
40.8583.67 × 10−21.86 × 10−240.8863.09 × 10−21.12 × 10−2
50.873.60 × 10−22.10 × 10−250.9087.66 × 10−31.28 × 10−2
EXP-NoVaSEXP-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9464.14 × 10−22.18 × 10−210.9523.97 × 10−22.67 × 10−2
20.9484.02 × 10−22.26 × 10−220.9444.22 × 10−22.80 × 10−2
30.964.78 × 10−22.19 × 10−230.9583.99 × 10−22.74 × 10−2
40.9584.16 × 10−22.06 × 10−240.9383.86 × 10−22.80 × 10−2
50.9564.27 × 10−22.07 × 10−250.9444.20 × 10−22.97 × 10−2
Simple-NoVaSSimple-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9544.21 × 10−22.87 × 10−210.9583.05 × 10−22.03 × 10−2
20.9483.98 × 10−22.81 × 10−220.9363.08 × 10−22.26 × 10−2
30.944.47 × 10−22.91 × 10−230.9363.45 × 10−22.13 × 10−2
40.9484.26 × 10−22.80 × 10−240.943.42 × 10−22.28 × 10−2
50.9464.32 × 10−22.93 × 10−250.9383.52 × 10−22.16 × 10−2
GS-NoVaSGS-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9544.21 × 10−22.87 × 10−210.9483.21 × 10−22.16 × 10−2
20.9483.98 × 10−22.81 × 10−220.9463.26 × 10−22.18 × 10−2
30.9424.75 × 10−22.84 × 10−230.9483.46 × 10−22.28 × 10−2
40.9484.26 × 10−22.80 × 10−240.9523.17 × 10−22.21 × 10−2
50.9464.32 × 10−22.93 × 10−250.9463.22 × 10−22.02 × 10−2
GE-NoVaSGE-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9544.30 × 10−22.08 × 10−210.9483.40 × 10−22.19 × 10−2
20.9523.93 × 10−22.09 × 10−220.9423.24 × 10−22.09 × 10−2
30.9484.36 × 10−22.03 × 10−230.9463.63 × 10−22.15 × 10−2
40.9484.20 × 10−22.05 × 10−240.953.08 × 10−22.39 × 10−2
50.954.29 × 10−22.07 × 10−250.9443.55 × 10−22.90 × 10−2
Table 34. Results of interval predictions for data generated from ST-GARCH(1,1).
Table 34. Results of interval predictions for data generated from ST-GARCH(1,1).
L2L1
GARCH(1,1)GARCH(1,1).
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.892.80 × 10−32.14 × 10−310.8942.53 × 10−33.80 × 10−3
20.8882.19 × 10−32.79 × 10−320.9022.88 × 10−34.28 × 10−3
30.9042.15 × 10−32.87 × 10−330.8842.48 × 10−33.20 × 10−3
40.9082.07 × 10−32.04 × 10−340.90122.73 × 10−34.21 × 10−3
50.8962.00 × 10−32.09 × 10−350.892.02 × 10−34.92 × 10−3
EXP-NoVaSEXP-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9563.35 × 10−32.57 × 10−310.9583.70 × 10−32.69 × 10−3
20.963.40 × 10−32.58 × 10−320.9443.60 × 10−32.54 × 10−3
30.9463.60 × 10−32.73 × 10−330.9663.64 × 10−32.65 × 10−3
40.9423.42 × 10−32.51 × 10−340.9563.52 × 10−32.46 × 10−3
50.9443.56 × 10−32.70 × 10−350.9623.66 × 10−32.53 × 10−3
Simple-NoVaSSimple-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9523.44 × 10−32.32 × 10−310.9443.44 × 10−32.75 × 10−3
20.9563.49 × 10−32.24 × 10−320.9383.65 × 10−32.75 × 10−3
30.9583.46 × 10−32.04 × 10−330.9343.67 × 10−32.81 × 10−3
40.9543.40 × 10−32.10 × 10−340.9463.60 × 10−32.72 × 10−3
50.9463.73 × 10−32.19 × 10−350.9383.51 × 10−32.51 × 10−3
GS-NoVaSGS-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9523.49 × 10−32.18 × 10−310.9463.73 × 10−32.16 × 10−3
20.9563.30 × 10−32.29 × 10−320.9563.86 × 10−32.31 × 10−3
30.9483.51 × 10−32.29 × 10−330.9543.84 × 10−32.27 × 10−3
40.9483.54 × 10−32.36 × 10−340.953.82 × 10−32.15 × 10−3
50.9443.56 × 10−32.70 × 10−350.9463.76 × 10−32.22 × 10−3
GE-NoVaSGE-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9463.37 × 10−32.10 × 10−310.9543.64 × 10−32.60 × 10−3
20.9443.51 × 10−32.76 × 10−320.9523.60 × 10−32.55 × 10−3
30.9463.60 × 10−32.73 × 10−330.9563.61 × 10−32.55 × 10−3
40.9483.44 × 10−32.49 × 10−340.953.76 × 10−32.68 × 10−3
50.9483.69 × 10−32.61 × 10−350.9443.46 × 10−32.31 × 10−3
Table 35. Results of interval predictions for data generated from SV-GARCH(1,1).
Table 35. Results of interval predictions for data generated from SV-GARCH(1,1).
L2L1
GARCH(1,1)GARCH(1,1)
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.8782.24 × 10−22.99 × 10−210.8841.22 × 10−21.93 × 10−2
20.8622.11 × 10−22.51 × 10−220.8821.28 × 10−22.10 × 10−2
30.8962.49 × 10−22.78 × 10−230.8721.17 × 10−21.60 × 10−2
40.872.14 × 10−22.18 × 10−240.8781.26 × 10−21.46 × 10−2
50.8922.32 × 10−22.33 × 10−250.8761.33 × 10−21.63 × 10−2
EXP-NoVaSEXP-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9542.71 × 10−22.59 × 10−210.9441.84 × 10−22.33 × 10−2
20.9522.80 × 10−22.67 × 10−220.932.01 × 10−22.46 × 10−2
30.9562.79 × 10−22.66 × 10−230.9541.96 × 10−22.45 × 10−2
40.952.84 × 10−22.65 × 10−240.9422.19 × 10−22.62 × 10−2
50.9683.07 × 10−22.84 × 10−250.931.88 × 10−22.43 × 10−2
Simple-NoVaSSimple-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9522.68 × 10−22.79 × 10−210.9422.08 × 10−22.40 × 10−2
20.952.77 × 10−22.92 × 10−220.9322.46 × 10−22.53 × 10−2
30.9462.78 × 10−22.74 × 10−230.952.63 × 10−22.84 × 10−2
40.9583.05 × 10−22.90 × 10−240.9262.25 × 10−22.65 × 10−2
50.9542.66 × 10−22.73 × 10−250.9342.07 × 10−22.41 × 10−2
GS-NoVaSGS-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9522.68 × 10−22.79 × 10−210.952.17 × 10−22.42 × 10−2
20.952.77 × 10−22.92 × 10−220.962.05 × 10−22.40 × 10−2
30.9522.84 × 10−22.96 × 10−230.952.18 × 10−22.44 × 10−2
40.9442.28 × 10−22.34 × 10−240.952.22 × 10−22.48 × 10−2
50.9462.34 × 10−22.18 × 10−250.9422.28 × 10−22.55 × 10−2
GE-NoVaSGE-NoVaS
STEPSCVRLENST.ERRSTEPSCVRLENST.ERR
10.9542.71 × 10−22.59 × 10−210.9522.02 × 10−22.35 × 10−2
20.9482.52 × 10−22.29 × 10−220.9482.18 × 10−22.47 × 10−2
30.9562.79 × 10−22.66 × 10−230.9481.95 × 10−22.29 × 10−2
40.952.84 × 10−22.65 × 10−240.9422.29 × 10−22.58 × 10−2
50.9462.74 × 10−22.46 × 10−250.9462.10 × 10−22.34 × 10−2

Share and Cite

MDPI and ACS Style

Chen, J.; Politis, D.N. Optimal Multi-Step-Ahead Prediction of ARCH/GARCH Models and NoVaS Transformation. Econometrics 2019, 7, 34. https://0-doi-org.brum.beds.ac.uk/10.3390/econometrics7030034

AMA Style

Chen J, Politis DN. Optimal Multi-Step-Ahead Prediction of ARCH/GARCH Models and NoVaS Transformation. Econometrics. 2019; 7(3):34. https://0-doi-org.brum.beds.ac.uk/10.3390/econometrics7030034

Chicago/Turabian Style

Chen, Jie, and Dimitris N. Politis. 2019. "Optimal Multi-Step-Ahead Prediction of ARCH/GARCH Models and NoVaS Transformation" Econometrics 7, no. 3: 34. https://0-doi-org.brum.beds.ac.uk/10.3390/econometrics7030034

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

Article Metrics

Back to TopTop