Next Article in Journal
Numerical Study of Entropy Generation in a Flowing Nanofluid Used in Micro- and Minichannels
Next Article in Special Issue
Moving Frames of Reference, Relativity and Invariance in Transfer Entropy and Information Dynamics
Previous Article in Journal
Machine Learning with Squared-Loss Mutual Information
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

The Relation between Granger Causality and Directed Information Theory: A Review

by
Pierre-Olivier Amblard
1,2,* and
Olivier J. J. Michel
1
1
GIPSAlab/CNRS UMR 5216/ BP46, 38402 Saint Martin d'Hères cedex, France
2
The University of Melbourne, Department of Mathematics and Statistics, Parkville, VIC, 3010, Australia
*
Author to whom correspondence should be addressed.
Submission received: 14 November 2012 / Revised: 19 December 2012 / Accepted: 19 December 2012 / Published: 28 December 2012
(This article belongs to the Special Issue Transfer Entropy)

Abstract

:
This report reviews the conceptual and theoretical links between Granger causality and directed information theory. We begin with a short historical tour of Granger causality, concentrating on its closeness to information theory. The definitions of Granger causality based on prediction are recalled, and the importance of the observation set is discussed. We present the definitions based on conditional independence. The notion of instantaneous coupling is included in the definitions. The concept of Granger causality graphs is discussed. We present directed information theory from the perspective of studies of causal influences between stochastic processes. Causal conditioning appears to be the cornerstone for the relation between information theory and Granger causality. In the bivariate case, the fundamental measure is the directed information, which decomposes as the sum of the transfer entropies and a term quantifying instantaneous coupling. We show the decomposition of the mutual information into the sums of the transfer entropies and the instantaneous coupling measure, a relation known for the linear Gaussian case. We study the multivariate case, showing that the useful decomposition is blurred by instantaneous coupling. The links are further developed by studying how measures based on directed information theory naturally emerge from Granger causality inference frameworks as hypothesis testing.

1. Introduction

This review deals with the analysis of influences that one system, be it physical, economical, biological or social, for example, can exert over another. In several scientific fields, the finding of the influence network between different systems is crucial. As examples, we can think of gene influence networks [1,2], relations between economical variables [3,4], communication between neurons or the flow of information between different brain regions [5], or the human influence on the Earth climate [6,7], and many others.
The context studied in this report is illustrated in Figure 1. For a given system, we have at disposal a number of different measurements. In neuroscience, these can be local field potentials recorded in the brain of an animal. In solar physics, these can be solar indices measured by sensors onboard some satellite. In the study of turbulent fluids, these can be the velocity measured at different scales in the fluid (or can be as in Figure 1, the wavelet analysis of the velocity at different scales). For these different examples, the aim is to find dependencies between the different measurements, and if possible, to give a direction to the dependence. In neuroscience, this will allow to understand how information flows between different areas of the brain. In solar physics, this will allow to understand the links between indices and their influence on the total solar irradiance received on Earth. In the study of turbulence, this can confirm the directional cascade of energy from large down to small scales.
Figure 1. Illustration of the problem of information flow in networks of stochastic processes. Each node of the network is associated to a signal. Edges between nodes stand for dependence (shared information) between the signals. The dependence can be directed or not. This framework can be applied to different situations such as solar physics, neuroscience or the study of turbulence in fluids, as illustrated by the three examples depicted here.
Figure 1. Illustration of the problem of information flow in networks of stochastic processes. Each node of the network is associated to a signal. Edges between nodes stand for dependence (shared information) between the signals. The dependence can be directed or not. This framework can be applied to different situations such as solar physics, neuroscience or the study of turbulence in fluids, as illustrated by the three examples depicted here.
Entropy 15 00113 g001
In a graphical modeling approach, each signal is associated to a particular node of a graph, and dependencies are represented by edges, directed if a directional dependence exists. The questions addressed in this paper concern the assessment of directional dependence between signals, and thus concern the inference problem of estimating the edge set in the graph of signals considered.
Climatology and neuroscience were already given as examples by Norbert Wiener in 1956 [8], a paper which inspired econometrist Clive Granger to develop what is now termed Granger causality [9]. Wiener proposed in this paper that a signal x causes another time series y, if the past of x has a strictly positive influence on the quality of prediction of y. Let us quote Wiener [8]:
“As an application of this, let us consider the case where f 1 ( α ) represents the temperature at 9 A.M. in Boston and f 2 ( α ) represents the temperature at the same time in Albany. We generally suppose that weather moves from west to east with the rotation of the earth; the two quantities 1 C and its correlate in the other direction will enable us to make a precise statement containing some if this content and then verify whether this statement is true or not. Or again, in the study of brain waves we may be able to obtain electroencephalograms more or less corresponding to electrical activity in different part of the brain. Here the study of coefficients of causality running both ways and of their analogues for sets of more than two functions f may be useful in determining what part of the brain is driving what other part of the brain in its normal activity.”
In a wide sense, Granger causality can be summed up as a theoretical framework based on conditional independence to assess directional dependencies between time series. It is interesting to note that Norbert Wiener influenced Granger causality, as well as another field dedicated to the analysis of dependencies: information theory. Information theory has led to the definition of quantities that measure the uncertainty of variables using probabilistic concepts. Furthermore, this has led to the definition of measures of dependence based on the decrease in uncertainty relating to one variable after observing another one. Usual information theory is, however, symmetrical. For example, the well-known mutual information rate between two stationary time series is symmetrical under an exchange of the two signals: the mutual information assesses the undirectional dependence. Directional dependence analysis viewed as an information-theoretic problem requires the breaking of the usual symmetry of information theory. This was realized in the 1960s and early 1970s by Hans Marko, a German professor of communication. He developed the bidirectional information theory in the Markov case [10]. This theory was later generalized by James Massey and Gerhard Kramer, to what we may now call directed information theory [11,12].
It is the aim of this report to review the conceptual and theoretical links between Granger causality and directed information theory.
Many information-theoretic tools have been designed for the practical implementation of Granger causality ideas. We will not show all of the different measures proposed, because they are almost always particular cases of the measures issued from directed information theory. Furthermore, some measures might have been proposed in different fields (and/or at different periods of time) and have received different names. We will only consider the well-accepted names. This is the case, for example, of “transfer entropy”, as coined by Schreiber in 2000 [13], but which appeared earlier under different names, in different fields, and might be considered under slightly different hypotheses. Prior to developing a unified view of the links between Granger causality and information theory, we will provide a survey of the literature, concentrating on studies where information theory and Granger causality are jointly presented.
Furthermore, we will not review any practical aspects, nor any detailed applications. In this spirit, this report is different from [14], which concentrated on the estimation of information quantities, and where the review is restricted to transfer entropy. For reviews on the analysis of dependencies between systems and for applications of Granger causality in neuroscience, we refer to [15,16]. We will mention however some important practical points in our conclusions, where we will also discuss some current and future directions of research in the field.

1.1. What Is, and What Is Not, Granger Causality

We will not debate the meaning of causality or causation. We instead refer to [17]. However, we must emphasize that Granger causality actually measures a statistical dependence between the past of a process and the present of another. In this respect, the word causality in Granger causality takes on the usual meaning that a cause occurs prior to its effect. However, nothing in the definitions that we will recall precludes that signal x can simultaneously be Granger caused by y and be a cause of y! This lies in the very close connection between Granger causality and the feedback between times series.
Granger causality is based on the usual concept of conditioning in probability theory, whereas approaches developed for example in [17,18] relied on causal calculus and the concept of intervention. In this spirit, intervention is closer to experimental sciences, where we imagine that we can really, for example, freeze some system and measure the influence of this action on another process. It is now well-known that causality in the sense of between random variables can be inferred unambiguously only in restricted cases, such as directed acyclic graph models [17,18,19,20]. In the Granger causality context, there is no such ambiguity and restriction.

1.2. A Historical Viewpoint

