VaR (Value-at-Risk) and CVaR (Conditional Value-at-Risk) have become the standard measures of market risk management. Their popularity has led to a large number of publications on this topic in recent years. Definition, description of the properties and comparative analysis of these risk measures can be found, for example, in    . Various methods for their evaluation and forecasting that represents different approaches are proposed. Most of the methods that provide explicit formulas for CVaR estimation are described in  . Optimization approach for CVaR evaluation is given in   . Non-parametric methods of estimation can be found, for example, in   . A large number of works devoted to the method of VaR and CVaR estimating based on the stochastic time series model. The basic ideas of the approach can be found for example in    . A significant number of works show the practical application of the approach for estimating and forecasting of stock indices, see for example     .
At the same time, during the global financial turmoil, the problem of constructing of new approaches for VaR and CVaR estimating and forecasting remains relevant. In this paper, we propose a new method for VaR and CVaR prediction for financial time series. The method takes into account the most statistically significant extreme values of data and the presence of the long-range dependence that is typical for financial time series   . For the convenience of practical application, the method is formulated as an incremental algorithm. At each step, the system of tests is proposed to evaluate the quality of the obtained results.
The proposed algorithm is used for forecasting VaR and CVaR for the time series of daily log return Nikkey225 Stock Index. The analysis of the obtained forecast estimates confirms their high quality. The formatter will need to create these components, incorporating the applicable criteria that follow.
2. Key Definitions
The continuously distributed random variable with finite mean defined on the probability space is considered. Here is the information set containing all available at the time information about the time series. Series is assumed to be stationary. It is accepted that the time series has the property of the long-range dependence  : there is and so that:
For a fixed confidence level dynamic risk measures and are defined as  :
denotes expectation with respect to .
The aim of the study is to construct a model for and where is an arbitrary moment of time. The forecasting values are determined by extrapolation of the values of this model (21 cm × 28.5 cm).
3. Forecast Methodology
In the article  , the most popular methods for dynamic and estimating are analyzed, their classification is given and the recommendations for their use are proposed. In accordance with the formulated in the article the structural scheme of selection of dynamic risk measures estimation the approach based on a stochastic time series model is chosen.
Suppose that the time series is a trajectory of stochastic process, that is:
where conditional mean and variation are defined on the information space , (independent, identically distributed random variables with a conditional distribution function ). Let is a random variable with the same distribution as any random variable from . Then    :
It is necessary to construct the forecast model for to determine its days forecast and to estimate and for a random variable . Then the forecasting values for dynamic risk measures can be found under the following formulas:
Hereinafter it is assumed that the trend, that defines , is absent (or removed from the data)  . Please do not revise any of the current designations.
4. An Algorithm for Constructing the Dynamic Risk Measures VaR and CVaR Forecast Taking into Account the Long-Range Dependence Presence
For the convenience of the practical application the proposed method for and forecasting is formulated as an incremental algorithm.
Step 1. For the time series a time series of variances (TSV) is constructed. General analysis of the studied time series and the TSV is carried out, the dependence of time series members (and their squares) from their previous values; the volatility and normality are analyzed.
Step 2. The TSV is tested on the long-range dependence. The Hurst parameter is estimated using five standard methods: the aggregated variance method, the method of absolute values of the aggregated series, the periodogram method, the method of residuals of regression, the R/S method  . Average value is chosen as the Hurst parameter estimation.
Step 3. The model for forecasting is estimated using the model and taking into account the long-range dependence of the WFD. The actualization of the model by reducing it to the model AR (∞) is performed. The method of smoothing of the autocorrelation function (ACF) proposed by the authors in  (the new method) is used. The least square method is used to determine the autoregression coefficients . So the problem is reduced to the infinite system of Yule-Walker equations  :
The regression equation for ACF based on the definition of the long-range dependence (1) is used to get estimates for , , . With the help of the optimization procedure  the Hurst parameter estimate and the estimates are corrected.
Using instead of the reduced system of normal Equations (5) is constructed and using the Holetskogo method the vector of assessments is found. As it is shown in  the solution of the reduced system converges to the exact solution.
The lag of the reduced AR model is determined using the information criterions: AIK (Akaike information criterion), HQC (Hannan-Quinn information criterion), SBIC (Bayesian information criterion)  . The lag value is chosen on the basis of minimum deviation.
The quality of the obtained AR model is checked. The variance ratio test  is used to test if the residuals of the model are iid (independent and identically distributed). The resulting model is used to obtain .
Step 4. The residuals of the model (2) are analyzed. Using (step 3) the implementations of a random variable are built. are analyzed on iid (the variance ratio test) and other properties. In accordance with the results using the classification scheme given in  , the method to get and estimates is chosen. The estimates , are obtained.
Step 5. With the results of steps 3 and 4 the model for dynamic risk measures estimating (3) is ready. After building the dynamic risk measures estimations и their quality is analyzed using the Kupiec test, the Kristoffersen test and the test   .
Step 6. The built dynamic risk measures model is used to get the forecast. Using the model from step 3 the -step forecast for is built by the formulas:
Using the estimates , (step 4) the -step forecast for dynamic risk measures and (4) is obtained. Then the
index of the time series is increased by and the procedure is repeated as many times as necessary. Thus in each cycle of the algorithm application the model is updated to take into account new data.
Step 7. Using the back testing procedure, the quality of the predicted values and (step 6) is checked, the prediction errors , , are calculated. For estimates the BPoE-test  is used.
Schematic description of the proposed method is shown in Figure 1. Please take note of the following items when proofreading spelling and grammar.
5. Numerical Testing of the Algorithm
To demonstrate the proposed algorithm a forecast for dynamic risk measures for the time series of log returns on a daily basis is built. Data are collected from the oldest and the most well-known index of Asian markets (the time series )―a composite index of the 225 largest companies publicly traded in for the period from 2005 to 2015. has a relatively low homogeneous volatility. The aim of this study is to forecast risk measures at a regular market behavior, so data without three time intervals with high volatility of the global financial system (01.07.2008-01.07.2009, 01.01.2011-01.07.2011, 01.02.2013-01.12.2013) are considered. Historical data of are not available online, but upon request.
Table 1 demonstrates the descriptive statistics for the time-series and the squared series .
Figure 1. Schematic description of the proposed method of dynamic risk measures VaR and CVaR forecasting.
Table 1. Basic descriptive statistics of the N225_RED.
Skewness is about 0 and kurtosis is about 3, so the distributions are close to normal. Ljung-Box test  results for confirm the dependence of data (and squared data) on their previous values (the values of Q-statistic are larger than critical value 12.017).
Consider the half of the general sample-843 values. The estimates of the Hurst parameter are and (step 2). These values confirm the long-range dependence of the time series.
Simulate (step 3) using the method (the designation ) and for comparison the standard methodology (the designation ). The standard methodology uses the model with the coefficients found by the maximum likelihood method (MLH). The lag of the reduced model is . The results of the variance ratio test (0.99 < 1.96 for the method and 0.69 < 1.96 for the standard method) confirm that the residuals of the models are iid.
For both models are found (step 4) and their analysis is carried out. The results of the variance ratio test (0.98 < 1.96 for the method and 0.97 < 1.96 for the standard method) show that the residuals of the model (2) are iid.
Estimates , are obtained using the following methods  : the historical simulation method , the explicit formulas under the assumption of a normal distribution with the maximum likelihood estimates of the parameters ,the explicit formulas using and functions with the maximum likelihood estimates of the parameters ( and respectively), the empirical method . Table 2 demonstrates the results of the estimating.
Using the results of steps 3, 4 estimates (3) of the dynamic and (step 5) are obtained. Figure 2 demonstrates the simulated dynamic and (first 836 values) where the explicit formulas under the assumption of a normal distribution were used for risk measures model residuals estimating.
Figure 2. Simulated and predicted values (last 5 values) of dynamic risk measures with new method.
Table 2. The estimates of the statics , .
Conduct the analysis of quality (step 5) for estimates using the Kupiec test ( values of statistics ), the Kristoffersen test ( values of statistics ) and their combination ( values of statistics ). The obtained estimates are reliable if values exceed the given significance level (0.1 in our case). To analyze the estimates the test with statistics , , is used. If the estimates are good the statistics, , , are close to zero. Table 3 presents the results of the analysis of the dynamic risk measures estimates for the method and Table 4 for the standard method.
The analysis of the results shows that the method (on the assumption of the normal distribution of residuals) gives the best estimates for both methods: values of statistics are essentially more than 0.1. This is consistent with the results of the basic analysis (Table 1) and is confirmed by the results of the Jarque-Bera test  conducted for (5.624 < 5.649 and 4.38 < 5.649 for both proposed and standard methods, respectively). At the same time, all estimates obtained with method show the poor quality. The popular historical simulated method gives quality estimates only with the method. In addition all values of statistics for the method are greater than the appropriate values for the standard method. V-test shows good results for both and standard methods.
The built models are used for dynamic risk measures forecasting. Forecasting procedure is performed on the window length equal to the half of the general sample power (843 values). 5-day forecast is built (see Figure 2). Thus, it is assumed that the parameters of the model are adequate for a period 5 (or more) days, the estimates of static risk measures at the forecast horizon are unchanged.
The forecasting procedure (steps 2 - 6) is repeated 168 times, and each time 5 new values (the accumulation window) are added. Table 5 presents minimum , maximum and average values of the Hurst parameter obtained for the windows.
Table 5 shows, the values of the Hurst parameter confirm the long-range dependence for all windows (all values are essentially greater than 0.5). The minimum range of values (about 7%) shows stability of this characteristic.
Table 3. The results of the analysis of the dynamic risk measures estimates (SACF method).
Table 4. The results of the analysis of the dynamic risk measures estimations (standard method).
Table 5. Hurst parameter estimates.
Figure 3 demonstrates the results of variance forecasting using (6) with SACF and standard methods.
Visual comparison of the predicted and real values shows that the proposed new method better describes the dynamic behavior of the time series. Extreme values obtained with the new method are much closer to real values. The new method also exhibits less lag in extreme values determination. This can be explained by the fact that the new method uses the ACF prediction and takes into account the property of the long-range dependence. It should also be noted that the optimization procedure in the determination of the Hurst parameter has significantly improved the forecast stability.
Figure 3. Real values of variance for real data and forecast estimates obtained by and standard methods.
Table 6. , estimates for different windows (SACF method).
Table 7. , estimates for different windows(standardmethod).
ues (max-min) for the method is less than the range of values for the standard method due to the fact that the proposed method explicitly uses the smoothing procedure of ACF and as a result the distribution function is more stable.
The obtained results are used to get the time series of dynamic risk measures estimates (3). As an example Figure 4 shows the forecast estimates for and , obtained with the use of the paramdistr method.
Table 8 and Table 9 show that the prediction errors obtained for the method is less than the prediction errors obtained for the standard method. This proves the advantage of the proposed method. In addition Table 8 shows that methods and gives the best estimates. This once again confirms the previously accepted hypothesis of data normal distribution.
The quality of built forecast estimates is analyzed with BPoE test. Table 10 shows the BPoE values obtained for the initial data and for the
Figure 4. Forecast estimates for and for 168 windows.
Table 8. The prediction errors of and (SACF method).
Table 9. The prediction errors of and (standard method).
Table 10. Results of BPoЕ-test for .
forecast estimates with the use of the new method and the standard method . These values are compared with the chosen level of risk measures . The results show the high quality of the forecast estimates obtained by the new method.
Table 10 may be used to liken the methods used for forecasting by comparing the value of confidence level for predicted and real risk measures values. The standard method based on demonstrates the best results. At the same time, the proposed method shows the best results with the historical simulation method . Table 10 shows that the deviation of for the new method is substantially less than the deviation for the standard method .
In the article, a multi-step procedure for constructing the dynamic risk measures VaR and CVaR forecast is proposed. The procedure is designed for volatile series with the long-range dependence and is based on the heteroscedastic time series model. The optimization procedure for constructing and forecasting of ACF is used to find the model parameters. For the convenience of practical application, the prediction procedure is formulated as an algorithm. To test the proposed algorithm, the risk measures forecast for the time series of daily log return is built. Different tests carried out at different stages of the algorithm confirm the good quality of the obtained estimates.