arxiv: v2 [q-fin.st] 9 Jul 2015

Similar documents
Asymmetric Dependence, Tail Dependence, and the. Time Interval over which the Variables Are Measured

A simple graphical method to explore tail-dependence in stock-return pairs

Behaviour of multivariate tail dependence coefficients

A simple nonparametric test for structural change in joint tail probabilities SFB 823. Discussion Paper. Walter Krämer, Maarten van Kampen

Copulas and dependence measurement

arxiv: v2 [q-fin.st] 10 May 2013

The Instability of Correlations: Measurement and the Implications for Market Risk

Tail Dependence of Multivariate Pareto Distributions

Trivariate copulas for characterisation of droughts

Simulation of Tail Dependence in Cot-copula

A Goodness-of-fit Test for Copulas

Dependence Patterns across Financial Markets: a Mixed Copula Approach

Risk Measures with Generalized Secant Hyperbolic Dependence. Paola Palmitesta. Working Paper n. 76, April 2008

Financial Econometrics and Volatility Models Copulas

Markov Switching Regular Vine Copulas

Quaderni di Dipartimento. Small Sample Properties of Copula-GARCH Modelling: A Monte Carlo Study. Carluccio Bianchi (Università di Pavia)

Copulas, a novel approach to model spatial and spatio-temporal dependence

EVANESCE Implementation in S-PLUS FinMetrics Module. July 2, Insightful Corp

Empirical properties of large covariance matrices in finance

Quasi-copulas and signed measures

2 (U 2 ), (0.1) 1 + u θ. where θ > 0 on the right. You can easily convince yourself that (0.3) is valid for both.

A measure of radial asymmetry for bivariate copulas based on Sobolev norm

Multivariate Asset Return Prediction with Mixture Models

ASSOCIATIVE n DIMENSIONAL COPULAS

When is a copula constant? A test for changing relationships

CONTAGION VERSUS FLIGHT TO QUALITY IN FINANCIAL MARKETS

Copulas. Mathematisches Seminar (Prof. Dr. D. Filipovic) Di Uhr in E

Multivariate Non-Normally Distributed Random Variables

Gaussian Slug Simple Nonlinearity Enhancement to the 1-Factor and Gaussian Copula Models in Finance, with Parametric Estimation and Goodness-of-Fit

arxiv: v2 [math.pr] 23 Jun 2014

Dependence and VaR Estimation:An Empirical Study of Chinese Stock Markets using Copula. Baoliang Li WISE, XMU Sep. 2009

Tail dependence in bivariate skew-normal and skew-t distributions

Determining and Forecasting High-Frequency Value-at-Risk by Using Lévy Processes

On the Estimation and Application of Max-Stable Processes

MAXIMUM ENTROPIES COPULAS

Simulation of High-Dimensional t-student Copulas with a Given Block Correlation Matrix

x. Figure 1: Examples of univariate Gaussian pdfs N (x; µ, σ 2 ).

arxiv:physics/ v2 [physics.soc-ph] 8 Dec 2006

When is a copula constant? A test for changing relationships

Marginal Specifications and a Gaussian Copula Estimation

Multivariate Distributions

Robustness of a semiparametric estimator of a copula

Asymptotic distribution of the sample average value-at-risk

Normal Probability Plot Probability Probability

How to select a good vine

Quantitative Methods in High-Frequency Financial Econometrics:Modeling Univariate and Multivariate Time Series

Notes on Random Variables, Expectations, Probability Densities, and Martingales

A Measure of Monotonicity of Two Random Variables

Framework for Analyzing Spatial Contagion between Financial Markets

Dependence. Practitioner Course: Portfolio Optimization. John Dodson. September 10, Dependence. John Dodson. Outline.

Convolution Based Unit Root Processes: a Simulation Approach

Modelling and Estimation of Stochastic Dependence

arxiv:physics/ v1 [physics.soc-ph] 18 Aug 2006

VaR vs. Expected Shortfall

Copula based Probabilistic Measures of Uncertainty with Applications

arxiv: v1 [physics.data-an] 7 Mar 2015

Tail Approximation of Value-at-Risk under Multivariate Regular Variation

Tail dependence coefficient of generalized hyperbolic distribution

Quantifying volatility clustering in financial time series

X

Solutions of the Financial Risk Management Examination

The Use of Copula Functions in Pricing Weather Contracts for the California Wine Industry

A Non-Parametric Approach of Heteroskedasticity Robust Estimation of Vector-Autoregressive (VAR) Models

Copula-based top-down approaches in financial risk aggregation

arxiv: v1 [q-fin.st] 5 Apr 2007

Correlation: Copulas and Conditioning

Imputation Algorithm Using Copulas

A copula goodness-of-t approach. conditional probability integral transform. Daniel Berg 1 Henrik Bakken 2

