Value at Risk and Self Similarity

Similar documents
Value at Risk and Self Similarity

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

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

STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL

Introduction to Algorithmic Trading Strategies Lecture 8

Case Study: Heavy-Tailed Distribution and Reinsurance Rate-making

BROWNIAN MOTION Antonella Basso, Martina Nardon

Homework Assignments

An Introduction to Stochastic Calculus

Modelling financial data with stochastic processes

Module 10:Application of stochastic processes in areas like finance Lecture 36:Black-Scholes Model. Stochastic Differential Equation.

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

Calculating VaR. There are several approaches for calculating the Value at Risk figure. The most popular are the

Financial Engineering. Craig Pirrong Spring, 2006

Math 416/516: Stochastic Simulation

CHOICE THEORY, UTILITY FUNCTIONS AND RISK AVERSION

IEOR E4602: Quantitative Risk Management

Recent Advances in Fractional Stochastic Volatility Models

Characterization of the Optimum

Probability. An intro for calculus students P= Figure 1: A normal integral

A new approach for scenario generation in risk management

Mathematics of Finance Final Preparation December 19. To be thoroughly prepared for the final exam, you should

Fractional Black - Scholes Equation

M5MF6. Advanced Methods in Derivatives Pricing

Advanced Financial Economics Homework 2 Due on April 14th before class

MEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL

Financial Risk Forecasting Chapter 9 Extreme Value Theory

S t d with probability (1 p), where

The Constant Expected Return Model

King s College London

AMH4 - ADVANCED OPTION PRICING. Contents

EE266 Homework 5 Solutions

Chapter 1. Bond Pricing (continued)

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

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

NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours

Measuring Financial Risk using Extreme Value Theory: evidence from Pakistan

BROWNIAN MOTION II. D.Majumdar

Bandit Problems with Lévy Payoff Processes

Financial Econometrics

Estimation of Value at Risk and ruin probability for diffusion processes with jumps

Market Risk: FROM VALUE AT RISK TO STRESS TESTING. Agenda. Agenda (Cont.) Traditional Measures of Market Risk

AN EXTREME VALUE APPROACH TO PRICING CREDIT RISK

Lecture 17. The model is parametrized by the time period, δt, and three fixed constant parameters, v, σ and the riskless rate r.

THE OPTIMAL ASSET ALLOCATION PROBLEMFOR AN INVESTOR THROUGH UTILITY MAXIMIZATION

Strategies for Improving the Efficiency of Monte-Carlo Methods

Pricing Dynamic Solvency Insurance and Investment Fund Protection

Weierstrass Institute for Applied Analysis and Stochastics Maximum likelihood estimation for jump diffusions

1 Geometric Brownian motion

Risk Measurement in Credit Portfolio Models

Analysis of truncated data with application to the operational risk estimation

Beyond the Black-Scholes-Merton model

Universität Regensburg Mathematik

Distortion operator of uncertainty claim pricing using weibull distortion operator

Week 2 Quantitative Analysis of Financial Markets Hypothesis Testing and Confidence Intervals

Hedging Under Jump Diffusions with Transaction Costs. Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo

Asset Allocation Model with Tail Risk Parity

Continuous-Time Pension-Fund Modelling

Comparison of Estimation For Conditional Value at Risk

Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies

Return dynamics of index-linked bond portfolios

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

Slides for Risk Management

King s College London

Risk Neutral Valuation

Liquidity and Risk Management

CONTINGENT CAPITAL WITH DISCRETE CONVERSION FROM DEBT TO EQUITY

A No-Arbitrage Theorem for Uncertain Stock Model

An Efficient Numerical Scheme for Simulation of Mean-reverting Square-root Diffusions

1 The continuous time limit

From Discrete Time to Continuous Time Modeling

Window Width Selection for L 2 Adjusted Quantile Regression

Asymptotic methods in risk management. Advances in Financial Mathematics

On modelling of electricity spot price

Dependence Modeling and Credit Risk

ELEMENTS OF MONTE CARLO SIMULATION

