Vladimir Spokoiny (joint with J.Polzehl) Varying coefficient GARCH versus local constant volatility modeling.

Size: px
Start display at page:

Download "Vladimir Spokoiny (joint with J.Polzehl) Varying coefficient GARCH versus local constant volatility modeling."

Transcription

1 W e ie rstra ß -In stitu t fü r A n g e w a n d te A n a ly sis u n d S to c h a stik STATDEP 2005 Vladimir Spokoiny (joint with J.Polzehl) Varying coefficient GARCH versus local constant volatility modeling. Mohrenstr Berlin spokoiny@wias-berlin.de

2 Outline Introduction: GARCH modeling. Varying coefficient GARCH. Estimation by AWS Semiparametric modeling. Simulated results and applications to some financial datasets Conclusion and Outlooks STATDEP

3 Volatility Estimation Problem Let S t an asset price process, and R t = log(s t /S t 1 ) log-returns. Conditional Heteroscedasticity Model: R t = σ t ε t, ε t N (0, 1). Typical Problem: estimate volatility θ T = σ 2 T from R 1,..., R T. STATDEP

4 GARCH(1,1) model Returns R 1,..., R T follow R t = σ t ε t with ε t N (0, 1). GARCH(1,1) model: X s = σs 2 fulfills X s = Ψ s θ = ω + αy s 1 + βx s 1. Here Ψ s = ( ) 1, Y s 1, X s 1 and θ = (ω, α, β). Maximum likelihood estimator: Define for a given θ = (ω, α, β) the process X s (θ) by X s (θ) = Ψ s (θ)θ = ω + αy s 1 + βx s 1 (θ). The MLE θ θ = argsup θ L(θ) = arginf θ s ( ) R 2 s X s (θ) + log X s(θ). If ε t N (0, 1), then θ a quasi MLE. STATDEP

5 Estimating equation The MLE θ fulfills the estimating equation dl(θ)/dθ = 0 leading to ( ) Rs 2 X s (θ)) X s (θ) 2 X s (θ) = 0. s For a given θ, the process X s (θ) is defined by the structural linear equation X s (θ) = Ψ s (θ)θ = ω + αy s 1 + βx s 1 (θ), s > t 0, X t0 = η. Here η, the initial value. Similarly for the derivatives X s (θ) = dx s (θ)/dθ and 2 X s (θ) = d 2 X s (θ)/dθ 2 : X s (θ) = Ψ s (θ) + β X s 1 (θ), s > t 0, X t0 = 0, 2 X s (θ) = Ψ s (θ) + Ψ s (θ) + β 2 X s 1 (θ), s > t 0, 2 X t0 = 0 with Ψ s (θ) = (0, 0, X s 1 (θ)). STATDEP

6 Numerical solution by Newton-Raphson Fix some initial value θ (0). Let, for k 1, θ (k 1) be the estimated parameter vector after step k 1. Compute X s (k) = X s (θ (k 1) ) and the derivatives X s (k) d 2 X s (θ (k 1) )/dθ 2. Define the update θ (k) as θ (k) = θ (k 1) + (B (k) ) 1 S (k) = dx s (θ (k 1) )/dθ and 2 X (k) s = with S (k) = s B (k) = s + s X (k) 2( ) Rs 2 X (k) s ( X (k) 2 X (k) s s s X (k) s X s (k), ) X (k) 2( ) ( Rs 2 X (k) 2 s s X (k) s X (k) s ( X (k) s ) ) 2 X (k) s. It is recommended to check that L(θ (k) ) < L(θ (k 1) ). Otherwise consider θ (k) = θ (k 1) + ρ(b (k) ) 1 S (k) for some ρ < 1, e.g. ρ = 1/2. STATDEP

7 Advantages and problems of the parametric (GARCH) modeling 1.Well developed algorithms 2. Nice asymptotic theory. Root-n consistence and asymptotic normality of the estimator θ under mild regularity assumptions. 3.Good in-sample properties. 4.Possibility to mimic the important stylized facts of the financial time series like volatility clustering, leptokurtic returns and excess kurtosis etc. Problems: 1. the parameter estimates show quite high variability. Especially estimation of β requires about 500 observations. 2. The parametric structure and stationarity of the process is crucial. If the GARCH assumption is violated, the MLE estimator θ is often completely misspecified. STATDEP

8 Example: Change point GARCH Estimates of ω Estimates of α Estimates of β ω 1.2 Median estimate.05/.95 quantiles.25/.75 quantiles true parameter α β time time time The true parameters (red) and the pointwise quantiles of the MLE s ω t, α t, β t obtained from the last 500 observations before t for the change point GARCH model. STATDEP

9 Models with time varying coefficients Returns R 1,..., R T follow R s = σ s ε s with ε s N (0, 1). Varying coefficient GARCH(1,1) model: X s = σ 2 s fulfills X s = Ψ s θ s = ω s + α s Y s 1 + β s X s 1. Here Y s = R 2 s and ω s, α s, β s some functions of time (predictable processes). Examples: Change point GARCH: θ s piecewise constant, Mikosch and Starica (2000). Smooth transition models: θ s, a smooth function of s, Fan, Jiang, Zhang, Zhou (2001). θ s piecewise smooth. We only assume that for every t there exists a time interval in which θ s θ t, Polzehl and Spokoiny (2000,2002,2003), Mercurio and Spokoiny (2004). STATDEP

10 Nonparametric modeling Define Θ = (θ t, t T ). The process X t (Θ) is described by X s (Θ) = Ψ s (Θ)θ s = ω s + α s Y s 1 + β s X s 1 (Θ) with Ψ s (Θ) = (1, Y s 1, X s 1 (Θ)). Leads to the log likelihood L(Θ) = t l(r t, X t (Θ)) with l(u, σ 2 ) = log φ(u/σ). Goal: maximize L(Θ) over the class of all admissible parameter processes: Θ = argmax Θ L(Θ). STATDEP

11 Local perturbation approach Problem: a local change of the parameter process Θ yields a global change of the latent process X(Θ). As a consequence, local likelihood approach does not apply directly. Idea: optimize L(Θ) by local perturbations of the process L(Θ). A local neighborhood of every time point t is described by weights w t,s [0, 1]. Important special cases: Kernel weights: w t,s = K( t s /h), where h, a bandwidth. Windowed weights: w t,s = 1(s I t ) for some time interval I t. Suppose that a preliminary estimate Θ = (θ t ) is fixed. Define for every θ, a process X t,s (θ) = X t,s (W t, θ; Θ ), s 1, by a recurrent equation X t,s (θ) = Ψ t,s (θ) (w t,s θ + (1 w t,s )θ s), s > t 0 where Ψ t,s (θ) = (1, Y s 1, X t,s 1 (θ)). This process describes the local perturbation at t of the process Θ. STATDEP

