Implied Systemic Risk Index (work in progress, still at an early stage)

Similar documents
Risk Aggregation with Dependence Uncertainty

Risk Aggregation with Dependence Uncertainty

Mean-Variance Optimal Portfolios in the Presence of a Benchmark with Applications to Fraud Detection

Financial Risk Management

Value-at-Risk bounds with variance constraints

All Investors are Risk-averse Expected Utility Maximizers. Carole Bernard (UW), Jit Seng Chen (GGY) and Steven Vanduffel (Vrije Universiteit Brussel)

Market Risk Analysis Volume II. Practical Financial Econometrics

All Investors are Risk-averse Expected Utility Maximizers

3.4 Copula approach for modeling default dependency. Two aspects of modeling the default times of several obligors

Correlation and Diversification in Integrated Risk Models

Modeling Co-movements and Tail Dependency in the International Stock Market via Copulae

Path-dependent inefficient strategies and how to make them efficient.

Market Risk Analysis Volume I

Lecture 1 of 4-part series. Spring School on Risk Management, Insurance and Finance European University at St. Petersburg, Russia.

Financial Models with Levy Processes and Volatility Clustering

CPSC 540: Machine Learning

Asset Allocation Model with Tail Risk Parity

Risk minimization and portfolio diversification

Equity correlations implied by index options: estimation and model uncertainty analysis

CPSC 540: Machine Learning

SOLVENCY AND CAPITAL ALLOCATION

BROWNIAN MOTION Antonella Basso, Martina Nardon

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

Asymptotic methods in risk management. Advances in Financial Mathematics

Operational Risk Aggregation

UPDATED IAA EDUCATION SYLLABUS

INTEREST RATES AND FX MODELS

Graduate School of Business, University of Chicago Business 41202, Spring Quarter 2007, Mr. Ruey S. Tsay. Solutions to Final Exam

INTERNATIONAL JOURNAL FOR INNOVATIVE RESEARCH IN MULTIDISCIPLINARY FIELD ISSN Volume - 3, Issue - 2, Feb

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 2010, Mr. Ruey S. Tsay Solutions to Final Exam

ESGs: Spoilt for choice or no alternatives?

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

Handbook of Financial Risk Management

FIN FINANCIAL INSTRUMENTS SPRING 2008

Approximating a multifactor di usion on a tree.

Master s in Financial Engineering Foundations of Buy-Side Finance: Quantitative Risk and Portfolio Management. > Teaching > Courses

Asymmetric Price Transmission: A Copula Approach

MEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL

ERM (Part 1) Measurement and Modeling of Depedencies in Economic Capital. PAK Study Manual

PORTFOLIO THEORY. Master in Finance INVESTMENTS. Szabolcs Sebestyén

Hedging Derivative Securities with VIX Derivatives: A Discrete-Time -Arbitrage Approach

Risk Measurement in Credit Portfolio Models

Market risk measurement in practice

Operational Risk Aggregation

Modeling of Price. Ximing Wu Texas A&M University

Multiple Objective Asset Allocation for Retirees Using Simulation

Valuation of performance-dependent options in a Black- Scholes framework

Lecture notes on risk management, public policy, and the financial system. Credit portfolios. Allan M. Malz. Columbia University

IEOR E4602: Quantitative Risk Management

A New Tool For Correlation Risk Management: The Market Implied Comonotonicity Gap

IEOR E4602: Quantitative Risk Management

Risk Quadrangle and Applications in Day-Trading of Equity Indices

Statistical Models and Methods for Financial Markets

IRC / stressed VaR : feedback from on-site examination

A Non-Parametric Technique of Option Pricing

ADVANCED OPERATIONAL RISK MODELLING IN BANKS AND INSURANCE COMPANIES

Internet Appendix for Asymmetry in Stock Comovements: An Entropy Approach

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

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

Mathematics in Finance

Downside Risk: Implications for Financial Management Robert Engle NYU Stern School of Business Carlos III, May 24,2004

Capital allocation: a guided tour

2017 IAA EDUCATION SYLLABUS

2. Copula Methods Background

Variable Annuities with fees tied to VIX

