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.

Popescuy i =K∑︁A k y i−k + Bu + Dw i (7)k=0The difference between this and Equation (6) is the presence of a 0-lag matrix A 0which, for easy tractability has zero diagonal entries and is sometimes present on theLHS. This 0-lag matrix is meant to model the sub-sampl<strong>in</strong>g <strong>in</strong>terval dynamic <strong>in</strong>teractionsamong observations, which appear <strong>in</strong>stantaneous, see Moneta et al. (2011) <strong>in</strong> thisvolume. Let us call this form zero lag SVAR. In electric- and magneto- encephalography(EEG/MEG) we often encounter the follow<strong>in</strong>g form:x i =K∑︁µA k x i−k + µ Bu + Dw i ,k=1y i = Cx i (8)Where C represents the observation matrix, or mix<strong>in</strong>g matrix and is determ<strong>in</strong>edby the conductivity/permeability of tissue, and accounts for the superposition of theelectromagnetic fields created by neural activity, which happens at nearly the speed oflight and therefore appears <strong>in</strong>stantaneous. Let us call this mixed output SVAR. F<strong>in</strong>ally,<strong>in</strong> certa<strong>in</strong> eng<strong>in</strong>eer<strong>in</strong>g applications we may see structured disturbances:y i =K∑︁θA k y i−k + θ Bu + D w w i (9)k=1Which we shall call covariate <strong>in</strong>novations SVAR (D w is a general nons<strong>in</strong>gular matrixunlike D which is diagonal). Another f<strong>in</strong>al SVAR form to consider would be one<strong>in</strong> which the 0-lag matrix ⊳A 0 is strictly upper triangular (upper triangular zero lagSVAR):y i = ⊳A 0 y i +K∑︁A k y i−k + ⊳Bu + Dw i (10)k=1F<strong>in</strong>ally, we may consider a upper or lower triangular co-variate <strong>in</strong>novations SVAR:y i =K∑︁A k y i−k + Bu + ⊳Dw i (11)k=0Where ⊳D is upper/lower triangular. The SVAR forms (6)-(10) may look different,and <strong>in</strong> fact each of them may uniquely represent physical processes and allow for direct<strong>in</strong>terpretation of parameters. From a statistical po<strong>in</strong>t of view, however, all four SVARDGPs <strong>in</strong>troduced above are equivalent s<strong>in</strong>ce they have identical cover.Lemma 3 The Gaussian covariate <strong>in</strong>novations SVAR DGP has the same cover as theGaussian mixed output SVAR DGP. Each of these sets has a redundancy of 2 N N! for44

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

Saved successfully!

Ooh no, something went wrong!