12 Local estimation Let Θ = (θ s) t t0 be a preliminary estimate and W t = ( w t,s )t t 0, a local model at t. Define the the local (weighted) pseudo MLE θ t of θ t as θ t = argsup L(W t, θ, Θ ) = argsup l ( Y s, X t,s (θ) ) θ θ for X t,s (θ) = Ψ t,s (θ) (w t,s θ + (1 w t,s )θ s) s > t 0. Remark: X t,s (θ) also depends on Θ. θ t solves the local estimating equation X t,s (θ) ( Y s X t,s (θ) ) X t,s (θ) 2 = 0. s Can be solved by the Newton-Raphson algorithm. s STATDEP

13 Estimation by Adaptive Weights Iterative adaptive procedure 1.Start with a local estimate θ (0) t for every t with a small bandwidth h 0. 2.Use the obtained local estimates new local models W (k) t θ (k 1) t to compute the weights w (k) t,s that describe = {w (k) t,s, s [t h k, t + h k ]} for a larger bandwidth h k 3.Recompute the local estimates θ (k) t for the models W (k) t. 4.Increase h k and repeat from step 2. STATDEP

14 Procedure (0) 1. Initialization: For every t, define β t = 0 and W (0) t = {w (0) t,s } with w (0) t,s θ (0) t = argmax θ L(W (0) t, θ) and k = Iteration: for every t = 1,..., T Calculate the adaptive weights: For every point s compute the penalty { s (k) t,s = λ 1 L(W (k 1) t,, Θ (k 1) ) L(W (k 1), Compute θ (k 1) t w (k) ( ) ( t,s = K loc t s 2 /h 2 (k)) (k) k Kst s t,s, W t t θ (k 1) s, Θ (k 1) ) = {w (k) t,s, s T }. = K l ( (t s)/h 0 2 ). Set Estimation of the parameter θ t : Define for every θ the process X (k) t,s (θ) by the recurrent formula ( ) X (k) t,s (θ) = Ψ (k) t,s (θ) w (k) t,s θ + (1 w (k) t,s ) θ (k 1) t, s > t 0, with Ψ (k) t,s (θ) = (1, Y s 1, X (k) t,s 1 θ (k) t = argsup θ Θ L (k) (W (k) t, θ, Θ (k 1) ) = arginf θ Θ (k) (θ)), and define the local MLE θ t as ( ) Rs/X 2 (k) t,s (θ) + log X (k) t,s (θ) s I t }., Θ(k) = ( θ (k) t ). 3. Stopping: Increase k by 1, set h k = ah k 1. If h k h max continue with step 2. Otherwise terminate. STATDEP

15 Semiparametric modeling The varying coefficient model σ 2 s = Ψ s θ s = ω s + α s Y s 1 + β s σ 2 s 1. may be too variable. Local estimation requires relatively large local neighborhoods, hence, low sensitivity to structural changes. Semiparametric approach: the parameters ω and α may depend on time, while β is a constant. Idea: iterate the steps of nonparametric estimation of ω s and α s and global estimation of β. STATDEP

16 Semiparametric modeling Denote γ t = (ω t, α t ) and Γ = (γ t ) t t0. For a fixed β, the process X s (Γ, β) and its derivatives are defined from X s (Γ, β) = ω s + α s R 2 s 1 + βx s 1 (Γ, β). With a fixed Γ, the estimating equation for β L(Γ, β)/ β = 0. With a fixed β, for a local model W t, an approximation Γ and the parameter vector γ, define the local perturbation X s (γ) = X s (W t, γ, Γ, β) using the recurrent formula X t,s (γ) = (1 w t,s ) ( ω s + α s R 2 s 1) + wt,s ( ω + αr 2 s 1 ) + βxt,s 1 (γ), s > t 0 and obtain the local estimate γ t by optimizing the corresponding L(γ) = L(W t, γ, Γ, β) w.r.t. γ. STATDEP

17 Local constant volatility modeling An important special case of the previous model with α s 0 and β s 0 : σ 2 s = ω s that is, σ 2 s coincides with the function (process) ω s. Local time homogeneity means that ω s is locally constant, see Mercurio and Spokoiny (2004). For the local constant case there exists an rigorous theory, see Mercurio and Spokoiny (2004) and Polzehl and Spokoiny (2002). Some properties: If σ 2 s σ 2, then the final estimate σ 2 s coincides with a high probability with the global MLE σ 2. The procedure delivers the optimal (in rate) sensitivity to the change points. If σ s is smooth in s, then we obtain the optimal nonparametric rate of estimation. STATDEP

18 Simulated results We use a set of six artificial examples to illustrate the predictive performance parametric, non- and semiparametric GARCH(1,1) models and the local constant volatility model from Polzehl and Spokoiny (2003). The sample size is set to T = Parameters are local constant and may change every 125 observations. 1.a parametric GARCH(1,1) model with common parameters ( β = 0.8 ), 2.a local constant constant volatility model ( β = 0.0 ), 3.semiparametric GARCH(1,1) models with small β, 4.entirely nonparametric GARCH with small β, 5.semiparametric GARCH(1,1) models with large β, 6.entirely nonparametric GARCH with large β. STATDEP

19 Parameters of simulated examples Ex 1: parametric GARCH(1,1) Ex 2: local const. volatility Ex 3: SP GARCH(1,1) parameter value ω t α t β t parameter value ω t α t β t parameter value ω t α t β t Index Index Index Ex 4: NP GARCH(1,1) Ex 5: NP GARCH(1,1) Ex 6: SP GARCH(1,1) parameter value ω t α t β t parameter value ω t α t β t parameter value ω t α t β t Index Index Index Parameters of simulated examples as functions of time. Procedure: All estimates are computed sequentially based on the observations from the past. For the parametric GARCH(1,1) model the last 250 observations are used. STATDEP

20 Criteria to compare the behaviour of the estimates: Mean estimated value of β t>250 A predictive likelihood risk PL(k) with horizon k = 10 ( ) 1 T k k PL(k) = log (T k 250)k X t+s t + X t+s X t+s t t=251 where X t+s t - the predicted volatility at t + s from R s for s t. Mean probability P EVaR (δ, k) of exceeding VaR at level δ and time horizon k ( 1 T k t+k ) P EVaR (δ, k) = P R s < VaR t (δ, k). T k 250 t=251 β t s=1 s=t+1 where Value at Risk (VaR) at level δ and time horizon k is defined as k VaR t (δ, k) = q δ s=1 X t+s t with q δ denoting the δ -quantile of the standard Gaussian distribution. a mean VaR at level δ and time horizon k MVaR(δ, k) = 1 T k 250 T k t=251 VaR t (δ, k) STATDEP

21 Results for the simulated examples Ex 1 Ex 2 Ex 3 Ex 4 Ex 5 Ex 6 Mean β Mean β Seq. GARCH Mean β Seq. NP-GARCH Mean β Seq. SP-GARCH PL(10) Seq. GARCH PL(10) Seq. SP-GARCH PL(10) Seq. NP-GARCH PL(10) Seq. loc.volat P EVaR (0.01, 10) Seq. GARCH P EVaR (0.01, 10) Seq. NP-GARCH P EVaR (0.01, 10) Seq. SP-GARCH P EVaR (0.01, 10) Seq. loc. volat MVaR(0.01, 10) Seq. GARCH MVaR(0.01, 10) Seq. NP-GARCH MVaR(0.01, 10) Seq. SP-GARCH MVaR(0.01, 10) Seq. loc. volat Simulation results for artificial examples 1-6. Simulation size 10. Mean estimated values of beta, mean predictive likelihood, probability of exceeding the Value at Risk and mean Value at Risk obtained for the sequential GARCH(1,1) model (last 250 observations), AWS for nonparametric GARCH(1,1), AWS for semiparametric GARCH(1,1) and using the local constant volatility AWS procedure. STATDEP