Operational risk Dependencies and the Determination of Risk Capital

P2.T6. Credit Risk Measurement & Management. Malz, Financial Risk Management: Models, History & Institutions

Cost-efficiency and Applications

Analytical Pricing of CDOs in a Multi-factor Setting. Setting by a Moment Matching Approach

Market Risk Analysis Volume IV. Value-at-Risk Models

2.1 Mathematical Basis: Risk-Neutral Pricing

Dynamic Copula Methods in Finance

ERM Sample Study Manual

HANDBOOK OF. Market Risk CHRISTIAN SZYLAR WILEY

TEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING

Lecture 4: Forecasting with option implied information

Quantitative Risk Management

CONCORDANCE MEASURES AND SECOND ORDER STOCHASTIC DOMINANCE PORTFOLIO EFFICIENCY ANALYSIS

Optimal construction of a fund of funds

Lecture 4 of 4-part series. Spring School on Risk Management, Insurance and Finance European University at St. Petersburg, Russia.

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

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

ESTIMATION OF UTILITY FUNCTIONS: MARKET VS. REPRESENTATIVE AGENT THEORY

LECTURE NOTES 10 ARIEL M. VIALE

Volatility Models and Their Applications

The Optimization Process: An example of portfolio optimization

Implementing Models in Quantitative Finance: Methods and Cases

Desirable properties for a good model of portfolio credit risk modelling

Some Simple Stochastic Models for Analyzing Investment Guarantees p. 1/36

Computer Exercise 2 Simulation

Monte Carlo Methods in Finance

Can we use kernel smoothing to estimate Value at Risk and Tail Value at Risk?

Dynamic Asset and Liability Management Models for Pension Systems

Key Words: emerging markets, copulas, tail dependence, Value-at-Risk JEL Classification: C51, C52, C14, G17

Log-Robust Portfolio Management

Portfolio Optimization using Conditional Sharpe Ratio

PORTFOLIO MODELLING USING THE THEORY OF COPULA IN LATVIAN AND AMERICAN EQUITY MARKET

Some characteristics of an equity security next-year impairment

Utility Indifference Pricing and Dynamic Programming Algorithm

Transcription:

Implied Systemic Risk Index (work in progress, still at an early stage) Carole Bernard, joint work with O. Bondarenko and S. Vanduffel IPAM, March 23-27, 2015: Workshop I: Systemic risk and financial networks Carole Bernard Implied Systemic Risk 1/35

Systemic Risk, Contagion, Dependence ˆ Contagion linked to coincidence of extreme returns. Studies on coexceedances using bank data reject the assumption that the modelling can be done using a multivariate Gaussian assumption. They find non linearities in the tail. Large common shocks are highly correlated compared to small shocks. ˆ Systemic risk measures: CoVaR, CoES, SRISK... are conditional risk measures. Existing studies on contagion and systemic risk measures are under the real-world measure. Our objective is to measure contagion/systemic risk using option prices (thus information on the risk-neutral probability). Carole Bernard Implied Systemic Risk 2/35

Outline 1 The CBOE implied correlation: use information on implied volatilities of at-the-money options 2 Toy examples with factor models to show that it may fail to capture changes in the dependence among assets. Importance of using the full marginal distributions 3 An algorithm to describe the set of possible dependence structures (copulas) that are consistent with the information ˆ marginal distribution of each individual asset return X i ˆ distribution of aggregated risk ω 1 X 1 +... + ω d X d 4 Algorithm useful to detect changes in implied dependence. forward looking measure of contagion/of systemic risk. to compute conditional risk measures (correlation in the tail), systemic risk measures under the risk neutral probability (model-free) Carole Bernard Implied Systemic Risk 3/35

CBOE implied correlation S T = i ω ix i,t with i ω i = 1 For observed at-the-money call option prices with maturity T, define σ S and σ i as follows C index,observed = BlackScholesCall(σ S, S 0 ) E [ [ ( (S T S 0 ) +] ) ] + = E S 0 e (r σ2 S 2 )T +σ S W T S 0 C Xi,observed = BlackScholesCall(σ i, X i,0 ) E [ [ ( (X i,t X i,0 ) +] ) ] + = E X i,0e (r σ2 i 2 )T +σ i W i,t X i,0 Use these implied volatilities to define an implied correlation index. Carole Bernard Implied Systemic Risk 4/35