In his Nobel Prize lecture in 2003, Clive W. Granger mentioned that in 1959, Denis Gabor pointed out the work of Wiener to him, as a hint to solve some of the difficulties he met in his work. Norbert Wiener’s paper is about the theory of prediction [8]. At the end of his paper, Wiener proposed that prediction theory could be used to define causality between time series. Granger further developed this idea, and came up with a definition of causality and testing procedures [3,21].
In these studies, the essential stones were laid. Granger’s causality states that a cause must occur before the effect, and that causality is relative to the knowledge that is available. This last statement deserves some comment. When testing for causality of one variable on another, it is assumed that the cause has information about the effect that is unique to it; i.e., this information is unknown to any other variable. Obviously, this cannot be verified for variables that are not known. Therefore, the conclusion drawn in a causal testing procedure is relative to the set of measurements that are available. A conclusion reached based on a set of measurements can be altered if new measurements are taken into account.
Mention of information theory is also present in the studies of Granger. In the restricted case of two Gaussian signals, Granger already noted the link between what he called the “causality indices” and the mutual information (Equation 5.4 in [21]). Furthermore, he already foresaw the generalization to the multivariate case, as he wrote in the same paper:
“In the case of q variables, similar equations exist if coherence is replaced by partial coherence, and a new concept of ‘partial information’ is introduced.”
Granger’s paper in 1969 does not contain much new information [3], but rather, it gives a refined presentation of the concepts.
During the 1970’s, some studies, e.g., [4,22,23], appeared that generalized along some of the directions Granger’s work, and related some of the applications to economics. In the early 1980’s, several studies were published that established the now accepted definitions of Granger causality [24,25,26,27]. These are natural extensions of the ideas built upon prediction, and they rely on conditional independence. Finally, the recent studies of Dalhaus and Eichler allowed the definitions of Granger causality graphs [28,29,30]. These studies provide a counterpart of graphical models of multivariate random variables to multivariable stochastic processes.
In two studies published in 1982 and 1984 [31,32], Geweke, another econometrician, set up a full treatment of Granger causality, which included the idea of feedback and instantaneous coupling. In [31], the study was restricted to the link between two time series (possibly multidimensional). In this study, Geweke defined an index of causality from x to y; It is the logarithm of the ratio of the asymptotic mean square error when predicting y from its past only, to the asymptotic mean square error when predicting y from its past and from the past of x. Geweke also defined the same kind of index for instantaneous coupling. When the innovation sequence is Gaussian, the mutual information rate between x and y decomposes as the sum of the indices of causality from x to y and from y to x with the index of instantaneous coupling. This decomposition was shown in the Gaussian case, and it remains valid in any case when the indices of causality are replaced by transfer entropy rates, and the instantaneous coupling index is replaced by an instantaneous information exchange rate. This link between Granger causality and directed information theory was further supported by [33,34,35] (without mention of instantaneous coupling in [34,35]), and the generalization to the non-Gaussian case by [36] (see also [37] for related results). However, prior to these recent studies, the generalization of Geweke’s idea to some general setting was reported in 1987, in econometry by Gouriéroux et al. [38], and in engineering by Rissannen&Wax [39]. Gouriéroux and his co-workers considered a joint Markovian representation of the signals, and worked in a decision-theoretic framework. They defined a sequence of nested hypotheses, whether causality was true or not, and whether instantaneous coupling was present or not. They then worked out the decision statistics using the Kullback approach to decision theory [40], in which discrepancies between hypotheses are measured according to the Kullback divergence between the probability measures under the hypotheses involved. In this setting, the decomposition obtained by Geweke in the Gaussian case was evidently generalised. In [39], the approach taken was closer to Geweke’s study, and it relied on system identification, in which the complexity of the model was taken into account. The probability measures were parameterized, and an information measure that jointly assessed the estimation procedure and the complexity of the model was used when predicting a signal. This allowed Geweke’s result to be extended to nonlinear modeling (and hence the non-Gaussian case), and provided an information-theoretic interpretation of the tests. Once again, the same kind of decomposition of dependence was obtained by these authors. We will see in Section 3 that the decomposition holds due to Kramers causal conditioning. These studies were limited to the bivariate case [38,39].
In the late 1990s, some studies began to develop in the physics community on influences between dynamical systems. A first route was taken that followed the ideas of dynamic system studies for the prediction of chaotic systems. To determine if one signal influenced another, the idea was to consider each of the signals as measured states of two different dynamic systems, and then to study the master-slave relationships between these two systems (for examples, see [41,42,43]). The dynamics of the systems was built using phase space reconstruction [44]. The influence of one system on another was then defined by making a prediction of the dynamics in the reconstructed phase space of one of the processes. To our knowledge, the setting was restricted to the bivariate case. A second route, which was also restricted to the bivariate case, was taken and relied on information-theoretic tools. The main contributions were from Paluš and Schreiber [13,45], with further developments appearing some years later [46,47,48,49]. In these studies, the influence of one process on the other was measured by the discrepancy between the probability measures under the hypotheses of influence or no influence. Naturally, the measures defined very much resembled the measures proposed by Gouriéroux et. al [38], and used the concept of conditional mutual information. The measure to assess whether one signal influences the other was termed transfer entropy by Schreiber. Its definition was proposed under a Markovian assumption, as was exactly done in [38]. The presentation by Paluš [45] was more direct and was not based on a decision-theoretic idea. The measure defined is, however, equivalent to the transfer entropy. Interestingly, Paluš noted in this 2001 paper the closeness of the approach to Granger causality, as per the quotation:
“the latter measure can also be understood as an information theoretic formulation of the Granger causality concept.”
Note that most of these studies considered bivariate analysis, with the notable exception of [46], in which the presence of side information (other measured time series) was explicitely considered.
In parallel with these studies, many others were dedicated to the implementation of Granger causality testing in fields as diverse as climatology (with applications to the controversial questions of global warming) and neuroscience; see [6,7,15,30,50,51,52,53,54], to cite but a few.
In a very different field, information theory, the problem of feedback has lead to many questions since the 1950s. We will not review or cite anything on the problem created by feedback in information theory as this is not within the scope of the present study, but some information can be found in [55]. Instead, we will concentrate on studies that are directly related to the subject of this review. A major breakthrough was achieved by James Massey in 1990 in a short conference paper [12]. Following the ideas of Marko on bidirectional information theory that were developed in the Markovian case [10], Massey re-examined the usual definition of what is called a discrete memoryless channel in information theory, and he showed that the usual definition based on some probabilistic assumptions prohibited the use of feedback. He then clarified the definition of memory and feedback in a communication channel. As a consequence, he showed that in a general channel used with feedback, the usual definition of capacity that relies on mutual information was not adequate. Instead, the right measure was shown to be directed information, an asymmetrical measure of the flow of information. These ideas were further examined by Kramer, who introduced the concept of causal conditioning, and who developed the first applications of directed information theory to communication in networks [11]. After some years, the importance of causal conditioning for the analysis of communication in systems with feedback was realized. Many studies were then dedicated to the analysis of the capacity of channels with feedback and the dual problem of rate-distortion theory [56,57,58,59]. Due to the rapid development in the study of networks (e.g., social networks, neural networks) and of the afferent connectivity problem, more recently many authors made connections between information theory and Granger causality [33,34,36,37,60,61,62]. Some of these studies were restricted to the Gaussian case, and to the bivariate case. Most of these studies did not tackle the problem of instantaneous coupling. Furthermore, several authors realized the importance of directed information theory to assess the circulation of information in networks [1,2,63,64].

1.3. Outline

Tools from directed information theory appear as natural measures to assess Granger causality. Although Granger causality can be considered as a powerful theoretical framework to study influences between signals mathematically, directed information theory provides the measures to test theoretical assertions practically. As already mentioned, these measures are transfer entropy (and its conditional versions), which assesses the dynamical part of Granger causality, and instantaneous information exchange (and its conditional versions), which assesses instantaneous coupling.
This review is structured here as follows. We will first give an overview of the definitions of Granger causality. These are presented in a multivariate setting. We go gradually from weak definitions based on prediction, to strong definitions based on conditional independence. The problem of instantaneous coupling is then discussed, and we show that there are two possible definitions for it. Causality graphs (after Eichler [28]) provide particular reasons to prefer one of these definitions. Section 3 introduces an analysis of Granger causality from an information-theoretic perspective. We insist on the concept of causal conditioning, which is at the root of the relationship studied. Section 4 then highlights the links. Here, we first restate the definitions of Granger causality using concepts from directed information theory. Then from a different point of view, we show how conceptual inference approaches lead to the measures defined in directed information theory. The review then closes with a discussion of some of the aspects that we do not present here intentionally, and on some lines for further research.

1.4. Notations

All of the random variables, vectors and signals considered here are defined in a common probability space ( Ω , B , P ) . They take values either in R or R d , d being some strictly positive integer, or they can even take discrete values. As we concentrate on conceptual aspects rather than technical aspects, we assume that the variables considered are “well behaved”. In particular, we assume finiteness of moments of sufficient order. We assume that continuously valued variables have a measure that is absolutely continuous with respect to the Lebesgue measure of the space considered. Hence, the existence of probability density functions is assumed. Limits are supposed to exist when needed. All of the processes considered in this report are assumed to be stationary.
We work with discrete time. A signal will generically be denoted as x ( k ) . This notation stands also for the value of the signal at time k. The collection of successive samples of the signal, x k , x k + 1 , , x k + n will be denoted as x k k + n . Often, an initial time will be assumed. This can be 0, 1, or . In any case, if we collect all of the sample of the signals from the initial time up to time n, we will suppress the lower index and write this collection as x n .
When dealing with multivariate signals, we use a graph-theoretic notation. This will simplify some connections with graphical modeling. Let V be an index set of finite cardinality | V | . x V = { x V ( k ) , k Z } is a d-dimensional discrete time stationary multivariate process for the probability space considered. For a V , x a is the corresponding component of x V . Likewise, for any subset A V , x A is the corresponding multivariate process ( x a 1 , , x | A | ) . We say that subsets A , B , C form a partition of V if they are disjoint and if A B C = V . The information obtained by observing x A up to time k is given by the filtration generated by { x A ( l ) , l k } . This is denoted as x A k . Furthermore, we will often identify x A with A in the discussion.
The probability density functions (p.d.f.) or probability mass functions (p.m.f) will be denoted by the same notation as p ( x A n ) . The conditional p.d.f. and p.m.f. are written as p ( x A n | x B m ) . The expected value is denoted as E [ . ] , E x [ . ] or E p [ . ] if we want to specify which variable is averaged, or under which probability measure the expected value is evaluated.
Independence between random variables and vectors x and y will be denoted as x y , while conditional independence given z will be written as x y z .

2. Granger’s Causality

The early definitions followed the ideas of Wiener: A signal x causes a signal y if the past of x helps in the prediction of y. Implementing this idea requires the performing of the prediction and the quantification of its quality. This leads to a weak, but operational, form of the definitions of Granger causality. The idea of improving a prediction is generalized by encoding it into conditional dependence or independence.

2.1. From Prediction-Based Definitions…

Consider a cost function g : R k R (k is some appropriate dimension), and the associated risk E [ g ( e ) ] , where e stands for an error term. Let a predictor of x B ( n ) be defined formally as x B ^ ( n + 1 ) = f ( x A n ) , where A and B are subsets of V, and f is a function between appropriate spaces, chosen to minimize the risk with e ( n ) : = x B ( n + 1 ) x B ^ ( n + 1 ) . Solvability may be granted if f is restricted to an element of a given class of functions, such as the set of linear functions. Let F be such a function class. Define:
R F B ( n + 1 ) | A n = inf f F E g x B ( n + 1 ) f ( x A n )
R F B ( n + 1 ) | A n is therefore the optimal risk when making a one-step-ahead prediction of the multivariate signal x B from the past samples of the multivariate signal x A . We are now ready to measure the influence of the past of a process on the prediction of another. To be relatively general and to prepare comments on the structure of the graph, this can be done for subsets of V. We thus choose A and B to be two disjoint subsets of V, and we define C : = V \ ( A B ) (we use ∖ to mean subtraction of a set). We study causality from x A to x B by measuring the decrease in the quality of the prediction of x B ( n ) when excluding the past of x A .
Let R F B ( n + 1 ) | V n be the optimal risk obtained for the prediction of x B from the past of all of the signals grouped in x V . This risk is compared with R F B ( n + 1 ) | ( V \ A ) n , where the past of x A is omitted. Then, for the usual costs functions, we have necessarily:
R F B ( n + 1 ) | V n R F B ( n + 1 ) | ( V \ A ) n
A natural first definition for Granger causality is:
Definition 1. 
x A Granger does not cause x B relative to V if and only if R F B ( n + 1 ) | V n = R F B ( n + 1 ) | ( V \ A ) n
This definition of Granger causality depends on the cost g chosen as well as on the class F of the functions considered. Usually, a quadratic cost function is chosen, for its simplicity and for its evident physical interpretation (a measure of the power of the error). The choice of the class of functions F is crucial. The result of the causality test in definition 1 can change when the class is changed. Consider the very simple example of x n + 1 = α x n + β y n 2 + ε n + 1 , where y n and ε n are zero-mean Gaussian independent and identically distributed (i.i.d.) sequences that are independent of each other. The covariance between x n + 1 and y n is zero, and using the quadratic loss and the class of linear functions, we conclude that y does not Granger cause x, because using a linear function of x n , y n to predict x would lead to the same minimal risk as using a linear function of x n only. However, y n obviously causes x n , but in a nonlinear setting.
The definition is given using the negative of the proposition. If by using the positive way, i.e., R F B ( n + 1 ) | V n < R F B ( n + 1 ) | ( V \ A ) n , Granger proposes to say that x A is a prima facie cause of x B relative to V, prima facie can be translated as “at a first glance”. This is used to insist that if V is enlarged by including other measurements, then the conclusion might be changed. This can be seen as redundant with the mention of the relativity to the observation set V, and we therefore do not use this terminology. However, a mention of the relativity to V must be used, as modification of this set can alter the conclusion. A very simple example of this situation is the chain x n y n z n , where, for example, x n is an i.i.d. sequence, y n + 1 = x n + ε n + 1 , z n + 1 = y n + η n + 1 , ε n , η n being independent i.i.d. sequences. Relative to V = { x , z } , x causes z if we use the quadratic loss and linear functions of the past samples of x (note here that the predictor z n + 1 must be a function of not only x n , but also of x n 1 ). However, if we include the past samples of y and V = { x , y , z } , then the quality of the prediction of z does not deteriorate if we do not use past samples of x. Therefore, x does not cause z relative to V = { x , y , z } .
The advantage of the prediction-based definition is that it leads to operational tests. If the quadratic loss is chosen, working in a parameterized class of functions, such as linear filters or Volterra filters, or even working in reproducing kernel Hilbert spaces, allows the implementation of the definition [65,66,67]. In such cases, the test can be evaluated efficiently from the data. From a theoretical point of view, the quadratic loss can be used to find the optimal function in a much wider class of functions: the measurable functions. In this class, the optimal function for the quadratic loss is widely known to be the conditional expectation [68]. When predicting x B from the whole observation set V, the optimal predictor is written as x B ^ ( n + 1 ) = E x B ( n + 1 ) | x V n . Likewise, elimination of A from V to study its influence on B leads to the predictor x B ^ ( n + 1 ) = E x B ( n + 1 ) | x B n , x C n , where V = C A B . These estimators are of little use, because they are too difficult, or even impossible, to compute. However, they highlight the importance of conditional distributions p ( x B ( n + 1 ) | x V n ) and p ( x B ( n + 1 ) | x B n , x C n ) in the problem of testing whether x A Granger causes x B relative to V or not.

