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.

Causal analysis of fMRIIn fact, Granger dist<strong>in</strong>guished true causal relations – only to be <strong>in</strong>ferred <strong>in</strong> the presenceof knowledge of the state of the whole universe – from “prima facie" causal relationsthat we refer to as “<strong>in</strong>fluence" <strong>in</strong> agreement with other authors (Commenges andGegout-Petit, 2009). Almost simultaneous with Grangers work, Akaike (Akaike, 1968),and Schweder (Schweder, 1970) <strong>in</strong>troduced similar concepts of <strong>in</strong>fluence, prompt<strong>in</strong>g(Valdes-Sosa et al., <strong>in</strong> press) to co<strong>in</strong> the term WAGS <strong>in</strong>fluence (for Wiener-Akaike-Granger-Schweder). This is a generalization of a proposal by placeAalen (Aalen, 1987;Aalen and Frigessi, 2007) who was among the first to po<strong>in</strong>t out the connections betweenGranger’s and Schweder’s <strong>in</strong>fluence concepts. With<strong>in</strong> this framework we candef<strong>in</strong>e several general types of WAGS <strong>in</strong>fluence, which are applicable to both Markovianand non-Markovian processes, <strong>in</strong> discrete or cont<strong>in</strong>uous time.For three vector time series X 1 (t), X 2 (t), X 3 (t) we wish to know if time series X 1 (t)is <strong>in</strong>fluenced by time series X 2 (t) conditional on X 3 (t). Here X 3 (t) can be consideredany set of relevant time series to be controlled for. Let X [a,b] = {X (t),t ∈ [a,b]} denotethe history of a time series <strong>in</strong> the discrete or cont<strong>in</strong>uous time <strong>in</strong>terval [a,b] The firstcategorical dist<strong>in</strong>ction is based on what part of the present or future of X 1 (t) can be predictedby the past or present of X 2 (τ 2 ) τ 2 ≤ t. This leads to the follow<strong>in</strong>g classification(Florens, 2003; Florens and Fougere, 1996):1. If X 2 (τ 2 ) : τ 2 < t , can <strong>in</strong>fluence any future value of X 1 (t) it is a global <strong>in</strong>fluence.2. If X 2 (τ 2 ) : τ 2 < t , can <strong>in</strong>fluence X 1 (t) at time t it is a local <strong>in</strong>fluence.3. If X 2 (τ 2 ) : τ 2 = t , can <strong>in</strong>fluence X 1 (t) it is a contemporaneous <strong>in</strong>fluence.A second dist<strong>in</strong>ction is based on predict<strong>in</strong>g the whole probability distribution (strong<strong>in</strong>fluence) or only given moments (weak <strong>in</strong>fluence). S<strong>in</strong>ce the most natural formal def<strong>in</strong>itionis one of <strong>in</strong>dependence, every <strong>in</strong>fluence type amounts to the negation of an<strong>in</strong>dependence statement. The two classifications give rise to six types of <strong>in</strong>dependenceand correspond<strong>in</strong>g <strong>in</strong>fluence as set out <strong>in</strong> Table 1.To illustrate, X 1 (t) is strongly, conditionally, and globally <strong>in</strong>dependent of X 2 (t)given X 3 (t), ifP( X 1 (∞,t]| X 1 (t,−∞], X 2 (t,−∞], X 3 (t,−∞]) = P( X 1 (∞,t]| X 1 (t,−∞], X 3 (t,−∞])That is: the probability distribution of the future values of X 1 does not depend on thepast of X 2 , given that the <strong>in</strong>fluence of the past of both X 1 and X 3 has been taken <strong>in</strong>toaccount. When this condition does not hold we say X 2 (t) strongly, conditionally, andglobally <strong>in</strong>fluences (SCGi) X 1 (t) given X 3 (t). Here we use a convention for <strong>in</strong>tervals[a,b) which <strong>in</strong>dicates that the left endpo<strong>in</strong>t is <strong>in</strong>cluded but not the right and thatb precedes a. Note that the whole future of X t is <strong>in</strong>cluded (hence the term “global").And the whole past of all time series is considered. This means these def<strong>in</strong>itions accommodatenon-Markovian processes (for Markovian processes, we only consider theprevious time po<strong>in</strong>t). Furthermore, these def<strong>in</strong>itions do not depend on an assumptionof l<strong>in</strong>earity or any given functional relationship between time series. Note also that83

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

Saved successfully!

Ooh no, something went wrong!