Eigenvalue spectra of time-lagged covariance matrices: Possibilities for arbitrage?
|
|
- Cassandra Morton
- 6 years ago
- Views:
Transcription
1 Eigenvalue spectra of time-lagged covariance matrices: Possibilities for arbitrage? Stefan Thurner London July 28
2 Foundations of theory of financial economics CAPM, Markowitz portfolio optimization, etc. Key: correlation matrices of timeseries of financial instruments efficient market hypothesis know: price timeseries are not GBM, BM,... if it were so: would there be hedge funds? London July 28 1
3 Market inefficiency and random matrices eigenvalue spectra of empirical equal-time covariance matrices compared to predictions of EV densities for Gaussian-randomness obtained from random matrix theory (RMT). eigenvalues which strongly depart from RMT spectrum contain information about market sectors largest eigenvalue identified as the market-mode EV cleaning of the original correlation matrices results in improved mean variance efficient frontier RMT provides full understanding why the Markowitz approach is close to useless in actual portfolio management. Reason: dominance of small eigenvalues in the noise regime London July 28 2
4 Eigenvalue spectrum of market data plot from Bochaud/Potters London July 28 3
5 Eigenvalue cleaning + mean variance frontier plot from Bochaud/Potters London July 28 4
6 How can this be used? RMT systematic search for non-random structures in large data sets no time lag: eigenvalue cleaning: Markowitz portfolio arbitrage with no time lag? maybe not arbitrage from time-shifted covariances: maybe yes London July 28 5
7 Ensembles of random matrices random matrix ensemble of N N matrices M with iid random variables distribution P (M) exp ( βn ) 2 Tr(MMT ) β different for different ensembles, i.e. if variables are complex or real Eigenvalue spectra and correlations of eigenvalues are known for: Gaussian orthogonal ensemble (GOE): real valued entries, symmetric random matrices, N Ginibre ensembles (GinOE, GinUE, GinSE): not symmetric (real, complex, quaternion). Eigenvalue densities available even for finite size London July 28 6
8 Random matrices from time-series Input: covariance matrices from N T data matrices X N assets (or instruments) at T observation points Matrix ensemble for N N covariance matrix C XX T is the Wishart ensemble (cornerstone of multivariate data analysis) uncorrelated Gaussian distributed data: exact EV spectrum of XX T the Marcenko-Pastur law is spectrum of time-lagged covariance matrix C τ T t ri tr j t τ, unknown. (results exist for symmetrized lagged correlation matrices problematic) London July 28 7
9 There is structure in lead-lag relations analysis of asymmetric time-lagged correlations forms a fundamental part of finance and econometrics. From practical point: Arbitraging asymmetric lead-lag relationships reported for U.S. stock market (lo & kinlay) lagged correlation function typically exhibits asymmetric peak Why? information adjustment asymmetries: lead-lag effects mainly explained by information adjustment asymmetry (brennan) some recent understanding of relation of strength of lagged correlations and time-shift τ (kertesz) London July 28 8
10 Eigenvalues of time-lagged covariance matrices N T data matrix X: N assets, T observation times, time-lag τ entries: log-return time-series (zero-mean unit-variance) of asset i at time t Time-lagged correlation function C τ (T ) (r i t r i t )(r j t τ r j t τ ) T = 1 T XD τx T r i t = ln S i t ln S i t 1 and D τ δ t,t+τ For τ, C τ is not symmetric Denote eigenvalues of C τ by λ i C τ random variables with a certain distribution,through specific construction, not purely random real asymmetric N N matrix with iid Gaussian entries. Can not expect a flat eigenvalue distribution as in the Ginibre-Girko case. London July 28 9
11 !"#.7,-$!! '+(!"$ Re(!) Im(!).5 "x(!), "y(!) )*'!( (b).6!!!"$ !!"#!!"#.5!!"$! %&'!(!"$ &*'!"#!.2 Q,1!.1.6 Re(!), Im(!).1 Re(!) (d) Im(!) "x(!), "y(!) ()&!' !.5!.5 $%&!'.5!.5 2 &e' Re(!), Im(!).15 *+1.5 Re(!) (f) Im(!) "x(!), "y(!) (m&!' 1.1.5!1!2!2 3!1 $e&!' '+( 1!2 2!1 Re(!)/Im(!) (h) Re(!) 2 Im(!).3 "x(!), "y(!) )*'!(!1 2.4 Q-./ !2!3!2 %&'!( 2!3!2!1 1 Re(!), Im(!) 2 3 London July 28 1
12 General arguments Idea: see distribution of EV in complex plane as a electrical charges. In present case EV density ρ(z) = ρ(x, y) is then given by Poisson equation ρ(x, y) = 1 φ(x, y) 4π with φ(x, y) = 1 N ln det ( (1z C T τ )(1z C τ ) ) c... c average over distribution of X, P (X) exp ( N 2 Tr(XXT ) ) Expand argument of determinant H = 1 z + C τ C T τ x(c τ + C T τ ) + iy(c τ C T τ ) any symmetric (anti-symmetric) contribution of Cτ ij real (imaginary) part of z only influences the London July 28 11
13 General arguments: consequences potential is a function of the radius r = x 2 + y 2 only, φ(x, y) = φ(r) in the limit N thus EV density is radial symmetric function ρ(x, y) = ρ(r) 1 2πr S dzρ(z) δ( z r) Checked validity to 4th order computation if ρ(r) is circular symmetric support S of the EV-spectrum is limited to a circle of r max, which can be computed. London July 28 12
14 Idea: determine ρ(r) based on its radial symmetry EV density of symmetric problem is obtained from well-known relation ρ S (x) = n δ(x x n ) = 1 π lim [ Im(G S (x iɛ)) ] ɛ Idea: use inverse Abel-transform to determine radial symmetric density ρ(r) ρ S ( 2x) = 2 x ρ(r)r r2 x 2dr Now reconstruct EV spectrumexactly for N via inverse Abel-transform and thus via the cuts of the Greens function of the symmetric problem ρ(r) = 1 π 2 r d dx lim [ ɛ Im(G S τ ( 2x iɛ)) ] dx x2 r 2 expect discrepancies in finite data due to finite-size effects London July 28 13
15 Application to lagged correlation matrices Greens function G(z) of the symmetrized problem C τ S = 1 2T X(D τ + D τ )X T is given by implicit solution of (burda) 1 Q 3z2 G 4 (z) 2 1 Q 2( 1 Q 1)zG3 (z) 1 Q (z2 ( 1 Q 1)2 )G 2 (z) + 2( 1 Q 1)zG(z) Q = Q T/N information-to-noise ratio remaining integral hard to solve in general but easy to solve numerically For Q = 1 exact solution is possible ρ Q=1 (r) = 1 6 π 5 r 3 " 243rΓ 3 ««5 5 Γ 4 4 Φ , 5 4, 3! 2, λ Γ 1 «Γ 2 4 «!# 7 Φ λ2,,, 4 2 London July 28 14
16 !(r) Im(")# Re(")# Im(")#Re(")!(r) Im(")# Re(")# Im(")#Re(")!(r) Im(")# Re(")# Im(")#Re(") r.1 (a) r.5 (b) r.5 (c) r London July 28 15
17 Summary theory used a well-known analogy to classical electrostatics to show radial symmetry of the potential for lagged correlation matrices introduced a novel method to calculate the exact radial eigenvalue-density via the inverse Abel-transform used existing results for symmetrized lagged correlation matrices as an input to our method to arrive at ρ Q (r) with the knowledge of the pure random spectrum: deviations indicate arbitragable structure in financial data London July 28 16
18 Application to markets financial data: 5 min data of the S&P5, Jan 2 22 Apr 2 24 empirical lagged correlation matrices C τ, τ = 5, 3 mins data cleaning X: N = 4 time-series at T = 4472 observations From X construct two surrogate data sets: Market mode removed data: X re : largest eigenvalue is market-mode. Market return (movement of the index) rt m = N j=1 v 1jrt, j at equal times, τ =. Regress as in CAPM, r t = α + βr t m + ɛ t and Xit res ɛ i t Randomly reshuffled data : X scr is a random permutation of all elements of X. This destroys correlation structures but keeps same distributions as in original data London July 28 17
19 Im(!) 1 (a)!1! Re(!).2 (b) Im(!) Im(!)!.2.2! Re(!) (c) Im(! res )!.2.2!.2! Re(!) (d)!1! Re(! res ) London July 28 18
20 Eigenvalue spectra! x ("),! y (") r#!(r) r theory Re(") Im(")!.4! Re("),Im(") Eigenvalues lying outside the random regime can be confidently associated with specific non-random structures London July 28 19
21 Interpretation of deviating eigenvalues Deviations from the theoretical pure-random prediction indicate correlation structure in data across time Which assets participate in a given eigenvector associated to a deviating eigenvalue? define inverse participation ratio for the eigenvectors u i IPR( u i ) N u ik 4 k=1 u ik entries in i-th EV IPR shows to what extent each of N = 4 assets contribute to EV u i small IPR: assets contribute equally large IPR: a few assets dominate the eigenvector London July 28 2
22 Indication for group structure IPR(u ik ).2.15 (a) IPR(u ik ) ! i ! i res IPR(u ik ) (b) res IPR(u ik ).1.1.2! res i ! res i largest EV relatively small IPR IPRs for residuals X res larger for the deviating eigenvalues group structure in the lagged-correlations London July 28 21
23 Sector organization in time-lagged data RMT with τ = : know eigenvectors u i of large eigenvalues are associated with sector organization of markets label sectors with s and define sk = { 1 if stock k belongs to sector s otherwise To visualize influence of each sector s to a given eigenvector i I si 1 N s N k=1 sk u ik 2 N s number of stocks in sector s London July 28 22
24 Sector GICS No. of Stocks N s Energy 1 22 Materials Industrials 2 44 Consumer Discretionary Consumer Staples 3 35 Healthcare 35 4 Financials 4 71 Information Technology Telecommunication 5 11 Utilities Table 1: Global Industry Classification Standard (GICS) for the 1 main sectors in the S&P5 from London July 28 23
25 contribution of sector 6 x 1! original! 1 market removed contribution of sector.2 res! GICS code GICS code.1! 2 =! 3 * 8 x 1!3! 2 res contribution of sector contribution of sector GICS code GICS code.2! 4 4 x 1!3! 3 res contribution of sector contribution of sector GICS code GICS code 8 x 1!3! 1.2! 4 res =!5 res contribution of sector contribution of sector London July GICS code GICS code
26 Lead-lag networks London July 28 25
27 Some observations Structures in network-view are associated with deviating eigenvalues via decomposition of the lagged correlation matrix in its right eigenvectors, C λi = UΛ i U 1 Λ i = diag(λ i ) diagonal matrix with only one entry at the respective position associated with eigenvalue λ i For all eigenvalues we found no indication for the leading stocks being the ones with the highest market capitalization (as would be implied by work of lo and kinlay) Close correspondence between the C λi and different sectors visible in the network-representation of the data confirms the validity of analysis Large negative eigenvalues are associated with time-lagged anticorrelations between various sectors London July 28 26
28 Time dependence of largest eigenvalues Q=1 Im(! 1 )= Im(! 1 )" Q=1 Im(! res 1 )= Im(! res 1 )" abs(! 1 ) n Q=1.25 abs(! res 1 ) n Q= (a) n 2 (b) n C 1 (T i ) for consecutive, non-overlapping time periods T i show clear deviations from the predicted support down to Q even though noise is drastically increased for low Q non-random structures present at very short time-scales similar for other ranks of EV London July 28 27
29 Predictive power? c(t 1,T 1+d ) c res (T 1,T 1+d ) 1.5 Q=1.25 Q= d d c(t n, T m ) = (Cij τ (T n ) C ij τ (T n ) ij )(C ij τ (T m ) C ij τ (T m ) ij ) ij σ Tn σ Tm correlation of matrix elements of lagged correlation matrices from nonoverlapping T n and T m. Average over all matrix-elements significance band 1/4 extremely significant, increase with Q London July 28 28
30 A strategy or two there are two ways: cleaning of the matrices in analogy to the method known for τ =. cleaned matrix at time T n allows for an strongly improved estimation of future lagged correlation matrices at times T m > T n fix timescales T n compute complex spectrum of data replace all EV within circle with a number rotate back: find significant correlations and identify temporarily stable ones use movement of leading instrument as trading signal at appropriate timescale use pseudo-sector information for non-standard stat arb strategy London July 28 29
31 Conclusion established random matrix theory to lagged cross-correlation matrices (Abel transformation) this opens maybe the most straight forward way for systematic quantitative search for lead-lag structures in a noisy world possibility for predictability of significant parts of lagged-correlation matrices based on measurements of past (non-overlapping) periods identified correlations can be used in a multitude of strategies as a by-product one gets a pseudo-sectorization if the trading universe (clustering in the lead-lag network of residuals) discussed issues of Q dependence (information to noise ratio) London July 28 3
Random Correlation Matrices, Top Eigenvalue with Heavy Tails and Financial Applications
Random Correlation Matrices, Top Eigenvalue with Heavy Tails and Financial Applications J.P Bouchaud with: M. Potters, G. Biroli, L. Laloux, M. A. Miceli http://www.cfm.fr Portfolio theory: Basics Portfolio
More informationEigenvector stability: Random Matrix Theory and Financial Applications
Eigenvector stability: Random Matrix Theory and Financial Applications J.P Bouchaud with: R. Allez, M. Potters http://www.cfm.fr Portfolio theory: Basics Portfolio weights w i, Asset returns X t i If expected/predicted
More informationEmpirical properties of large covariance matrices in finance
Empirical properties of large covariance matrices in finance Ex: RiskMetrics Group, Geneva Since 2010: Swissquote, Gland December 2009 Covariance and large random matrices Many problems in finance require
More informationARANDOM-MATRIX-THEORY-BASEDANALYSISOF STOCKS OF MARKETS FROM DIFFERENT COUNTRIES
Advances in Complex Systems, Vol. 11, No. 5 (2008) 655 668 c World Scientific Publishing Company ARANDOM-MATRIX-THEORY-BASEDANALYSISOF STOCKS OF MARKETS FROM DIFFERENT COUNTRIES RICARDO COELHO,PETERRICHMONDandSTEFANHUTZLER
More informationApplications of Random Matrix Theory to Economics, Finance and Political Science
Outline Applications of Random Matrix Theory to Economics, Finance and Political Science Matthew C. 1 1 Department of Economics, MIT Institute for Quantitative Social Science, Harvard University SEA 06
More informationPrincipal Component Analysis (PCA) Our starting point consists of T observations from N variables, which will be arranged in an T N matrix R,
Principal Component Analysis (PCA) PCA is a widely used statistical tool for dimension reduction. The objective of PCA is to find common factors, the so called principal components, in form of linear combinations
More informationarxiv:cond-mat/ v2 [cond-mat.stat-mech] 3 Feb 2004
arxiv:cond-mat/0312496v2 [cond-mat.stat-mech] 3 Feb 2004 Signal and Noise in Financial Correlation Matrices Zdzis law Burda, Jerzy Jurkiewicz 1 M. Smoluchowski Institute of Physics, Jagiellonian University,
More informationRandom Matrix Theory Lecture 1 Introduction, Ensembles and Basic Laws. Symeon Chatzinotas February 11, 2013 Luxembourg
Random Matrix Theory Lecture 1 Introduction, Ensembles and Basic Laws Symeon Chatzinotas February 11, 2013 Luxembourg Outline 1. Random Matrix Theory 1. Definition 2. Applications 3. Asymptotics 2. Ensembles
More informationMathematics for Economics and Finance
Mathematics for Economics and Finance Michael Harrison and Patrick Waldron B 375482 Routledge Taylor & Francis Croup LONDON AND NEW YORK Contents List of figures ix List of tables xi Foreword xiii Preface
More informationR = µ + Bf Arbitrage Pricing Model, APM
4.2 Arbitrage Pricing Model, APM Empirical evidence indicates that the CAPM beta does not completely explain the cross section of expected asset returns. This suggests that additional factors may be required.
More informationLectures 6 7 : Marchenko-Pastur Law
Fall 2009 MATH 833 Random Matrices B. Valkó Lectures 6 7 : Marchenko-Pastur Law Notes prepared by: A. Ganguly We will now turn our attention to rectangular matrices. Let X = (X 1, X 2,..., X n ) R p n
More informationRegression: Ordinary Least Squares
Regression: Ordinary Least Squares Mark Hendricks Autumn 2017 FINM Intro: Regression Outline Regression OLS Mathematics Linear Projection Hendricks, Autumn 2017 FINM Intro: Regression: Lecture 2/32 Regression
More informationFactor Mimicking Portfolios
Factor Mimicking Portfolios Bernt Arne Ødegaard 25 October 2018 Contents 1 Factor Mimicking Portfolios 1 1.1 Economic Tracking Portfolios.................................................. 1 1.2 Example.............................................................
More informationRoss (1976) introduced the Arbitrage Pricing Theory (APT) as an alternative to the CAPM.
4.2 Arbitrage Pricing Model, APM Empirical evidence indicates that the CAPM beta does not completely explain the cross section of expected asset returns. This suggests that additional factors may be required.
More information1 Intro to RMT (Gene)
M705 Spring 2013 Summary for Week 2 1 Intro to RMT (Gene) (Also see the Anderson - Guionnet - Zeitouni book, pp.6-11(?) ) We start with two independent families of R.V.s, {Z i,j } 1 i
More informationEconophysics III: Financial Correlations and Portfolio Optimization
FAKULTÄT FÜR PHYSIK Econophysics III: Financial Correlations and Portfolio Optimization Thomas Guhr Let s Face Chaos through Nonlinear Dynamics, Maribor 21 Outline Portfolio optimization is a key issue
More informationSignal and noise in financial correlation matrices
Physica A 344 (2004) 67 72 www.elsevier.com/locate/physa Signal and noise in financial correlation matrices Zdzis"aw Burda, Jerzy Jurkiewicz 1 M. Smoluchowski Institute of Physics, Jagiellonian University,
More informationThe Slow Convergence of OLS Estimators of α, β and Portfolio. β and Portfolio Weights under Long Memory Stochastic Volatility
The Slow Convergence of OLS Estimators of α, β and Portfolio Weights under Long Memory Stochastic Volatility New York University Stern School of Business June 21, 2018 Introduction Bivariate long memory
More informationThreading Rotational Dynamics, Crystal Optics, and Quantum Mechanics to Risky Asset Selection. A Physics & Pizza on Wednesdays presentation
Threading Rotational Dynamics, Crystal Optics, and Quantum Mechanics to Risky Asset Selection A Physics & Pizza on Wednesdays presentation M. Hossein Partovi Department of Physics & Astronomy, California
More informationCorrelation Matrices and the Perron-Frobenius Theorem
Wilfrid Laurier University July 14 2014 Acknowledgements Thanks to David Melkuev, Johnew Zhang, Bill Zhang, Ronnnie Feng Jeyan Thangaraj for research assistance. Outline The - theorem extensions to negative
More informationQuantum Chaos and Nonunitary Dynamics
Quantum Chaos and Nonunitary Dynamics Karol Życzkowski in collaboration with W. Bruzda, V. Cappellini, H.-J. Sommers, M. Smaczyński Phys. Lett. A 373, 320 (2009) Institute of Physics, Jagiellonian University,
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 1 Aug 2001
A Random Matrix Approach to Cross-Correlations in Financial Data arxiv:cond-mat/0108023v1 [cond-mat.stat-mech] 1 Aug 2001 Vasiliki Plerou 1,2, Parameswaran Gopikrishnan 1, Bernd Rosenow 1,3, Luís A. Nunes
More informationUniversity of Cambridge Engineering Part IIB Module 3F3: Signal and Pattern Processing Handout 2:. The Multivariate Gaussian & Decision Boundaries
University of Cambridge Engineering Part IIB Module 3F3: Signal and Pattern Processing Handout :. The Multivariate Gaussian & Decision Boundaries..15.1.5 1 8 6 6 8 1 Mark Gales mjfg@eng.cam.ac.uk Lent
More informationarxiv: v2 [cond-mat.stat-mech] 16 Jul 2012
EIGENVECTOR DYNAMICS: GENERAL THEORY AND SOME APPLICATIONS arxiv:1203.6228v2 [cond-mat.stat-mech] 16 Jul 2012 ROMAIN ALLEZ 1,2 AND JEAN-PHILIPPE BOUCHAUD 1 Abstract. We propose a general framework to study
More informationVast Volatility Matrix Estimation for High Frequency Data
Vast Volatility Matrix Estimation for High Frequency Data Yazhen Wang National Science Foundation Yale Workshop, May 14-17, 2009 Disclaimer: My opinion, not the views of NSF Y. Wang (at NSF) 1 / 36 Outline
More informationPrior Information: Shrinkage and Black-Litterman. Prof. Daniel P. Palomar
Prior Information: Shrinkage and Black-Litterman Prof. Daniel P. Palomar The Hong Kong University of Science and Technology (HKUST) MAFS6010R- Portfolio Optimization with R MSc in Financial Mathematics
More informationGI07/COMPM012: Mathematical Programming and Research Methods (Part 2) 2. Least Squares and Principal Components Analysis. Massimiliano Pontil
GI07/COMPM012: Mathematical Programming and Research Methods (Part 2) 2. Least Squares and Principal Components Analysis Massimiliano Pontil 1 Today s plan SVD and principal component analysis (PCA) Connection
More informationAssessing the dependence of high-dimensional time series via sample autocovariances and correlations
Assessing the dependence of high-dimensional time series via sample autocovariances and correlations Johannes Heiny University of Aarhus Joint work with Thomas Mikosch (Copenhagen), Richard Davis (Columbia),
More informationRobustní monitorování stability v modelu CAPM
Robustní monitorování stability v modelu CAPM Ondřej Chochola, Marie Hušková, Zuzana Prášková (MFF UK) Josef Steinebach (University of Cologne) ROBUST 2012, Němčičky, 10.-14.9. 2012 Contents Introduction
More informationProducts of Rectangular Gaussian Matrices
Products of Rectangular Gaussian Matrices Jesper R. Ipsen Department of Physics Bielefeld University in collaboration with Gernot Akemann Mario Kieburg Summer School on Randomness in Physics & Mathematics
More informationASSET PRICING MODELS
ASSE PRICING MODELS [1] CAPM (1) Some notation: R it = (gross) return on asset i at time t. R mt = (gross) return on the market portfolio at time t. R ft = return on risk-free asset at time t. X it = R
More informationMultivariate Time Series Analysis and Its Applications [Tsay (2005), chapter 8]
1 Multivariate Time Series Analysis and Its Applications [Tsay (2005), chapter 8] Insights: Price movements in one market can spread easily and instantly to another market [economic globalization and internet
More informationCleaning correlation matrices, Random Matrix Theory & HCIZ integrals
Cleaning correlation matrices, Random Matrix Theory & HCIZ integrals J.P Bouchaud with: M. Potters, L. Laloux, R. Allez, J. Bun, S. Majumdar http://www.cfm.fr Portfolio theory: Basics Portfolio weights
More informationFactor Models for Asset Returns. Prof. Daniel P. Palomar
Factor Models for Asset Returns Prof. Daniel P. Palomar The Hong Kong University of Science and Technology (HKUST) MAFS6010R- Portfolio Optimization with R MSc in Financial Mathematics Fall 2018-19, HKUST,
More informationNotes on Random Variables, Expectations, Probability Densities, and Martingales
Eco 315.2 Spring 2006 C.Sims Notes on Random Variables, Expectations, Probability Densities, and Martingales Includes Exercise Due Tuesday, April 4. For many or most of you, parts of these notes will be
More informationUniversality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium
Universality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium SEA 06@MIT, Workshop on Stochastic Eigen-Analysis and its Applications, MIT, Cambridge,
More informationIntro VEC and BEKK Example Factor Models Cond Var and Cor Application Ref 4. MGARCH
ntro VEC and BEKK Example Factor Models Cond Var and Cor Application Ref 4. MGARCH JEM 140: Quantitative Multivariate Finance ES, Charles University, Prague Summer 2018 JEM 140 () 4. MGARCH Summer 2018
More informationIdentifying Financial Risk Factors
Identifying Financial Risk Factors with a Low-Rank Sparse Decomposition Lisa Goldberg Alex Shkolnik Berkeley Columbia Meeting in Engineering and Statistics 24 March 2016 Outline 1 A Brief History of Factor
More information8.1 Concentration inequality for Gaussian random matrix (cont d)
MGMT 69: Topics in High-dimensional Data Analysis Falll 26 Lecture 8: Spectral clustering and Laplacian matrices Lecturer: Jiaming Xu Scribe: Hyun-Ju Oh and Taotao He, October 4, 26 Outline Concentration
More informationarxiv:physics/ v1 [physics.soc-ph] 15 May 2006
arxiv:physics/0605115v1 [physics.soc-ph] 15 May 2006 Asymmetric matrices in an analysis of financial correlations J. Kwapień 1, S. Drożdż 1,2, A.Z. Górski 1, P. Oświȩcimka 1 1 Institute of Nuclear Physics,
More informationMarkowitz Efficient Portfolio Frontier as Least-Norm Analytic Solution to Underdetermined Equations
Markowitz Efficient Portfolio Frontier as Least-Norm Analytic Solution to Underdetermined Equations Sahand Rabbani Introduction Modern portfolio theory deals in part with the efficient allocation of investments
More informationarxiv:physics/ v3 [physics.data-an] 29 Nov 2006
arxiv:physics/060604v3 [physics.data-an] 29 Nov 2006 Correlation matrix decomposition of WIG20 intraday fluctuations R.Rak a S. Drożdż ab J. Kwapień b P. Oświȩcimka b a Institute of Physics University
More informationNumerical Solutions to the General Marcenko Pastur Equation
Numerical Solutions to the General Marcenko Pastur Equation Miriam Huntley May 2, 23 Motivation: Data Analysis A frequent goal in real world data analysis is estimation of the covariance matrix of some
More informationTime Series Modeling of Financial Data. Prof. Daniel P. Palomar
Time Series Modeling of Financial Data Prof. Daniel P. Palomar The Hong Kong University of Science and Technology (HKUST) MAFS6010R- Portfolio Optimization with R MSc in Financial Mathematics Fall 2018-19,
More informationEstimation of the Global Minimum Variance Portfolio in High Dimensions
Estimation of the Global Minimum Variance Portfolio in High Dimensions Taras Bodnar, Nestor Parolya and Wolfgang Schmid 07.FEBRUARY 2014 1 / 25 Outline Introduction Random Matrix Theory: Preliminary Results
More informationSignal and Noise in Correlation Matrix
arxiv:cond-mat/0305627v2 [cond-mat.stat-mech] 3 Feb 2004 Signal and Noise in Correlation Matrix Z. Burda, A. Görlich, A. Jarosz and J. Jurkiewicz M. Smoluchowski Institute of Physics, Jagellonian University,
More informationSingular Value Decomposition and Principal Component Analysis (PCA) I
Singular Value Decomposition and Principal Component Analysis (PCA) I Prof Ned Wingreen MOL 40/50 Microarray review Data per array: 0000 genes, I (green) i,i (red) i 000 000+ data points! The expression
More informationConcentration Inequalities for Random Matrices
Concentration Inequalities for Random Matrices M. Ledoux Institut de Mathématiques de Toulouse, France exponential tail inequalities classical theme in probability and statistics quantify the asymptotic
More informationRandom Matrix: From Wigner to Quantum Chaos
Random Matrix: From Wigner to Quantum Chaos Horng-Tzer Yau Harvard University Joint work with P. Bourgade, L. Erdős, B. Schlein and J. Yin 1 Perhaps I am now too courageous when I try to guess the distribution
More informationEigenvalue PDFs. Peter Forrester, M&S, University of Melbourne
Outline Eigenvalue PDFs Peter Forrester, M&S, University of Melbourne Hermitian matrices with real, complex or real quaternion elements Circular ensembles and classical groups Products of random matrices
More informationIntroduction to Machine Learning
10-701 Introduction to Machine Learning PCA Slides based on 18-661 Fall 2018 PCA Raw data can be Complex, High-dimensional To understand a phenomenon we measure various related quantities If we knew what
More informationTHE UNIVERSITY OF CHICAGO Graduate School of Business Business 41912, Spring Quarter 2008, Mr. Ruey S. Tsay. Solutions to Final Exam
THE UNIVERSITY OF CHICAGO Graduate School of Business Business 41912, Spring Quarter 2008, Mr. Ruey S. Tsay Solutions to Final Exam 1. (13 pts) Consider the monthly log returns, in percentages, of five
More informationDepartment of Economics, UCSB UC Santa Barbara
Department of Economics, UCSB UC Santa Barbara Title: Past trend versus future expectation: test of exchange rate volatility Author: Sengupta, Jati K., University of California, Santa Barbara Sfeir, Raymond,
More informationEigenvalues and Singular Values of Random Matrices: A Tutorial Introduction
Random Matrix Theory and its applications to Statistics and Wireless Communications Eigenvalues and Singular Values of Random Matrices: A Tutorial Introduction Sergio Verdú Princeton University National
More informationErrata for Campbell, Financial Decisions and Markets, 01/02/2019.
Errata for Campbell, Financial Decisions and Markets, 01/02/2019. Page xi, section title for Section 11.4.3 should be Endogenous Margin Requirements. Page 20, equation 1.49), expectations operator E should
More informationStatistics 351 Probability I Fall 2006 (200630) Final Exam Solutions. θ α β Γ(α)Γ(β) (uv)α 1 (v uv) β 1 exp v }
Statistics 35 Probability I Fall 6 (63 Final Exam Solutions Instructor: Michael Kozdron (a Solving for X and Y gives X UV and Y V UV, so that the Jacobian of this transformation is x x u v J y y v u v
More informationA note on a Marčenko-Pastur type theorem for time series. Jianfeng. Yao
A note on a Marčenko-Pastur type theorem for time series Jianfeng Yao Workshop on High-dimensional statistics The University of Hong Kong, October 2011 Overview 1 High-dimensional data and the sample covariance
More informationMultivariate Statistical Analysis
Multivariate Statistical Analysis Fall 2011 C. L. Williams, Ph.D. Lecture 4 for Applied Multivariate Analysis Outline 1 Eigen values and eigen vectors Characteristic equation Some properties of eigendecompositions
More informationRegularization of Portfolio Allocation
Regularization of Portfolio Allocation Benjamin Bruder Quantitative Research Lyxor Asset Management, Paris benjamin.bruder@lyxor.com Jean-Charles Richard Quantitative Research Lyxor Asset Management, Paris
More informationCointegrated VAR s. Eduardo Rossi University of Pavia. November Rossi Cointegrated VAR s Financial Econometrics / 56
Cointegrated VAR s Eduardo Rossi University of Pavia November 2013 Rossi Cointegrated VAR s Financial Econometrics - 2013 1 / 56 VAR y t = (y 1t,..., y nt ) is (n 1) vector. y t VAR(p): Φ(L)y t = ɛ t The
More informationThe problem is to infer on the underlying probability distribution that gives rise to the data S.
Basic Problem of Statistical Inference Assume that we have a set of observations S = { x 1, x 2,..., x N }, xj R n. The problem is to infer on the underlying probability distribution that gives rise to
More informationNetwork Connectivity and Systematic Risk
Network Connectivity and Systematic Risk Monica Billio 1 Massimiliano Caporin 2 Roberto Panzica 3 Loriana Pelizzon 1,3 1 University Ca Foscari Venezia (Italy) 2 University of Padova (Italy) 3 Goethe University
More informationEcon671 Factor Models: Principal Components
Econ671 Factor Models: Principal Components Jun YU April 8, 2016 Jun YU () Econ671 Factor Models: Principal Components April 8, 2016 1 / 59 Factor Models: Principal Components Learning Objectives 1. Show
More informationFinancial Econometrics Return Predictability
Financial Econometrics Return Predictability Eric Zivot March 30, 2011 Lecture Outline Market Efficiency The Forms of the Random Walk Hypothesis Testing the Random Walk Hypothesis Reading FMUND, chapter
More informationComparison Method in Random Matrix Theory
Comparison Method in Random Matrix Theory Jun Yin UW-Madison Valparaíso, Chile, July - 2015 Joint work with A. Knowles. 1 Some random matrices Wigner Matrix: H is N N square matrix, H : H ij = H ji, EH
More informationFinancial Econometrics
Financial Econometrics Nonlinear time series analysis Gerald P. Dwyer Trinity College, Dublin January 2016 Outline 1 Nonlinearity Does nonlinearity matter? Nonlinear models Tests for nonlinearity Forecasting
More informationMultivariate Time Series: VAR(p) Processes and Models
Multivariate Time Series: VAR(p) Processes and Models A VAR(p) model, for p > 0 is X t = φ 0 + Φ 1 X t 1 + + Φ p X t p + A t, where X t, φ 0, and X t i are k-vectors, Φ 1,..., Φ p are k k matrices, with
More informationFinancial Econometrics Short Course Lecture 3 Multifactor Pricing Model
Financial Econometrics Short Course Lecture 3 Multifactor Pricing Model Oliver Linton obl20@cam.ac.uk Renmin University Financial Econometrics Short Course Lecture 3 MultifactorRenmin Pricing University
More informationSTATE COUNCIL OF EDUCATIONAL RESEARCH AND TRAINING TNCF DRAFT SYLLABUS
STATE COUNCIL OF EDUCATIONAL RESEARCH AND TRAINING TNCF 2017 - DRAFT SYLLABUS Subject :Business Maths Class : XI Unit 1 : TOPIC Matrices and Determinants CONTENT Determinants - Minors; Cofactors; Evaluation
More informationRandom Matrices and Multivariate Statistical Analysis
Random Matrices and Multivariate Statistical Analysis Iain Johnstone, Statistics, Stanford imj@stanford.edu SEA 06@MIT p.1 Agenda Classical multivariate techniques Principal Component Analysis Canonical
More informationOn the distinguishability of random quantum states
1 1 Department of Computer Science University of Bristol Bristol, UK quant-ph/0607011 Distinguishing quantum states Distinguishing quantum states This talk Question Consider a known ensemble E of n quantum
More informationRobust Portfolio Risk Minimization Using the Graphical Lasso
Robust Portfolio Risk Minimization Using the Graphical Lasso Tristan Millington & Mahesan Niranjan Department of Electronics and Computer Science University of Southampton Highfield SO17 1BJ, Southampton,
More informationLarge sample covariance matrices and the T 2 statistic
Large sample covariance matrices and the T 2 statistic EURANDOM, the Netherlands Joint work with W. Zhou Outline 1 2 Basic setting Let {X ij }, i, j =, be i.i.d. r.v. Write n s j = (X 1j,, X pj ) T and
More informationMotivation Non-linear Rational Expectations The Permanent Income Hypothesis The Log of Gravity Non-linear IV Estimation Summary.
Econometrics I Department of Economics Universidad Carlos III de Madrid Master in Industrial Economics and Markets Outline Motivation 1 Motivation 2 3 4 5 Motivation Hansen's contributions GMM was developed
More informationPeter S. Karlsson Jönköping University SE Jönköping Sweden
Peter S. Karlsson Jönköping University SE-55 Jönköping Sweden Abstract: The Arbitrage Pricing Theory provides a theory to quantify risk and the reward for taking it. While the theory itself is sound from
More informationThe Multivariate Gaussian Distribution
The Multivariate Gaussian Distribution Chuong B. Do October, 8 A vector-valued random variable X = T X X n is said to have a multivariate normal or Gaussian) distribution with mean µ R n and covariance
More informationFinancial Econometrics
Financial Econometrics Multivariate Time Series Analysis: VAR Gerald P. Dwyer Trinity College, Dublin January 2013 GPD (TCD) VAR 01/13 1 / 25 Structural equations Suppose have simultaneous system for supply
More informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationUnsupervised Learning: Dimensionality Reduction
Unsupervised Learning: Dimensionality Reduction CMPSCI 689 Fall 2015 Sridhar Mahadevan Lecture 3 Outline In this lecture, we set about to solve the problem posed in the previous lecture Given a dataset,
More informationFinancial Econometrics Lecture 6: Testing the CAPM model
Financial Econometrics Lecture 6: Testing the CAPM model Richard G. Pierse 1 Introduction The capital asset pricing model has some strong implications which are testable. The restrictions that can be tested
More informationMultivariate Random Variable
Multivariate Random Variable Author: Author: Andrés Hincapié and Linyi Cao This Version: August 7, 2016 Multivariate Random Variable 3 Now we consider models with more than one r.v. These are called multivariate
More informationProbabilities & Statistics Revision
Probabilities & Statistics Revision Christopher Ting Christopher Ting http://www.mysmu.edu/faculty/christophert/ : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 January 6, 2017 Christopher Ting QF
More informationMA Advanced Macroeconomics: Solving Models with Rational Expectations
MA Advanced Macroeconomics: Solving Models with Rational Expectations Karl Whelan School of Economics, UCD February 6, 2009 Karl Whelan (UCD) Models with Rational Expectations February 6, 2009 1 / 32 Moving
More informationECON4515 Finance theory 1 Diderik Lund, 5 May Perold: The CAPM
Perold: The CAPM Perold starts with a historical background, the development of portfolio theory and the CAPM. Points out that until 1950 there was no theory to describe the equilibrium determination of
More informationThe Matrix Dyson Equation in random matrix theory
The Matrix Dyson Equation in random matrix theory László Erdős IST, Austria Mathematical Physics seminar University of Bristol, Feb 3, 207 Joint work with O. Ajanki, T. Krüger Partially supported by ERC
More informationOptimal Investment Strategies: A Constrained Optimization Approach
Optimal Investment Strategies: A Constrained Optimization Approach Janet L Waldrop Mississippi State University jlc3@ramsstateedu Faculty Advisor: Michael Pearson Pearson@mathmsstateedu Contents Introduction
More informationOperator norm convergence for sequence of matrices and application to QIT
Operator norm convergence for sequence of matrices and application to QIT Benoît Collins University of Ottawa & AIMR, Tohoku University Cambridge, INI, October 15, 2013 Overview Overview Plan: 1. Norm
More informationVector Space Models. wine_spectral.r
Vector Space Models 137 wine_spectral.r Latent Semantic Analysis Problem with words Even a small vocabulary as in wine example is challenging LSA Reduce number of columns of DTM by principal components
More information1 Tridiagonal matrices
Lecture Notes: β-ensembles Bálint Virág Notes with Diane Holcomb 1 Tridiagonal matrices Definition 1. Suppose you have a symmetric matrix A, we can define its spectral measure (at the first coordinate
More informationEstimation of large dimensional sparse covariance matrices
Estimation of large dimensional sparse covariance matrices Department of Statistics UC, Berkeley May 5, 2009 Sample covariance matrix and its eigenvalues Data: n p matrix X n (independent identically distributed)
More informationVAR Model. (k-variate) VAR(p) model (in the Reduced Form): Y t-2. Y t-1 = A + B 1. Y t + B 2. Y t-p. + ε t. + + B p. where:
VAR Model (k-variate VAR(p model (in the Reduced Form: where: Y t = A + B 1 Y t-1 + B 2 Y t-2 + + B p Y t-p + ε t Y t = (y 1t, y 2t,, y kt : a (k x 1 vector of time series variables A: a (k x 1 vector
More informationNon white sample covariance matrices.
Non white sample covariance matrices. S. Péché, Université Grenoble 1, joint work with O. Ledoit, Uni. Zurich 17-21/05/2010, Université Marne la Vallée Workshop Probability and Geometry in High Dimensions
More informationThe Hilbert Space of Random Variables
The Hilbert Space of Random Variables Electrical Engineering 126 (UC Berkeley) Spring 2018 1 Outline Fix a probability space and consider the set H := {X : X is a real-valued random variable with E[X 2
More information3.1. The probabilistic view of the principal component analysis.
301 Chapter 3 Principal Components and Statistical Factor Models This chapter of introduces the principal component analysis (PCA), briefly reviews statistical factor models PCA is among the most popular
More informationSparse PCA with applications in finance
Sparse PCA with applications in finance A. d Aspremont, L. El Ghaoui, M. Jordan, G. Lanckriet ORFE, Princeton University & EECS, U.C. Berkeley Available online at www.princeton.edu/~aspremon 1 Introduction
More informationConvergence of spectral measures and eigenvalue rigidity
Convergence of spectral measures and eigenvalue rigidity Elizabeth Meckes Case Western Reserve University ICERM, March 1, 2018 Macroscopic scale: the empirical spectral measure Macroscopic scale: the empirical
More information1: PROBABILITY REVIEW
1: PROBABILITY REVIEW Marek Rutkowski School of Mathematics and Statistics University of Sydney Semester 2, 2016 M. Rutkowski (USydney) Slides 1: Probability Review 1 / 56 Outline We will review the following
More informationStatistical Machine Learning
Statistical Machine Learning Christoph Lampert Spring Semester 2015/2016 // Lecture 12 1 / 36 Unsupervised Learning Dimensionality Reduction 2 / 36 Dimensionality Reduction Given: data X = {x 1,..., x
More informationRandom Matrix Theory and its Applications to Econometrics
Random Matrix Theory and its Applications to Econometrics Hyungsik Roger Moon University of Southern California Conference to Celebrate Peter Phillips 40 Years at Yale, October 2018 Spectral Analysis of
More information