2.2. …To a Probabilistic Definition

The optimal predictors studied above are equal if the conditional probability distributions p ( x B ( n + 1 ) | x V n ) and p ( x B ( n + 1 ) | x B n , x C n ) are equal. These distributions are identical if and only if x B ( n + 1 ) and x A n are independent conditionally to x B n , x C n . A natural extension of definition 1 relies on the use of conditional independence. Once again, let A B C be a partition of V.
Definition 2. 
x A does not Granger cause x B relative to V if and only if x B ( n + 1 ) x A n x B n , x C n , n Z
This definition means that conditionally to the past of x C , the past of x A does not bring more information about x B ( n + 1 ) than is contained in the past of x B .
Definition 2 is far more general than definition 1. If x A does not Granger cause x B relatively to V in the sense of definition 1, it also does not in the sense of definition 2. Then, definition 2 does not rely on any function class and on any cost function. However, it lacks an inherent operational character: the tools to evaluate conditional independence remain to be defined. The assessment of conditional independence can be achieved using measures of conditional independence, and some of these measures will be the cornerstone to link directed information theory and Granger causality.
Note also that the concept of causality in this definition is again a relative concept, and that adding or deleting data from the observation set V might modify the conclusions.

2.3. Instantaneous Coupling

The definitions given so far concern the influence of the past of one process on the present of another. This is one reason that justifies the use of the term “causality”, when the definitions are actually based on statistical dependence. For an extensive discussion on the differences between causality and statistical dependence, we refer to [17].
There is another influence between the processes that is not taken into account by definitions 1 and 2. This influence is referred to as “instantaneous causality” [21,27]. However, we will prefer the term “instantaneous coupling”, specifically to insist that it is not equivalent to a causal link per se, but actually a statistical dependence relationship. The term “contemporaneous conditional independence” that is used in [28] could also be chosen.
Instantaneous coupling measures the common information between x A ( n + 1 ) and x B ( n + 1 ) that is not shared with their past. A definition of instantaneous coupling might then be that x A ( n + 1 ) and x B ( n + 1 ) are not instantaneously coupled if x A ( n + 1 ) x B ( n + 1 ) x A n , x B n , n . This definition makes perfect sense if the observation set is reduced to A and B, a situation we refer to as the bivariate case. However, in general, there is also side information C, and the definition must include this knowledge. However, this presence of side information then leads to two possible definitions of instantaneous coupling.
Definition 3. 
x A and x B are not conditionally instantaneously coupled relative to V if and only if x A ( n + 1 ) x B ( n + 1 ) x A n , x B n , x C n + 1 , n Z , where A B C is a partition of V.
The second possibility is the following:
Definition 4. 
x A and x B are not instantaneously coupled relative to V if and only if x A ( n + 1 ) x B ( n + 1 ) x A n , x B n , x C n , n Z
Note that definitions 3 and 4 are symmetrical in A and B (the application of Bayes theorem). The difference between definitions 3 and 4 resides in the conditioning on x C n + 1 instead of x C n . If the side information up to time n is considered only as in definition 4, the instantaneous dependence or independence is not conditional on the presence of the remaining nodes in C. Thus, this coupling is a bivariate instantaneous coupling: it does measure instantaneous dependence (or independence between A and B) without considering the possible instantaneous coupling between either A and C or B and C. Thus, instantaneous coupling found with definition 4 between A and B does not preclude the possibility that the coupling is actually due to couplings between A and C and/or B and C.
Inclusion of all of the information up to time n + 1 in the conditioning variables allows the dependence or independence to be tested between x A ( n + 1 ) and x B ( n + 1 ) conditionally to x C ( n + 1 ) .
We end up here with the same differences as those between correlation and partial correlation, or dependence and conditional independence for random variables. In graphical modeling, the usual graphs are based on conditional independence between variables [19,20]. These conditional independence graphs are preferred to independence graphs because of their geometrical properties ( e.g., d-separation [17]), which match the Markov properties possibly present in the multivariate distribution they represent. From a physical point of view, conditional independence might be preferable, specifically to eliminate “false” coupling due to third parties. In this respect, conditional independence is not the panacea, as independent variables can be conditionally dependent. The well-known example is the conditional coupling of independent x and y by their addition. Indeed, even if independent, x and y are conditionally dependent to z = x + y .

2.4. More on Graphs

Granger causality graphs were defined and studied in [28]. A causality graph is a mixed graph ( V , E d , E u ) that encodes Granger causality relationships between the components of x V . The vertex set V stores the indices of the components of x V . E d is a set of directed edges between vertices. A directed edge from a to b is equivalent to “ x a Granger causes x b relatively to V”. E u is a set of undirected edges. An undirected edge between x a and x b is equivalent to “ x a and x b are (conditionally if def.4 adopted) instantaneously coupled”. Interestingly, a Granger causality graph may have Markov properties (as in usual graphical models) reflecting a particular (spatial) structure of the joint probability distribution of the whole process { x V t } [28]. A taxonomy of Markov properties (local, global, block recursive) is studied in [28], and equivalence between these properties is put forward. More interestingly, these properties are linked with topological properties of the graph. Therefore, structural properties of the graphs are equivalent to a particular factorization of the joint probability of the multivariate process. We will not continue on this subject here, but this must be known since it paves the way to more efficient inference methods for Granger graphical modeling of multivariate processes (see first steps in this direction in [69,70]).

3. Directed Information Theory and Directional Dependence

Directed information theory is a recent extension of information theory, even if its roots go back to the 1960s and 1970s and the studies of Marko [10]. The developments began in the late 1990s, after the impetus given by James Massey in 1990 [12]. The basic theory was then extended by Gerhard Kramer [11], and then further developed by many authors [56,57,58,59,71] to cite a few. We provide here a short review of the essentials of directed information theory. We will, moreover, adopt a presentation close to the spirit of Granger causality to highlight the links between Granger causality and information theory. We begin by recalling some basics from information theory. Then, we describe the information-theoretic approach to study directional dependence between stochastic processes, first in the bivariate case, and then, from Section 3.5, for networks, i.e., the multivariate case.

3.1. Notation and Basics

Let H ( x A n ) = E [ log p ( x A n ) ] be the entropy of a random vector x A n , the density of which is p. Let the conditional entropy be defined as H ( x A n | x B n ) = E [ log p ( x A n | x B n ) ] . The mutual information I ( x A n ; y B n ) between x A n and x B n is defined as [55]:
I ( x A n ; x B n ) = H ( x B n ) H ( x B n | x A n ) = D K L p ( x A n , x B n ) p ( x A n ) p ( x B n )
where D K L ( p | | q ) = E p [ log p ( x ) / q ( x ) ] is the Kulback–Leibler divergence. D K L ( p | | q ) is 0 if and only if p = q , and it is positive otherwise. The mutual information effectively measures independence since it is 0 if and only if x A n and x B n are independent random vectors. As I ( x A n ; x B n ) = I ( x B n ; x A n ) , mutual information cannot handle directional dependence.
Let x C n be a third time series. It might be a multivariate process that accounts for side information (all of the available observations, but x A n and x B n ). To account for x C n , the conditional mutual information is introduced:
I ( x A n ; x B n | x C n ) = E D K L p ( x A n , x B n | x C n ) | | p ( x A n | x C n ) p ( x B m | x C n )
= D K L p ( x A n , x B n , x C n ) | | p ( x A n | x C n ) p ( x B n | x C n ) p ( x C n )
I ( x A n ; x B n | x C n ) is zero if and only if x A n and x B n are independent conditionally to x C n . Stated differently, conditional mutual information measures the divergence between the actual observations and those that would be observed under the Markov assumption ( x z y ) . Arrows can be misleading here, as by reversibility of Markov chains, the equality above holds also for ( y z x ) . This emphasizes how mutual information cannot provide answers to the information flow directivity problem.

3.2. Directional Dependence between Stochastic Processes; Causal Conditioning