Multivariate Normal-Laplace Distribution and Processes

Songklanakarin Journal of Science and Technology SJST R1 Sukparungsee

Asymmetric CAPM dependence for large dimensions: the Canonical Vine Autoregressive Model

Three-Stage Semi-parametric Estimation of T-Copulas: Asymptotics, Finite-Samples Properties and Computational Aspects

Overview of Extreme Value Theory. Dr. Sawsan Hilal space

ECONOMICS 7200 MODERN TIME SERIES ANALYSIS Econometric Theory and Applications

Multivariate Skew Distributions Based on the GT-Copula *

Accounting for extreme-value dependence in multivariate data

Efficient estimation of a semiparametric dynamic copula model

Generalized Autoregressive Score Models

LOCAL TIMES OF RANKED CONTINUOUS SEMIMARTINGALES

Asymmetry in Tail Dependence of Equity Portfolios

Asymptotic behaviour of multivariate default probabilities and default correlations under stress

Estimation of direction of increase of gold mineralisation using pair-copulas

A New Generalized Gumbel Copula for Multivariate Distributions

Price asymmetry between different pork cuts in the USA: a copula approach

arxiv: v3 [q-fin.gn] 22 Nov 2010

ARANDOM-MATRIX-THEORY-BASEDANALYSISOF STOCKS OF MARKETS FROM DIFFERENT COUNTRIES

Probabilities & Statistics Revision

arxiv:physics/ v1 [physics.soc-ph] 2 Aug 2006

Copulas and Measures of Dependence

Shape of the return probability density function and extreme value statistics

arxiv: v3 [q-fin.tr] 3 Apr 2018

Reducing Model Risk With Goodness-of-fit Victory Idowu London School of Economics

Estimating Global Bank Network Connectedness

VaR bounds in models with partial dependence information on subgroups

Modeling Ultra-High-Frequency Multivariate Financial Data by Monte Carlo Simulation Methods

Eigenvalue spectra of time-lagged covariance matrices: Possibilities for arbitrage?

arxiv:physics/ v2 [physics.soc-ph] 22 Apr 2007

TESTING FOR STRUCTURAL CHANGES IN EXCHANGE RATES DEPENDENCE BEYOND LINEAR CORRELATION. July 16, 2007

Dependence modelling of the joint extremes in a portfolio using Archimedean copulas : application to MSCI indices

An Econometric Study of Vine Copulas

Transcription:

Dependence structure of market states arxiv:153.94v [q-fin.st] 9 Jul 15 Desislava Chetalova, Marcel Wollschläger and Rudi Schäfer Fakultät für Physik, Universität Duisburg Essen, D 4748 Duisburg, Germany E-mail: desislava.chetalova@uni-due.de, marcel.wollschlaeger@uni-due.de, rudi.schaefer@uni-due.de Abstract. We study the dependence structure of market states by estimating empirical pairwise copulas of daily stock returns. We consider both original returns, which exhibit time-varying trends and volatilities, as well as locally normalized ones, where the non-stationarity has been removed. The empirical pairwise copula for each state is compared with a bivariate K-copula. This copula arises from a recently introduced random matrix model, in which non-stationary correlations between returns are modeled by an ensemble of random matrices. The comparison reveals overall good agreement between empirical and analytical copulas, especially for locally normalized returns. Still, there are some deviations in the tails. Furthermore, we find an asymmetry in the dependence structure of market states. The empirical pairwise copulas exhibit a stronger lower tail dependence, particularly in times of crisis.