22 Applications to financial time series We now apply our methodology to two time series, the German DAX index (August 1991 to July 2003) and the US$-GBP exchange rate (January 1990 to December 2000). Criteria: 1.the predictive empirical likelihood risk PEL(k) at different time horizons k ranging from 2 weeks to half a year. ( ) 1 T k k PEL(k) = log (T k 250)k X t+s t + R2 t+s X t+s t t=251 where X t+s t means the predicted volatility at t + s from R s for s t. 2.the excess probability P EVaR 3.the mean Value at Risk MVaR 4.tail index estimate of the standardized returns s=1 STATDEP

23 Analysis of German DAX DAX : Logarithmic returns Volatility rdax SQRT of volatility process Parametric GARCH(1,1) SQRT of volatility process AWS for Nonparametric GARCH(1,1) Volatility Volatility SQRT of volatility process AWS for Semiparametric GARCH(1,1) SQRT of volatility process local constant AWS (Volatility) Volatility DAX: Logarithmic returns (top) and estimated volatility processes. Given are global estimates (red) and sequential estimates (obtained from the last 250 observations, green) by parametric GARCH(1,1), AWS for nonparametric GARCH(1,1), AWS for semiparametric GARCH(1,1) and the local constant volatility model (from top to bottom) STATDEP

24 Cointegration in DAX: fact or artifact? Mean value of the estimate β t of the parameter β t using the parametric GARCH(1,1) model and its non- and semiparametric generalizations: GARCH(1,1) NP-GARCH(1,1) SP-GARCH(1,1) Global Sequential STATDEP

25 DAX Applications: Analysis of standardized returns Model assumption R t = σ t ε t with ε t N (0, 1). Define ε t = R t / σ t. DAX: of R t 2 DAX: of R 2 ^ t X t (Seq. GARCH(1,1)) DAX: of R 2 ^ t X t (Seq. NP GARCH(1,1)) DAX: of R 2 ^ t X t (Seq. SP GARCH(1,1)) DAX: of R 2 ^ t X t (Seq. local const. AWS) 2 DAX: of R t DAX: of R 2 ^ t X t (GARCH(1,1)) DAX: of R 2 ^ t X t (NP GARCH(1,1)) DAX: of R 2 ^ t X t (SP GARCH(1,1)) DAX: of R 2 ^ t X t (local const. AWS) of squared log returns and squared standardized residuals obtained for the four methods using sequential (top) and global (bottom) volatility estimates, respectively. Dotted straight line - the 95% significant level. STATDEP

26 Heavy tail behaviour The estimated tail index for the returns and standardized returns ε t : log residuals residuals residuals local const. returns GARCH NP-GARCH SP-GARCH Volatility Global Sequential Tail index of absolute logarithmic returns and standardized residuals. Critical values for Gaussian distributions with same sample size: (.95), (.99). STATDEP

27 Out-of-sample performance Predictive empirical likelihood PEL(k) for four different time horizons ranging from two weeks to half a year ( ) 1 T k k PEL(k) = log (T k 250)k X t+s t + R2 t+s. X t+s t t=251 Method two weeks one month three months six months GARCH(1,1) NP-GARCH(1,1) SP-GARCH(1,1) local const Mean predictive likelihood for different forecast horizons. The best result for each time horizon in blue. s=1 STATDEP

28 DAX: Value-at-Risk performance 1. The excess probability P EVaR Level GARCH(1,1) NP-GARCH(1,1) SP-GARCH(1,1) local const. volatility Probability to exceed the Value at Risk at 10 trading days. The best result in blue. 2. Mean VaR at level δ = 0.01 and δ = 0.05 and time horizon k = 10 MVaR(δ, k) = 1 T k 250 T k t=251 VaR t (δ, k) Level GARCH(1,1) NP-GARCH(1,1) SP-GARCH(1,1) local const. volatility Value at Risk at 10 trading days. The best result in blue. STATDEP

29 Analysis of the US$-GBP exchange rate US$ GBP Exchange rate : Logarithmic returns rusdgbp SQRT of volatility process Parametric GARCH(1,1) Volatility Volatility Volatility estimated sequential (n=250) estimated sequential (n=250) estimated sequential (n=250) SQRT of volatility process AWS for Nonparametric GARCH(1,1) SQRT of volatility process AWS for Semiparametric GARCH(1,1) SQRT of volatility process local constant AWS (Volatility) Volatility estimated sequential (n=250) DAX: Logarithmic returns (top) and estimated volatility processes. Given are global estimates (red) and sequential estimates (obtained from the last 250 observations, green) by parametric GARCH(1,1), AWS for nonparametric GARCH(1,1), AWS for semiparametric GARCH(1,1) and the local constant volatility model (from top to bottom) STATDEP

30 Cointegration in US$-GBP exchange rate: fact or artifact? Mean value of the estimate β t of the parameter β t using the parametric GARCH(1,1) model and its non- and semiparametric generalizations: GARCH(1,1) NP-GARCH(1,1) SP-GARCH(1,1) Global Sequential STATDEP