The dependence between the components of the stochastic process x V is encoded in the full generality by the joint probability distributions p ( x V n ) . If V is partitioned into subsets A , B , C , studying dependencies between A and B then requires that p ( x V n ) is factorized into terms where x A and x B appear. For example, as p ( x V n ) = p ( x A n , x B n , x C n ) , we can factorize the probability distribution as p ( x B n | x A n , x C n ) p ( x A n , x C n ) , which appears to emphasize a link from A to B. Two problems appear, however: first, the presence of C perturbs the analysis (more than this, A and C have a symmetrical role here); secondly, the factorization does not take into account the arrow of time, as the conditioning is considered over the whole observations up to time n.
Marginalizing x C out makes it possible to work directly on p ( x A n , x B n ) . However, this eliminates all of the dependence between A and B that might exist via C, and therefore this might lead to an incorrect assessment of the dependence. As for Granger causality, this means that dependence analysis is relative to the observation set. Restricting the study to A and B is what we referred to as the bivariate case, and this allows the basic ideas to be studied. We will therefore present directed information first in the bivariate case, and then turn to the full multivariate case.
The second problem is at the root of the measure of directional dependence between stochastic processes. Assuming that x A ( n ) and x B ( n ) are linked by some physical (e.g., biological, economical) system, it is natural to postulate that their dependence is constrained by causality: if A B , then an event occurring at some time in A will influence B later on. Let us come back to the simple factorization above for the bivariate case. We have p ( x A n , x B n ) = p ( x B n | x A n ) p ( x A n ) , and furthermore (We implicitly choose 1 here as the initial time):
p ( x B n | x A n ) = i = 1 n p x B ( i ) | x B i 1 , x A n
where for i = 1 , the first term is p ( x B ( 1 ) | x A ( 1 ) ) . The conditional distribution quantifies a directional dependence from A to B, but it lacks the causality property mentioned above, as p x B ( i ) | x B i 1 , x A n quantifies the influence of the whole observation x A n (past and future of i) on the present x B ( i ) knowing its past x B i 1 . The causality principle would require the restriction of the prior time i to the past of A only. Kramer defined “causal conditioning” precisely in this sense [11]. Modifying Equation (6) accordingly, we end up with the definition of the causal conditional probability distribution:
p ( x B n x A n ) : = i = 1 n p x B ( i ) | x B i 1 , x A i
Remarkably this provides an alternative factorization of the joint probability. As noted by Massey [12], p ( x A n , y B n ) can then be factorized as ( x B n 1 stands for the delayed collections of samples of x B . If the time origin is finite, 0 or 1, the first element of the list x B n 1 should be understood as a wild card , which does not influence the conditioning.):
p ( x A n , x B n ) = p ( x B n x A n ) p ( x A n x B n 1 )
Assuming that x A is the input of a system that creates x B , p ( x A n x B n 1 ) = i p ( x A ( i ) | x A i 1 , x B i 1 ) characterizes the feedback in the system: each of the factors controls the probability of the input x A at time i conditionally to its past and to the past values of the output x B . Likewise, the term p ( x B n x A n ) = i p ( x B ( i ) | x B i 1 , x A i ) characterizes the direct (or feedforward) link in the system.
Several interesting simple cases occur:
  • In the absence of feedback in the link from A to B, there is the following:
    p ( x A ( i ) | x A i 1 , x B i 1 ) = p ( x A ( i ) | x A i 1 ) , i 2
    or equivalently, in terms of entropies,
    H ( x A ( i ) | x A i 1 , x B i 1 ) = H ( x A ( i ) | x A i 1 ) , i 2
    and as a consequence:
    p ( x A n x B n 1 ) = p ( x A n )
  • Likewise, if there is only a feedback term, then p ( x B ( i ) | x B i 1 , x A i ) = p ( x B ( i ) | x B i 1 ) and then:
    p ( x B n x A n ) = p ( x B n )
  • If the link is memoryless, i.e., the output x B does not depend on the past, then:
    p ( x B ( i ) | x A i , y B i 1 ) = p ( x B ( i ) | x A ( i ) ) i 1
These results allow the question of whether x A influences x B to be addressed. If it does, then the joint distribution has the factorization of Equation (8). However, if x A does not influence x B , then p ( x B n x A n ) = p ( x B n ) , and the factorization of the joint probability distribution simplifies to p ( x A n x B n 1 ) p ( x B n ) . Kullback divergence between the probability distributions for each case generalizes the definition of mutual information to the directional mutual information:
I ( x A n x B n ) = D K L p ( x A n , x B n ) p ( x A n x B n 1 ) p ( x B n )
This quantity measures the loss of information when it is incorrectly assumed that x A does not influence x B . This was called directed information by Massey [12]. Expanding the Kullback divergence allows different forms for the directed information to be obtained:
I ( x A n x B n ) = i = 1 n I x A i ; x B ( i ) | x B i 1
= H ( x B n ) H x B n x A n
where we define the “causal conditional entropy”:
H x B n x A n = E log p x B n x A n
= i = 1 n H x B ( i ) | x B i 1 , x A i
Note that causal conditioning might involve more than one process. This leads to define the causal conditional directed information as:
I ( x A n x B n x C n ) : = H x B n x C n H x B n x A n , x C n = i = 1 n I x A i ; x B ( i ) | x B i 1 , x C i
The basic properties of the directed information were studied by Massey and Kramer [11,12,72], and some are recalled below. As a Kullback divergence, the directed information is always positive or zero. Then, simple algebraic manipulation allows the decomposition to be obtained:
I ( x A n x B n ) + I ( x B n 1 x A n ) = I ( x A n ; x B n )
Equation (20) is fundamental, as it shows how mutual information splits into the sum of a feedforward information flow I ( x A n x B n ) and a feedback information flow I ( x B n 1 x A n ) . In the absence of feedback, p ( x A n x B n 1 ) = p ( x A n ) and I ( x A n ; x B n ) = I ( x A n x B n ) . Equation (20) allows the conclusion that the mutual information is always greater than the directed information, as I ( x B n 1 x A n ) is always positive or zero (as directed information). It is zero if and only if:
I ( x A ( i ) ; x B i 1 | x A i 1 ) = 0   i = 2 , , n
or equivalently:
H ( x A ( i ) | x A i 1 , x B i 1 ) = H ( x A ( i ) | x A i 1 ) i = 2 , , n
This situation corresponds to the absence of feedback in the link A B , whence the fundamental result that the directed information and the mutual information are equal if the channel is free of feedback. This result implies that mutual information over-estimates the directed information between two processes in the presence of feedback. This was thoroughly studied in [11,57,58,59], in a communication-theoretic framework.
The decomposition of Equation (20) is surprising, as it shows that the mutual information is not the sum of the directed information flowing in both directions. Instead, the following decomposition holds:
I ( x A n x B n ) + I ( x B n x A n ) = I ( x A n ; x B n ) + I ( x A n x B n x A n 1 )
where:
I ( x A n x B n | | x A n 1 ) = i I ( x A i ; x B ( i ) | x B i 1 , x A i 1 ) = i I ( x A ( i ) ; x B ( i ) | x B i 1 , x A i 1 )
This demonstrates that I ( x A n x B n ) + I ( x B n x A n ) is symmetrical, but is in general not equal to the mutual information, except if and only if I ( x A ( i ) ; x B ( i ) | x B i 1 , x A i 1 ) = 0 , i = 1 , , n . As the term in the sum is the mutual information between the present samples of the two processes conditioned on their joint past values, this measure is a measure of instantaneous dependence. It is indeed symmetrical in A and B. The term I ( x A n x B n | | x A n 1 ) = I ( x B n x A n | | x B n 1 ) will thus be named the instantaneous information exchange between x A and x B , and will hereafter be denoted as I ( x A n x B n ) . Like directed information, conditional forms of the instantaneous information exchange can be defined, as for example:
I ( x A n x B n x C n ) : = I ( x A n x B n | | x A n 1 , x C n )
which quantifies an instantaneous information exchange between A and B causally conditionally to C.

3.3. Directed Information Rates

Entropy and mutual information in general increase linearly with the length n of the recorded time series. Shannon’s information rate for stochastic processes compensates for the linear growth by considering A ( x ) = lim n + A ( x n ) / n ( if the limit exists), where A ( x n ) denotes any information measure on the sample x n of length n. For the important class of stationary processes (see e.g., [55]), the entropy rate turns out to be the limit of the conditional entropy:
lim n + 1 n H ( x A n ) = lim n + H ( x A ( n ) | x A n 1 )
Kramer generalized this result for causal conditional entropies [11], thus defining the directed information rate for stationary processes as:
I ( x A x B ) = lim n + 1 n i = 1 n I ( x A i ; x B ( i ) | x B i 1 ) = lim n + I ( x A n ; x B ( n ) | x B n 1 )
This result holds also for the instantaneous information exchange rate. Note that the proof of the result relies on the positivity of the entropy for discrete valued stochastic processes. For continuously valued processes, for which the entropy can be negative, the proof is more involved and requires the methods developed in [73,74,75], and see also [58].

3.4. Transfer Entropy and Instantaneous Information Exchange

As introduced by Schreiber in [13,47], transfer entropy evaluates the deviation of the observed data from a model, assuming the following joint Markov property:
p ( x B ( n ) | x B n k + 1 n 1 , x A n l + 1 n 1 ) = p ( x B ( n ) | x B n k + 1 n 1 )
This leads to the following definition:
T ( x A n l + 1 n 1 x B n k + 1 n ) = E log p ( x B ( n ) | x B n k + 1 n 1 , x A n l + 1 n 1 ) p ( x B ( n ) | x B n k + 1 n 1 )
Then T ( x A n l + 1 n 1 x B n k + 1 n ) = 0 if and only if Equation (28) is satisfied. Although in the original definition, the past of x in the conditioning might begin at a different time m n , for practical reasons m = n is considered. Actually, no a priori information is available about possible delays, and setting m = n allows the transfer entropy to be compared with the directed information.
By expressing the transfer entropy as a difference of conditional entropies, we get:
T ( x A n l + 1 n 1 x B n k + 1 n ) = H ( x B ( n ) | x B n k + 1 n 1 ) H ( x B ( n ) | x B n k + 1 n 1 , x A n l + 1 n 1 ) = I ( x A n l + 1 n 1 ; x B ( n ) | x B n k + 1 n 1 )
For l = n = k and choosing 1 as the time origin, the identity I ( x , y ; z | w ) = I ( x ; z | w ) + I ( y ; z | x , w ) leads to:
I ( x A n ; x B ( n ) | x B n 1 ) = I ( x A n 1 ; x B ( n ) | x B n 1 ) + I ( x A ( n ) ; x B ( n ) | x A n 1 , x B n 1 ) = T ( x A n 1 x B n ) + I ( x A ( n ) ; x B ( n ) | x A n 1 , x B n 1 )
For stationary processes, letting n and provided the limits exist, for the rates, we obtain:
I ( x A x B ) = T ( x A x B ) + I ( x A x B )
Transfer entropy is the part of the directed information that measures the influence of the past of x A on the present of x B . However it does not take into account the possible instantaneous dependence of one time series on another, which is handled by directed information.
Moreover, as defined by Schreiber in [13,47], only I ( x A i 1 ; x B ( i ) | x B i 1 ) is considered in T, instead of its sum over i in the directed information. Thus stationarity is implicitly assumed and the transfer entropy has the same meaning as a rate. A sum over delays was considered by Paluš as a means of reducing errors when estimating the measure [48]. Summing over n in Equation (31), the following decomposition of the directed information is obtained:
I ( x A n x B n ) = I ( x A n 1 x B n ) + I ( x A n x B n )
Equation (33) establishes that the influence of one process on another can be decomposed into two terms that account for the past and for the instantaneous contributions. Moreover, this explains the presence of the term I ( x A n x B n ) in the r.h.s. of Equation (23): Instantaneous information exchange is counted twice in the l.h.s. terms I ( x A n x B n ) + I ( x B n x A n ) , but only once in the mutual information I ( x A n ; x B n ) . This allows Equation (23) to be written in a slightly different form, as:
I ( x A n 1 x B n ) + I ( x B n 1 x A n ) + I ( x A n x B n ) = I ( x A n ; x B n )
which is very appealing, as it shows how dependence as measured by mutual information decomposes as the sum of the measures of directional dependences and the measure of instantaneous coupling.

3.5. Accounting for Side Information