1. Introduction The concept of copulas was introduced by Sklar in 1959 [1, ] to study the linkage between multivariate distribution functions and their univariate marginals. Since then, copulas have gained growing importance as a tool for modeling statistical dependence of random variables in many fields. In finance, the usage of copulas is relatively new but it has already found application in risk management, see e.g., [3, 4, 5, 6, 7, 8, 9], derivative pricing, see e.g., [1, 11, 1, 13, 14, 15], and portfolio optimization, see e.g., [16, 17, 18, 19]. For an overview of the literature on applications of copulas in finance the reader is referred to [, 1]. Copulas allow to separate the dependence structure of random variables from their marginal distributions. This is sometimes useful in statistical applications as the dependence structure and the marginal distributions can be modeled separately and joined together resulting in new multivariate distributions with different behavior. For a discussion on difficulties in the application of copulas the reader is referred to [, 3, 4]. Here, we simply view the copulas as providing a standardized way for the study of statistical dependences. The marginal distributions are mapped to the uniform distribution; the statistical dependence is considered in terms of the marginal cumulative distribution functions. Recently, we identified market states as clusters of similar correlation matrices and studied their corresponding correlation structures [5]. The correlation structure, however, does not capture the full statistical dependence between financial return time series. Here, we choose a copula approach to study the dependence structure of market states [6]. To this end, we estimate empirical copulas for many stock return pairs and average over all of them to obtain an empirical pairwise copula for each market state. We stress that the identification of market states relies on the correlation matrices, the copulas are used only for analyzing the states and not defining them. To estimate the empirical copulas we use both original and locally normalized returns. The original return time series exhibit time-varying trends and volatilities [7, 8, 9]. These have to be taken into account when transforming the marginal distributions to uniform ones. To this end, we apply the method of local normalization [3], which leads to stationary time series while preserving the correlations between them. The resulting empirical copulas provide different information. The copulas for the original returns describe the dependence structure on a global scale, i.e., for the full time horizon, whereas the copulas for the locally normalized returns describe the dependence structure on a local scale. The empirical pairwise copulas for each market state are compared with a bivariate K-copula, which arises from a random matrix approach introduced in [31, 3]. It models the non-stationarity of true correlations by an ensemble of random matrices. The model yields a multivariate return distribution in terms of a modified Bessel function of the second kind, a so-called K-distribution. In [5] the K-distribution was found to provide a good description of the heavy-tailed empirical return distributions for each market state. Here, we aim to arrive at a consistent description within the random matrix model studying the agreement between K-copula and empirical dependence structure for each market state. In addition, our study provides further evidence for asymmetric dependencies between financial returns [33, 34, 35]. We find an asymmetry in the tail dependence of empirical pairwise copulas which we study in more detail. The paper is organized as follows. In section we review the basic concepts of copulas, stating the main result in the copula theory, the Sklar s theorem, which

we use to derive a K-copula. In section 3 we present the data set and recapitulate the identification of market states for the Nasdaq Composite market in the period 199 13 as done in [5]. In section 4 we study the empirical copula densities for each market state and compare them with the K-copula. We conclude our findings in section 5.. Copula We begin with a short introduction to the concept of copulas in.1. For more details with an emphasis on the statistical and mathematical foundations of copulas the reader is referred to the textbooks of Joe [36] and Nelsen [37]. In. we present the K-copula which plays a central role in this study. We restrict ourselves to the bivariate case, since later on we study empirical pairwise copulas..1. Basic Concepts Consider two random variables X and Y. The joint distribution of X and Y contains all the statistical information about them. It can be expressed either in terms of the joint probability density function (pdf) f X,Y (x, y) or in terms of the joint cumulative distribution function (cdf) F X,Y (x, y), where F X,Y (x, y) = x dx y dy f X,Y (x, y ). (1) From the joint pdf f X,Y (x, y) one can extract the individual distributions of X and Y as follows f X (x) = dy f X,Y (x, y), () and analogously for Y. The densities f X (x) and f Y (y), called marginal probability density functions, and the corresponding marginal cumulative distribution functions F X (x) and F Y (y) describe the individual statistical behavior of the random variables. When dealing with correlated random variables, one is interested in their statistical dependence. The Pearson correlation coefficient is commonly used as a measure of dependence. It is defined as C X,Y = Cov(X, Y ) σ X σ Y, (3) where Cov(X, Y ) is the covariance of both random variables and σ X and σ Y are the respective standard deviations. However, the correlation coefficient only measures the linear dependence between the random variables. Copulas provide a natural way to study the statistical dependence of random variables. The transformation U i = F i (i) i = X, Y (4) leads to new random variables with uniform distributions on the unit interval, called the rank of X and Y, respectively. Their joint distribution is called a copula. It 3

describes the dependence structure of the random variables X and Y separated from their marginal distributions. A central result in the copula theory is Sklar s theorem which enables us to separate any multivariate distribution function into two components: the marginal distributions of each random variable and their copula F X,Y (x, y) = Cop X,Y (F X (x), F Y (y)). (5) If the marginal distribution functions are continuous, the copula that satisfies equation (5) is given by Cop X,Y (u, v) = F X,Y (F 1 1 (u), F (v)), (6) X where F 1 1 X and FY represent the inverse cumulative distribution functions, the so called quantile functions. This equation allows to extract the dependence structure directly from the joint distribution function. From the copula (6) one can compute the copula density as follows.. K-copula cop X,Y (u, v) = Y u v Cop X,Y (u, v). (7) The K-copula arises from a random matrix model introduced in [31] to model timevarying correlations between financial time series [38, 39]. It was first used in [6], where it was found to describe the empirical dependencies in financial data much better than a Gaussian copula. Consider a market consisting of K stocks. At each time t (t = 1,..., T ) we assume that the return vector r(t) = (r 1 (t),..., r K (t)) is drawn from a multivariate normal distribution with a covariance matrix Σ t g(r Σ t ) = 1 det πσt exp ( 1 r Σ 1 t r ) 4, (8) where we suppress the argument t of r to simplify our notation. We now model the time-dependent covariance matrix Σ t by a Wishart random matrix AA, where the K N model matrix A is drawn from a Gaussian distribution with the pdf w(a Σ, N) = KN N π 1 det Σ N exp ( N tr A Σ 1 A ). (9) Here, Σ represents the average covariance matrix estimated over the sample of all r(t), t = 1,..., T. Hence, the time-dependent covariance matrices are modeled by an ensemble of Wishart matrices AA which fluctuate around the sample average Σ. Averaging the multivariate normal distribution (8) with the random covariance matrix AA over the Wishart ensemble leads to a K-distribution for the multivariate returns g (r Σ, N) = d[a] w(a Σ, N) g(r AA ) 1 = (π) K Γ(N/) det Σ N K N = Γ(N/) det(πσ) K K N dz z N πn 1 e z z ( Nr Σ r) 1 Nr Σ 1 r K N, K ( exp N ) 4z r Σ 1 r (1)