Are stylized facts irrelevant in option-pricing?

Risk management. Introduction to the modeling of assets. Christian Groll

arxiv:cond-mat/ v2 [cond-mat.str-el] 5 Nov 2002

Long-Term Risk Management

Mathematics in Finance

From default probabilities to credit spreads: Credit risk models do explain market prices

Modeling via Stochastic Processes in Finance

STOCHASTIC VOLATILITY AND OPTION PRICING

Business Statistics 41000: Probability 3

Lévy models in finance

1 Rare event simulation and importance sampling

Calibration of Interest Rates

Fat tails and 4th Moments: Practical Problems of Variance Estimation

SADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1. By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD

Market Risk Prediction under Long Memory: When VaR is Higher than Expected

ROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices

INTEREST RATES AND FX MODELS

Version A. Problem 1. Let X be the continuous random variable defined by the following pdf: 1 x/2 when 0 x 2, f(x) = 0 otherwise.

An Application of Extreme Value Theory for Measuring Financial Risk in the Uruguayan Pension Fund 1

Lecture 5: Fundamentals of Statistical Analysis and Distributions Derived from Normal Distributions

Fast Convergence of Regress-later Series Estimators

Lecture 7: Computation of Greeks

American Option Pricing Formula for Uncertain Financial Market

- 1 - **** d(lns) = (µ (1/2)σ 2 )dt + σdw t

Transcription:

Value at Risk and Self Similarity by Olaf Menkens School of Mathematical Sciences Dublin City University (DCU) St. Andrews, March 17 th, 2009

Value at Risk and Self Similarity 1 1 Introduction The concept of Value at Risk (VaR) measures the risk of a portfolio. More precisely, it is a statement of the following form: With probability 1 q the potential loss will not exceed the Value at Risk figure. Speaking in mathematical terms, this is simply the q quantile of the distribution of the change of value for a given portfolio P. More specifically, ( V ar q P d ) = F 1 P d (1 q), where P d is the change of value for a given portfolio over d days (the d day return) and F P d is the distribution function of P d.

Value at Risk and Self Similarity 2 In the daily use it is common to derive the Value at Risk figure of d days from the one of one day by multiplying the Value at Risk figure of one day by d, that is ( V ar q P d) = d V ar q (P 1). (1) Even banking supervisors recommend this procedure (see the Basel Committee on Banking Supervision). On the other hand, if the changes of the value of the considered portfolio P are self similar with Hurst coefficient H (with H 2 1 ), equation (1) has to be modified in the following way: ( V ar q P d) = d H V ar q (P 1). (2)

Value at Risk and Self Similarity 3 To verify this equation, let us first recall the definition of self similarity. Definition 1.1 A real valued process (X(t)) t R is self similar with index H > 0 (H ss) if for all a > 0, the finite dimensional distributions of (X(at)) ( t R are identical to the finite dimensional distributions of a H X(t) ), i.e., if for any a > 0 t R (X(at)) t R d = ( a H X(t) ) t R. This implies F X(at) (x) = F a H X(t) (x) for all a > 0 and t R = P ( a H X(t) < x ) = P ( X(t) < a H x ) = F X(t) ( a H x ).

Value at Risk and Self Similarity 4 2 Possible Size of Risk Underestimation Value at Risk Underestimation for Various Days Given H = 0.55 and H = 0.6 Days H = 0.55 H = 0.6 d d 0.55 d 0.55 d 1 2 Relative Difference in Percent d 0.6 d 0.6 d 1 2 Relative Difference in Percent 1 1 0 0 1 0 0 2 1.46 0.05 3.53 1.52 0.1 7.18 5 2.42 0.19 8.38 2.63 0.39 17.46 10 3.55 0.39 12.2 3.98 0.82 25.89 30 6.49 1.02 18.54 7.7 2.22 40.51 50 8.6 1.53 21.6 10.46 3.39 47.88 100 12.59 2.59 25.89 15.85 5.85 58.49 250 20.84 5.03 31.79 27.46 11.65 73.7 This table shows d H, the difference between d H and d, and the relative difference dh d for d various days d and for H = 0.55 and H = 0.6.

