Financial Time Series Analysis: Part II

Transcription

1 Department of Mathematics and Statistics, University of Vaasa, Finland Spring 2017

2 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

3 Background 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

4 Background Example 1 Consider the following daily close-to-close SP500 values. S&P 500 Daily Closing Prices and Returns [Jan to Feb ] Return (%) Index 200 day rolling mean Return 22 day rolling mean Index Date

5 Background For example a rolling k = 22 days ( month of trading days) mean m t = 1 k t u=t k+1 r u, (1) and volatility (annualized standard deviation) s t = 252 t (r u m t ) k 2, (2) u=t k+1 t = k,... T give visual idea of volatility dynamics of the returns.

6 Background S&P 500 Daily Closing Prices and Volatility [Jan to Feb ] Index 200 day rolling average Absolute return 22 day rolling volatility Index Volatitliy (%, p.a) Date

7 Background S&P 500 Daily Returns and Volatility [Jan to Feb ] Return 22 day rolling average Absolute return 22 day rolling volatility Return (%) Volatitliy (%, p.a) Date

8 Background Because squared observations are the building blocks of the variance of the series, the results suggest that the variation (volatility) of the series is time dependent. This leads to the so called ARCH-family of models. 3 Note: Volatility not directly observable!! Methods: a) Implied volatility b) Realized volatility c) Econometric modeling (stochastic volatility, ARCH) 3 The inventor of this modeling approach is Robert F. Engle (1982). Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation. Econometrica, 50,

9 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

10 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

11 ARCH-modles Consider a time series y t with mean µ = E[y t ] and denote by u t = y t µ the deviation of y t from its mean, t = 1,..., T. The conditional distribution of u t given information up to time point t 1 is denoted as u t F t 1 D(0, σ 2 t ), (3) where F t is the information available at time t (usually the past values of u t ; u 1,..., u t 1 ), D is an appropriated (conditional) distribution (e.g., normal or t-distribution). Then u t follows an ARCH(q) process if its (conditional) variance σ 2 t is of the form σ 2 t = var[u t F t 1 ] = ω + α 1 u 2 t 1 + α 2 u 2 t α q u 2 t q. (4)

12 ARCH-modles Furthermore, it is assumed that ω > 0, α i 0 for all i and α α q < 1. For short it is denoted u t ARCH(q). This reminds essentially an AR(q) process for the squared residuals, because defining ν t = u 2 t σ 2 t, we can write u 2 t = ω + α 1 u 2 t 1 + α 2 u 2 t α q u 2 t q + ν t. (5) Nevertheless, the error term ν t is time heteroscedastic, which implies that the conventional estimation procedure used in AR-estimation does not produce optimal results here.

13 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