where K ν is the modified Bessel function of the second kind of order ν = (K N)/. It depends only on the average covariance matrix Σ = σcσ estimated over the whole sample, where C is the average correlation matrix and σ = diag (σ 1,..., σ K ) the diagonal matrix of the standard deviations, and a free parameter N, which governs the variance of the Wishart ensemble var([aa ] kl ) = Σ kl + Σ kkσ ll N 5, (11) where Σ kl is the kl-th element of the average covariance matrix Σ. Thus, it characterizes the strength of fluctuations around the average covariance matrix in the considered sample. The larger N, the smaller the fluctuations around Σ, eventually vanishing in the limit N. The K-copula is the dependence structure which arises for the K-distribution (1). For the bivariate case K =, the pdf of the vector r = (r 1, r ) reads f c,n (r 1, r ) = g (r Σ, N) = 1 Γ(N/) dz zn/ 1 e z 1 c Here, we used the covariance matrix ( ) σ Σ = 1 σ 1 σ c σ 1 σ c σ = ( N 4πz exp N 4z ( ) 1 c c 1 r1 cr 1 r + r ) 1 c. (1), (13) where c denotes the average correlation coefficient, estimated over the whole sample. We chose the standard deviations one, σ 1 = σ = 1, since the copula is independent of the marginal distributions. Then, the marginal distribution densities are identical, f 1 (r 1 ) = f (r ), where f 1 (r 1 ) = dr f c,n (r 1, r ) = 1 Γ(N/) According to equation (6) the bivariate K-copula is given by ( N dz z N/ 1 e z 4πz exp N ) 4z r 1. (14) Cop c,n (u, v) = F c,n (F 1 1 (u), F 1 (v)), (15) where c and N are the parameters of the copula, F c,n denotes the cumulative distribution function of the bivariate distribution (1), and F 1 the inverse distribution function of the marginal cdf given by F c,n (r 1, r ) = = = r 1 r 1 r 1 dξ dξ dξ r r r dζ f c,n (ξ, ζ) dζ dz Γ(N/) z N/ 1 e z 1 c N dζ N N(ξ cξζ+ζ ) 1 c π Γ(N/) N 1 c ( N 4πz exp N 4z K N ξ cξζ + ζ ) 1 c ( ) N(ξ cξζ + ζ ) 1 c (16)

6 and F 1 (r 1 ) = = r 1 r 1 dξ f 1 (ξ) = dξ r 1 N Nξ N 1 1 dξ Γ(N/) π Γ(N/) N 1 K 1 N N dz z N/ 1 e ( z 4πz exp N4z ) ξ ( Nξ ). (17) The K-copula density can be obtained from the K-copula (15) by differentiation cop c,n (u, v) = Cop c,n (u, v) u v = f c,n(f1 1 (u), F 1 (v)) f 1 (F1 1 (u))f (F 1 (v)). (18) It depends only on the average correlation c and the free parameter N, which characterizes the fluctuations around Σ. Figure 1 shows the K-copula density for different parameter values. The stronger the average correlation c and the lower the parameter N, the higher is the probability for extreme co-movements. Furthermore, we note that the K-copula is a symmetric copula. It is based on the elliptical distribution (1) and thus it belongs to the class of elliptical copulas. 3. Identification of market states We now present the data set and identify market states as clusters of correlation matrices with similar correlation structures, as done in [5]. These market states will be the object of our empirical study in the next section. We consider K = 58 stocks of the Nasdaq Composite index traded in the -year period from January 199 to December 13 [4] which corresponds to 554 trading days. For each stock k we calculate the return time series r k (t) = S k(t + t) S k (t) S k (t), k = 1,..., K, (19) where S k (t) is the price of the k-th stock at time t and t is the return interval which we chose to be one trading day. Empirical return time series exhibit time-varying drift and volatilities. The correlation coefficient (3) averages over these time-dependent parameters which results in an estimation error of the correlations. In order to eliminate this kind of error we employ the method of local normalization [3]. For each return time series k we subtract the local mean and divide by the local standard deviation ˆr k (t) = r(t) r(t) n r (t) n r(t) n, () where... n denotes the local average which runs over the n most recent sampling points. For daily data we use n = 13 as discussed in [3]. The local normalization removes the local trends and variable volatilities while preserving the correlations between the time series. Using the locally normalized daily returns we now obtain a set of 131 correlation matrices estimated on disjoint two-month intervals of the -year observation period.