31 US$-GBP exchange rate applications: Analysis of standardized returns Model assumption Y t = σ t ε t with ε t N (0, 1). Define ε t = Y t / σ t. DAX: of R t 2 USD GBP: of R 2 ^ t X t (Seq. GARCH(1,1)) USD GBP: of R 2 ^ t X t (Seq. NP GARCH(1,1 USD GBP: of R 2 ^ t X t (Seq. SP GARCH(1,1 USD GBP: of R 2 ^ t X t (Seq. local const. AWS 2 DAX: of R t USD GBP: of R 2 ^ t X t (GARCH(1,1)) USD GBP: of R 2 ^ t X t (NP GARCH(1,1)) USD GBP: of R 2 ^ t X t (SP GARCH(1,1)) USD GBP: of R 2 ^ t X t (local const. AWS) of squared log returns and squared standardized residuals obtained for the four methods using sequential (top) and global (bottom) volatility estimates, respectively. Dotted straight line - the 95% significant level. STATDEP

32 US$-GBP exchange rate: Heavy tail behaviour The estimated tail index for the returns and standardized returns ε t : log residuals residuals residuals local const. returns GARCH NP-GARCH SP-GARCH Volatility Global Sequential Tail index of absolute logarithmic returns and standardized residuals. Critical values for Gaussian distributions with same sample size: (.95), (.99). STATDEP

33 US$-GBP exchange rate: Out-of-sample performance Predictive empirical likelihood PEL(k) for four different time horizons ranging from two weeks to half a year ( ) 1 T k k PEL(k) = log (T k 250)k X t+s t + R2 t+s. X t+s t t=251 Method two weeks one month three months six months GARCH(1,1) NP-GARCH(1,1) SP-GARCH(1,1) local const. Volatility Mean predictive likelihood for different forecast horizons. The best result for each time horizon in blue. s=1 STATDEP

34 US$-GBP exchange rate: Value-at-Risk performance 1. The excess probability P EVaR Level GARCH(1,1) NP-GARCH(1,1) SP-GARCH(1,1) local const. volatility Probability to exceed the Value at Risk at 10 trading days. The best result in blue. 2. Mean VaR at level δ = 0.01 and δ = 0.05 and time horizon k = 10 MVaR(δ, k) = 1 T k 250 T k t=251 VaR t (δ, k) Level GARCH(1,1) NP-GARCH(1,1) SP-GARCH(1,1) local const. volatility Value at Risk at 10 trading days. The best result in blue. STATDEP

35 Conclusion and Outlook The new approach can be applied and studied in a unified way for a wide class of different models presents a consistent way of selecting the tuning parameters demonstrates a very reasonable numerical performance the simplest local constant modeling is slightly preferable as far as the in sample properties or short time ahead forecasting is concerned. STATDEP

36 Extension: Parametric generalized linear modeling Model: R t = σ t ε t and the process X t = g(σt 2 ) satisfies the linear equation: X t = ω + α 1 Y t α p Y t p + β 1 X t β q X t q where Y t = h(r t ) for some fixed functions g, h. Example 1. ARCH(p) X t = σt 2 and Y t = Rt 2 : σt 2 = ω + α 1 Rt α p Rt p. 2 Example 2. GARCH(p,q) X t = σ 2 t and Y t = R 2 t X t = ω + α 1 R 2 t α p R 2 t p + β 1 σ 2 t β q σ 2 t q Example 3. Exponential (G)ARCH(p) X t = log σ 2 t and Y t = h(r t ) where h(u) = log(u 2 + ɛ) or h(u) = log log(e u + e + ɛ) for some ɛ 0 : log σ 2 t = ω + α 1 h(r t 1 ) α p h(y t p ) + β 1 log σ 2 t β q log σ 2 t q STATDEP

Course information FN3142 Quantitative finance

Course information FN3142 Quantitative finance Course information 015 16 FN314 Quantitative finance This course is aimed at students interested in obtaining a thorough grounding in market finance and related empirical methods. Prerequisite If taken

More information

Amath 546/Econ 589 Univariate GARCH Models

Amath 546/Econ 589 Univariate GARCH Models Amath 546/Econ 589 Univariate GARCH Models Eric Zivot April 24, 2013 Lecture Outline Conditional vs. Unconditional Risk Measures Empirical regularities of asset returns Engle s ARCH model Testing for ARCH

More information

Statistical Inference and Methods

Statistical Inference and Methods Department of Mathematics Imperial College London d.stephens@imperial.ac.uk http://stats.ma.ic.ac.uk/ das01/ 14th February 2006 Part VII Session 7: Volatility Modelling Session 7: Volatility Modelling

More information

Modelling financial data with stochastic processes

Modelling financial data with stochastic processes Modelling financial data with stochastic processes Vlad Ardelean, Fabian Tinkl 01.08.2012 Chair of statistics and econometrics FAU Erlangen-Nuremberg Outline Introduction Stochastic processes Volatility

More information

12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006.

12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006. 12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006. References for this Lecture: Robert F. Engle. Autoregressive Conditional Heteroscedasticity with Estimates of Variance

More information

Conditional Heteroscedasticity

Conditional Heteroscedasticity 1 Conditional Heteroscedasticity May 30, 2010 Junhui Qian 1 Introduction ARMA(p,q) models dictate that the conditional mean of a time series depends on past observations of the time series and the past

More information

Risk Management and Time Series

Risk Management and Time Series IEOR E4602: Quantitative Risk Management Spring 2016 c 2016 by Martin Haugh Risk Management and Time Series Time series models are often employed in risk management applications. They can be used to estimate

More information

ARCH and GARCH models

ARCH and GARCH models ARCH and GARCH models Fulvio Corsi SNS Pisa 5 Dic 2011 Fulvio Corsi ARCH and () GARCH models SNS Pisa 5 Dic 2011 1 / 21 Asset prices S&P 500 index from 1982 to 2009 1600 1400 1200 1000 800 600 400 200

More information

Booth School of Business, University of Chicago Business 41202, Spring Quarter 2014, Mr. Ruey S. Tsay. Solutions to Midterm

Booth School of Business, University of Chicago Business 41202, Spring Quarter 2014, Mr. Ruey S. Tsay. Solutions to Midterm Booth School of Business, University of Chicago Business 41202, Spring Quarter 2014, Mr. Ruey S. Tsay Solutions to Midterm Problem A: (30 pts) Answer briefly the following questions. Each question has

More information

Time series: Variance modelling

Time series: Variance modelling Time series: Variance modelling Bernt Arne Ødegaard 5 October 018 Contents 1 Motivation 1 1.1 Variance clustering.......................... 1 1. Relation to heteroskedasticity.................... 3 1.3

More information

Financial Risk Management

Financial Risk Management Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #4 1 Correlation and copulas 1. The bivariate Gaussian copula is given

More information

Financial Econometrics Jeffrey R. Russell. Midterm 2014 Suggested Solutions. TA: B. B. Deng

Financial Econometrics Jeffrey R. Russell. Midterm 2014 Suggested Solutions. TA: B. B. Deng Financial Econometrics Jeffrey R. Russell Midterm 2014 Suggested Solutions TA: B. B. Deng Unless otherwise stated, e t is iid N(0,s 2 ) 1. (12 points) Consider the three series y1, y2, y3, and y4. Match

More information

Chapter 6 Forecasting Volatility using Stochastic Volatility Model

Chapter 6 Forecasting Volatility using Stochastic Volatility Model Chapter 6 Forecasting Volatility using Stochastic Volatility Model Chapter 6 Forecasting Volatility using SV Model In this chapter, the empirical performance of GARCH(1,1), GARCH-KF and SV models from

More information

Supplemental Online Appendix to Han and Hong, Understanding In-House Transactions in the Real Estate Brokerage Industry

Supplemental Online Appendix to Han and Hong, Understanding In-House Transactions in the Real Estate Brokerage Industry Supplemental Online Appendix to Han and Hong, Understanding In-House Transactions in the Real Estate Brokerage Industry Appendix A: An Agent-Intermediated Search Model Our motivating theoretical framework

More information

Adaptive Interest Rate Modelling

Adaptive Interest Rate Modelling Modelling Mengmeng Guo Wolfgang Karl Härdle Ladislaus von Bortkiewicz Chair of Statistics C.A.S.E. - Center for Applied Statistics and Economics Humboldt-Universität zu Berlin http://lvb.wiwi.hu-berlin.de

More information

Booth School of Business, University of Chicago Business 41202, Spring Quarter 2010, Mr. Ruey S. Tsay. Solutions to Midterm

Booth School of Business, University of Chicago Business 41202, Spring Quarter 2010, Mr. Ruey S. Tsay. Solutions to Midterm Booth School of Business, University of Chicago Business 41202, Spring Quarter 2010, Mr. Ruey S. Tsay Solutions to Midterm Problem A: (30 pts) Answer briefly the following questions. Each question has

More information

ESTIMATION OF UTILITY FUNCTIONS: MARKET VS. REPRESENTATIVE AGENT THEORY

ESTIMATION OF UTILITY FUNCTIONS: MARKET VS. REPRESENTATIVE AGENT THEORY ESTIMATION OF UTILITY FUNCTIONS: MARKET VS. REPRESENTATIVE AGENT THEORY Kai Detlefsen Wolfgang K. Härdle Rouslan A. Moro, Deutsches Institut für Wirtschaftsforschung (DIW) Center for Applied Statistics

More information

Financial Econometrics

Financial Econometrics Financial Econometrics Volatility Gerald P. Dwyer Trinity College, Dublin January 2013 GPD (TCD) Volatility 01/13 1 / 37 Squared log returns for CRSP daily GPD (TCD) Volatility 01/13 2 / 37 Absolute value

More information

Financial Risk Forecasting Chapter 9 Extreme Value Theory

Financial Risk Forecasting Chapter 9 Extreme Value Theory Financial Risk Forecasting Chapter 9 Extreme Value Theory Jon Danielsson 2017 London School of Economics To accompany Financial Risk Forecasting www.financialriskforecasting.com Published by Wiley 2011

More information

On modelling of electricity spot price

On modelling of electricity spot price , Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction

More information

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2009, Mr. Ruey S. Tsay. Solutions to Final Exam

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2009, Mr. Ruey S. Tsay. Solutions to Final Exam The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2009, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (42 pts) Answer briefly the following questions. 1. Questions

More information

1. You are given the following information about a stationary AR(2) model:

1. You are given the following information about a stationary AR(2) model: Fall 2003 Society of Actuaries **BEGINNING OF EXAMINATION** 1. You are given the following information about a stationary AR(2) model: (i) ρ 1 = 05. (ii) ρ 2 = 01. Determine φ 2. (A) 0.2 (B) 0.1 (C) 0.4

More information

A Hidden Markov Model Approach to Information-Based Trading: Theory and Applications

A Hidden Markov Model Approach to Information-Based Trading: Theory and Applications A Hidden Markov Model Approach to Information-Based Trading: Theory and Applications Online Supplementary Appendix Xiangkang Yin and Jing Zhao La Trobe University Corresponding author, Department of Finance,

More information

Lecture 5a: ARCH Models

Lecture 5a: ARCH Models Lecture 5a: ARCH Models 1 2 Big Picture 1. We use ARMA model for the conditional mean 2. We use ARCH model for the conditional variance 3. ARMA and ARCH model can be used together to describe both conditional

More information

Lecture 17: More on Markov Decision Processes. Reinforcement learning

Lecture 17: More on Markov Decision Processes. Reinforcement learning Lecture 17: More on Markov Decision Processes. Reinforcement learning Learning a model: maximum likelihood Learning a value function directly Monte Carlo Temporal-difference (TD) learning COMP-424, Lecture

More information

Market risk measurement in practice

Market risk measurement in practice Lecture notes on risk management, public policy, and the financial system Allan M. Malz Columbia University 2018 Allan M. Malz Last updated: October 23, 2018 2/32 Outline Nonlinearity in market risk Market

More information

Chapter 5 Univariate time-series analysis. () Chapter 5 Univariate time-series analysis 1 / 29

Chapter 5 Univariate time-series analysis. () Chapter 5 Univariate time-series analysis 1 / 29 Chapter 5 Univariate time-series analysis () Chapter 5 Univariate time-series analysis 1 / 29 Time-Series Time-series is a sequence fx 1, x 2,..., x T g or fx t g, t = 1,..., T, where t is an index denoting

More information

Monetary Economics Final Exam

Monetary Economics Final Exam 316-466 Monetary Economics Final Exam 1. Flexible-price monetary economics (90 marks). Consider a stochastic flexibleprice money in the utility function model. Time is discrete and denoted t =0, 1,...

More information

IMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS

IMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS IMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS DENIS BELOMESTNY AND MARKUS REISS 1. Introduction The aim of this report is to describe more precisely how the spectral calibration method

More information

4 Reinforcement Learning Basic Algorithms

4 Reinforcement Learning Basic Algorithms Learning in Complex Systems Spring 2011 Lecture Notes Nahum Shimkin 4 Reinforcement Learning Basic Algorithms 4.1 Introduction RL methods essentially deal with the solution of (optimal) control problems

More information

Model Construction & Forecast Based Portfolio Allocation:

Model Construction & Forecast Based Portfolio Allocation: QBUS6830 Financial Time Series and Forecasting Model Construction & Forecast Based Portfolio Allocation: Is Quantitative Method Worth It? Members: Bowei Li (303083) Wenjian Xu (308077237) Xiaoyun Lu (3295347)

More information

Web-based Supplementary Materials for. A space-time conditional intensity model. for invasive meningococcal disease occurence

Web-based Supplementary Materials for. A space-time conditional intensity model. for invasive meningococcal disease occurence Web-based Supplementary Materials for A space-time conditional intensity model for invasive meningococcal disease occurence by Sebastian Meyer 1,2, Johannes Elias 3, and Michael Höhle 4,2 1 Department

More information

Chapter 7: Estimation Sections

Chapter 7: Estimation Sections 1 / 40 Chapter 7: Estimation Sections 7.1 Statistical Inference Bayesian Methods: Chapter 7 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions 7.4 Bayes Estimators Frequentist Methods:

More information

Estimating Pricing Kernel via Series Methods

Estimating Pricing Kernel via Series Methods Estimating Pricing Kernel via Series Methods Maria Grith Wolfgang Karl Härdle Melanie Schienle Ladislaus von Bortkiewicz Chair of Statistics Chair of Econometrics C.A.S.E. Center for Applied Statistics

More information

Analysis of truncated data with application to the operational risk estimation

Analysis of truncated data with application to the operational risk estimation Analysis of truncated data with application to the operational risk estimation Petr Volf 1 Abstract. Researchers interested in the estimation of operational risk often face problems arising from the structure

More information

Extend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty

Extend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty Extend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty George Photiou Lincoln College University of Oxford A dissertation submitted in partial fulfilment for

More information

Nonparametric Risk Management with Generalized Hyperbolic Distributions

Nonparametric Risk Management with Generalized Hyperbolic Distributions Nonparametric Risk Management with Generalized Hyperbolic Distributions Ying Chen Wolfgang Härdle Center for Applied Statistics and Economics Institut für Statistik and Ökonometrie Humboldt-Universität

More information

SADDLE POINT APPROXIMATION AND VOLATILITY ESTIMATION OF VALUE-AT-RISK

SADDLE POINT APPROXIMATION AND VOLATILITY ESTIMATION OF VALUE-AT-RISK Statistica Sinica 20 (2010), 1239-1256 SADDLE POINT APPROXIMATION AND VOLATILITY ESTIMATION OF VALUE-AT-RISK Maozai Tian and Ngai Hang Chan Renmin University of China and Chinese University of Hong Kong

More information

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay. Solutions to Final Exam

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay. Solutions to Final Exam The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (40 points) Answer briefly the following questions. 1. Consider

More information

A gentle introduction to the RM 2006 methodology

A gentle introduction to the RM 2006 methodology A gentle introduction to the RM 2006 methodology Gilles Zumbach RiskMetrics Group Av. des Morgines 12 1213 Petit-Lancy Geneva, Switzerland gilles.zumbach@riskmetrics.com Initial version: August 2006 This

More information

Occasional Paper. Risk Measurement Illiquidity Distortions. Jiaqi Chen and Michael L. Tindall

Occasional Paper. Risk Measurement Illiquidity Distortions. Jiaqi Chen and Michael L. Tindall DALLASFED Occasional Paper Risk Measurement Illiquidity Distortions Jiaqi Chen and Michael L. Tindall Federal Reserve Bank of Dallas Financial Industry Studies Department Occasional Paper 12-2 December

More information

Research Article The Volatility of the Index of Shanghai Stock Market Research Based on ARCH and Its Extended Forms

Research Article The Volatility of the Index of Shanghai Stock Market Research Based on ARCH and Its Extended Forms Discrete Dynamics in Nature and Society Volume 2009, Article ID 743685, 9 pages doi:10.1155/2009/743685 Research Article The Volatility of the Index of Shanghai Stock Market Research Based on ARCH and

More information

Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS001) p approach

Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS001) p approach Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS001) p.5901 What drives short rate dynamics? approach A functional gradient descent Audrino, Francesco University

More information

CS340 Machine learning Bayesian model selection

CS340 Machine learning Bayesian model selection CS340 Machine learning Bayesian model selection Bayesian model selection Suppose we have several models, each with potentially different numbers of parameters. Example: M0 = constant, M1 = straight line,

More information

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay. Solutions to Final Exam.

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay. Solutions to Final Exam. The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (32 pts) Answer briefly the following questions. 1. Suppose

More information

A Multivariate Analysis of Intercompany Loss Triangles

A Multivariate Analysis of Intercompany Loss Triangles A Multivariate Analysis of Intercompany Loss Triangles Peng Shi School of Business University of Wisconsin-Madison ASTIN Colloquium May 21-24, 2013 Peng Shi (Wisconsin School of Business) Intercompany

More information

Analyzing Oil Futures with a Dynamic Nelson-Siegel Model

Analyzing Oil Futures with a Dynamic Nelson-Siegel Model Analyzing Oil Futures with a Dynamic Nelson-Siegel Model NIELS STRANGE HANSEN & ASGER LUNDE DEPARTMENT OF ECONOMICS AND BUSINESS, BUSINESS AND SOCIAL SCIENCES, AARHUS UNIVERSITY AND CENTER FOR RESEARCH

More information

Internet Appendix for Asymmetry in Stock Comovements: An Entropy Approach

Internet Appendix for Asymmetry in Stock Comovements: An Entropy Approach Internet Appendix for Asymmetry in Stock Comovements: An Entropy Approach Lei Jiang Tsinghua University Ke Wu Renmin University of China Guofu Zhou Washington University in St. Louis August 2017 Jiang,

More information

Pakes (1986): Patents as Options: Some Estimates of the Value of Holding European Patent Stocks

Pakes (1986): Patents as Options: Some Estimates of the Value of Holding European Patent Stocks Pakes (1986): Patents as Options: Some Estimates of the Value of Holding European Patent Stocks Spring 2009 Main question: How much are patents worth? Answering this question is important, because it helps

More information

Forecasting jumps in conditional volatility The GARCH-IE model

Forecasting jumps in conditional volatility The GARCH-IE model Forecasting jumps in conditional volatility The GARCH-IE model Philip Hans Franses and Marco van der Leij Econometric Institute Erasmus University Rotterdam e-mail: franses@few.eur.nl 1 Outline of presentation

More information

Modelling of Long-Term Risk

Modelling of Long-Term Risk Modelling of Long-Term Risk Roger Kaufmann Swiss Life roger.kaufmann@swisslife.ch 15th International AFIR Colloquium 6-9 September 2005, Zurich c 2005 (R. Kaufmann, Swiss Life) Contents A. Basel II B.

More information

6. Genetics examples: Hardy-Weinberg Equilibrium

6. Genetics examples: Hardy-Weinberg Equilibrium PBCB 206 (Fall 2006) Instructor: Fei Zou email: fzou@bios.unc.edu office: 3107D McGavran-Greenberg Hall Lecture 4 Topics for Lecture 4 1. Parametric models and estimating parameters from data 2. Method

More information

Problem set 5. Asset pricing. Markus Roth. Chair for Macroeconomics Johannes Gutenberg Universität Mainz. Juli 5, 2010

Problem set 5. Asset pricing. Markus Roth. Chair for Macroeconomics Johannes Gutenberg Universität Mainz. Juli 5, 2010 Problem set 5 Asset pricing Markus Roth Chair for Macroeconomics Johannes Gutenberg Universität Mainz Juli 5, 200 Markus Roth (Macroeconomics 2) Problem set 5 Juli 5, 200 / 40 Contents Problem 5 of problem

More information

Two hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER

Two hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER Two hours MATH20802 To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER STATISTICAL METHODS Answer any FOUR of the SIX questions.

More information

Dependence Structure and Extreme Comovements in International Equity and Bond Markets

Dependence Structure and Extreme Comovements in International Equity and Bond Markets Dependence Structure and Extreme Comovements in International Equity and Bond Markets René Garcia Edhec Business School, Université de Montréal, CIRANO and CIREQ Georges Tsafack Suffolk University Measuring

More information

Value at Risk with Stable Distributions

Value at Risk with Stable Distributions Value at Risk with Stable Distributions Tecnológico de Monterrey, Guadalajara Ramona Serrano B Introduction The core activity of financial institutions is risk management. Calculate capital reserves given

More information

Window Width Selection for L 2 Adjusted Quantile Regression

Window Width Selection for L 2 Adjusted Quantile Regression Window Width Selection for L 2 Adjusted Quantile Regression Yoonsuh Jung, The Ohio State University Steven N. MacEachern, The Ohio State University Yoonkyung Lee, The Ohio State University Technical Report

More information

The Use of Penultimate Approximations in Risk Management

The Use of Penultimate Approximations in Risk Management The Use of Penultimate Approximations in Risk Management www.math.ethz.ch/ degen (joint work with P. Embrechts) 6th International Conference on Extreme Value Analysis Fort Collins CO, June 26, 2009 Penultimate

More information

Assicurazioni Generali: An Option Pricing Case with NAGARCH

Assicurazioni Generali: An Option Pricing Case with NAGARCH Assicurazioni Generali: An Option Pricing Case with NAGARCH Assicurazioni Generali: Business Snapshot Find our latest analyses and trade ideas on bsic.it Assicurazioni Generali SpA is an Italy-based insurance

More information

DYNAMIC ECONOMETRIC MODELS Vol. 8 Nicolaus Copernicus University Toruń Mateusz Pipień Cracow University of Economics

DYNAMIC ECONOMETRIC MODELS Vol. 8 Nicolaus Copernicus University Toruń Mateusz Pipień Cracow University of Economics DYNAMIC ECONOMETRIC MODELS Vol. 8 Nicolaus Copernicus University Toruń 2008 Mateusz Pipień Cracow University of Economics On the Use of the Family of Beta Distributions in Testing Tradeoff Between Risk

More information

Portfolio Optimization. Prof. Daniel P. Palomar

Portfolio Optimization. Prof. Daniel P. Palomar Portfolio Optimization Prof. Daniel P. Palomar The Hong Kong University of Science and Technology (HKUST) MAFS6010R- Portfolio Optimization with R MSc in Financial Mathematics Fall 2018-19, HKUST, Hong

More information

Modelling volatility - ARCH and GARCH models

Modelling volatility - ARCH and GARCH models Modelling volatility - ARCH and GARCH models Beáta Stehlíková Time series analysis Modelling volatility- ARCH and GARCH models p.1/33 Stock prices Weekly stock prices (library quantmod) Continuous returns:

More information

Strategies for High Frequency FX Trading

Strategies for High Frequency FX Trading Strategies for High Frequency FX Trading - The choice of bucket size Malin Lunsjö and Malin Riddarström Department of Mathematical Statistics Faculty of Engineering at Lund University June 2017 Abstract

More information

IEOR E4602: Quantitative Risk Management

IEOR E4602: Quantitative Risk Management IEOR E4602: Quantitative Risk Management Basic Concepts and Techniques of Risk Management Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com

More information

Backtesting Trading Book Models

Backtesting Trading Book Models Backtesting Trading Book Models Using Estimates of VaR Expected Shortfall and Realized p-values Alexander J. McNeil 1 1 Heriot-Watt University Edinburgh ETH Risk Day 11 September 2015 AJM (HWU) Backtesting

More information

Value at risk might underestimate risk when risk bites. Just bootstrap it!

Value at risk might underestimate risk when risk bites. Just bootstrap it! 23 September 215 by Zhili Cao Research & Investment Strategy at risk might underestimate risk when risk bites. Just bootstrap it! Key points at Risk (VaR) is one of the most widely used statistical tools

More information

Variance clustering. Two motivations, volatility clustering, and implied volatility

Variance clustering. Two motivations, volatility clustering, and implied volatility Variance modelling The simplest assumption for time series is that variance is constant. Unfortunately that assumption is often violated in actual data. In this lecture we look at the implications of time

More information

Practice Exam 1. Loss Amount Number of Losses

Practice Exam 1. Loss Amount Number of Losses Practice Exam 1 1. You are given the following data on loss sizes: An ogive is used as a model for loss sizes. Determine the fitted median. Loss Amount Number of Losses 0 1000 5 1000 5000 4 5000 10000

More information

SUPPLEMENT TO EQUILIBRIA IN HEALTH EXCHANGES: ADVERSE SELECTION VERSUS RECLASSIFICATION RISK (Econometrica, Vol. 83, No. 4, July 2015, )

SUPPLEMENT TO EQUILIBRIA IN HEALTH EXCHANGES: ADVERSE SELECTION VERSUS RECLASSIFICATION RISK (Econometrica, Vol. 83, No. 4, July 2015, ) Econometrica Supplementary Material SUPPLEMENT TO EQUILIBRIA IN HEALTH EXCHANGES: ADVERSE SELECTION VERSUS RECLASSIFICATION RISK (Econometrica, Vol. 83, No. 4, July 2015, 1261 1313) BY BEN HANDEL, IGAL

More information

Recursive estimation of piecewise constant volatilities 1

Recursive estimation of piecewise constant volatilities 1 Recursive estimation of piecewise constant volatilities 1 by Christian Höhenrieder Deutsche Bundesbank, Berliner Allee 14 D-401 Düsseldorf, Germany, Laurie Davies Fakultät Mathematik, Universität Duisburg-Essen

More information

Financial Econometrics Lecture 5: Modelling Volatility and Correlation

Financial Econometrics Lecture 5: Modelling Volatility and Correlation Financial Econometrics Lecture 5: Modelling Volatility and Correlation Dayong Zhang Research Institute of Economics and Management Autumn, 2011 Learning Outcomes Discuss the special features of financial

More information

Homework Problems Stat 479

Homework Problems Stat 479 Chapter 2 1. Model 1 is a uniform distribution from 0 to 100. Determine the table entries for a generalized uniform distribution covering the range from a to b where a < b. 2. Let X be a discrete random

More information

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay. Solutions to Final Exam

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay. Solutions to Final Exam The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (40 points) Answer briefly the following questions. 1. Describe

More information

U n i ve rs i t y of He idelberg

U n i ve rs i t y of He idelberg U n i ve rs i t y of He idelberg Department of Economics Discussion Paper Series No. 613 On the statistical properties of multiplicative GARCH models Christian Conrad and Onno Kleen March 2016 On the statistical

More information

Data-analytic Approaches to the Estimation of Value-at-Risk

Data-analytic Approaches to the Estimation of Value-at-Risk Data-analytic Approaches to the Estimation of Value-at-Risk Jianqing Fan, and Juan Gu Department of Statistics, Chinese University of Hong Kong Shatin,N.T., Hong Kong GF Securities Co., Ltd,Guangzhou,Guangdong,China

More information

Journal of Economics and Financial Analysis, Vol:1, No:1 (2017) 1-13

Journal of Economics and Financial Analysis, Vol:1, No:1 (2017) 1-13 Journal of Economics and Financial Analysis, Vol:1, No:1 (2017) 1-13 Journal of Economics and Financial Analysis Type: Double Blind Peer Reviewed Scientific Journal Printed ISSN: 2521-6627 Online ISSN:

More information

Financial Time Series and Volatility Prediction using NoVaS Transformations

Financial Time Series and Volatility Prediction using NoVaS Transformations Financial Time Series and Volatility Prediction using NoVaS Transformations Dimitris N. Politis Dimitrios D. Thomakos December 21, 2006 Abstract We extend earlier work on the NoVaS transformation approach

More information

Forecasting Volatility of USD/MUR Exchange Rate using a GARCH (1,1) model with GED and Student s-t errors

Forecasting Volatility of USD/MUR Exchange Rate using a GARCH (1,1) model with GED and Student s-t errors UNIVERSITY OF MAURITIUS RESEARCH JOURNAL Volume 17 2011 University of Mauritius, Réduit, Mauritius Research Week 2009/2010 Forecasting Volatility of USD/MUR Exchange Rate using a GARCH (1,1) model with

More information

Modeling Capital Market with Financial Signal Processing

Modeling Capital Market with Financial Signal Processing Modeling Capital Market with Financial Signal Processing Jenher Jeng Ph.D., Statistics, U.C. Berkeley Founder & CTO of Harmonic Financial Engineering, www.harmonicfinance.com Outline Theory and Techniques

More information

Asymptotic results discrete time martingales and stochastic algorithms

Asymptotic results discrete time martingales and stochastic algorithms Asymptotic results discrete time martingales and stochastic algorithms Bernard Bercu Bordeaux University, France IFCAM Summer School Bangalore, India, July 2015 Bernard Bercu Asymptotic results for discrete

More information

John Hull, Risk Management and Financial Institutions, 4th Edition

John Hull, Risk Management and Financial Institutions, 4th Edition P1.T2. Quantitative Analysis John Hull, Risk Management and Financial Institutions, 4th Edition Bionic Turtle FRM Video Tutorials By David Harper, CFA FRM 1 Chapter 10: Volatility (Learning objectives)

More information

The Effects of Monetary Policy on Asset Price Bubbles: Some Evidence

The Effects of Monetary Policy on Asset Price Bubbles: Some Evidence The Effects of Monetary Policy on Asset Price Bubbles: Some Evidence Jordi Galí Luca Gambetti September 2013 Jordi Galí, Luca Gambetti () Monetary Policy and Bubbles September 2013 1 / 17 Monetary Policy

More information

Adaptive Monitoring of Intraday Market Data

Adaptive Monitoring of Intraday Market Data Enzo Giacomini Nikolaus Hautsch Vladimir Spokoiny CASE - Center for Applied Statistics and Economics Humboldt-Universität zu Berlin Motivation 1-2 Ultra-High Frequency Data Ultra-high frequency, Engle

More information

Two-step conditional α-quantile estimation via additive models of location and scale 1

Two-step conditional α-quantile estimation via additive models of location and scale 1 Two-step conditional α-quantile estimation via additive models of location and scale 1 Carlos Martins-Filho Department of Economics IFPRI University of Colorado 2033 K Street NW Boulder, CO 80309-0256,

More information

Short-selling constraints and stock-return volatility: empirical evidence from the German stock market

Short-selling constraints and stock-return volatility: empirical evidence from the German stock market Short-selling constraints and stock-return volatility: empirical evidence from the German stock market Martin Bohl, Gerrit Reher, Bernd Wilfling Westfälische Wilhelms-Universität Münster Contents 1. Introduction

More information

Modeling dynamic diurnal patterns in high frequency financial data

Modeling dynamic diurnal patterns in high frequency financial data Modeling dynamic diurnal patterns in high frequency financial data Ryoko Ito 1 Faculty of Economics, Cambridge University Email: ri239@cam.ac.uk Website: www.itoryoko.com This paper: Cambridge Working

More information

SYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data

SYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data SYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu September 5, 2015

More information

MAFS Computational Methods for Pricing Structured Products

MAFS Computational Methods for Pricing Structured Products MAFS550 - Computational Methods for Pricing Structured Products Solution to Homework Two Course instructor: Prof YK Kwok 1 Expand f(x 0 ) and f(x 0 x) at x 0 into Taylor series, where f(x 0 ) = f(x 0 )

More information

The method of Maximum Likelihood.

The method of Maximum Likelihood. Maximum Likelihood The method of Maximum Likelihood. In developing the least squares estimator - no mention of probabilities. Minimize the distance between the predicted linear regression and the observed

More information

An Online Algorithm for Multi-Strategy Trading Utilizing Market Regimes

An Online Algorithm for Multi-Strategy Trading Utilizing Market Regimes An Online Algorithm for Multi-Strategy Trading Utilizing Market Regimes Hynek Mlnařík 1 Subramanian Ramamoorthy 2 Rahul Savani 1 1 Warwick Institute for Financial Computing Department of Computer Science

More information

Lecture 7: Bayesian approach to MAB - Gittins index

Lecture 7: Bayesian approach to MAB - Gittins index Advanced Topics in Machine Learning and Algorithmic Game Theory Lecture 7: Bayesian approach to MAB - Gittins index Lecturer: Yishay Mansour Scribe: Mariano Schain 7.1 Introduction In the Bayesian approach

More information

Definition 9.1 A point estimate is any function T (X 1,..., X n ) of a random sample. We often write an estimator of the parameter θ as ˆθ.

Definition 9.1 A point estimate is any function T (X 1,..., X n ) of a random sample. We often write an estimator of the parameter θ as ˆθ. 9 Point estimation 9.1 Rationale behind point estimation When sampling from a population described by a pdf f(x θ) or probability function P [X = x θ] knowledge of θ gives knowledge of the entire population.

More information

Scenario reduction and scenario tree construction for power management problems

Scenario reduction and scenario tree construction for power management problems Scenario reduction and scenario tree construction for power management problems N. Gröwe-Kuska, H. Heitsch and W. Römisch Humboldt-University Berlin Institute of Mathematics Page 1 of 20 IEEE Bologna POWER

More information

Chapter 7: Estimation Sections

Chapter 7: Estimation Sections 1 / 31 : Estimation Sections 7.1 Statistical Inference Bayesian Methods: 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions 7.4 Bayes Estimators Frequentist Methods: 7.5 Maximum Likelihood

More information

An overview of some financial models using BSDE with enlarged filtrations

An overview of some financial models using BSDE with enlarged filtrations An overview of some financial models using BSDE with enlarged filtrations Anne EYRAUD-LOISEL Workshop : Enlargement of Filtrations and Applications to Finance and Insurance May 31st - June 4th, 2010, Jena

More information

Effects of Outliers and Parameter Uncertainties in Portfolio Selection

Effects of Outliers and Parameter Uncertainties in Portfolio Selection Effects of Outliers and Parameter Uncertainties in Portfolio Selection Luiz Hotta 1 Carlos Trucíos 2 Esther Ruiz 3 1 Department of Statistics, University of Campinas. 2 EESP-FGV (postdoctoral). 3 Department

More information

Cross-Sectional Distribution of GARCH Coefficients across S&P 500 Constituents : Time-Variation over the Period

Cross-Sectional Distribution of GARCH Coefficients across S&P 500 Constituents : Time-Variation over the Period Cahier de recherche/working Paper 13-13 Cross-Sectional Distribution of GARCH Coefficients across S&P 500 Constituents : Time-Variation over the Period 2000-2012 David Ardia Lennart F. Hoogerheide Mai/May

More information

Implementing an Agent-Based General Equilibrium Model

Implementing an Agent-Based General Equilibrium Model Implementing an Agent-Based General Equilibrium Model 1 2 3 Pure Exchange General Equilibrium We shall take N dividend processes δ n (t) as exogenous with a distribution which is known to all agents There

More information

Chapter 7: Estimation Sections

Chapter 7: Estimation Sections Chapter 7: Estimation Sections 7.1 Statistical Inference Bayesian Methods: 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions Frequentist Methods: 7.5 Maximum Likelihood Estimators

More information

A potentially useful approach to model nonlinearities in time series is to assume different behavior (structural break) in different subsamples

A potentially useful approach to model nonlinearities in time series is to assume different behavior (structural break) in different subsamples 1.3 Regime switching models A potentially useful approach to model nonlinearities in time series is to assume different behavior (structural break) in different subsamples (or regimes). If the dates, the

More information