Value at Risk and Self Similarity 5 3 Estimation of the Hurst Exponent via Quantiles There are different methods for estimating the Hurst exponent: the p th moment method (with p N), possibly together with a quantile plot (the so called QQ plot), the R or the R/S statistics (by Hurst or Lo) which captures only the long range dependencies, and the Hill estimator or the peaks over threshold method (POT), which are both in widespread use in the extreme value theory (and which does not account for long range dependencies), just to mention the most popular ones.

Value at Risk and Self Similarity 6 3.1 Theoretical Foundations In order to derive an alternative estimation of the Hurst exponent, let us recall, that V ar q ( P d ) = d H V ar q ( P 1 ), if ( P d) is H ss. Given this, it is easy to derive that ( log (V ar q P d)) ( = H log (d) + log (V ar q P 1)). (3) Thus the Hurst exponent can be derived from the gradient of a linear regression in a log log plot. In particular, it is not necessary that the higher moments of the underlying process exist for this approach. it is possible to observe the evolution of the estimation of the Hurst coefficient along the various quantiles.

Value at Risk and Self Similarity 7 Various Value at Risk Figures for the SAP Daily Closing Stock Prices, based on Commercial Return 0.12 VaR for the Quantiles 0.07 to 0.3 0.07 0.3 0.1 0.08 VaR 0.06 0.04 0.02 0 0 2 4 6 8 10 12 14 16 k Day Return This graphic shows various Value at Risk figures based on the daily closing stock prices for the SAP stock, beginning on June 18 th, 1990 and ending on January 13 th, 2000. The Value at Risk has been calculated for the k day commercial return (with k = 1,...,16) and for various quantiles, ranging from 7 percent to 30 percent.

Value at Risk and Self Similarity 8 Log Log plot for the SAP Stock, Based on Commercial Return log log Plot of the VaR over the k Day Returns for the Quantiles 0.07 to 0.3 3 3.5 4 Logarithm to the Basis 2 of the VaR 4.5 5 5.5 6 6.5 7 7.5 0 0.5 1 1.5 2 2.5 3 3.5 4 Logarithm to the Basis 2 of the k Day Return This graphic shows the logarithm to the basis 2 of various Value at Risk figures based on the daily closing stock prices for the SAP stock, beginning on June 18 th, 1990 and ending on January 13 th, 2000. The Value at Risk has been calculated for the 2 j day commercial return (with j = 0,...,4) and for various quantiles, ranging from 7 percent to 30 percent (the black dashed lines) from top to bottom, respectively. The solid blue lines show the linear regression for the various quantiles.

Value at Risk and Self Similarity 9 There are several approaches for calculating the Value at Risk figure. The most popular are the variance covariance approach, historical simulation, Monte Carlo simulation, and extreme value theory.

Value at Risk and Self Similarity 10 In practice banks often estimate the Value at Risk via order statistics, which is the focus of this talk. Let G j:n (x) be the distribution function of the j th order statistics. Since the probability, that exactly j observations (of a total of n observations) are less or equal x, is given by n! j! (n j)! F(x)j (1 F(x)) n j, it can be verified, that G j:n (x) = n k=j n! k! (n k)! F(x)k (1 F(x)) n k. (4) This is the probability, that at least j observations are less or equal x given a total of n observations.

Value at Risk and Self Similarity 11 Equation (4) implies, that the self similarity holds also for the distribution function of the j th order statistics of a self similar random variable. In this case, one has G j:n,p d(x) = G j:n,p 1 ( d H x ). It is important, that one has n observations for ( P 1) as well as for ( P d ), otherwise the equation does not hold. This shows that the j th order statistics preserves and therefore shows the self similarity of a self similar process. Thus the j th order statistics can be used to estimate the Hurst exponent as it will be done in this talk.