7 Figure 1: K-copula densities copc,n (u, v) for different parameter values. Top: N = 5, (left) c = and (right) c =.5, bottom: c =., (left) N = 3 and (right) N = 3. To identify market states we perform a clustering analysis based on the partitioning around medoids algorithm [41], where the number of clusters is estimated via the gap statistic [4]. The clustering analysis separates the set of 131 correlation matrices into six groups based on the similarity of their correlation structures. Each group is associated with a market state. We point out that the identification of market states is performed ex-post, the clustering algorithm has the correlation matrices of all times. Figure shows the time evolution of the six market states. We observe that the market switches back and forth between states. Sometimes it remains in a state for a long time, sometimes it jumps briefly to another state and returns back or evolves further. On longer time scales the market evolves towards new states, whereas previous states die out. In the next section we will study the statistical dependence for each market state. To obtain the return time series for each state, we proceed as follows: We take the complete return time series, r(t) or r (t), and divide it into a sequence of disjoint two month intervals. We merge all intervals belonging to a given state according to the cluster analysis. We note that the return time series for the six market states differ

8 market state 6 5 4 3 1 199 1996 4 8 1 Figure : Time evolution of the market in the observation period 199 13. Each point represents a correlation matrix measured over a two-month time window. in length. 4. Empirical results We present the empirical pairwise copula densities for each market states in 4.1 and compare them with the K-copula densities in 4.. In 4.3 we study the asymmetry of the tail dependence of the empirical copula densities in more detail. 4.1. Empirical pairwise copulas for each market state To estimate the empirical pairwise copula of two return time series r k (t) and r l (t), we first have to transform them into uniformly distributed time series. To achieve this, we employ the empirical distribution function u k (t) = F k (r k (t)) = 1 T T τ=1 1{r k (τ) r k (t)} 1 T, (1) where 1 is the indicator function, T denotes the length of the time series and the factor 1/ ensures that the values of the transformed time series u k (t) lie in the interval (, 1). The empirical copula density of the time series r k (t) and r l (t) is then the two-dimensional histogram of the transformed time series u = u k (t) and v = u l (t). An accurate estimation of the copula density requires a large amount of data. Thus, for each state we compute the copula densities of all K(K 1)/ stock pairs as two-dimensional histograms of the transformed time series and then average over all pairs cop (i) (u, v) = K 1 K(K 1) K k=1 l=k+1 cop (i) (u, v), i = 1,..., 6, () where the upper index i denotes the state number. For the bin size of the histograms we choose u = v =.5. We note that as the length of the time series for each market state is rather short, we are not able to study the full copulas for each stock pair k, l separately. Hence, we cannot make any direct statements about similarity and

dispersion regarding the full dependence structure. Our results only yield statements about the empirical dependence structure on average. Figure 3 shows the empirical pairwise copula densities for the original returns (19) for each of the six market states. We observe a variation of the dependence structure from state to state, particularly visible in the tails. In state 1, which covers the period from 199 to roughly, we find a rather flat copula density, indicating low dependence between return pairs. In state we observe deviations from the flat copula density particularly in the tails, which become more and more pronounced in state 3 and 4. State 5, first appearing during the financial crisis in 8, exhibits the strongest dependence. The dependence decreases again in state 6. Figure 4 shows the empirical pairwise copula densities for the locally normalized returns (). The dependence structures of the six states are mostly preserved after performing local normalization. Deviations are observed in the lower-left and the upper-right corners where the copula densities for the original returns exhibit higher peaks. The reason for that is the time-varying volatility. During periods of high volatility stocks tend to have large returns which contribute to the corners of the copula density. On the other hand, in periods of high volatility the correlations between stocks become stronger. This leads to higher peaks in the corners of the dependence structure for the original returns. It is important to note that the copulas for the original and the locally normalized returns contain different statements. The copulas for the original returns describe the dependence structure for the full time horizon. On the other hand, the copulas for the locally normalized returns provide information about the statistical dependence on a local scale. Furthermore, we observe that the empirical copula densities are asymmetric with respect to opposite corners. We find stronger dependence in the lower tail than in the upper one, that is, the dependence between large negative returns is stronger than the dependence between large positive ones. This asymmetry is an important feature of empirical copula densities and thus we will discuss it in more details in section 4.3. 4.. Comparison with the K-copula We now compare the empirical pairwise copula densities for each market state with the K-copula densities. Again, we consider both the original and the locally normalized returns. The K-copula density is obtained in the following way: We calculate the K-copula according to equation (15) where the integrals are performed numerically. The K-copula density for each bin of size u = v =.5 is then estimated by 9 cop c,n (u, v) = Cop c,n (u, v) Cop c,n (u, v v) Cop c,n (u u, v) + Cop c,n (u u, v v). (3) For comparison we compute the difference between empirical and analytical copula density for each state cop (i) (u, v) cop (i) c,n (u, v), i = 1,..., 6, (4) where c is the average correlation coefficient of all K(K 1)/ stock pairs for the considered state. The free parameter N is estimated by a fit which minimizes the mean squared difference between empirical and analytical copula density. The parameter values for each market state are summarized in table 1. The differences between the

