11.07.2015 Views

Causality in Time Series - ClopiNet

Causality in Time Series - ClopiNet

Causality in Time Series - ClopiNet

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Robust Statistics for <strong>Causality</strong><strong>in</strong>stances <strong>in</strong> which the matrices D w is the product of and unitary and diagonal matrices,the matrix C is a unitary matrix and the matrix A 0 is a permutation of an uppertriangular matrix.Proof Star<strong>in</strong>g with the def<strong>in</strong>ition of covariate <strong>in</strong>novations SVAR <strong>in</strong> Equation (9) weuse the variable transformation y = D w x and obta<strong>in</strong> the mixed-output form (trivial).The set of Guassian random variables is closed under scalar multiplication (and hencesign change) and addition. This means that the variance if the <strong>in</strong>novations term <strong>in</strong>Equation (9) can be written as:Σ w = D T wD w = D T wU T UD wWhere U is a unitary (orthogonal, unit 2-norm) matrix. S<strong>in</strong>ce all <strong>in</strong>novations termelements are zero mean, the covariance matrix is the sole descriptor of the Gaussian<strong>in</strong>novations term. This <strong>in</strong> turn means that any other matrix D ′ w = D T wU T substituted <strong>in</strong>tothe DGP described <strong>in</strong> Equation (9) amounts to a stochastically equivalent DGP. Thematrix D ′ w can belong to a number of general sets of matrices, one of which is the set ofnons<strong>in</strong>gular upper triangular matrices (the transformation is achievable through the QRdecomposition of Σ w ). Another such set is lower triangular matrix set. Both are subsetsof the set of matrices sometimes named ‘psychologically upper triangular’, mean<strong>in</strong>g apermutation of an upper triangular matrix.If we constra<strong>in</strong> D w to be of the form D w = UD, i.e. such that (by polar decomposition)it is the product of a unitary and a diagonal positive def<strong>in</strong>ite matrix, the onlystochastically equivalent transformations of D w are a symmetry preserv<strong>in</strong>g permutationof its rows/columns and a sign change <strong>in</strong> one of the columns (this is a property oforthogonal matrices such as U). There are N! such permutations and 2 N possible signchanges. For the general case, <strong>in</strong> which the <strong>in</strong>put u has no special properties, there areno other redundancies <strong>in</strong> the SVAR model (s<strong>in</strong>ce chang<strong>in</strong>g any parameter <strong>in</strong> A and Bwill otherwise change the output). Without loss of generality then, we can write thetransformation from covariate <strong>in</strong>novations to mixed output SVAR form as:x i =y i =K∑︁θA k y i−k + θ Bu + UD w w ik=1K∑︁U T ( θ A k )Ux i−k + U T ( θ B)u + D w w ik=1y i = U T x iS<strong>in</strong>ce the transformation U is one to one and <strong>in</strong>vertible, and s<strong>in</strong>ce this transformationis what allows a (restricted) a covariate noise SVAR to map, one to one, onto amixed output SVAR, the card<strong>in</strong>ality of both covers is the same.Now consider the zero-lag SVAR form:45

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!