Value at Risk and Self Similarity 12 3.2 Error of the Quantile Estimation Obviously, (3) can only be applied, if V ar q ( P d ) 0. Moreover, close to zero, a possible error in the quantile estimation will lead to an error in (3), which is much larger than the original error from the quantile estimation. Let l be the number, which represents the q th quantile of the order statistics with n observations. With this, x l is the q th quantile of a given time series X, which consists of n observations and with q = n l. Let X be a stochastic process with a differentiable density function f > 0. Moreover the variance of x l is σ 2 x l = q (1 q) n (f(x l )) 2, where f is the density function of X and f must be strictly greater than zero.

Value at Risk and Self Similarity 13 The propagation of errors are calculated by the total differential. Thus, the propagation of this error in (3) is given by σ log(xl ) = 1 x l = ( )1 q (1 q) 2 n (f(x l )) 2 q (1 q) n 1 x l f(x l ).

Value at Risk and Self Similarity 14 1. Example: The normal distribution For example, if X N(0, σ 2 ) the propagation of the error can be written as q (1 q) 2π σ σ log(x) = ( ) n x exp x2 2σ 2 q (1 q) 2π = ( ), n y exp y2 2 where the substitution σ y = x has been used. This shows, that the error is independent of the variance of the underlying process, if this underlying process is normal distributed.

Value at Risk and Self Similarity 15 Error Function for the Normal Distribution 1.8 1.6 1.4 Error of the logarithm of the quantile Error of the quantile estimation 1.2 Error 1 0.8 0.6 0.4 0.2 0 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 Quantile This figure shows the error curves of the quantile estimation (red dash dotted line) and of the logarithm of the quantile estimation (blue solid line) for the normal distribution. This error function is based on the normal density function with standard deviation 1 and with n = 5000.

Value at Risk and Self Similarity 16 Error Function for the Normal Distribution, a Zoom In 0.065 0.06 Error of the logarithm of the quantile Error of the quantile estimation 0.055 0.05 Error 0.045 0.04 0.035 0.03 0.025 0.02 0.015 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 Quantile This figure shows the error curves of the quantile estimation (red dash dotted line) and of the logarithm of the quantile estimation (blue solid line) for the normal distribution. This error function is based on the normal density function with standard deviation 1 and with n = 5000.

Value at Risk and Self Similarity 17 2. Example: The Cauchy distribution Similarly, if X is Cauchy with mean zero, the propagation of the error can be shown to be q (1 q) σ log(x) = = n q (1 q) n π (x2 + σ 2 ) σ x π (y2 + 1) y where the substitution σ y = x has also been used.,

Value at Risk and Self Similarity 18 Error Function for the Cauchy Distribution, a Zoom In 0.08 0.075 Error of the logarithm of the quantile Error of the quantile estimation 0.07 0.065 Error 0.06 0.055 0.05 0.045 0.04 0.035 0.03 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 Quantile This figure shows the error curves of the quantile estimation (red dash dotted line) and of the logarithm of the quantile estimation (blue solid line) for the Cauchy distribution. This error function is based on the Cauchy density function with standard deviation 1 and with n = 5000.

Value at Risk and Self Similarity 19 4 Estimating the Scaling Law for Some Stocks A self similar process with Hurst exponent H can not have a drift. Since it is recognized that financial time series do have a drift, they can not be self similar. Because of this the wording scaling law instead of Hurst exponent will be used when talking about financial time series, which have not been detrended. Since the scaling law is more relevant in practice than in theory,only those figures are depicted which are based on commercial returns.

Value at Risk and Self Similarity 20 Estimation of the Scaling Law for 24 DAX Stocks, Left Side 0.65 Hurst Exponent Estimation for the Quantiles 0.01 to 0.3 for 24 Time Series Hurst Exponent 0.6 0.55 0.5 0.45 0.4 0.35 dcx 0.3 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 vow hvm kar meo bas veb hoe bay The underlying regression uses 5 points, starting with the 1 day and ending with the 16 day returns. This plot is based on commercial returns and order statistics. Quantile Shown are the estimation for the scaling law for 24 DAX stocks. The underlying time series is a commercial return. Shown are the lower (left) quantiles. The following stocks are denoted explicitly: DaimlerChysler (dcx), Karstadt (kar), Volkswagen (vow), Metro (meo), Hypovereinsbank (hvm), BASF (bas), Veba (veb), Hoechst (hoe), and Bayer (bay).