Figure 3: Empirical pairwise copula density cop (i) (u, v) for the original returns. 1

Figure 4: Empirical pairwise copula density cop (i) (u, v) for the locally normalized returns. 11

1 returns state 1 state state 3 state 4 state 5 state 6 original c.46.13.17.5.4. N 41.7 11.7 8.4 5.6.8 1. loc. normalized c.48.13.19.8.43.5 N 7.7 8.6 9.8 15.8 7.4.4 Table 1: Parameters of the K-copula density for the original and the locally normalized returns. returns state 1 state state 3 state 4 state 5 state 6 original.11.41 1..47 5.53.84 loc. normalized.47.38.7 1.43.94 1.5 Table : Mean squared differences between empirical and K-copula densities. empirical copula density and the K-copula density for each state are presented in figure 5 for the original and in figure 6 for the locally normalized returns. Overall, we find a good agreement between empirical and analytical copula densities. The K-copula seems to capture the dependence structure of the first three states very well. Small deviations from the K-copula density are observed for state 4 and 6. Only the dependence structure of state 5 cannot be captured by the K-copula. For the locally normalized returns we find a better agreement, which is reflected in the smaller mean squared differences, see table. Deviations are observed mainly in the corners of the copula densities. The empirical copula densities exhibit stronger dependence in the tails, i.e., the probability for extreme co-movements is underestimated by the K-copula. It is important to note that the K-copula captures the empirical dependence structure rather well with only one free parameter. However, due to its symmetric nature it cannot account for the asymmetry observed in the data. The skewed Student s t-copula is an alternative proposed by [43] which was found to account for asymmetric dependencies in financial data [44, 45]. It captures the empirical dependence structure of the original returns better than the K-copula due to the presence of an additional parameter which accounts for the asymmetry [6]. Nevertheless, here we confine ourselves to the comparison of the empirical copulas with the K-copula as we aim to arrive at a consistent picture within the random matrix model. We discuss the asymmetry in more detail in the following. 4.3. Asymmetry of the tail dependence Asymmetric dependence between returns has been reported by several authors, see e.g., [33, 34, 35]. Our study provides further evidence revealing a stronger lower tail dependence in the empirical copula densities of market states. We now take a closer look at this asymmetry. To this end, we estimate the tail dependence in the four

13 Figure 5: Difference between the empirical copula density and the K-copula density cop (i) (u, v) cop (i) c,n (u, v) for the original returns.

14 Figure 6: Difference between the empirical copula density and the K-copula density cop (i) (u, v) cop (i) c,n (u, v) for the locally normalized returns.

15 corners for all K(K 1)/ empirical pairwise copulas according to. LL (i) =.8 du. 1 UL (i) = du 1 UU (i) =.8. LU (i) = du du. 1.8 1.8 dv cop (i) (u, v), dv cop (i) (u, v), dv cop (i) (u, v), dv cop (i) (u, v), (5) where the upper index i = 1,..., 6 denotes the state number and the lower indices represent a stock pair k, l. Here, LL and UU refer to the lower-lower and upperupper corners, respectively, and represent the positive tail dependence, whereas UL and LU refer to the upper-lower and lower-upper corners, respectively, and represent the negative tail dependence. The asymmetry in the tail dependence can now be quantified by the differences α (i) = UU(i) LL(i), β (i) = LU(i) UL(i), i = 1,..., 6, (6) where α (i) captures the asymmetry of the positive and β (i) the asymmetry of the negative tail dependence for each stock pair. Figure 7 shows the histograms of the asymmetry values exemplarily for state 5. For the returns we find a negative offset for the values of α (5), whereas the values of β(5) are centered around zero. This indicates, on average, an asymmetry in the positive tail dependence, i.e., simultaneous large negative returns are more likely to occur than simultaneous large positive returns. On the other hand, we do not find such an asymmetry in the negative tail dependence. For the locally normalized returns, we find much weaker asymmetry in the positive tail dependence and once again no asymmetry for the negative tail dependence. We note that the means of the asymmetry values are the relevant quantities, the asymmetry values for each pair are distributed around the mean due to statistical fluctuations. The standard deviation for α (5) is.1 and for β (5).7. Indeed, the asymmetry effect is very small. Still, it is clearly visible, see figure 3. In the following, we study the asymmetry values for each market state. Figure 8 shows the mean asymmetry values ᾱ (i) for each market state, obtained by averaging over all α (i) for a given state. On a local scale the asymmetry in the positive tail dependence is much weaker. Only state 5 still exhibits a certain amount of asymmetry. On the other hand, the asymmetry in the negative tail dependence is negligibly small for both original and locally normalized returns.