The preceding developments aimed at the proposing of definitions of the information flow between x A and x B ; however, whenever A and B are connected to other parts of the network, the flow of information between A and B might be mediated by other members of the network. Time series observed on nodes other than A and B are hereafter referred to as side information. The available side information at time n is denoted as x C n , with A , B , C forming a partition of V. Then, depending on the type of conditioning (usual or causal) two approaches are possible. Usual conditioning considers directed information from A to B that is conditioned on the whole observation x C n . However, this leads to the consideration of causal flows from A to B that possibly include a flow that goes from A to B via C in the future! Thus, an alternate definition for conditioning is required. This is given by the definition of Equation (19) of the causal conditional directed information:
I ( x A n x B n x C n ) : = H x B n x C n H x B n x A n , x C n = i = 1 n I x A i ; x B ( i ) | x B i 1 , x C i
Does the causal conditional directed information decompose as the sum of a causal conditional transfer entropy and a causal conditional instantaneous information exchange, as it does in the bivariate case? Applying twice the chain rule for conditional mutual information, we obtain:
I ( x A n x B n x C n ) = I ( x A n 1 x B n x C n 1 ) + I ( x A n x B n x C n ) + Δ I ( x C n x B n )
In this equation, I ( x A n 1 x B n x C n 1 ) is termed the “causal conditional transfer entropy”. This measures the flow of information from A to B by taking into account a possible route via C. If the flow of information from A to B is entirely relayed by C, the “causal conditional transfer entropy” is zero. In this situation, the usual transfer entropy is not zero, indicating the existence of a flow from A to B. Conditioning on C allows the examination of whether the route goes through C. The term:
I ( x A n x B n x C n ) : = I ( x A n x B n x A n 1 , x C n )
= i = 1 n I ( x A ( i ) ; x B ( i ) | x B i 1 , x A i 1 , x C i )
is the “causal conditional information exchange”. It measures the conditional instantaneous coupling between A and B. The term Δ I ( x C n x B n ) emphasizes the difference between the bivariate and the multivariate cases. This extra term measures an instantaneous coupling and is defined by:
Δ I ( x C n x B n ) = I ( x C n x B n x A n 1 ) I ( x C n x B n )
An alternate decomposition to Equation (36) is:
I ( x A n x B n x C n ) = I ( x A n 1 x B n x C n ) + I ( x A n x B n x C n )
which emphasizes that the extra term comes from:
I ( x A n 1 x B n x C n ) = I ( x A n 1 x B n x C n 1 ) + Δ I ( x C n x B n )
This demonstrates that the definition of the conditional transfer entropy requires conditioning on the past of C. If not, the extra term appears and accounts for instantaneous information exchanges between C and B, due to the addition of the term x C ( i ) in the conditioning. This extra term highlights the difference between the two different natures of instantaneous coupling. The first term,
I ( x C n x B n x A n 1 ) = i I ( x C ( i ) ; x B ( i ) | x A i 1 , x B i 1 , x C i 1 )
describes the intrinsic coupling in the sense that it does not depend on parties other than C and B. The second coupling term,
I ( x C n x B n ) = i I ( x C ( i ) ; x B ( i ) | x B i 1 , x C i 1 )
is relative to the extrinsic coupling, as it measures the instantaneous coupling at time i that is created by variables other than B and C.
As discussed in Section 2.3, the second definition for instantaneous coupling considers conditioning on the past of the side information only. Causally conditioning on x C n 1 does not modify the results of the bivariate case. In particular, we still get the elegant decomposition:
I ( x A n x B n x C n 1 ) = I ( x A n 1 x B n x C n 1 ) + I ( x A n x B n x C n 1 )
and therefore, the decomposition of Equation (34) is generalized to:
I ( x A n 1 x B n x C n 1 ) + I ( x B n 1 x A n x C n 1 ) + I ( x A n x B n x C n 1 ) = I ( x A n ; x B n x C n 1 )
where:
I ( x A n ; x B n x C n 1 ) = i I x A n ; x B ( i ) | x B i 1 , x C i 1
is the causally conditioned mutual information.
Finally, let us consider that for jointly stationary time series, the causal directed information rate is defined similarly to the bivariate case as:
I ( x A x B x C ) = lim n + 1 n i = 1 n I x A i ; x B ( i ) | x B i 1 , x C i
= lim n + I x A n ; x B ( n ) | x B n 1 , x C n
In this section we have emphasized on Kramer’s causal conditioning, both for the definition of directed information and for taking into account side information. We have also shown that Schreiber’s transfer entropy is the part of the directed information that is dedicated to the strict sense of causal information flow (not accounting for simultaneous coupling). The next section more explicitly revisits the links between Granger causality and directed information theory.

4. Inferring Granger Causality and Instantaneous Coupling

Granger causality in its probabilistic form is not operational. In practical situations, for assessing Granger causality between time series, we cannot use the definition directly. We have to define dedicated tools to assess the conditional independence. We use this inference framework to show the links between information theory and Granger causality. We begin by re-expressing Granger causality definitions in terms of some measures that arise from directed information theory. Therefore, in an inference problem, these measures can be used as tools for inference. However, we show in the following sections that these measures naturally emerge from the more usual statistical inference strategies. In the following, and as above, we use the same partitioning of V into the union of disjoint subsets of A, B and C.

4.1. Information-theoretic Measures and Granger Causality

As anticipated in the presentation of directed information, there are profound links between Granger causality and directed information measures. Granger causality relies on conditional independence, and it can also be defined using measures of conditional independence. Information-theoretic measures appear as natural candidates. Recall that two random elements are independent if and only if their mutual information is zero. Moreover, two random elements are independent conditionally to a third one if and only if the conditional mutual information is zero. We can reconsider definitions 2, 3 and 4 and recast them in terms of information-theoretic measures.
Definition 2 stated that x A does not Granger cause x B relative to V if and only if x B ( n + 1 ) x A n x B n , x C n , n 1 . This can be alternatively rephrased into:
Definition 5. 
x A does not Granger cause x B relative to V if and only if I ( x A n 1 x B n x C n 1 ) = 0 n 1 since x B ( i ) x A i x A i 1 , x C i 1 , 1 i n is equivalent to I ( x B ( i ) ; x A i x A i 1 , x C i 1 ) = 0 1 i n .
Otherwise stated, the transfer entropy from A to B causally conditioned on C is zero if and only if A does not Granger cause B relative to V. This shows that causal conditional transfer entropy can be used to assess Granger causality.
Likewise, we can give alternative definitions of instantaneous coupling.
Definition 6. 
x A and x B are not conditionally instantaneously coupled relative to V if and only if I ( x A n x B n x C n ) n 1 ,
or if and only if the instantaneous information exchange causally conditioned on C is zero. The second possible definition of instantaneous coupling is equivalent to:
Definition 7. 
x A and x B are not instantaneously coupled relative to V if and only if I ( x A n x B n x C n 1 ) n 1 ,
or if and only if the instantaneous information exchange causally conditioned on the past of C is zero.
Note that in the bivariate case only (when C is not taken into account), the directed information I ( x A n x B n ) summarizes both the Granger causality and the coupling, as it decomposes as the sum of the transfer entropy I ( x A n 1 x B n ) and the instantaneous information exchange I ( x A n 1 x B n ) .

4.2. Granger Causality Inference

We consider the practical problem of inferring the graph of dependence between the components of a multivariate process. Let us assume that we have measured a multivariate process x V ( n ) for n T . We want to study the dependence between each pair of components (Granger causality and instantaneous coupling between any pair of components relative to V).
We can use the result of the preceding section to evaluate the directed information measures on the data. When studying the influence from any subset A to any subset B, if the measures are zero, then there is no causality (or no coupling); if they are strictly positive, then A Granger causes B relative to V (or A and B are coupled relative to V). This point of view has been adopted in many of the studies that we have already referred to (e.g., [14,16,37,47,76]), and it relies on estimating the measures from the data. We will not review the estimation problem here.
However, it is interesting to examine more traditional frameworks for testing Granger causality, and to examine how directed information theory naturally emerges from these frameworks. To begin with, we show how the measures defined emerge from a binary hypothesis-testing view of Granger causality inference. We then turn to prediction and model-based approaches. We will review how Geweke’s measures of Granger causality in the Gaussian case are equivalent to directed information measures. We will then present a more general case adopted by [37,38,39,77,78,79] and based on a model of the data.

4.2.1. Directed Information Emerges from a Hypotheses-testing Framework

In the inference problem, we want to determine whether or not x A Granger causes (is coupled with) x B relative to V. This can be formulated as a binary hypothesis testing problem. For inferring dependencies between A and B relative to V, we can state the problem as follows.
Assume we observe x V ( n ) , n T . Then, we want to test: “ x A does not Granger cause x B ”, against “ x A causes x B ”; and “ x A and x B are instantaneously coupled” against “ x A are x B not instantaneously coupled”. We will refer to the first test as the Granger causality test, and to the second one, as the instantaneous coupling test.
In the bivariate case, for which the Granger causality test indicates:
H 0 : p 0 ( x B ( i ) x A i 1 , x B i 1 ) = p ( x B ( i ) x B i 1 ) , i T H 1 : p 1 ( x B ( i ) x A i 1 , x B i 1 ) = p ( x B ( i ) x A i 1 , x B i 1 ) , i T
this leads to the testing of different functional forms of the conditional densities of x B ( i ) given the past of x A . The likelihood of the observation under H 1 is the full joint probability p ( x A T , x B T ) = p ( x A T x B T ) p ( x B T x A T 1 ) . Under H 0 we have p ( x B T x A T 1 ) = p ( x B T ) and the likelihood reduces to p ( x A T x B T ) p ( x B T x A T 1 ) = p ( x A T x B T ) p ( x B T ) . The log likelihood ratio for the test is:
l ( x A T , x B T ) : = log p ( x A T , x B T H 1 ) p ( x A T , x B T H 0 ) = log p ( x B T x A T 1 ) p ( x B T )
= i = 1 T log p ( x B ( i ) x A i 1 , x B i 1 ) p ( x B ( i ) x B i 1 )
For example, in the case where the multivariate process is a positive Harris recurrent Markov chain [80], the law of large numbers applies and we have under hypothesis H 1 :
1 T l ( x A T , x B T ) T + T ( x A x B ) a.s.
where T ( x A x B ) is the transfer entropy rate. Thus from a practical point of view, as the amount of data increases, we expect the log likelihood ratio to be close to the transfer entropy rate (under H 1 ). Turning the point of view, this can justify the use of an estimated transfer entropy to assess Granger causality. Under H 0 , 1 T l ( x A T , x B T ) converges to lim T + ( 1 / T ) D K L p ( x A T x B T ) p ( x B T ) p ( x A T x B T ) p ( x B T x A T 1 ) , which can be termed “the Lautum transfer entropy rate” that extends the “Lautum directed information” defined in [71]. Directed information can be viewed as a measure of the loss of information when assuming x A does not causally influence x B when it actually does. Likewise, “Lautum directed information” measures the loss of information when assuming x A does causally influence x B , when actually it does not.
For testing instantaneous coupling, we will use the following:
H 0 : p 0 ( x A ( i ) , x B ( i ) x A i 1 , x B i 1 ) = p ( x A ( i ) x A i 1 , x B i 1 ) p ( x B ( i ) x A i 1 , x B i 1 ) , i T H 1 : p 1 ( x A ( i ) , x B ( i ) x A i 1 , x B i 1 ) = p ( x A ( i ) , x B ( i ) x A i 1 , x B i 1 ) , i T
where under H 0 , there is no coupling. Then, under H 1 and some hypothesis on the data, the likelihood ratio converges almost surely to the information exchange rate I ( x A x B ) .
A related encouraging result due to [71] is the emergence of the directed information in the false-alarm probability error rate. Merging the two tests Equations (48) and (52), i.e., testing both for causality and coupling, or neither, the test is written as:
H 0 : p 0 ( x B ( i ) x A i , x B i 1 ) = p ( x B ( i ) x B i 1 ) , i T H 1 : p 1 ( x B ( i ) x A i , x B i 1 ) = p ( x B ( i ) x A i , x B i 1 ) , i T
Among the tests with a probability of miss P M that is lower than some positive value ε > 0 , the best probability of false alarm P F A follows exp T I ( x A x B ) when T is large. For the case studied here, this is the so-called Stein lemma [55]. In the multivariate case, there is no such result in the literature. An extension is proposed here. However, this is restricted to the case of instantaneously uncoupled time series. Thus, we assume for the end of this subsection that:
p ( x A ( i ) , x B ( i ) , x C ( i ) x A i 1 , x B i 1 , x C i 1 ) = α = A , B , C p ( x α ( i ) x A i 1 , x B i 1 , x C i 1 ) , i T
which means that there is no instantaneous exchange of information between the three subsets that form a partition of V. This assumption has held in most of the recent studies that have applied Granger causality tests. It is, however, unrealistic in applications where the dynamics of the processes involved are faster than the sampling period adopted (see [27] for a discussion in econometry). Consider now the problem of testing Granger causality of A on B relative to V. The binary hypothesis test is given by:
H 0 : p 0 ( x B ( i ) x A i 1 , x B i 1 , x C i 1 ) = p ( x B ( i ) x B i 1 , x C i 1 ) , i T H 1 : p 1 ( x B ( i ) x A i 1 , x B i 1 , x C i 1 ) = p ( x B ( i ) x A i 1 , x B i 1 , x C i 1 ) , i T
The log likelihood ratio reads as:
l ( x A T , x B T , x C T ) = i = 1 T log p ( x B ( i ) x A i 1 , x B i 1 , x C i 1 ) p ( x B ( i ) x B i 1 , x C i 1 )
Again, by assuming that the law of large numbers applies, we can conclude that under H 1
1 T l ( x A T , x B T , x C T ) T + T ( x A x B x C ) a.s.
This means that the causal conditional transfer entropy rate is the limit of the log likelihood ratio as the amount of data increases.