Value at Risk and Self Similarity 21 Estimation of the Scaling Law for 24 DAX Stocks, Left Side 0.65 Quantitative Characteristics of the Hurst Exponent Estimation for the Quantiles 0.01 to 0.3 for 24 Time Series Hurst Exponent 0.6 0.55 0.5 0.45 0.4 0.35 Mean Mean+/ Stddeviation Maximum Minimum Actual Hurst Exponent 0.3 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 Quantile The green solid line is the mean, the magenta dash dotted lines are the mean plus/minus the standard deviation, the red triangles are the minimum, and the blue upside down triangles are the maximum of the estimation for the scaling law, which are based on 24 DAX stocks. The underlying time series is a commercial return. Shown are the lower (left) quantiles. The underlying regression uses 5 points, starting with the 1 day and ending with the 16 day returns. This plot is based on commercial returns and order statistics.

Value at Risk and Self Similarity 22 Error of the Estimation of the Scaling Law for 24 DAX Stocks Error of the Hurst Exponent Estimation Quantitative Characteristics of the Error of the Regression for the Hurst Exponent Estimation for the Quantiles 0.01 to 0.3 for 24 Time Series 0.08 0.07 0.06 0.05 0.04 0.03 0.02 0.01 Mean Mean+/ Stddeviation Maximum Minimum 0 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 Quantile The green solid line is the mean, the magenta dash dotted lines are the mean plus/minus the standard deviation, the red triangles are the minimum, and the blue upside down triangles are the maximum of the error curves of the estimation for the scaling law, which are based on 24 DAX stocks. The underlying time series is a commercial return. The underlying regression uses 5 points, starting with the 1 day and ending with the 16 day returns. This plot is based on commercial returns and order statistics.

Value at Risk and Self Similarity 23 Estimation of the Scaling Law for 24 DAX Stocks, Right Side Hurst Exponent 0.75 0.7 0.65 0.6 0.55 0.5 0.45 Quantitative Characteristics of the Hurst Exponent Estimation for the Quantiles 0.7 to 0.99 for 24 DAX Stocks Mean Mean+/ Stddeviation Maximum Minimum Actual Hurst Exponent 0.4 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 Quantile The green solid line is the mean, the magenta dash dotted lines are the mean plus/minus the standard deviation, the red triangles are the minimum, and the blue upside down triangles are the maximum of the estimation for the scaling law, which are based on 24 DAX stocks. The underlying time series is a commercial return. Shown are the upper (right) quantiles. The underlying regression uses 5 points, starting with the 2 0 day and ending with the 2 4 day returns. This plot is based on commercial returns and order statistics. These figures are based on the absolute VaR.

Value at Risk and Self Similarity 24 Error of the Estimation of the Scaling Law for 24 DAX Stocks Error of the Hurst Exponent Estimation Quantitative Characteristics of the Error of the Regression for the Hurst Exponent Estimation for the Quantiles 0.99 to 0.7 for 24 DAX Stocks 0.03 0.025 0.02 0.015 0.01 0.005 Mean Mean+/ Stddeviation Maximum Minimum 0 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 Quantile The green solid line is the mean, the magenta dash dotted lines are the mean plus/minus the standard deviation, the red triangles are the minimum, and the blue upside down triangles are the maximum of the error curves of the estimation for the scaling law, which are based on 24 DAX stocks. The underlying time series is a commercial return. The underlying regression uses 5 points, starting with the 2 0 day and ending with the 2 4 day returns. This plot is based on commercial returns and order statistics. These figures are based on the absolute VaR.