CBOE implied correlation The CBOE correlation index is defined by ρ cboe = σ 2 S d 2 d 1 i=1 i=1 ω2 i σ2 i j>i ω iω j σ i σ j... assuming that the index implied volatility σ S and the individual implied volatilities σ i for i = 1,..., d are such that σ 2 S = d d 1 ωi 2 σi 2 + 2 ω i ω j σ i σ j ρ ij i=1 i=1 where ρ ij is constant equal to ρ cboe Given this assumption, d 1 i=1 j>i ρ cboe = ω iω j σ i σ j ρ ij d 1 i=1 j>i ω iω j σ i σ j j>i Carole Bernard Implied Systemic Risk 5/35

Comments on the CBOE implied correlation not always between -1 and 1. not always a feasible correlation parameter. The CBOE implied correlation index can be very far from the weighted average of pairwise correlations. No setting in which the assumption of a constant pairwise correlation ρ among assets logreturns allows us to find that the CBOE implied correlation is equal to ρ exactly. But it works well in a multivariate Black Scholes with constant pairwise correlation ρ ij. It is affected by changes in marginal distributions and not only in the dependence. It does not give any information on the dependence in the tail (global measure). Carole Bernard Implied Systemic Risk 6/35

Proposal ˆ The CBOE implied correlation index makes use of the implied volatilities of at-the-money option prices only ˆ Use all strikes to get the full marginal distribution of X i and S and infer the dependence structures that are compatible with this information. Method based on the Rearrangement Algorithm. Using this approach, we can compute for example the average pairwise Pearson correlation ρ := d 1 i=1 d 1 i=1 j>i ω iω j σ i σ j ρ ij j>i ω iω j σ i σ j Then we compare ρ with the CBOE index and with the true correlation in toy examples... Carole Bernard Implied Systemic Risk 7/35

In a multidimensional Black Scholes model with heterogeneous volatilities (from 20% to 70%) with logreturns that are multivariate Gaussian with homogeneous correlation matrix ρ ij = ρ for all i j. d = 10 assets, weights are all equal to 1/d 0.8 0.7 CBOE index with T=6mths Average Pearson correlation by rearrangement with T=6mths 0.6 Dependence index 0.5 0.4 0.3 0.2 0.1 0-0.1-0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 ; The CBOE index is roughly equal to ρ and ρ. Carole Bernard Implied Systemic Risk 8/35

Factor model with 2 assets returns Define X i = 100e r v 2 i 2 +v i W i (1), Z = 100e r σ2 Z 2 +σ Z W Z (1) W 1, W 2 and W Z are Brownian motions. W Z is independent of W 1 and W 2 and W 1 and W 2 have correlation ρ 12. Carole Bernard Implied Systemic Risk 9/35

Factor model with 2 assets returns Define X i = 100e r v 2 i 2 +v i W i (1), Z = 100e r σ2 Z 2 +σ Z W Z (1) W 1, W 2 and W Z are Brownian motions. W Z is independent of W 1 and W 2 and W 1 and W 2 have correlation ρ 12. Let I be a variable indicating in which regime we are. S 1 = (1 I)X 1 + IZ S 2 = (1 I)X 2 + IZ In one regime (when I = 1), S 1 and S 2 are perfectly dependent (equal here) and in the other regime, 2-dimensional Black-Scholes. Carole Bernard Implied Systemic Risk 9/35

Define Factor model with 2 assets returns X i = 100e r v 2 i 2 +v i W i (1), Z = 100e r σ2 Z 2 +σ Z W Z (1) W 1, W 2 and W Z are Brownian motions. W Z is independent of W 1 and W 2 and W 1 and W 2 have correlation ρ 12. Let I be a variable indicating in which regime we are. S 1 = (1 I)X 1 + IZ S 2 = (1 I)X 2 + IZ In one regime (when I = 1), S 1 and S 2 are perfectly dependent (equal here) and in the other regime, 2-dimensional Black-Scholes. Index: = S 1 2 + S 2 2. I = 1 Z<zq where z q is the Value-at-Risk at level q of Z. Carole Bernard Implied Systemic Risk 9/35