4.2.2. Linear Prediction based Approach and the Gaussian Case

Following definition 1 and focusing on linear models and the quadratic risk R ( e ) = E [ e 2 ] , Geweke introduced the following indices for the study of stationary processes [31,32]:
F x A x B = lim n + R ( x B ( n ) | x B n 1 , x A n 1 ) R ( x B ( n ) | x B n 1 , x A n )
F x A x B x C = lim n + R ( x B ( n ) | x B n 1 , x A n 1 , x C n ) R ( x B ( n ) | x B n 1 , x A n 1 , x C n )
F x A x B = lim n + R ( x B ( n ) | x B n 1 ) R ( x B ( n ) | x B n 1 , x A n 1 )
F x A x B x C = lim n + R ( x B ( n ) | x B n 1 , x C n 1 ) R ( x B ( n ) | x B n 1 , x A n 1 , x C n 1 )
Geweke demonstrated the efficiency of these indices for testing Granger causality and instantaneous coupling (bivariate and multivariate cases). In the particular Gaussian and bivariate case, he gave explicit results for the statistics of the tests, and furthermore he showed that:
F x A x B + F x B x A + F x A x B = I ( x A ; x B )
where I ( x A ; x B ) is the mutual information rate. This relationship, which was already sketched out in [21], is nothing but Equation (34). Indeed, in the Gaussian case, F x A x B = I ( x A x B ) and F x A x B = I ( x A x B ) stem from the knowledge that the entropy rate of a Gaussian stationary process is the logarithm of the asymptotic power of the one-step-ahead prediction [55]. Likewise, we can show that F x A x B x C = I ( x A x B x C ) and F x A x B x C = I ( x A x B x C ) holds.
In the multivariate case, conditioning on the past of the side information, i.e., x C n 1 , in the definition of F x A x B x C , a decomposition analogous to Equation (62) holds and is exactly that given by Equation (44).

4.2.3. The Model-based Approach

In a more general framework, we examine how a model-based approach can be used to test for Granger causality, and how directed information comes into play.
Let us consider a rather general model in which x V ( t ) is a multivariate Markovian process that statisfies:
x V ( t ) = f θ x V t k t 1 + w V ( t )
where f θ : R k | V | R | V | is a function belonging to some functional class F , and where w V is a multivariate i.i.d. sequence, the components of which are not necessarily mutually independent. Function f θ might (or might not) depend on θ, a multidimensional parameter. This general model considers each signal as an AR model (linear or not) with exogeneous inputs; f θ can also stand for a function belonging to some reproducing kernel Hilbert space, which can be estimated from the data [66,67,81]. Using the partition A , B , C , this model can be written equivalently as:
x A ( t ) = f A , θ A x A t k t 1 , x B t k t 1 , x C t k t 1 + w A ( t ) x B ( t ) = f B , θ B x A t k t 1 , x B t k t 1 , x C t k t 1 + w B ( t ) x C ( t ) = f C , θ C x A t k t 1 , x B t k t 1 , x C t k t 1 + w C ( t )
where the functions f . , θ . are the corresponding components of f θ . This relation can be used for inference in a parametric setting: the functional form is assumed to be known and the determination of the function is replaced by the estimation of the parameters θ A , B , C . This can also be used in a nonparametric setting, in which case the function f is searched for in an appropriate functional space, such as an rkHs associated to a kernel [81].
In any case, for studying the influence of x A to x B relative to V, two models are required for x B : one in which x B explicitly depends on x A , and the other one in which x B does not depend on x A . In the parametric setting, the two models can be merged into a single model, in such a way that some components of the parameter θ B are, or not, zero, which depends on whether A causes B or not. The procedure then consists of testing nullity (or not) of these components. In the linear Gaussian case, this leads to the Geweke indices discussed above. In the nonlinear (non-Gaussian) case, the Geweke indices can be used to evaluate the prediction in some classes of nonlinear models (in the minimum mean square error sense). In this latter case, the decomposition of the mutual information, Equation (62), has no reason to remain valid.
Another approach base relies on directly modeling the probability measures. This approach has been used recently to model spiking neurons and to infer Granger causality between several neurons working in the class of generalized linear models [37,79]. Interestingly, the approach has been used either to estimate the directed information [37,77] or to design a likelihood ratio test [38,79]. Suppose we wish to test whether “ x A Granger causes x B relative to V” as a binary hypothesis problem (as in Section 4.2.1.). Forgetting the problem of instantaneous coupling, the problem is then to choose between the hypotheses:
H 0 : p 0 ( x B ( i ) x V i 1 ) = p ( x B ( i ) x V i 1 ; θ 0 ) , i T H 1 : p 1 ( x B ( i ) x V i 1 ) = p ( x B ( i ) x V i 1 ; θ 1 ) , i T
where the existence of causality is entirely reflected into the parameter θ. To be more precise, θ 0 should be seen as a restriction of θ 1 when its components linked to x A are set to zero. As a simple example using the model approach discussed above, consider the simple linear Gaussian model
x B ( t ) = i > 0 θ A ( i ) x A ( t i ) + i > 0 θ B ( i ) x B ( t i ) + i > 0 θ C ( i ) x C ( t i ) + w B ( t )
where w B ( t ) is an i.i.d. Gaussian sequence, and θ A , θ B , θ C are multivariate impulse responses of appropriate dimensions. Define θ 1 = ( θ A , θ B , θ C ) and θ 0 = ( 0 , θ B , θ C ) . Testing for Granger causality is then equivalent to testing θ = θ 1 ; furthermore, the likelihood ratio can be implemented due to the Gaussian assumption. The example developed in [37,79] assumes that the probability that neuron b ( b A C = V ) sends a message at time t ( x b ( t ) = 1 ) to its connected neighbors is given by the conditional probability
Pr x b ( t ) = 1 | x V t ; θ = U ( i > 0 θ A ( i ) x A ( t i ) + i > 0 θ b ( i ) x b ( t i ) + i > 0 θ E b ( i ) x E b ( t i ) + w b ( t ) )
where U is some decision function, the output of which belongs to [ 0 ; 1 ] , A represents the subset of neurons that can send information to b, and E b represents external inputs to b. Defining this probability for all b V completely specifies the behavior of the neural network V.
The problem is a composite hypothesis testing problem, in which parameters defining the likelihoods have to be estimated. It is known that there is no definitive answer to this problem [82]. An approach that relies on an estimation of the parameters using maximum likelihood can be used. Letting Ω be the space where parameter θ is searched for and Ω 0 the subspace where θ 0 lives, then the generalized log likelihood ratio test reads:
l ( x A T , x B T ) : = log sup θ Ω p ( x V T ; θ ) sup θ Ω 0 p ( x V T ; θ ) = log p ( x V T ; θ 1 T ^ ) p ( x V T ; θ 0 T ^ )
where θ i T ^ denotes the maximum likelihood estimator of θ under hypothesis i. In the linear Gaussian case, we will recover exactly the measures developed by Geweke. In a more general case, and as illustrated in Section 4.2.1, as the maximum likelihood estimates are efficient, we can conjecture that the generalized log likelihood ratio will converge to the causal conditional transfer entropy rate if sufficiently relevant conditions are imposed on the models (e.g., Markov processes with recurrent properties). This approach was described in [38] in the bivariate case.

5. Discussion and Extensions