16 8 8 6 6 pdf 4 pdf 4.6.4....4 8.6.4....4 8, 6 6 pdf 4 pdf 4.6.4....4.6.4....4 Figure 7: Histograms of the asymmetry values for all stock pairs, exemplarily for state 5. Left: asymmetry for the positive tail dependence α (5), right: asymmetry for the negative tail dependence: β (5). Top: for returns, bottom: for locally normalized returns. The dashed red lines represent the corresponding mean values. 5. Conclusion We studied the dependence structure of market states by means of a copula approach. To this end, we estimated the empirical pairwise copulas for each state and compared them with the bivariate K-copula. The bivariate K-copula arises from a random matrix model where the non-stationarity of correlations is taken into account by an ensemble of random matrices. It is a symmetric, elliptical copula, which depends on two parameters: the average correlation coefficient, estimated over the considered sample of returns, and a free parameter which characterizes the fluctuations around the average correlation in this sample. We estimated the empirical pairwise copulas for both original and locally normalized returns. The local normalization removes the time-varying trends and volatilities while preserving the dependence structure. The corresponding copula describes the dependence structure on a local scale, whereas the copula of the original returns provides information about the dependence structure on a global scale, i.e., for the full time horizon. Overall, the K-copula captures the empirical dependence structure of market states. We found a good agreement, in particular for the copulas estimated on a local scale. Thus, we obtain a consistent description within the random matrix model: The K-distribution describes the heavytailed return distributions while the K-copula captures the corresponding dependence structure. However, we found an asymmetry in the positive tail dependence, i.e.,

17..4.. Α i...4..6 1 3 4 5 6.4 1 3 4 5 6..4.. Α i...4..6 1 3 4 5 6 market state.4 1 3 4 5 6 market state Figure 8: The mean asymmetry values for each state. Left: asymmetry for the positive tail dependence ᾱ (i), right: asymmetry for the negative tail dependence: β(i). Top: for returns, bottom: for locally normalized returns. a stronger lower tail dependence, indicating a larger probability for simultaneous extreme negative returns. This asymmetry cannot be captured by our model. It is more pronounced on a global scale. On a local scale we find a much weaker asymmetry. However, in times of crisis the asymmetry is still clearly present. References [1] A. Sklar. Fonctions de répartition à n dimensions et leurs marges. Publ. Inst. Statist. Univ. Paris, 8:9 31, 1959. [] A. Sklar. Random variables, joint distribution functions, and copulas. Kybernetika, 9:449 46, 1973. [3] P. Embrechts, A. McNeil, and D. Straumann. Correlation and dependency in risk management: Properties and pitfalls. In Risk Management: Value at Risk and Beyond, Cambridge,. Cambridge University Press. [4] P. Embrechts, F. Lindskog, and A. McNeil. Modelling dependence with copulas and application to risk management. In Handbook of Heavy Tailed Distributions in Finance, Amsterdam, 3. Elsevier. [5] A. McNeil, R. Frey, and P. Embrechts. Quantitative Risk Management: Concepts, Techniques and Tools. Princeton University Press, Princeton, NJ, 5. [6] J. V. Rosenberg and T. Schuermann. A general approach to integrated risk management with skewed fat-tailed risks. J. Finan. Econ., 79:569 614, 6. [7] E. Kole, K. Koedijk, and M. Verbeek. Selecting copulas for risk management. J. Bank. Finance, 31:45 43, 7. [8] D. Brigo, A. Pallavicini, and R. Torresetti. Credit Models and the Crisis: A Journey into CDOs, Copulas, Correlations and Dynamic Models. Wiley Finance, Chichester, 1. [9] A. Meucci. A new breed of copulas for risk and portfolio management. Risk, 4:1 16, 11.