Value at Risk and Self Similarity 25 Estimation of the Scaling Law for the DJI and its Stocks, LQ 0.6 Quantitative Characteristics of the Hurst Exponent Estimation for the Quantiles 0.01 to 0.3 for 31 Time Series Hurst Exponent 0.55 0.5 0.45 0.4 0.35 0.3 0.25 0.2 0.15 Mean Mean+/ Stddeviation Maximum Minimum Actual Hurst Exponent 0.1 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 Quantile The green solid line is the mean, the magenta dash dotted lines are the mean plus/minus the standard deviation, the red triangles are the minimum, and the blue upside down triangles are the maximum of the estimation for the scaling law, which are based on the DJI and its 30 stocks. The underlying time series is a commercial return. Shown are the lower (left) quantiles. The underlying regression uses 5 points, starting with the 1 day and ending with the 16 day returns. This plot is based on commercial returns and order statistics.

Value at Risk and Self Similarity 26 Estimation of the Scaling Law for the DJI and its Stocks, RQ Hurst Exponent 0.7 0.65 0.6 0.55 0.5 0.45 0.4 Quantitative Characteristics of the Hurst Exponent Estimation for the Quantiles 0.7 to 0.99 for 31 DJI Stocks Mean Mean+/ Stddeviation Maximum Minimum Actual Hurst Exponent 0.35 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 Quantile The green solid line is the mean, the magenta dash dotted lines are the mean plus/minus the standard deviation, the red triangles are the minimum, and the blue upside down triangles are the maximum of the estimation for the scaling law, which are based on the DJI and its 30 stocks. The underlying time series is a commercial return. Shown are the upper (right) quantiles. The underlying regression uses 5 points, starting with the 2 0 day and ending with the 2 4 day returns. This plot is based on commercial returns and order statistics. These figures are based on the absolute VaR.

Value at Risk and Self Similarity 27 5 Determining the Hurst Exponent It has been already stated, that the financial time series can not be self similar. However, it is possible that the detrended financial time series are self similar with Hurst exponent H. This will be scrutinized in the following where the financial time series have been detrended. Since the Hurst exponent is more relevant in theory than in practice, only those figures are shown which are based on logarithmic returns.

Value at Risk and Self Similarity 28 Hurst Exponent Estimation for 24 DAX Stocks, Left Quantile 0.62 Quantitative Characteristics of the Hurst Exponent Estimation for the Quantiles 0.01 to 0.3 for 24 Time Series Hurst Exponent 0.6 0.58 0.56 0.54 0.52 0.5 0.48 Mean Mean+/ Stddeviation Maximum Minimum Actual Hurst Exponent 0.46 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 Quantile The green solid line is the mean, the magenta dash dotted lines are the mean plus/minus the standard deviation, the red triangles are the minimum, and the blue upside down triangles are the maximum of the estimation for the Hurst exponent, which are based on 24 DAX stocks. The underlying time series are logarithmic returns, which have been detrended. Shown are the lower (left) quantiles. The underlying regression uses 5 points, starting with the 1 day and ending with the 16 day returns. This plot is based on logarithmic returns and order statistics. The underlying log processes have been detrended.

Value at Risk and Self Similarity 29 Hurst Exponent Estimation for 24 DAX Stocks, Right Quantile Hurst Exponent 0.7 0.65 0.6 0.55 0.5 0.45 0.4 Quantitative Characteristics of the Hurst Exponent Estimation for the Quantiles 0.7 to 0.99 for 24 DAX Stocks Mean Mean+/ Stddeviation Maximum Minimum Actual Hurst Exponent 0.35 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 Quantile The underlying regression uses 5 points, starting with the 2 0 day and ending with the 2 4 day returns. This plot is based on logarithmic returns and order statistics. These figures are based on the absolute VaR. The underlying processes have been detrended The green solid line is the mean, the magenta dash dotted lines are the mean plus/minus the standard deviation, the red triangles are the minimum, and the blue upside down triangles are the maximum of the estimation for the Hurst exponent, which are based on 24 DAX stocks. The underlying time series are logarithmic returns, which have been detrended. Shown are the upper (right) quantiles.