Granger causality was developed originally in econometrics, and it is now transdisciplinary, with the literature on the subject being widely dispersed. We have tried here to sum up the profound links that exist between Granger causality and directed information theory. The key ingredients to build these links are conditional independence and the recently introduced causal conditioning.
We have eluded the important question of the practical use of the definitions and measures presented here. Some of the measures can be used and implemented easily, especially in the linear Gaussian case. In a more general case, different approaches can be taken. The information-theoretic measures can be estimated, or the prediction can be explicitly carried out and the residuals used to assess causality.
Many studies have been carried out over the last 20 years on the problem of estimation of information-theoretic measures. We refer to [83,84,85,86,87] for information on the different ways to estimate information measures. Recent studies into the estimation of entropy and/or information measures are [88,89,90]. The recent report by [76] extensively details and applies transfer entropy in neuroscience using k-nearest neighbors type of estimators. Concerning the applications, important reviews include [14,16], where some of the ideas discussed here are also mentioned, and where practicalities such as the use of surrogate data, for example, are extensively discussed. Applications for neuroscience are discussed in [15,30,50,51,79].
Information-theoretic measures of conditional independence based on Kullback divergence were chosen here to illustrate the links between Granger causality and (usual) directed information theory. Other type of divergence could have been chosen (see e.g., [91,92]); metrics in probability space could also be useful in the assessing of conditional independence. As an illustration, we refer to the study of Fukumizu and co-workers [93], where conditional independence was evaluated using the Hilbert–Schmidt norm of an operator between reproducing kernel Hilbert spaces. The operator generalizes the partial covariance between two random vectors given a third one, and is called the conditional covariance operator. Furthermore, the Hilbert–Schmidt norm of conditional covariance operator can be efficiently estimated from data. A related approach is also detailed in [94].
Many important directions can be followed. An issue is in the time horizon over which the side information is considered in definition 2. As done for instantaneous coupling, we could have chosen to condition by x C n + 1 instead of x C n . This proposition made recently in [35,95] allows in certain circumstances to eliminate the effect of common inputs to A, B and C. It is denoted as partial Granger causality. As noted in [35] this is particularly useful when the common inputs are very powerful and distributed equally likely among all the nodes. If this definition is adopted, then according to Equation (40), the directed information I ( x A n x B n x C n ) decomposes as the sum of instantaneous information exchange I ( x A n x B n x C n ) with the adequate formulation of the transfer entropy for this definition I ( x A n 1 x B n x C n ) . Despite this nice result, a definitive interpretation remains unclear within the probabilistic description presented here. Even in the usual linear setting as developed in [32] this definition leads to some difficulties. Indeed, Geweke’s analysis relies on the possibility to invert the Wold decomposition of the time series, representing the times series as a possibly infinite autoregression with the innovation sequence as input. All the existing dynamical structure (finite order autoregression and moving average input representing exogeneous inputs) is then captured by Geweke’s approach. The analysis in [35,95] assumes that residuals may not be white, and identifiability issues may then arise in this case. Other important issues are the following. Causality between nonstationary processes has rarely been considered (see however [76] for an ad-hoc approach in neuroscience). A very promising methodology is to adopt a graphical modeling way of thinking. The result of [28] on the structural properties of Markov–Granger causality graphs can be used to identify such graphs from real datasets. First steps in this direction were proposed by [69,70]. Assuming that the network under study is a network of sparsely connected nodes and that some Markov properties hold, efficient estimation procedures can be designed, as is the case in usual graphical modeling.

Acknowledgements

P.O.A. is supported by a Marie Curie International Outgoing Fellowship from the European Community. We thank the anonymous reviewers for pointing out to interesting references that were missing in the early draft of the paper.