[1] J. V. Rosenberg. Nonparametric pricing of multivariate contingent claims. J. Derivatives, 1:9 6, 3. [11] M. N. Bennett and J. E. Kennedy. Quanto pricing with copulas. J. Derivatives, 1:6 45, 4. [1] U. Cherubini, E. Luciano, and W. Vecchiato. Copula Methods in Finance. Wiley, 4. [13] R. W. J. van den Goorbergh, C. Genest, and B. J. M. Werker. Multivariate option pricing using dynamic copula models. Insurance Math. Econom., 37:11 114, 5. [14] J. Hull and A. White. Valuing credit derivatives using an implied copula approach. J. Derivatives, 14:8 8, 6. [15] M. Hofert and M. Scherer. CDO pricing with nested Archimedean copulas. Quant. Finance, 11:775 787, 11. [16] D. A. Hennessy and H. E. Lapan. The use of Archimedean copulas to model portfolio allocations. Math. Finance, 1:143 154,. [17] A. J. Patton. On the out-of-sample importance of skewness and asymmetric dependence for asset allocation. J. Finan. Econ., :13 168, 4. [18] A. Di Clemente and C. Romano. Measuring and optimizing portfolio credit risk: a copula-based approach. Econ. Notes, 33:35 367, 4. [19] H. Boubaker and N. Sghaier. Portfolio optimization in the presence of dependent financial returns with long memory: A copula based approach. J. Bank. Finance, 37:361 377, 13. [] C. Genest, M. Gendron, and M. Bourdeau-Brien. The advent of copulas in finance. Europ. J. Finance, 15:69 618, 9. [1] A. J. Patton. Handbook of Economic Forecasting. In Copula-Based Models for Financial Time Series, Berlin, 1. Springer. [] T. Mikosch. Copulas: tales and facts. Extremes, 9:3, 6. [3] C. Genest, and Rémillard. Discussion of Copulas: tales and facts. Extremes, 9:7 36, 6. [4] H. Joe. Discussion of Copulas: tales and facts. Extremes, 9:37 41, 6. [5] D. Chetalova, R. Schäfer, and T. Guhr. Zooming into market states. J. Stat. Mech., P19, 15. [6] M. Wollschläger and R. Schäfer. Impact of non-stationarity on estimating and modeling empirical copulas of daily stock returns. arxiv:156.854, 15. [7] F. Black. Studies of stock price volatility changes. In Meetings of the American Statistical Association, Business and Economical Statistics Section, Washington DC, 1976. American Statistical Association. [8] A. A. Christie. The stochastic behavior of common stock variances : Value, leverage and interest rate effects. J. Finan. Econ., 1(4):47 43, 198. [9] G. W. Schwert. Why does stock market volatility change over time? J. Finance, 44:1115 1153, 1989. [3] R. Schäfer and T. Guhr. Local normalization: Uncovering correlations in non-stationary financial time series. Physica A, 389:3856, 1. [31] T. A. Schmitt, D. Chetalova, R. Schäfer, and T. Guhr. Non-stationarity in financial time series: Generic features and tail behavior. Europhys. Lett., 13:583, 13. [3] D. Chetalova, T. A. Schmitt, R. Schäfer, and T. Guhr. Portfolio return distributions: Sample statistics with stochastic correlations. Int. J. Theor. Appl. Finance, 18:1551, 15. [33] F. M. Longin and B. Solnik. Extreme correlation of international equity markets. J. Finance, 56:649 676, 1. [34] A. Ang and J. Chen. Asymmetric correlations of equity portfolios. J. Finan. Econ., 63:443 494,. [35] Y. Hong, J. Tu, and G. Zhou. Asymmetries in stock returns: statistical tests and economic evaluation. Rev. Finan. Stud., :1547 1581, 7. [36] H. Joe. Multivariate Models and Multivariate Dependence Concepts. Chapman & Hall, London, 1997. [37] R. B. Nelsen. An Introduction to Copulas. Springer, New York, 6. [38] D. J. Fenn, M. A. Porter, S. Williams, M. McDonald, N. F. Johnson, and N. S. Jones. Temporal evolution of financial-market correlations. Phys. Rev. E, 84:619, 11. [39] M. C. Münnix, T. Shimada, R. Schäfer, F. Leyvraz, T. H. Seligman, T. Guhr, and H. E. Stanley. Identifying states of a financial market. Sci. Rep., 644, 1. [4] NASDAQ Composite data obtained from http://finance.yahoo.com/. [41] L. Kaufman and P. J. Rousseeuw. Finding Groups in Data: An Introduction to Cluster Analysis. Wiley, New York, 199. [4] R. Tibshirani, G. Walther, and T. Hastie. Estimating the number of clusters in a data set via the gap statistic. J. R. Statist. Soc. B, 63:411 43, 1. [43] S. Demarta and A. McNeil. The t copula and related copulas. Int. Stat. Rev., 73:111 19, 5. 18

[44] W. Sun, S. V. Stoyanov, S. Rachev, and F. J. Fabozzi. Multivariate skewed Students s t copula in the analysis of nonlinear and asymmetric dependence in the German equity market. Stud. Nonlinear Dynam. Econometrics, 1:157 157, 8. [45] M. Ammann and S. Süss. Asymmetric dependence patterns in financial time series. Europ. J. Finance, 15:73 719, 9. 19