Value at Risk and Self Similarity 30 Hurst Exponent Estimation for the DJI and its Stocks, LQ 0.65 Quantitative Characteristics of the Hurst Exponent Estimation for the Quantiles 0.01 to 0.3 for 31 Time Series Hurst Exponent 0.6 0.55 0.5 0.45 0.4 0.35 Mean Mean+/ Stddeviation Maximum Minimum Actual Hurst Exponent 0.3 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 Quantile The green solid line is the mean, the magenta dash dotted lines are the mean plus/minus the standard deviation, the red triangles are the minimum, and the blue upside down triangles are the maximum of the estimation for the Hurst exponent, which are based on the DJI and its 30 stocks. The underlying time series are logarithmic returns, which have been detrended. Shown are the lower (left) quantiles. The underlying regression uses 5 points, starting with the 1 day and ending with the 16 day returns. This plot is based on logarithmic returns and order statistics. The underlying log processes have been detrended.

Value at Risk and Self Similarity 31 Hurst Exponent Estimation for the DJI and its Stocks, RQ Hurst Exponent 0.7 0.65 0.6 0.55 0.5 0.45 0.4 0.35 Quantitative Characteristics of the Hurst Exponent Estimation for the Quantiles 0.7 to 0.99 for 31 DJI Stocks Mean Mean+/ Stddeviation Maximum Minimum Actual Hurst Exponent 0.3 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 Quantile The underlying regression uses 5 points, starting with the 2 0 day and ending with the 2 4 day returns. This plot is based on logarithmic returns and order statistics. These figures are based on the absolute VaR. The underlying processes have been detrended The green solid line is the mean, the magenta dash dotted lines are the mean plus/minus the standard deviation, the red triangles are the minimum, and the blue upside down triangles are the maximum of the estimation for the Hurst exponent, which are based on the DJI and its 30 stocks. The underlying time series are logarithmic returns, which have been detrended. Shown are the upper (right) quantiles.

Value at Risk and Self Similarity 32 6 Interpretation of the Hurst Exponent for Financial Time Series There are different possibilities of which we discuss the two most relevant ones: Fractional Brownian Motion: H describes the persistence of the process. 0 < H < 2 1 means that the process is anti persistent. < H < 1 means that the process is persistent. 1 2 α stable Lévy Process: H = α 1 for α [1,2] determines how much the process is heavy tailed.

Value at Risk and Self Similarity 33 The bigger H is, the more heavy tailed the process is.

Value at Risk and Self Similarity 34 Fractional Brownian Motion Definition 6.1 Let 0 < H 1. A real valued Gaussian process (B H (t)) t 0 is called fractional Brownian motion if E[B H (t)] = 0 and E[B H (t) B H (s)] = 1 2 { t 2H + s 2H t s 2H} E [ B H (1) 2]. Note that the distribution of a Gaussian process is determined by its mean and covariance structure. Hence, the two conditions given in the above definition specify a unique Gaussian process. is a Brownian motion up to a multiplicative con- ( B 1/2 (t) ) t 0 stant.

Value at Risk and Self Similarity 35 A Sample Path of the Fractional Brownian Motion with H = 0.1 7 6 5 4 3 2 1 0 1 2 3 0 100 200 300 400 500 600 700 800 900 1000 Time Sample Path of the fractional Brownian motion One sample path of the fractional Brownian motion with parameters σ = 1, H = 0.1, N = 1000.

Value at Risk and Self Similarity 36 A Sample Path of the Fractional Brownian Motion with H = 0.4 20 15 10 5 0 5 10 15 0 100 200 300 400 500 600 700 800 900 1000 Time Sample Path of the fractional Brownian motion One sample path of the fractional Brownian motion with parameters σ = 1, H = 0.4, N = 1000.