14 Properties of ARCH-processes Consider (for the sake of simplicity) ARCH(1) process σ 2 t = ω + αu 2 t 1 (6) with ω > 0 and 0 α < 1 and u t u t 1 N(0, σ 2 t ). (a) u t is white noise: (i) Constant mean (zero): E[u t ] = E[ E t 1 [u t ] ] = E[0] = 0. (7) }{{} =0 Note E t 1 [u t ] = E[u t F t 1 ], the conditional expectation given information up to time t The law of iterated expectations: Consider time points t1 < t 2 such that F t1 F t2, then for any t > t 2 ] [ ] ] ] E t1 [E t2 [u t ] = E E [u t F t2 F t1 = E [u t F t1 = E t1 [u t ]. (8)

15 Properties of ARCH-processes (ii) Constant variance: Using again the law of iterated expectations, we get var[u t ] = E [ ] [ [ ]] ut 2 = E Et 1 u 2 t = E [ ] [ ] σt 2 = E ω + αu 2 t 1 = ω + αe [ ] ut 1 2. = ω(1 + α + α α n ) + α n+1 E [ ut n 1 2 ] }{{} 0, as n = ω ( n lim n i=0 αi) (9) ω = 1 α.

16 Properties of ARCH-processes (iii) Autocovariances: Exercise, show that autocovariances are zero, i.e., E[u t u t+k ] = 0 for all k 0. (Hint: use the law of iterated expectations.) (b) The unconditional distribution of u t is symmetric, but nonnormal: (i) Skewness: Exercise, show that E [ u 3 t ] = 0. (ii) Kurtosis: Exercise, show that under the assumption u t u t 1 N(0, σ 2 t ), and that α < 1/3, the kurtosis E [ ut 4 ] ω 2 = 3 (1 α) 2 1 α 2 1 3α 2. (10) Hint: If X N(0, σ 2 ) then E [ (X µ) 4] = 3(σ 2 ) 2 = 3σ 4.

17 Properties of ARCH-processes Because (1 α 2 )/(1 3α 2 ) > 1 we have that E [ ut 4 ] ω 2 > 3 (1 α) 2 = 3 [var[u t]] 2, (11) we find that the kurtosis of the unconditional distribution exceed that what it would be, if u t were normally distributed. Thus the unconditional distribution of u t is nonnormal and has fatter tails than a normal distribution with variance equal to var[u t ] = ω/(1 α).

18 Properties of ARCH-processes (c) Standardized variables: Write z t = u t σ 2 t (12) then z t NID(0, 1), i.e., normally and independently distributed. Thus we can always write u t = z t σ 2 t, (13) where z t independent standard normal random variables (strict white noise). This gives us a useful device to check after fitting an ARCH model the adequacy of the specification: Check the autocorrelations of the squared standardized series (see, Econometrics I).

19 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

20 Estimation Instead of the normal distribution, more popular conditional distributions in ARCH-modeling are t-distribution and the generalized error distribution (ged), of which the normal distribution is a special case. The unknown parameters are estimated usually by the method of the maximum likelihood (ML).

21 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

22 Generalized ARCH (GARCH) In practice the ARCH needs fairly many lags. Usually far less lags are needed by modifying the model to σt 2 = ω + αut βσt 1, 2 (14) with ω > 0, α > 0, β 0, and α + β < 1. The model is called the Generalized ARCH (GARCH) model. Usually the above GARCH(1,1) is adequate in practice. Econometric packages call α (coefficient of ut 1 2 ) the ARCH parameter and β (coefficient of σt 1 2 ) the GARCH parameter.

23 Generalized ARCH model (GARCH) Applying backward substitution, one easily gets σ 2 t = ω 1 β + α β j 1 ut j 2 (15) j=1 an ARCH( ) process. Thus the GARCH term captures all the history from t 2 backwards of the shocks u t.

24 Generalized ARCH models (GARCH) Example 2 MA(1)-GARCH(1,1) model of SP500 returns estimated with conditional normal, u t F t 1 N(0, σ 2 t ), and GED, u t F t 1 GED(0, σ 2 t ) The model is r t = µ + u t + θu t 1 σ 2 t = ω + αu 2 t 1 + βσ2 t 1. (16)

25 Generalized ARCH models (GARCH) S&P 500 Returns [Jan to Feb 8, 2017]

26 Generalized ARCH models (GARCH) Using R-program with package fgarch (see garchfit(formula = ret ~ arma(2, 0) + garch(1, 1), data = sp1, trace = FALSE) Conditional Distribution: norm Coefficient(s): mu ar1 ar2 omega alpha1 beta Std. Errors: based on Hessian Error Analysis: Estimate Std. Error t value Pr(> t ) mu e-05 *** ar *** ar omega e-10 *** alpha < 2e-16 *** beta < 2e-16 *** --- Signif. codes: 0 *** ** 0.01 * Log Likelihood: normalized: AIC: BIC:

27 Generalized ARCH models (GARCH) garchfit(formula = ret ~ arma(2, 0) + garch(1, 1), data = sp1, cond.dist = "ged", trace = FALSE) Conditional Distribution: ged Coefficient(s): mu ar1 ar2 omega alpha1 beta1 shape Std. Errors: based on Hessian Error Analysis: Estimate Std. Error t value Pr(> t ) mu e-09 *** ar e-05 *** ar * omega e-06 *** alpha < 2e-16 *** beta < 2e-16 *** shape < 2e-16 *** --- Signif. codes: 0 *** ** 0.01 * Log Likelihood: normalized: AIC: BIC:

28 Generalized ARCH models (GARCH) Goodness of fit (AIC, BIC) is marginally better for GED. The shape parameter estimate is < 2 (statistically significantly) indicated fat tails. Otherwise the coefficient estimates are about the same.

29 Generalized ARCH models (GARCH) Residual Diagnostics: Autocorrelations of squared standardized residuals: Normal: GED: Lag LB p-value LB p-value JB = , df = 2, p-value = (dropping the outlier JB = 43.96)

30 Generalized ARCH models (GARCH) SP500 Return residuals [MA(1)-GED-GARCH(1,1)] Conditional Volatility Volatility [sqrt(252*h_t)] Days

31 Generalized ARCH models (GARCH) Autocorrelations of the squared standardized residuals pass (approximately) the white noise test. Nevertheless, the normality of the standardized residuals is strongly rejected. Usually this affects mostly to stantard errors. Common practice is to use some sort of robust standard errors (e.g. White).

32 Generalized ARCH models (GARCH) The variance function can be extended by including regressors (exogenous or predetermined variables), x t, in it σ 2 t = ω + αu 2 t 1 + βσ 2 t 1 + πx t. (17) Note that if x t can assume negative values, it may be desirable to introduce absolute values x t in place of x t in the conditional variance function. For example, with daily data a Monday dummy could be introduced into the model to capture the non-trading over the weekends in the volatility.

33 Generalized ARCH models (GARCH) Example 3 Monday effect in SP500 returns and/or volatility? Pulse (additive) effect on mean and innovative effect on volatility y t = φ 0 + φ m M t + φ(y t 1 φ m M t 1 ) + u t σt 2 = ω + πm t + αut βσ2 t 1, (18) where M t = 1 if Monday, zero otherwise.

34 Generalized ARCH models (GARCH) SAS: conditional t-distribution proc autoreg data = tmp; model sp500 = mon/ nlag = 1 garch = (p=1, q=1) dist = t; hetero mon; run;... The AUTOREG Procedure GARCH Estimates SSE Observations 2321 MSE Uncond Var Log Likelihood Total R-Square SBC AIC MAE AICC MAPE Normality Test Pr > ChiSq <.0001 Standard Approx Variable DF Estimate Error t Value Pr > t Variable Label Intercept mon AR ARCH ARCH <.0001 GARCH <.0001 TDFI <.0001 Inverse of t DF HET E E

35 Generalized ARCH models (GARCH) Degrees of freedom estimate: 1/ No empirical evidence of a Monday effect in returns or volatility.

36 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

37 ARCH-M Model The regression equation may be extended by introducing the variance function into the equation y t = x tθ + λg(σ 2 t ) + u t, (19) where u t GARCH, and g is a suitable function (usually square root or logarithm). This is called the ARCH in Mean (ARCH-M) model [Engle, Lilien and Robbins (1987) 5 ]. The ARCH-M model is often used in finance, where the expected return on an asset is related to the expected asset risk. The coefficient λ reflects the risk-return tradeoff. 5 Econometrica, 55,

38 ARCH-M Model Example 4 Does the daily mean return of SP500 depend on the volatility level? Model AR(1)-GARCH(1, 1)-M with conditional t-distribution y t = φ 0 + φ 1 y t 1 + λ σ 2 t + u t σ 2 t = ω + αu 2 t 1 + βσ2 t 1, (20) where u t F t 1 t ν (t-distribution with ν degrees of freedom, to be estimated). proc autoreg data = tmp; model sp500 = / nlag = 1 garch = (p=1, q=1, mean = sqrt) dist = t; /* t-dist */ run;

39 ARCH-M Model The AUTOREG Procedure GARCH Estimates SSE Observations 2321 MSE Uncond Var. Log Likelihood Total R-Square SBC AIC MAE AICC MAPE Normality Test Pr > ChiSq <.0001 Standard Approx Variable DF Estimate Error t Value Pr > t Intercept AR ARCH ARCH <.0001 GARCH <.0001 DELTA TDFI <.0001 Inverse of t DF

40 ARCH-M Model The volatility term in the mean equation (DELTA) is not statistically significant neither is the constant, indicating no discernible drift in SP500 index (again ˆν = 1/ ).

41 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

42 Asymmetric ARCH: TARCH, EGARCH, PARCH A stylized fact in stock markets is that downward movements are followed by higher volatility than upward movements. A rough view of this can be obtained from the cross-autocorrelations of z t and z 2 t, where z t defined in (12).

43 Asymmetric ARCH: TARCH, EGARCH, PARCH Example 5 Cross-autocorrelations of z t and z 2 t from Ex 4 Cross-correlations of z and z^2 ACF Lag Some autocorrelations signify possible leverage effect.

44 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

45 The TARCH model Threshold ARCH, TARCH (Zakoian 1994, Journal of Economic Dynamics and Control, , Glosten, Jagannathan and Runkle 1993, Journal of Finance, ) is given by [TARCH(1,1)] σ 2 t = ω + αu 2 t 1 + γu 2 t 1d t 1 + βσ 2 t 1, (21) where ω > 0, α, β 0, α γ + β < 1, and d t = 1, if u t < 0 (bad news) and zero otherwise. The impact of good news is α and bad news α + γ. Thus, γ 0 implies asymmetry. Leverage exists if γ > 0.

46 The TARCH model Example 6 SP500 returns, MA(1)-TARCH model (EViews,

47 The TARCH model Dependent Variable: SP500 Method: ML - ARCH (Marquardt) - Generalized error distribution (GED) Sample (adjusted): 1/04/2000 3/27/2009 Included observations: 2321 after adjustments Convergence achieved after 28 iterations MA Backcast: 1/03/2000 Presample variance: backcast (parameter = 0.7) GARCH = C(3) + C(4)*RESID(-1)^2 + C(5)*RESID(-1)^2*(RESID(-1)<0) + C(6)*GARCH(-1) ============================================================== Variable Coefficient Std. Error z-stat Prob C MA(1) ============================================================== Variance Equation ============================================================== C RESID(-1)^ RES(-1)^2*(RES(-1)<0) GARCH(-1) ============================================================== R-squared Mean dependent var Adjusted R-squared S.D. dependent var S.E. of regression Akaike info criterion Sum squared resid Schwarz criterion Log likelihood Hannan-Quinn criter F-statistic Durbin-Watson stat Prob(F-statistic) ==============================================================

48 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

49 The EGARCH model Nelson (1991) (Econometrica, ) proposed the Exponential GARCH (EGARCH) model for the variance function of the form (EGARCH(1,1)) log σ 2 t = ω + β log σ 2 t 1 + α z t 1 + γz t 1, (22) where z t = u t / σ 2 t is the standardized shock. Again, the impact is asymmetric if γ 0, and leverage is present if γ < 0.

50 The EGARCH model Example 7 MA(1)-EGARCH(1,1)-M estimation results.

51 The EGARCH model Dependent Variable: SP500 Method: ML - ARCH (Marquardt) - Generalized error distribution (GED) Sample (adjusted): 1/04/2000 3/27/2009 Included observations: 2321 after adjustments Convergence achieved after 25 iterations MA Backcast: 1/03/2000 Presample variance: backcast (parameter = 0.7) LOG(GARCH) = C(3) + C(4)*ABS(RESID(-1)/@SQRT(GARCH(-1))) + C(5) *RESID(-1)/@SQRT(GARCH(-1)) + C(6)*LOG(GARCH(-1)) ============================================================= Variable Coefficient Std. Error z-statistic Prob C MA(1) ============================================================= Variance Equation ============================================================= C(3) C(4) C(5) C(6) ============================================================= GED PARAM =============================================================... ============================================================= Inverted MA Roots.07 =============================================================

52 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

53 Ding, Granger, and Engle (1993). A long memory property of stock market returns and a new model. Journal of Empirical Finance. PARC(1,1) σ δ t = ω + βσ δ t 1 + α( u t 1 γu t 1 ) δ, (23) where γ is the leverage parameter. Again γ > 0 implies leverage.

54 Power ARCH (PARCH) Example 8 R: fgarch::garchfit MA(1)-APARCH(1,1) results for SP500 returns. R parametrization: MA(1), mu and theta gfa <- fgarch::garchfit(sp500r~arma(0,1) + aparch(1,1), data = sp500r, cond.dist = "ged", trace=f) Mean and Variance Equation: data ~ arma(0, 1) + aparch(1, 1) Conditional Distribution: ged Error Analysis: Estimate Std. Error t value Pr(> t ) mu ma *** omega *** alpha e-11 *** gamma < 2e-16 *** beta < 2e-16 *** delta e-09 *** shape < 2e-16 *** --- Signif. codes: 0 *** ** 0.01 * Log Likelihood: normalized:

55 Power ARCH (PARCH) The leverage parameter ( gamma1 ) estimates to unity. ˆδ = ( ) does not deviate significantly from unity (standard deviation process).

56 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

57 Integrated GARCH (IGARCH) Often in GARCH ˆα + ˆβ 1. Engle and Bollerslev (19896). Modelling the persistence of conditional variances, Econometrics Reviews 5, 1 50, introduce integrated GARCH with α + β = 1. σ 2 t = ω + αu 2 t 1 + (1 α)σ 2 t 1. (24) Close to the EWMA (Exponentially Weighted Moving Average) specification σ 2 t = αu 2 t 1 + (1 α)σ 2 t 1 (25) favored often by practitioners (e.g. RiskMetrics). Unconditional variance does not exist [more details, see Nelson (1990). Stationarity and persistence in in the GARCH(1,1) model. Econometric Theory 6, ].

58 Predicting Volatility 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

59 Predicting Volatility Predicting Volatility Predicting with the GARCH models is straightforward. Generally a k-period forward prediction is of the form σt k 2 = E [ ] [ ] t u 2 t+k = E u 2 t+k F t (26) k = 1, 2,....

60 Predicting Volatility Predicting Volatility Because u t = σ t z t, (27) where generally z t i.i.d(0, 1). Thus in (26) E t [ u 2 t+k ] [ = E t σ 2 t+k zt+k 2 ] [ ] [ ] = E t σ 2 t+k Et z 2 t+k (z t are i.i.d(0, 1)) (28) = E t [ σ 2 t+k ]. This can be utilize to derive explicit prediction formulas in most cases.

61 Predicting Volatility Predicting Volatility ARCH(1): σ 2 t+1 = ω + αu 2 t (29) σ 2 t k = ω(1 αk 1 ) 1 α + α k 1 σ 2 t+1 = σ 2 + α k 1 (σ 2 t+1 σ2 ), where σ 2 = var[u t ] = ω 1 α Recursive fromula: σ σt k 2 = t+1 2 for k = 1 ω + ασt k 1 2 for k > 1 (30) (31) (32)

62 Predicting Volatility Predicting Volatility GARCH(1,1): σ 2 t+1 = ω + αu 2 t + βσ 2 t. (33) where σ 2 t k = ω(1 (α+β)k 1 ) 1 (α+β) + (α + β) k 1 σ 2 t+1 = σ 2 + (α + β) k 1 (σ 2 t+1 σ2 ), σ 2 = (34) ω 1 α β. (35) Recursive fromula: σ σt k 2 = t+1 2 for k = 1 ω + (α + β)σt k 1 2 for k > 1 (36)

63 Predicting Volatility Predicting Volatility IGARCH: σ 2 t+1 = ω + αu 2 t + (1 α)σ 2 t. (37) σ 2 t k = (k 1)ω + σ2 t+1. (38) Recursive fromula: σ σt k 2 = t+1 2 for k = 1 ω + σt k 1 2 for k > 1 (39)

64 Predicting Volatility Predicting Volatility TGARCH: σ 2 t+1 = ω + αu 2 t + γu 2 t d t + βσ 2 t. (40) with σ 2 t k = ω(1 (α+ 1 2 γ+β)k 1 ) 1 (α+ 1 2 γ+β) + (α γ + β)k 1 σ 2 t+1 = σ 2 + (α γ + β)k 1 (σ 2 t+1 σ2 ) σ 2 = (41) ω 1 α 1 2 γ β. (42) Recursive fromula: σ σt k 2 = t+1 2 for k = 1 ω + (α γ + β)σ2 t k 1 for k > 1 (43)

65 Predicting Volatility Predicting Volatility EGARCH and APARCH prediction equations are a bit more involved. Recursive formulas are more appropriate in these cases. The volatility forecasts are applied for example in Value At Risk computations.

66 Predicting Volatility Predicting Volatility Evaluation of predictions unfortunately not that straightforward! (see, Andersen and Bollerslev (1998) International Economic Review.)

67 Realized volatility 1 Volatility Models Background ARCH-models Properties of ARCH-processes Estimation of ARCH models Generalized ARCH models (GARCH) ARCH-M Model Asymmetric ARCH: TARCH, EGARCH, PARCH The TARCH model The EGARCH model Power ARCH (PARCH) Integrated GARCH (IGARCH) Predicting Volatility Realized volatility

68 Realized volatility Realized volatility Under certain assumptions, volatility during a period of time can be estimated more and more precisely as the frequency of the returns increases. Daily (log) returns r t are sums of intraday returns (e.g. returns calculated at 30 minutes interval) r t = m r t (h) (44) h=1 where r t (h) = log P t (h) log P t (h 1) is the day s t intraday return in time interval [h 1, h], h = 1,..., m, P t (h) is the price at time point h within the day t, P t (0) is the opening price and P t (m) is the closing price.

69 Realized volatility Realized volatility The realized variance for day t is defined as m ˆσ t,m 2 = rt 2 (h) (45) and the realized volatility is ˆσ t,m = h=1 ˆσ 2 t,m which is typically presented in percentages per annum (i.e., scaled by the square root of the number of trading days and presented in percentages). Under certain conditions it can be shown that ˆσ t,m σ t as m, i.e., when the intraday return interval 0. For a resent survey on RV, see: Andersen, T.G. and L. Benzoni (2009). Realized volatility. In Handbook of Financial Time Series, T.G. Andersen, R.A. Davis, J-P. Kreiss and T. Mikosh (eds), Springer, New York, pp SSRN version is available at: id= &rec=1&srcabs=903659