Within this toy model Factor model with 2 assets returns ˆ By simulation, get prices for at-the-money calls on S and X i. ˆ Estimate implied volatilities σ i and σ S. ˆ Compute CBOE index from implied volatilities. ˆ For our approach (that I will describe later) we need to specify the marginal distributions (and not just the implied volatilities) 1 Use lognormal distribution for X 1, X 2 and S with logmean r σ2 i 2 and r σ2 S 2 respectively and logvariance σi 2 and σs 2. 2 Use empirical distributions obtained by simulation (correct margins). Carole Bernard Implied Systemic Risk 10/35

Change of regime driven by Z I = 1 Z<zq where z q is the Value-at-Risk at level q of Z. 1 0.8 S i = (1-I)X i +I Z with I=1 Z<quantileq CBOE index Pearson correlation between logreturns of S 1 and S 2 in the model Average Pearson correlation by rearrangement (LN margins) 0.6 Dependence index 0.4 0.2 0-0.2-0.4 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 probability q of I=1 in the mixture Carole Bernard Implied Systemic Risk 11/35

How to explain this graph? Change of regime driven by Z Carole Bernard Implied Systemic Risk 12/35

Change of regime driven by Z How to explain this graph? S 1, S 2, S 1 + S 2 are far from lognormally distributed... q = 0.25 q = 0.5 q = 0.75 0.03 0.035 0.025 Probability Distribution 0.025 0.02 0.015 0.01 0.005 density of S 1 density of S 2 density of index Probability Distribution 0.03 0.025 0.02 0.015 0.01 0.005 density of S 1 density of S 2 density of index Probability Distribution 0.02 0.015 0.01 0.005 density of S 1 density of S 2 density of index 0 0 50 100 150 200 250 S 0 0 50 100 150 200 250 S 0 0 50 100 150 200 250 S Carole Bernard Implied Systemic Risk 12/35

Change of regime driven by Z - Correct margins I = 1 Z<zq where z q is the Value-at-Risk at level q of Z. We apply our method with the correct margins (and not with Lognormal). Better than CBOE! 1 0.8 S i = (1-I)X i +I Z with I=1 Z<quantileq CBOE index Pearson correlation between logreturns of S and S in the model 1 2 Average Pearson correlation by rearrangement (correct margins) 0.6 Dependence index 0.4 0.2 0-0.2-0.4 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 probability q of I=1 in the mixture Other example with I independent ( indep ) Carole Bernard Implied Systemic Risk 13/35

Consequences Marginal distributions matter a lot The CBOE implied correlation is model-free, but it is roughly equal to the average pairwise correlation assuming ˆ logreturns are normal ˆ Gaussian dependence Our approach allows us to compute the average pairwise correlation with the information about margins of the index components and of the index. In fact, our approach finds the set of dependence structures consistent with margin informations (i.e. full joint distribution of (X 1, X 2,..., X d )). We can thus compute anything... Let us explain how? Carole Bernard Implied Systemic Risk 14/35

Algorithm to infer dependence Inputs ˆ Distributions of X i for i = 1, 2,..., d (discretized) ˆ Distribution of the index S (discretized) Output The joint distribution of (X 1, X 2,..., X d ) Carole Bernard Implied Systemic Risk 15/35

Algorithm to infer dependence N = 4 observations of d = 3 variables: X 1, X 2, X 3 M = 1 1 2 0 6 3 4 0 0 6 3 4 Each column: marginal distribution Interaction among columns: dependence Rearrange the order of the elements per column Same margins but effect on the sum! Find the right rearrangement. Carole Bernard Implied Systemic Risk 16/35

Use of the Rearrangement Algorithm first used to minimize var(x 1 + X 2 +... + X d ) Why do we need an algorithm? Carole Bernard Implied Systemic Risk 17/35