Value at Risk and Self Similarity 37 A Sample Path of the Fractional Brownian Motion with H = 0.6 40 30 20 10 0 10 20 0 100 200 300 400 500 600 700 800 900 1000 Time Sample Path of the fractional Brownian motion One sample path of the fractional Brownian motion with parameters σ = 1, H = 0.6, N = 1000.

Value at Risk and Self Similarity 38 A Sample Path of the Fractional Brownian Motion with H = 0.9 0 20 40 60 80 100 120 140 160 180 0 100 200 300 400 500 600 700 800 900 1000 Time Sample Path of the fractional Brownian motion One sample path of the fractional Brownian motion with parameters σ = 1, H = 0.9, N = 1000.

Value at Risk and Self Similarity 39 Stable Lévy Processes Definition 6.2 A (cadlag) stochastic process (X t ) t 0 on (Ω, F, P) with values in R d such that X 0 = 0 is called a Lévy process if it possesses the following properties: 1. Independent increments: for every increasing sequence of times t 0... t n, the random variables X t0, X t1 X t0,..., X tn X tn 1 are independent. 2. Stationary increments: the distribution of X t+h X t does not depend on t. 3. Stochastic continuity: ǫ > 0, lim h 0 P ( X t+h X t ǫ ) = 0. Definition 6.3 A probability measure µ on R d is called (strictly) stable, if for any a > 0, there exits b > 0 such that ˆµ(θ) a = ˆµ(b θ) for all θ R d.

Value at Risk and Self Similarity 40 Theorem 6.4 If µ on R d is stable, there exits a unique α (0,2] such that b = a 1/α. Such a µ is referred to as α stable. When α = 2, µ is a mean zero Gaussian probability measure. Theorem 6.5 Suppose (X t ) t 0 is a Lévy process. Then L(X(1)) is stable if and only if (X t ) is self similar. The index α of stability and the exponent H of self similarity satisfy α = 1 H. If Z α is an R d valued random variable with a α stable distribution, 0 < α < 2, then for any γ (0, α), E[ Z α γ ] <, but E[ Z α α ] =. Note that being heavy tailed implies that the variance is infinite.

Value at Risk and Self Similarity 41 0.35 0.3 0.25 0.2 0.15 0.1 0.05 Cauchy Density versus Normal Density Cauchy density with location δ = 0 and scale γ = 1. Moreover, normal density with µ = 0 and σ = 1.4142. 0 5 4 3 2 1 0 1 2 3 4 5 x Cauchy density normal density

Value at Risk and Self Similarity 42 7 Conclusion The main results are that the scaling coefficient 0.5 has to be used very carefully for financial time series and there are substantial doubts about the self similarity of the underlying processes of financial time series.

Value at Risk and Self Similarity 43 Concerning the scaling law, it is better to use a scaling law of 0.55 for the left quantile and a scaling law of 0.6 for the right quantile (the short positions) in order to be on the safe side. It is important to keep in mind that these figures are only based on the (highly traded) DAX and Dow Jones Index stocks. Considering low traded stocks might yield even higher maximal scaling laws. These numbers should be set by the market supervision institutions like the SEC. However, it is possible for the banks to reduce their Value at Risk figures if they use the correct scaling law numbers. For instance, the Value at Risk figure of a well diversified portfolio of Dow Jones Index stocks would be reduced in this way about 12%, since it would have a scaling law of approximately 0.44.

Value at Risk and Self Similarity 44 8 Extensions Using bootstrap methods in order to overcome the lack of data. Developing a test on self similarity on the quantiles which overcomes the phenomena already described by Granger and Newbold.

Value at Risk and Self Similarity 45 References 1. Paul Embrechts and Makoto Maejima, Selfsimilar Processes, Princeton University Press, Princeton, 2002. 2. Olaf Menkens, Value at Risk and Self Similarity, In Numerical Methods for Finance, John Miller, John Appleby and David Edelman (Eds.), CRC Press/Chapman & Hall, 2007. 3. Gennardy Samorodnitsky and Murad S. Taqqu, Stable Non Gaussian Processes, Chapman & Hall, London, 1994.