References

  1. Rao, A.; Hero, A.O.; States, D.J.; Engel, J.D. Inference of Biologically Relevant Gene Influence Networks Using the Directed Information Criterion. In Proceedings of the ICASSP, Toulouse, France, 15–19 May, 2006.
  2. Rao, A.; Hero, A.O.; States, D.J.; Engel, J.D. Motif discovery in tissue-specific regulatory sequences using directed information. EURASIP J. Bioinf. Syst. Biol. 2007, 13853. [Google Scholar] [CrossRef] [PubMed]
  3. Granger, C.W.J. Investigating causal relations by econometrics models and cross-spectral methods. Econometrica 1969, 37, 424–438. [Google Scholar] [CrossRef]
  4. Sims, C.A. Money, income and causality. Am. Econ. Rev. 1972, 62, 540–552. [Google Scholar]
  5. Sporns, O. The Networks of the Brain; MIT Press: Cambridge, MA, USA, 2010. [Google Scholar]
  6. Kaufmann, R.K.; Stern, D.I. Evidence for human influence on climat from hemispheric temperature relations. Nature 1997, 388, 39–44. [Google Scholar] [CrossRef]
  7. Triacca, U. On the use of Granger causality to investigate the human influence on climate. Theor. Appl. Clim. 2001, 69, 137–138. [Google Scholar] [CrossRef]
  8. Wiener, N. The theory of prediction. In Modern Mathematics for the Engineer; MacGrawHill: New York, NY, USA, 1956; pp. 165–190. [Google Scholar]
  9. Granger, C.W.J. Times Series Anlaysis, Cointegration and Applications. In Nobel Lecture; Stockholm, Sweden, 2003. [Google Scholar]
  10. Marko, H. The bidirectional communication theory–a generalization of information theory. IEEE Trans. Commun. 1973, 21, 1345–1351. [Google Scholar] [CrossRef]
  11. Kramer, G. Directed Information for Channels with Feedback. PhD thesis, Swiss Federal Institute of Technology Zurich, 1998. [Google Scholar]
  12. Massey, J.L. Causality, Feedback and Directed Information. In Proceedings of the International Symposium on Information Theory and its Applications, Waikiki, HI, USA, November 1990.
  13. Schreiber, T. Measuring information transfer. Phys. Rev. Lett. 2000, 85, 461–465. [Google Scholar] [CrossRef] [PubMed]
  14. Hlavackova-Schindler, K.; Palus, M.; Vejmelka, M.; Bhattacharya, J. Causality detection based on information-theoretic approaches in time series analysis. Phys. Rep. 2007, 441, 1–46. [Google Scholar] [CrossRef]
  15. Gourévitch, B.; Bouquin-Jeannès, R.L.; Faucon, G. Linear and nonlinear causality between signals: Methods, example and neurophysiological applications. Biol. Cybern. 2006, 95, 349–369. [Google Scholar] [CrossRef] [PubMed]
  16. Palus, M. From nonlinearity to causality: Statistical testing and inference of physical mechanisms underlying complex dynamics. Contemp. Phys. 2007, 48, 307–348. [Google Scholar] [CrossRef]
  17. Pearl, J. Causality: Models, Reasoning and Inference; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  18. Lauritzen, S. Chapter 2: Causal inference from graphical models. In Complex Stochastic Systems; Barndroff-Nielsen, O., Cox, D.R., Kluppelberg, C., Eds.; Chapman and Hall: London, UK, 2001; pp. 63–108. [Google Scholar]
  19. Lauritzen, S. Graphical Models; Oxford University Press: Oxford, UK, 1996. [Google Scholar]
  20. Whittaker, J. Graphical Models in Applied Multivariate Statistics; Wiley& Sons: Weinheim, Germany, 1989. [Google Scholar]
  21. Granger, C.W.J. Economic processes involving feedback. Inf. Control 1963, 6, 28–48. [Google Scholar] [CrossRef]
  22. Caines, P.E.; Chan, C.W. Feedback between stationary stochastic processes. IEEE Trans. Autom. Control 1975, 20, 498–508. [Google Scholar] [CrossRef]
  23. Hosoya, Y. On the granger condition for non-causality. Econometrica 1977, 45, 1735–1736. [Google Scholar] [CrossRef]
  24. Chamberlain, G. The general equivalence of granger and sims causality. Econometrica 1982, 50, 569–581. [Google Scholar] [CrossRef]
  25. Florens, J.P.; Mouchart, M. A note on noncausality. Econometrica 1982, 50, 583–591. [Google Scholar] [CrossRef]
  26. Granger, C.W.J. Some recent developments in a concept of causality. J. Econ. 1988, 39, 199–211. [Google Scholar] [CrossRef]
  27. Granger, C.W.J. Testing for causality: A personal viewpoint. J. Econ. Dyn. Control 1980, 2, 329–352. [Google Scholar] [CrossRef]
  28. Eichler, M. Graphical modeling of multivariate time series. Proba. Theory Relat. Fields 2011. [Google Scholar] [CrossRef]
  29. Dahlaus, R.; Eichler, M. Causality and graphical models in time series analysis. In Highly Structured Stochastic Systems; Green, P., Hjort, N., Richardson, S., Eds.; Oxford University Press: Oxford, UK, 2003. [Google Scholar]
  30. Eichler, M. On the evaluation of information flow in multivariate systems by the directed transfer function. Biol. Cybern. 2006, 94, 469–482. [Google Scholar] [CrossRef] [PubMed]
  31. Geweke, J. Measurement of linear dependence and feedback between multiple time series. J. Am. Stat. Assoc. 1982, 77, 304–313. [Google Scholar] [CrossRef]
  32. Geweke, J. Measures of conditional linear dependence and feedback between times series. J. Am. Stat. Assoc. 1984, 79, 907–915. [Google Scholar] [CrossRef]
  33. Amblard, P.O.; Michel, O.J.J. Sur Différentes Mesures de Dépendance Causales Entre Signaux Alé Atoires (On Different Measures of Causal Dependencies between Random Signals). In Proceedings of the GRETSI, Dijon, France, September 8-11, 2009.
  34. Barnett, L.; Barrett, A.B.; Seth, A.K. Granger causality and transfer entropy are equivalent for Gaussian variables. Phys. Rev. Lett. 2009, 103, 238707. [Google Scholar] [CrossRef]
  35. Barrett, A.B.; Barnett, L.; Seth, A.K. Multivariate Granger causality and generalized variance. Phys. Rev. E 2010, 81, 041907. [Google Scholar] [CrossRef]
  36. Amblard, P.O.; Michel, O.J.J. On directed information theory and Granger causality graphs. J. Comput. Neurosci. 2011, 30, 7–16. [Google Scholar] [CrossRef] [PubMed]
  37. Quinn, C.J.; Coleman, T.P.; Kiyavash, N.; Hastopoulos, N.G. Estimating the directed information to infer causal relationships in ensemble neural spike train recordings. J. Comput. Neurosci. 2011, 30, 17–44. [Google Scholar] [CrossRef] [PubMed]
  38. Gouriéroux, C.; Monfort, A.; Renault, E. Kullback causality measures. Ann. Econ. Stat. 1987, (6–7), 369–410. [Google Scholar]
  39. Rissanen, J.; Wax, M. Measures of mutual and causal dependence between two time series. IEEE Trans. Inf. Theory 1987, 33, 598–601. [Google Scholar] [CrossRef]
  40. Kullback, S. Information Theory and Statistics; Dover: NY, USA, 1968. [Google Scholar]
  41. Quian Quiroga, R.; Arnhold, J.; Grassberger, P. Learning driver-response relashionship from synchronisation patterns. Phys. Rev. E 2000, 61, 5142–5148. [Google Scholar] [CrossRef]
  42. Lashermes, B.; Michel, O.J.J.; Abry, P. Measuring Directional Dependences of Information Flow between Signal and Systems. In Proceedings of the PSIP’03, Grenoble, France, 29–31 January, 2003.
  43. Le Van Quyen, M.; Martinerie, J.; Adam, C.; Varela, F. Nonlinear analyses of interictal eeg map the brain interdependences in human focal epilepsy. Physica D 1999, 127, 250–266. [Google Scholar] [CrossRef]
  44. Kantz, H.; Schreiber, T. Nonlinear Time Series Analysis, 2nd ed.; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  45. Palus, M.; Komarek, V.; Hrncir, Z.; Sterbova, K. Synchronisation as adjustment of information rates: Detection from bivariate time series. Phys. Rev. E 2001, 63, 046211:1–6. [Google Scholar] [CrossRef]
  46. Frenzel, S.; Pompe, B. Partial mutual information for coupling analysis of multivariate time series. Phys. Rev. Lett. 2007, 99, 204101. [Google Scholar] [PubMed]
  47. Kaiser, A.; Schreiber, T. Information transfer in continuous processes. Physica D 2002, 166, 43–62. [Google Scholar] [CrossRef]
  48. Palus, M.; Vejmelka, M. Directionality of coupling from bivariate time series: How to avoid false causalities and missed connections. Phys. Rev. E 2007, 75, 056211:2–056211:14. [Google Scholar] [CrossRef]
  49. Vejmelka, M.; Palus, M. Inferring the directionality of coupling with conditional mutual information. Phys. Rev. E 2008, 77, 026214. [Google Scholar] [CrossRef]
  50. Eichler, M. A graphical approach for evaluating effective connectivity in neural systems. Phil. Trans. R. Soc. B 2005, 360, 953–967. [Google Scholar] [PubMed]
  51. Kaminski, M.; Ding, M.; Truccolo, W.; Bressler, S. Evaluating causal relations in neural systems: Granger causality, directed transfer functions and statistical assessment of significance. Biol. Cybern. 2001, 85, 145–157. [Google Scholar] [PubMed]
  52. Lungarella, M.; Sporns, O. Mapping information flow in sensorimotor networks. PLOS Comput. Biol. 2006, 2, 1301–1312. [Google Scholar] [CrossRef] [PubMed]
  53. Mosedale, T.J.; Stephenson, D.B.; Collins, M.; Mills, T.C. Granger causality of coupled climate processes: Ocean feedback on the north Atlantic oscillation. J. Clim. 2006, 19, 1182–1194. [Google Scholar] [CrossRef]
  54. Saito, Y.; Harashima, H. Recent Advances in EEG and EMG Data Processing; chapter Tracking of information within multichannel EEG record-causal analysis in EEG; Elsevier: Amsterdam, The Netherlands, 1981; pp. 133–146. [Google Scholar]
  55. Cover, J.; Thomas, B. Elements of Information Theory, 2nd ed.; Wiley: Weinheim, Germany, 2006. [Google Scholar]
  56. Kim, Y.H. A coding theorem for a class of stationary channel with feedback. IEEE Trans. Inf. Theory 2008, 54, 1488–1499. [Google Scholar] [CrossRef]
  57. Tatikonda, S.C. Control Under Communication Constraints . In PhD thesis; MIT: Cambridge, MA, USA, 2000. [Google Scholar]
  58. Tatikonda, S.; Mitter, S. The capacity of channels with feedback. IEEE Trans. Inf. Theory 2009, 55, 323–349. [Google Scholar] [CrossRef]
  59. Venkataramanan, R.; Pradhan, S.S. Source coding with feed-forward: Rate-distortion theorems and error exponents for a general source. IEEE Trans Inf. Theory 2007, 53, 2154–2179. [Google Scholar] [CrossRef]
  60. Amblard, P.O.; Michel, O.J.J. Information Flow through Scales. In Proceedings of the IMA Conference on Maths and Signal processing, Cirencester, UK, 16–18 December, 2008; p. 78.
  61. Amblard, P.O.; Michel, O.J.J. Measuring information flow in networks of stochastic processes. 2009; arXiv:0911.2873. [Google Scholar]
  62. Solo, V. On Causality and Mutual Information. In Proceedings of the 47th IEEE conference on Decision and Control, Cancun, Mexico, 9–11 December, 2008.
  63. Kamitake, T.; Harashima, H.; Miyakawa, H.; Saito, Y. A time-series analysis method based on the directed transinformation. Electron. Commun. Jpn. 1984, 67, 1–9. [Google Scholar] [CrossRef]
  64. Al-Khassaweneh, M.; Aviyente, S. The relashionship between two directed information measures. IEEE Sig. Proc. Lett. 2008, 15, 801–804. [Google Scholar] [CrossRef]
  65. Amblard, P.-O.; Michel, O.J.J.; Richard, C.; Honeine, P. A Gaussian Process Regression Approach for Testing Granger Causality between Time Series Data. In Proceedings of the ICASSP, Osaka, Japan, 25–30 March, 2012.
  66. Amblard, P.O.; Vincent, R.; Michel, O.J.J.; Richard, C. Kernelizing Geweke’s Measure of Granger Causality. In Proceedings of the IEEE Workshop on MLSP, Santander, Spain, 23–26 September, 2012.
  67. Marinazzo, D.; Pellicoro, M.; Stramaglia, S. Kernel-Granger causality and the analysis of dynamical networks. Phys. Rev. E 2008, 77, 056215. [Google Scholar] [CrossRef]
  68. Lehmann, E.L.; Casella, G. Theory of Point Estimation, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 1998. [Google Scholar]
  69. Quinn, C.J.; Kiyavas, N.; Coleman, T.P. Equivalence between Minimal Generative Model Graphs and Directed Information Graph. In Proceeding of the ISIT, St. Petersburg, Russia, 31 July–5 August, 2011.
  70. Runge, J.; Heitzig, J.; Petoukhov, V.; Kurths, J. Escping the curse of dimensionality in estimating multivariate transfer entropy. Phys. Rev. Lett. 2012, 108, 258701. [Google Scholar] [CrossRef] [PubMed]
  71. Permuter, H.H.; Kim, Y.-H.; Weissman, T. Interpretations of directed information in portfolio theory, data compression, and hypothesis testing. IEEE Trans. Inf. Theory 2011, 57, 3248–3259. [Google Scholar] [CrossRef]
  72. Massey, J.L.; Massey, P.C. Conservation of Mutual and Directed Information. In Proceedings of the International Symposium on Information Theory and its Applications, Adelalaïde, Australia, 4–7 September, 2005.
  73. Gray, R.M.; Kieffer, J.C. Mutual information rate, distorsion and quantization in metric spaces. IEEE Trans. Inf. Theory 1980, 26, 412–422. [Google Scholar] [CrossRef]
  74. Gray, R.M. Entropy and Information Theory; Springer-Verlag: Berlin/Heidelberg, Germany, 1990. [Google Scholar]
  75. Pinsker, M.S. Information and Information Stability of Random Variables; Holden Day: San Francisco, USA, 1964. [Google Scholar]
  76. Vicente, R.; Wibral, M.; Lindner, M.; Pipa, G. Transfer entropy–a model-free measure of effective connectivity for the neurosciences. J. Comput. Neurosci. 2011, 30, 45–67. [Google Scholar] [CrossRef] [PubMed]
  77. Barnett, L.; Bossomaier, T. Transfer entropy as log-likelihood ratio. Phys. Rev. Lett. 2012, 109, 138105. [Google Scholar] [CrossRef] [PubMed]
  78. Kim, S.; Brown, E.N. A General Statistical Framework for Assessing Granger Causality. In Proceedings of IEEE Icassp, Prague, Czech Republic, 22–27 May, 2010; pp. 2222–2225.
  79. Kim, S.; Putrino, D.; Ghosh, S.; Brown, E.N. A Granger causality measure for point process models of ensembled neural spiking activity. PLOS Comput. Biol. 2011. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Meyn, S.; Tweedie, R.L. Markov Chains and Stochastic Stability, 2nd ed.; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar]
  81. Schölkopf, B.; Smola, A.J. Learning with Kernels; MIT Press: Cambridge, MA, USA, 2002. [Google Scholar]
  82. Lehmann, E.L.; Romano, J.P. Testing Statistical Hypotheses, 3rd ed.; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  83. Beirlant, J.; Dudewicz, E.J.; Gyorfi, L.; van Der Meulen, E.C. Nonparametric entropy estimation: An overview. Int. J. Math. Stat. Sci. 1997, 6, 17–39. [Google Scholar]
  84. Goria, M.N.; Leonenko, N.N.; Mergell, V.V.; Novi Invardi, P.L. A new class of random vector entropy estimators and its applications in testing statistical hypotheses. J. Nonparam. Stat. 2005, 17, 277–297. [Google Scholar] [CrossRef]
  85. Kraskov, A.; Stogbauer, H.; Grassberger, P. Estimating mutual information. Phys. Rev. E 2004, 69, 066138. [Google Scholar] [CrossRef]
  86. Kozachenko, L.F.; Leonenko, N.N. Sample estimate of the entropy of a random vector. Problems Inf. Trans. 1987, 23, 95–101. [Google Scholar]
  87. Paninski, L. Estimation of entropy and mutual information. Neural Comput. 2003, 15, 1191–1253. [Google Scholar] [CrossRef]
  88. Leonenko, N.N.; Pronzato, L.; Savani, V. A class of Rényi information estimators for multidimensional densities. Ann. Stat. 2008, 36, 2153–2182. [Google Scholar] [CrossRef] [Green Version]
  89. Sricharan, K.; Raich, R.; Hero, A.O. Estimation of non-linear functionals of densities with confidence. IEEE Trans. Inf. Theory 2012, 58, 4135–4159. [Google Scholar] [CrossRef] [Green Version]
  90. Wang, Q.; Kulkarni, S.; Verdu, S. Divergence estimation for multidimensional densities via-nearest-neighbor distances. IEEE Trans. Inf. Theory 2009, 55, 2392–2405. [Google Scholar] [CrossRef]
  91. Basseville, M. Divergence measures for statistical data processing—an annotated bibliography. Signal Process. 2012, in press. [Google Scholar] [CrossRef]
  92. Bercher, J.F. Escort entropies and divergences and related canonical distribution. Phys. Lett. A 2011, 375, 2969–2973. [Google Scholar] [CrossRef] [Green Version]
  93. FukumizuU, K.; Gretton, A.; Sun, X.; Scholkopf, B. Kernel Measures of Conditional Dependence. In NIPS, Vancouver, Canada, 3–8 December, 2007.
  94. Seth, S.; Príncipe, J.C. Assessing Granger non-causality using nonparametric measure of conditional independence. IEEE Trans. Neural Netw. Learn. Syst. 2012, 23, 47–59. [Google Scholar] [CrossRef] [PubMed]
  95. Guo, S.; Seth, A.K.; Kendrick, K.M.; Zhou, C.; Feng, J. Partial Granger causality–eliminating exogeneous inputs and latent variables. J. Neurosci. Methods 2008, 172, 79–93. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Amblard, P.-O.; Michel, O.J.J. The Relation between Granger Causality and Directed Information Theory: A Review. Entropy 2013, 15, 113-143. https://doi.org/10.3390/e15010113

AMA Style

Amblard P-O, Michel OJJ. The Relation between Granger Causality and Directed Information Theory: A Review. Entropy. 2013; 15(1):113-143. https://doi.org/10.3390/e15010113

Chicago/Turabian Style

Amblard, Pierre-Olivier, and Olivier J. J. Michel. 2013. "The Relation between Granger Causality and Directed Information Theory: A Review" Entropy 15, no. 1: 113-143. https://doi.org/10.3390/e15010113

Article Metrics

Back to TopTop