Use of the Rearrangement Algorithm first used to minimize var(x 1 + X 2 +... + X d ) Why do we need an algorithm? When d = 2, then the minimum variance is the lower Fréchet-Hoeffding bound or extreme negative dependence (antimonotonic) var(f 1 1 (U) + F 1 2 (1 U)) var(x 1 + X 2 ) Carole Bernard Implied Systemic Risk 17/35

Use of the Rearrangement Algorithm first used to minimize var(x 1 + X 2 +... + X d ) Why do we need an algorithm? When d = 2, then the minimum variance is the lower Fréchet-Hoeffding bound or extreme negative dependence (antimonotonic) var(f 1 1 (U) + F 1 2 (1 U)) var(x 1 + X 2 ) When d 2, the Fréchet lower bound does not exist: ˆ Wang and Wang (2011) study complete mixability (X 1 F 1,..., X d F d are completely mixable if there exists a dependence structure between X 1,... X d such that X 1 + X 2 +... + X d = cst) ˆ Puccetti and Rüschendorf (2012): algorithm (RA) useful to approximate the minimum variance. Carole Bernard Implied Systemic Risk 17/35

Inputs: ˆ X 1 F 1,..., X d F d Solving for the minimum variance ˆ Goal: look for copulas such that min var(x 1 + X 2 +... + X d ) It s a NP complete problem: there are no efficient algorithms but we develop an heuristic that performs very well in practice. Carole Bernard Implied Systemic Risk 18/35

Rearrangement Algorithm to solve the minimum variance N = 4 observations of d = 3 variables: X 1, X 2, X 3 M = 1 1 2 0 6 3 4 0 0 6 3 4 Each column: marginal distribution Interaction among columns: dependence Carole Bernard Implied Systemic Risk 19/35

Rearrangement Algorithm: Sum with Minimum Variance minimum variance with d = 2 risks X 1 and X 2 Antimonotonicity: var(x a 1 + X 2) var(x 1 + X 2 ) How about in d dimensions? Carole Bernard Implied Systemic Risk 20/35

Rearrangement Algorithm: Sum with Minimum Variance minimum variance with d = 2 risks X 1 and X 2 Antimonotonicity: var(x a 1 + X 2) var(x 1 + X 2 ) How about in d dimensions? Use of the rearrangement algorithm on the original matrix M. Aggregate Risk with Minimum Variance Columns of M are rearranged such that they become anti-monotonic with the sum of all other columns. k {1, 2,..., d}, X a k antimonotonic with j k X j After each step, var ( X a k + ) ( j k X j var X k + ) j k X j where X a k is antimonotonic with j k X j Carole Bernard Implied Systemic Risk 20/35

Bernard, Department of Statistics and Actuarial Science at the University of Waterl uwaterloo.ca). Carole Bernard Implied Systemic Risk 21/35 Outline Implied Correlation Toy examples Algorithm Implied Index Toy example Conclusions 6 6 4 4 3 3 1 1 2 S N = Aggregate risk with minimum variance 0 0Step 0 1: First column0 16 9 3 X 2 + X 3 6 6 4 10 4 3 2 5 1 1 1 2 0 0 0 0 becomes 0 6 4 1 3 2 4 1 1 6 0 0 set...

Aggregate risk with minimum variance X 2 + X 3 6 6 4 10 4 3 2 5 1 1 1 2 0 0 0 0 X 1 + X 3 4 3 5 6 0 6 4 1 3 2 4 1 1 6 0 0 0 3 4 1 6 2 4 1 1 6 0 0 X 1 + X 2 3 7 5 6 becomes becomes becomes 0 6 4 1 3 2 4 1 1 6 0 0 0 3 4 1 6 2 4 1 1 6 0 0 0 3 4 1 6 0 4 1 2 6 0 1 lumns are antimonotonic with the sum of the others: X 2 + X 3 4 7 X 1 + X 3 0 3 4 4 0 3 4 Carole Bernard Implied Systemic Risk 22/35

1 6 2 7 becomes 1 6 0 (5) 4 1 1 5 4 1 2 6 0 0 6 6 0 1 Aggregate risk with minimum variance All columns are antimonotonic with the sum of the others: Each column is antimonotonic with the sum of the others: X 2 + X 3 X 1 + X 3 X 1 + X 2,, Outline Implied Correlation Toy examples Algorithm Implied Index Toy example Conclusions 0 3 4 1 6 0 4 1 2 6 0 1 7 6 3 1 0 3 4 1 6 0 4 1 2 6 0 1 4 1 6 7 0 3 4 1 6 0 4 1 2 6 0 1 3 7 5 6 Minimum variance sum 0 3 4 1 6 0 4 1 2 6 0 1 X 1 + X 2 + X 3 7 7 7 7 S N = (6) Carole Bernard Implied Systemic Risk 23/35

1 6 2 7 becomes 1 6 0 (5) e antimonotonic 4with 1 1 the 5 sum of 4the 1 2others: 6 0 0 6 6 0 1 Aggregate risk with minimum variance 2 + X 3 All columns are antimonotonic X 1 + X 3 7 with the sum of the others: Each column is antimonotonic 0 3 with 4 the sum of 4 the others: 6, X 2 + X 3 1 6 0 X 1 + X 3 1, X 1 + X 2 3 0 3 4 7 4 10 324 4 6 0 3 4 3 1 6 0 6, 1 6 0 1, 1 6 0 7 1 4 1 2 6 3 0 1 4 1 2 7 6 4 1 2 5 6 0 1 1 6 0 1 7 6 0 1 6 ance sum Minimum variance sum X 1 + X 1 + X 2 X 2 + X 3 X 3 0 3 40 3 4 7 7 1 6 0 4 1 2 N = 7 1 6 0 (6) 4 1 2 S N = 7 7 7 6 0 1 7 6 0 1 7 Outline Implied Correlation Toy examples Algorithm Implied Index Toy example Conclusions 0 1 4 6 The minimum variance of the sum is equal to 0! (ideal case of a constant sum (complete mixability, see Wang and Wang (2011)) Carole Bernard Implied Systemic Risk 23/35

Block Rearrangement Algorithm With more than 3 variables, we can improve the standard algorithm (which proceeds column by column) by proceeding by block: Split d columns into two subsets Π and Π. and make sure that i Π X i is in reverse order with i Π X i (working paper with D. McLeish) 1 In general, many local minima for the variance of the sum: 2 By starting with a random initial matrix, and reproducing the experience several times, we are able to approximate the set of all copulas that minimize the variance of the sum. Carole Bernard Implied Systemic Risk 24/35

Inputs: Using the Block RA to infer the dependence ˆ X 1 F 1,... X d F d ˆ the cdf of ω 1 X 1 +... + ω d X d G is known for some ω i R Question Describe the set of possible dependence structures (copulas) that are consistent with this information. Carole Bernard Implied Systemic Risk 25/35

Inputs: Method: Block RA to infer the dependence ˆ X 1 F 1,... X d F d ˆ X 1 +... + X d G Method: ˆ Matrix of n rows (for discretization step) by d + 1 columns. ˆ In each of the first d columns F 1 j ( i n + 1 ), i = 1, 2,..., n ˆ In the last column ( ) i G 1, i = 1, 2,..., n n + 1 ˆ Apply the Block RA on the full matrix Output: Extract the d first columns, and they describe a discrete copula that is consistent with the information on the cdfs of the risks and of their sum. Carole Bernard Implied Systemic Risk 26/35

A proposed global dependence measure Instead of Pearson correlation, we can use Spearman s rho ϱ ij := Spearman s rho(x i, X j ) = ρ(f i (X i ), F j (X j )) (correlation between the ranks) It is not affected by changes in marginal distributions (and thus not sensitive to changes in the volatility parameter) We can consider the average pairwise Spearman s rho ϱ := d 1 i=1 d 1 i=1 j>i ω iω j ϱ ij j>i ω iω j Compared to the CBOE, it is not affected by changes in the volatilities of the individual components of the index. Carole Bernard Implied Systemic Risk 27/35

Empirical work (coming soon) ˆ From option prices on Dow Jones 30 (use all strikes) to estimate the marginal distribution of index ˆ From option prices on components of Dow Jones 30 to estimate their marginal distributions ˆ Compare this proposal with the CBOE index. 1 Not affected by changes in volatility 2 Use full information from option prices 3 Same empirical conclusions? Similar to the VIX, implied correlation exhibits a tendency to increase when the S&P 500 decreases. Carole Bernard Implied Systemic Risk 28/35

An Implied Systemic Risk Measure Of interest to go beyond a global measure of dependence. Systemic risk measurement is closely related to ˆ contagion effects in the tail ˆ extreme events / coexceedances (tail dependence) Our approach allows to study the dependence in the tail. A natural measure to study is for example a pairwise average of ϱ tail ij := Spearman s rho(x i, X j scenarios) These scenarios can be driven by the aggregate risk or some sector, or some individual institutions. Carole Bernard Implied Systemic Risk 29/35

More empirical work (coming soon) ˆ From option prices on DJ 30 and on its components (use all strikes) estimate marginal distributions ˆ Study systemic risk contribution of each of the 30 institutions within the DJ 30 ] E Q [Xi X i < quantile i and check whether the order is consistent with what is found under the real-world probability measure (same spirit as SES of Acharya et al. (2010), SRisk of Brownlees and Engle (2014)...) Carole Bernard Implied Systemic Risk 30/35

Using the Block RA to infer the dependence Example: start with a situation for which we know the dependence, and see if we can recover this information. A one-period financial market ˆ with maturity T. ˆ with two assets LogNormally distributed (as in Black-Scholes) with r = 0.01, σ 1 = 15% and σ 2 = 40%. S 1 0 = S 2 0 = 1. ˆ a Clayton copula with parameter 3. ˆ obtain the distribution of the sum G by simulation Carole Bernard Implied Systemic Risk 31/35

S 2 2.5 2 1.5 1 0.5 Information obtained by simulation 0.5 1 1.5 2 2.5 S 1 F 2 (S 2 ) (rank) 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 F 1 (S 1 ) (rank) Carole Bernard Implied Systemic Risk 32/35

S 2 2.5 2 1.5 1 0.5 Information obtained by simulation 0.5 1 1.5 2 2.5 S 1 F 2 (S 2 ) (rank) 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 F 1 (S 1 ) (rank) Pearson correlation = ρ (X 1, X 2 ) 0.78 Denote q S α = Quantile α (X 1 + X 2 ) ρ {X1,X 2 X 1+X 2 q S 25%} 0.81 ; ρ {X 1,X 2 X 1+X 2 q S 75%} 0.15 ρ {X1,X 2 X 1+X 2 [q S ]} 0.26 25%,qS 75% Carole Bernard Implied Systemic Risk 32/35

Information that we obtain using the information on F 1, F 2 and G ONLY and the BRA ran 500 times Pearson correlation = ρ (X 1, X 2 ) [0.7800, 0.7801] q S α = Quantile α (X 1 + X 2 ) ρ {X1,X 2 X 1 +X 2 q S 25%} ρ {X1,X 2 X 1 +X 2 q S 75%} [0.813, 0.818] [ 0.15, 0.14] ρ {X1,X 2 X 1 +X 2 [q25% S,qS ]} [0.24, 0.27] 75% Carole Bernard Implied Systemic Risk 33/35

but The method works in higher dimensions ˆ Not able to reproduce a single pairwise correlation, especially if X 1,... X d have same marginal distributions. ˆ But able to reproduce an average correlation, an average tail correlation... ˆ Intervals are wider in higher dimensions than in two dimensions because of uncertainty on the copula even if one knows the distribution of the sum. Carole Bernard Implied Systemic Risk 34/35

Conclusions & Research Directions Develop an efficient algorithm for inferring the dependence among variables for which we know the marginal distributions and the distribution of a weighted sum. Use it to develop new indicators of implied dependence among assets, richer than the implied correlation from CBOE. We hope to find an indicator that is forward looking and can have some predictive power... Carole Bernard Implied Systemic Risk 35/35

Conclusions & Research Directions Develop an efficient algorithm for inferring the dependence among variables for which we know the marginal distributions and the distribution of a weighted sum. Use it to develop new indicators of implied dependence among assets, richer than the implied correlation from CBOE. We hope to find an indicator that is forward looking and can have some predictive power... Thank You Carole Bernard Implied Systemic Risk 35/35

References Aït-Sahalia, Yacine, and Andrew W. Lo. Nonparametric estimation of state?price densities implicit in financial asset prices. Journal of Finance 53.2 (1998): 499-547. Bernard, C., X. Jiang, and R. Wang (2014): Risk Aggregation with Dependence Uncertainty, Insurance: Mathematics and Economics. Bernard, C., McLeish D. (2014): Algorithms for Finding Copulas Minimizing the Variance of Sums, Working Paper. Bernard, C., L. Rüschendorf, and S. Vanduffel (2014): VaR Bounds with a Variance Constraint, Working Paper. Bondarenko, Oleg. Estimation of risk-neutral densities using positive convolution approximation. Journal of Econometrics 116.1 (2003): 85-112. Bernard, C., Vanduffel S. (2014): A new approach to assessing model risk in high dimensions, available on SSRN. Embrechts, P., G. Puccetti, and L. Rüschendorf (2013): Model uncertainty and VaR aggregation, Journal of Banking & Finance. Haus U.-U. (2014): Bounding Stochastic Dependence, Complete Mixability of Matrices, and Multidimensional Bottleneck Assignment Problems, Working Paper. Puccetti, G., and L. Rüschendorf (2012): Computation of sharp bounds on the distribution of a function of dependent risks, Journal of Computational and Applied Mathematics, 236(7), 1833 1840. Wang, B., and R. Wang (2011): The complete mixability and convex minimization problems with monotone marginal densities, Journal of Multivariate Analysis, 102(10), 1344 1360. Carole Bernard Implied Systemic Risk 36/35

Acknowledgements ˆ GRI in Financial Services and the Louis Bachelier Institute ˆ BNP Paribas Fortis Chair in Banking. ˆ Research project on Risk Aggregation and Diversification with Steven Vanduffel for the Canadian Institute of Actuaries. ˆ Humboldt Research Foundation. ˆ Project on Systemic Risk funded by the Global Risk Institute in Financial Services. ˆ Natural Sciences and Engineering Research Council of Canada ˆ Society of Actuaries Center of Actuarial Excellence Research Grant Carole Bernard Implied Systemic Risk 37/35

Change of regime independent from X 1, X 2 and Z I takes the value 1 with probability q (independent from X 1, X 2 and Z), σ 1 = 0.1, σ 2 = 0.2, σ Z = 0.4, ρ 12 = 0.2. 1 S i = (1-I)X i +I Z 0.9 0.8 Dependence index 0.7 0.6 0.5 0.4 0.3 CBOE index Pearson correlation between logreturns of S 1 and S 2 in the model Average Pearson correlation by rearrangement (LN margins) 0.2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 probability q of I=1 in the mixture Back Carole Bernard Implied Systemic Risk 38/35

Change of regime independent from X 1, X 2 and Z How to explain the discrepancy between the CBOE index and the actual correlation in the model? Carole Bernard Implied Systemic Risk 39/35

Change of regime independent from X 1, X 2 and Z How to explain the discrepancy between the CBOE index and the actual correlation in the model? S 1, S 2, S 1 + S 2 are not lognormally distributed... q = 0.25 q = 0.5 q = 0.75 0.035 0.025 0.018 Probability Distribution 0.03 0.025 0.02 0.015 0.01 0.005 density of S 1 density of S 2 density of index Probability Distribution 0.02 0.015 0.01 0.005 density of S 1 density of S 2 density of index Probability Distribution 0.016 0.014 0.012 0.01 0.008 0.006 0.004 0.002 density of S 1 density of S 2 density of index 0 0 50 100 150 200 250 300 S 0 0 50 100 150 200 250 300 S 0 0 50 100 150 200 250 300 S Back Carole Bernard Implied Systemic Risk 39/35