Publisher of the most cited open-access journals in their fields.

Original Research ARTICLE

Front. Robot. AI, 08 January 2016 | https://doi.org/10.3389/frobt.2015.00035

Hierarchical Quantification of Synergy in Channels

  • 1Max Planck Institute for Mathematics in the Sciences, Leipzig, Germany
  • 2Faculty of Mathematics and Computer Science, University of Leipzig, Leipzig, Germany
  • 3Santa Fe Institute, Santa Fe, NM, USA

The decomposition of channel information into synergies of different order is an open, active problem in the theory of complex systems. Most approaches to the problem are based on information theory and propose decompositions of mutual information between inputs and outputs in several ways, none of which is generally accepted yet. We propose a new point of view on the topic. We model a multi-input channel as a Markov kernel. We can project the channel onto a series of exponential families, which form a hierarchical structure. This is carried out with tools from information geometry in a way analogous to the projections of probability distributions introduced by Amari. A Pythagorean relation leads naturally to a decomposition of the mutual information between inputs and outputs into terms, which represent single node information, pairwise interactions, and in general n-node interactions. The synergy measures introduced in this paper can be easily evaluated by an iterative scaling algorithm, which is a standard procedure in information geometry.

1. Introduction

In complex systems like biological networks, for example neural networks, a basic principle is that their functioning is based on the correlation and interaction of their different parts. While correlation between two sources is well understood, and can be quantified by Shannon’s mutual information (see, for example, Kakihara (1999)), there is still no generally accepted theory for interactions of three nodes or more. If we label one of the nodes as “output,” the problem is equivalent to determine how much two (or more) input nodes interact to yield the output. This concept is known in common language as “synergy,” which means “working together,” or performing a task that would not be feasible by the single parts separately.

There are a number of important works which address the topic, but the problem is still considered open. The first generalization of mutual information was interaction information (introduced in McGill (1954)), defined for three nodes in terms of the joint and marginal entropies:

I(X:Y:Z)=H(X,Y,Z)+H(X,Y)+H(X,Z)+H(Y,Z)+H(X)H(Y)H(Z).

Interaction information is defined symmetrically on the joint distribution, but most approaches interpret it by looking at a channel, rather than a joint distribution, (X,Y) → Z. For example, we can rewrite equation (1) equivalently in terms of mutual information (choosing Z as “output”):

I(X:Y:Z)=I(X,Y:Z)I(X:Z)I(Y:Z),

where we see that it can mean intuitively “how much the whole (X,Y) gives more (or less) information about Z than the sum of the parts separately.” Another expression, again equivalent, is

I(X:Y:Z)=I(X:Y|Z)I(X:Y),

which we can interpret as “how much conditioning over Z changes the correlation between X and Y ” (see Schneidmann et al. (2003a)). Unlike mutual information, interaction information carries a sign:

I > 0: synergy. Conditioning on one node increases the correlation between the remaining nodes. Or, the whole gives more information than the sum of the parts. Example: XOR function.

I < 0: redundancy. Conditioning on one node decreases, or explains away the correlation between the remaining nodes. Or, the whole gives less information than the sum of the parts. Example: X = Y = Z.

I = 0: 3-independence. Conditioning on one node has no effect on the correlation between the remaining nodes. Or, the whole gives the same amount of information as the parts separately. The nodes can nevertheless still be conditionally dependent. Example: independent nodes.1

As argued in Schneidmann et al. (2003b), Williams and Beer (2010), and Griffith and Koch (2014), however, this is not the whole picture. There are systems which exhibit both synergetic and redundant behavior, and interaction information only quantifies the difference of synergy and redundancy, with a priori no way to tell the two apart. In a system with highly correlated inputs, for example, the synergy would remain unseen, as it would be canceled by the redundancy. Moreover, this picture breaks down for more than three nodes. Another problem pointed out in Schneidmann et al. (2003b) and Amari (2001) is that redundancy (as for example in X = Y = Z) can be described in terms of pairwise interactions, not triple, while synergy (as in the XOR function) is purely threewise. Therefore, I compares and mixes information quantities of different nature.

A detailed explanation of the problem for two inputs is presented in Williams and Beer (2010) and it yields a decomposition (“Partial Information Decomposition,” PID) as follows: there exist two non-negative quantities, Synergy and Redundancy, such that

I(X,Y:Z)=I(X:Z)+I(Y:Z)+SynRed,

or equivalently:

I(X:Y:Z)=SynRed.

Moreover, they define unique information for the inputs X and Y as

UI(X)=I(X:Z)Red,
UI(Y)=I(Y:Z)Red,

so that the total mutual information is decomposed positively:

I(X,Y:Z)=UI(X)+UI(Y)+Red+Syn.

What these quantities intuitively mean is

• Redundancy – information available in both inputs;

• Unique information – information available only in one of the inputs;

• Synergy – information available only when both inputs are present, arising purely from their interaction.

In this formulation, if one finds a measure of synergy, one can automatically define compatible measures of redundancy and unique information (and vice versa), provided that the measure of synergy is always larger or equal to I(X:Y :Z), and that the resulting measure of redundancy is less or equal than I(X:Z) and I(Y :Z). Synergy, redundancy, and unique information are defined on a channel, and choosing a different channel with the same joint distribution (e.g., (Y, Z) → X) may yield a different decomposition.

In Griffith and Koch (2014), an overview of (previous) measures of synergy and their shortcomings in standard examples is presented. In the same paper is then presented a newer measure for synergy, defined equivalently in Bertschinger et al. (2014) as

CI(X,Y;Z):=I(X,Y:Z)minpIp(X,Y:Z),

where ∧ is the space of distributions with prescribed marginals:

={qP(X,Y,Z)|q(X,Z)=p(X,Z),q(Y,Z)=p(Y,Z)}.

This measure satisfies interesting properties (proven in Griffith and Koch (2014) and Bertschinger et al. (2014)), which make it compatible with Williams and Beer’s PID, and with the intuition in most examples. However, it was proven in Rauh et al. (2015) that such an approach can not work in the desired way for more than three nodes (two inputs).

Our approach uses information geometry (Amari and Nagaoka, 2000), extending previous work on hierarchical decompositions (Amari, 2001) and complexity (Ay, 2015). (Compare the related approach on information decomposition pursued in Harder et al. (2013).) The main tools of the present paper are KL-projections and the Pythagorean relation that they satisfy. This allows (as in Amari (2001)) to form hierarchies of interactions of different orders in a geometrical way. In the present problem, we decompose mutual information between inputs and outputs of a channel k, for two inputs, as

I(X,Y:Z)=d1(k)+d2(k),

where d2 quantifies synergy (as in Equation (8)), and d1 integrates all the lower order terms (UI, Red), quantifying the so-called union information (see Griffith and Koch (2014)). One may want to use this measure of synergy to form a complete decomposition analogous to Equation (8), but this does not work, as in general it is not true that d2I(X:Y :Z). For this reason, we keep the decomposition more coarse, and we do not divide union information into unique and redundant.

For more inputs X1, … , XN, the decomposition generalizes to

I(X1,,XN:Z)=d1(k)++dN(k)=i=1Ndi(k),

where higher orders of synergy appear.

Until now, there seems to be no way of rewriting the decomposition of Griffith and Koch (2014) and Bertschinger et al. (2014) in a way consistent with information geometry, and more in general, Williams and Beer’s PID seems hard to write as a geometric decomposition. A comparison between d2 and the measure CI of Griffith and Koch (2014) and Bertschinger et al. (2014) is presented in Section 5. There we show that d2CI, and we argue, with a numerical example, that CI overestimates synergy at least in one case.

For a small number of inputs (≲ 5), our quantities are easily computable with the standard algorithms of information geometry (like iterative scaling (Csiszár and Shields, 2004)). This allowed us to get precise quantities for all the examples considered.

1.1. Technical Definitions

We consider a set of N input nodes V = {1, … , N}, taking values in the sets X1, … , XN, and an output node, taking values in the set Y. We write the input globally as X: = X1 × … × XN. For example, in biology, Y can be the phenotype and X can be a collection of genes determining Y. We denote by F(Y) the set of real functions on Y and with P(X) the set of probability measures on X.

We can model the channel from X to Y as a Markov kernel (called also stochastic kernel, transition kernel, or stochastic map) k, that assigns to each xX a probability measure on Y (for a detailed treatment, see Kakihara (1999)). Here, we will consider only finite systems, so we can think of a channel simply as a transition matrix (or stochastic matrix) whose rows sum to one.

k(x;y)0x,y;yk(x;y)=1x.

The space of channels from X to Y will be denoted by K(X;Y). We will denote by X and Y also the corresponding random variables, whenever this does not lead to confusion.

Conditional probabilities define channels: if p(X,Y) ∈ P(X,Y) and the marginal p(X) is strictly positive, then p(Y |X) ∈ K(X;Y) is a well-defined channel. Vice versa, if kK(X;Y), given pP(X) we can form a well-defined joint probability:

pk(x,y):=p(x)k(x;y)x,y.

An “input distribution” pP(X) is crucial also to extend the notion of divergence from probability distributions to channels. The most natural way of doing it is the following.

Definition 1: Let pP(X), let k, mK(X;Y). Then

Dp(k||m):=x,yp(x)k(x;y)logk(x;y)m(x;y).

Defined this way, Dp is affine in p. Moreover, it has an important compatibility property. Let p,q be joint probability distributions on X × Y, and let D be the KL-divergence. Then

D(p(X,Y)||q(X,Y))=D(p(X)||q(X))+Dp(X)(p(Y|X)||q(Y|X)).

We will now illustrate our geometric ideas in channels with one, two, and three input nodes, and then we present some examples. The general case will be addressed in Section 4.

2. Geometric Idea of Synergy

2.1. Mutual Information as Motivation

It is a well-known fact in information theory that Shannon’s mutual information can be written as a KL-divergence:

Ip(X:Y)=D(p(X,Y)||p(X)p(Y)).

From the point of view of information geometry, this can be interpreted as a “distance” between the real distribution and a product distribution that has exactly the same marginals, but maximal entropy (see Figure 1). In other words, we have

Ip(X:Y)=infqP(X)rP(Y)D(p(X,Y)||q(X)r(Y)).
FIGURE 1
www.frontiersin.org

Figure 1. For two binary nodes, the family of product distributions is a surface in a 3-dimensional simplex.

The distribution given by p(X)p(Y) is optimal in the sense that

p(X)p(Y)=arg minqP(X)rP(Y)D(p(X,Y)||q(X)r(Y)).

The divergence between p and a submanifold is, as usual in geometry, the “distance” between p and the “closest point” on that submanifold, which in our case is the geodesic projection w.r.t. the mixture connection.

2.2. Extension to Channels

We can use the same insight with channels. Instead of a joint distribution on N nodes, we consider a channel from an input X to an output Y. Suppose we have a family E of channels, and a channel k that may not be in E. Then, just as in geometry, we can define the “distance” between k and E.

Definition 2: Let p be an input distribution. The divergence between a channel k and a family of channels E is given by

Dp(k||E):=infmEDp(k||m).

If the minimum is uniquely realized, we call the channel

πEk:=arg minmEDp(k||m)

the KL-projection of k on E (and simply “a” KL-projection if it is not unique).

We will always work with compact families, so the minima will always be realized, and for strictly positive p they will be unique (see Section 4 for the details).

We will consider families E for which the KL-divergence satisfies a Pythagorean equality (see Figure 2 below for some intuition):

Dp(k||m)=Dp(k||πEk)+Dp(πEk||m)

for every mE. These families (technically, closures of exponential families) are defined in Section 4.

FIGURE 2
www.frontiersin.org

Figure 2. Illustration of the Pythagoras theorem for projections.

2.3. One Input

Consider first one-input node X, with input distribution p(X), and one output node Y. A constant channel k in K(X;Y) is a channel whose entries do not depend on X (more precisely: k(x;y) = k(x′;y) for any x,x′, y). This denomination is motivated by the following properties:

• They correspond to channels that do not use any information from the input to generate the output.

• The output distribution given by k is a probability distribution on Y which does not depend on X.

• Deterministic constant channels are precisely constant functions.

We call E0 the family of constant channels. Take now any channel kK(X;Y). If we want to quantify the dependence in k of Y on X, we can then look at the divergence of k from the constant channels:

d1(k):=Dp(k||E0).

The minimum is realized in πE0k. We have that

d1(k)=Dp(k||πE0k)=x,yp(x)k(x;y)logk(x;y)πE0k(y)
=HpπE0k(Y)Hpk(Y|X)=Ipk(X:Y),

so that consistently with our intuition, the dependence of Y on X is just the mutual information. From the channel point of view, it is simply the divergence from the constant channels. (A rigorous calculation is done in Section 4.)

2.4. Two Inputs

Consider now two input nodes with input probability p and one output node. We can again define the family E0 of constant channels, and the same calculations give

Dp(k||E0)=Ipk(X1,X2:Y).

This time, though, we can say a lot more: the quantity above can be decomposed. In analogy with the independence definition for probability distributions, we would like to define a split channel as a product channel of its parts: p(y|x1, x2) = p(y|x1)p(y|x2). Unfortunately, the term on the right would be in general not normalized, so we replace the condition by a weaker one. We call the channel k(X1, X2; Y) split if it can be written as

k(x1,x2;y)=ϕ0(x1,x2)ϕ1(x1;y)ϕ2(x2;y)

for some functions ϕ0, ϕ1, and ϕ2, which in general are not themselves channels (in particular, ϕi(xi; y) ≠ p(y|xi)). We call E1 the family of split channels. This family corresponds to those channels that do not have any synergy. This is a special case of an exponential family, analogous to the family of product distributions of Figure 1. The examples “single node” and “split channel” in the next section belong exactly to this family. Take now any channel k(X1, X2; Y). In analogy with mutual information, we call synergy the divergence:

d2(k):=Dp(k||E1).

Simply speaking, our synergy is quantified as the deviation of the channel from the set E1 of channels without synergy.

We can now project k first to E1, and then to E0. Since E0 is a subfamily of E1, the following Pythagoras relation holds from Equation (22):

Dp(k||πE0k)=Dp(k||πE1k)+Dp(πE1k||πE0k).

If, in analogy with the one-input case, we call the last quantity d1, we get from Equations (26) and (28):

Ipk(X1,X2:Y)=d2(k)+d1(k).

The term d1 measures how much information comes from single nodes (but it does not tell which nodes). The term d2 measures how much information comes from the synergy of X1 and X2 in the channel. The example “XOR” in the next section will show this.

If we call E2 the whole K(X;Y), we get E0E1E2 and

di(k):=Dp(πEik||πEi1k).

2.5. Three Inputs

Consider now three nodes X1, X2, and X3 with input probability p and a channel k. We have again

Dp(k||E0)=Ipk(X1,X2,X3:Y).

This time we can decompose the mutual information in different ways. We can, for example, look at split channels, i.e., in the form:

k(x1,x2,x3;y)=ϕ0(x)ϕ1(x1;y)ϕ2(x2;y)ϕ3(x3;y)

for some ϕ0, ϕ1, ϕ2, and ϕ3. As in the previous case, we call this family E1. Or we can look at more interesting channels, the ones in the form:

k(x1,x2,x3;y)=ϕ0(x)ϕ12(x1,x2;y)ϕ13(x1,x3;y)ϕ23(x2,x3;y)

for some ϕ0, ϕ12, ϕ13, and ϕ23. We call this family E2, and it is easy to see that

E0E1E2E3,

where E0 denotes again the constant channels, and E3 denotes the whole K(X;Y). We define again

di(k):=Dp(πEik||πEi1k).

This time, the Pythagorean relation can be nested, and it gives us

Ipk(X1,X2,X3:Y)=d3(k)+d2(k)+d1(k).

The difference between pairwise synergy and threewise synergy is shown in the “XOR” example in the next section.

Now that we have introduced the measure for a small number of inputs, we can study the examples from the literature (Griffith and Koch, 2014) and show that our measure is consistent with the intuition. The general case will be more in rigor introduced in Section 4.

3. Examples

Here, we present some examples of decomposition for well-known channels. All the quantities have been computed using an algorithm analogous to iterative scaling (as in Csiszár and Shields (2004)).

3.1. Single Node Channel

The easiest example is considering a channel which only depends on X1, i.e.,

I(X:Y)=I(X1:Y).

For example, consider 3 binary input nodes X1, X2, and X3 with constant input probability and one binary output node Y which is an exact copy of X1.

Then, we have exactly one bit of single node information and no higher order terms (see Figure 3). Geometrically, k lies in E1, so the only non-zero divergence in equation (37) is d1(k). As one would expect, d2(k) and d3(k) vanish, as there is no synergy of orders 2 and 3.

FIGURE 3
www.frontiersin.org

Figure 3. Synergies of different orders for the single-node channel, Example 3.1.

3.2. Split Channel

The second easiest example is a more general channel which obeys equation (33). In particular, consider 3 binary input nodes X1, X2, and X3 with constant input probability (so, the xi are independent), and output Y = X1 × X2 × X3. As channel, we simply take the identity map (x1, x2, x3) ↦ (x1, x2, x3) ∈ Y. In this particular case:

I(X:Y)=iI(Xi;Y).

We have 3 bits of mutual information (see Figure 4), which are all single node (but from different nodes). Since

k(x1,x2,x3;y)=ϕ1(x1;y1)ϕ2(x2;y2)ϕ3(x3;y3),

which is a special case of Equation (33), kE1, and so d2(k) and d3(k) in Equation (37) are again zero.

FIGURE 4
www.frontiersin.org

Figure 4. Synergies of different orders for the split channel, Example 3.2.

3.3. Correlated Inputs

Consider 3 perfectly correlated binary nodes, each one with uniform marginal probability. As output take a perfect copy of one (hence, all) of the inputs. We have again one bit of mutual information, which could come from any of the nodes, but no synergy, as no two nodes are interacting in the channel. The input distribution has correlation, but this has no effect on the channel, since the channel is simply copying the value of X1 (or X2 or X3, equivalently). Therefore, again kE1. Of the terms in Equation (37), again the only non-zero is d1(k) (see Figure 5).

FIGURE 5
www.frontiersin.org

Figure 5. Synergies of different orders for correlated inputs, Example 3.3.

This example in the literature is used to motivate the notion of redundancy. A “redundant channel” is in our decomposition exactly equivalent to a single node channel, since it contains exactly the same amount of information.

3.4. Parity (XOR)

The standard example of synergy is given by the XOR function and more generally by the parity function between two or more nodes.

For example, consider 3 binary input nodes X1, X2, and X3 with constant input probability, and one binary output node Y which is given by X1X2. We have 1 bit of mutual information, which is purely arising from a pairwise synergy (of X1 and X2), so this time kE2. The function XOR is pure synergy, so d2(k) is the only non-zero term in Equation (37) (see Figure 6).

FIGURE 6
www.frontiersin.org

Figure 6. Synergies of different orders for the binary XOR function, Example 3.4.

If instead Y is given by the threewise parity function, or X1X2X3, we have again 1 bit of mutual information, which now is purely arising from a threewise synergy, so here kE3, and the only non-zero term in Equation (37) is d3(k) (see Figure 7).

FIGURE 7
www.frontiersin.org

Figure 7. Synergies of different orders for the three wise parity function, Example 3.4.

In these examples, there are no terms of lower order synergy, but the generic elements of E2 and E3 usually do contain a non-zero lower part. Consider, for instance, the next example.

3.5. AND and OR

The other two standard logic gates, AND and OR, share the same decomposition. Consider two binary nodes X1, and X2 with uniform probability, and let Y be X1X2 (or X1X2). There is again one bit of mutual information, which comes mostly from single nodes, but also from synergy (see Figure 8).

FIGURE 8
www.frontiersin.org

Figure 8. Synergies of different orders for the AND (and OR) function, Example 3.5.

Geometrically, this means that AND and OR are channels in E2, which lie close to the submanifold E1.

3.6. XorLoses

Here, we present a slightly more complicated example, coming from Griffith and Koch (2014). We have three binary nodes X1, X2, and X3, where X1 and X2 have uniform probabilities, and an output node Y = X1X2, just like in the “XOR” example. Now we take X3 to be perfectly correlated with Y = X1X2, so that Y could get the information either from X3 or from the synergy between X1 and X2. We have one bit of mutual information, which can be seen as entirely coming from X3, and so the synergy between X1 and X2 is not adding anything (see Figure 9).

FIGURE 9
www.frontiersin.org

Figure 9. Synergies of different orders for the XorLoses channel, Example 3.6.

3.7. XorDuplicate

Again from Griffith and Koch (2014), we have 3 binary nodes X1, X2, and X3, where X1 and X2 have uniform probabilities, while X3 = X1. The output is X1X2 = X3X2, so it could get the information either from the synergy between X1 and X2 or X2 and X3. There is one bit of mutual information, which is coming from a pairwise interaction (see Figure 10). Again, it does not matter between whom.

FIGURE 10
www.frontiersin.org

Figure 10. Synergies of different orders for the XorDuplicate channel, Example 3.7.

It should be clear from the examples here that decomposing only by order, and not by the specific subsets, is crucial. For example, in the “input correlation” example, there is no natural way to decide from which single node the information comes, even if it is clear that the interaction is of order 1.

4. General Case

Here, we try to give a general formulation, for N inputs, of the quantities defined in Section 2. As in the Section “Introduction,” we call the set of input nodes V of cardinality N, and we consider a subset I of the nodes. We denote the joint random variable (Xi, iI) by XI, and we denote the complement of I in V by Ic. The case N = 3 in Section 2 should motivate the following definition.

Definition 3: Let IV. We call FI the space of functions that only depend on XI and Y :

FI:={fF(X,Y)|f(xI,xIc;y)=f(xI,xIc;y)xIc,xIc}.

Let 0 ≤ iN. We call Fi the space of channels which can be written as a product of functions of FI with the order of I at most i:

Ei:=clkK(X;Y)|ϕIFI,ϕ0F(X)|k=ϕ0IϕI;|I|i,

where cl denotes the closure in K(X;Y). Intuitively, this means that Ei does not only contain terms in the form given in the curly brackets, but also limits of such terms. Stated differently, the closure of a set includes not only the set itself, but also its boundary. This is important because when we project to a family, the projection may lie on the boundary. In order for the result to exist, the boundary must then be included.

This way:

E0 is the space of constant channels;

E N is the whole K(X;Y);

EiEj if and only if ij;

• For N ≤ 3 we recover exactly the quantities of Section 2.

The family Ei is also the closure of the family in the form:

1Z(X)expIfI(X;Y)|fIFI;|I|i,

where

Z(x):=yexpIfI(x;y).

Such families are known in the literature as exponential families (see, for example, Amari and Nagaoka (2000)). In particular, it is compact (for finite N), so that the infimum of any function on Ei is always a minimum. This means that for a channel k and an input distribution p:

Dp(k||Ei):=infmEiDp(k||m)=minmEiDp(k||m)

always exists. If it is unique, for example if p is strictly positive, we define the unique KL-projection as

πEik:=arg minmEiDp(k||m).

πEik has the property that it defines the same output probability on Y.

Definition 4: Let kK(X;Y), let 1 ≤ i ≤ N. Then the i-wise synergy of k is (if the KL-projections are unique):

di(k):=Dp(πEik||πEi1k).

For more clarity, we call the 1-wise synergy “single node information” or “single-node dependence.”

For kK(X;Y) = EN, we can look at its divergence from E0. If we denote πE0k by k0:

Dp(k||E0)=Dp(k||k0)=x,yp(x)k(x;y)logk(x;y)k0(y).

If k is not strictly positive, we take the convention 0 log(0/0) = 0, and we discard the zero terms from the sum. Since the output distributions are the same, but k0 is constant in x, it can not happen that for some (x;y), k0(x;y) = 0 but k(x;y) ≠ 0. (The very same is true for all KL-projections πEik, since Dp(πEik||k0)Dp(k||k0).) For all other terms, Equation (48) becomes:

Dp(k||E0)=x,yp(x)k(x;y)log k(x;y)x,yp(x)k(x;y)logk0(y)
=Hpk(Y|X)ykp(y)logk0(y)
=Hpk(Y|X)yk0p(y)logk0(y)
=Hpk(Y|X)+Hk0p(Y)=Hpk(Y|X)+Hkp(Y)
=Ipk(X:Y).

On the other hand, the Pythagorean relation (22) implies:

Dp(k||k0)=Dp(k||πEN1k)+Dp(πEN1k||k0),

and iterating:

Dp(k||k0)=Dp(k||πEN1k)+Dp(πEN1k||πEN2k)++Dp(πE1k||k0).

In the end, we get

I(X:Y)=i=1NDp(πEik||πEi1k)=i=1Ndi(k).

This decomposition is always non-negative, and it depends on the input distribution. The terms in Equation (56) can be in general difficult to compute exactly. Nevertheless, they can be approximated with iterative procedures.

5. Comparison with Two Recent Approaches

The measure of synergy, or respectively complementary information, defined in Griffith and Koch (2014) and Bertschinger et al. (2014) is

CIp(Y:X1,X2):=Ip(Y:X1,X2)minpIp(Y:X1,X2),

where ∧ is the space of prescribed marginals:

={qP(X1,X2,Y)|q(X1,Y)=p(X1,Y),q(X2,Y)=p(Y2,Y)}.

Our measure of synergy can be written, for two inputs, in a similar form:

d2(k)=Dp(k||πE1k)=Ip(Y:X1,X2)minpIp(Y:X1,X2),

where △, in addition to the constraints of ∧, prescribes also the input:

={qP(X1,X2,Y)|q(X1,Y)=p(X1,Y),q(X2,Y)=p(Y2,Y),q(X1,X2)=p(X1,X2)}.

Clearly △ ⊆ ∧, so

minpIp(Y:X1,X2)minpIp(Y:X1,X2),

which implies that

d2(k)CIp(Y:X1,X2).

We argue that not prescribing the input leads to overestimating synergy, because the subtraction in Equation (57) includes a possible difference in the correlation of the input distributions.

For example, consider X1, X2, Y binary and correlated, but not perfectly correlated. (For perfectly correlated nodes, as in Section 3, △ = ∧, so there is no difference between the two measures.) In detail, consider the channel:

kβ(x1,x2;y):=expβy(x1+x2)yexpβy(x1+x2),

and the input distribution:

pα(x1,x2):=expαx1x2x1,x2expαx1x2.

For α, β → ∞, the correlation becomes perfect, and the two measures of synergy are both zero. For 0 < α, β < ∞, our measure d2(kβ) is zero, as clearly kβE1. CI is more difficult to compute, but we can give a (non-zero) lower bound in the following way. First, we fix two values β = β0, α = α0. We consider the joint distribution pα0kβ0, and look at the marginals:

pα0kβ0(X1,Y),pα0kβ0(X2,Y).

We define the family ∧ as the set of joint probabilities which have exactly these marginals. If we increase β, we can always find an α such that the marginals do not change:

pαkβ(X1,Y)=pα0kβ0(X1,Y),pαkβ(X2,Y)=pα0kβ0(X2,Y),

i.e., such that pα kβ ∈ ∧. Now we can look at the mutual information of pα kβ and of pα0kβ0. If they differ, and (for example) the former is larger, then

Ipαkβ(Y:X1,X2)Ipα0kβ0(Y:X1,X2)Ipαkβ(Y:X1,X2)minpIp(Y:X1,X2)=CIpαkβ

is a well-defined lower bound for CIpαkβ. With a numerical simulation, we can show graphically that the mutual information is indeed not constant within the families ∧ (see Figure 11).

FIGURE 11
www.frontiersin.org

Figure 11. Mutual information and fixed marginals. The shades of blue represent the amount of Ip(Y :X1,X2) as a function of α,β (brighter is higher). Each red line represents a family ∧ of fixed marginals. While the lines of fixed mutual information and the families of fixed marginals look qualitatively similar, they do not coincide exactly, which means that Ip varies within the ∧.

From the picture, we can see that the red lines (families ∧ for different initial values) approximate well the lines of constant mutual information, at least qualitatively, but they are not exactly equal. This means that for most points p of ∧, the quantity:

CIp(Y:X1,X2):=Ip(Y:X1,X2)minpIp(Y:X1,X2)

will be non-zero. More explicitly, we can plot the increase in mutual information as p varies in ∧, for example, as a function of β (see Figure 12). This is always larger than or equal to the difference between the mutual information and its minimum in ∧ (i.e., CI). We can see that it is positive, which implies that CIp is also positive.

FIGURE 12
www.frontiersin.org

Figure 12. Lower bound for CI versus β. For each β ∈ [0.7, 3] we can find an α such that the joint pα kβ lies in ∧. The increase in mutual information as β varies is a lower bound for CI, which is then in general non-zero.

We can see in Figure 11 that especially for very large or very small values of α and β (i.e., very strong or very weak correlation), CI captures the behavior of mutual information quite well. These limits are precisely deterministic and constant kernels, for which most approaches in quantifying synergy coincide. This is the reason why the examples studied in Griffith and Koch (2014) give quite a satisfying result for CI (in their notation, SVK). For the less studied (and computationally more complex) intermediate values, like 1 < α, β < 2, the approximation is instead far from accurate, and in that interval (see Figure 12), there is a sharp increase in I, which leads to overestimating synergy.

6. Conclusion

Using information geometry, we have defined a non-negative decomposition of the mutual information between inputs and outputs of a channel.

The decomposition divides the mutual information into contributions of the different orders of synergy in a unique way. It does not, however, divide the mutual information into contributions of the different subsets of input nodes as Williams and Beer’s PID (Williams and Beer, 2010) would require.

For two inputs, our measure of synergy is closely related to the well-received quantification of synergy in Griffith and Koch (2014) and Bertschinger et al. (2014). Our measure though works in the desired way for an arbitrary (finite) number of inputs. Differently from Griffith and Koch (2014) and Bertschinger et al. (2014), anyway, we do not define a measure for redundant or “shared” information, nor unique information of the single inputs or subsets.

The decomposition depends on the choice of an input distribution. In case of input correlation, redundant information is counted automatically only once. This way, there is no need to quantify redundancy separately.

In general, there is no way to compute our quantities in closed form, but they can be approximated by an iterative-scaling algorithm (see, for example, Csiszár and Shields (2004)). The results are consistent with the intuitive properties of synergy, outlined in Williams and Beer (2010) and Griffith and Koch (2014).

Author Contributions

The research was initiated by NA. It was carried out by both authors, with main contribution by PP, who also wrote the article. Both authors read and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The handling editor Mikhail Prokopenko declares that despite having collaborated with author Nihat Ay, the review process was handled objectively and no conflict of interest exists.

Footnote

  1. ^For an example in which I = 0 but the nodes are not independent, see Williams and Beer (2010).

References

Amari, S. (2001). Information geometry on a hierarchy of probability distributions. IEEE Trans. Inf. Theory 47, 1701–1709. doi: 10.1109/18.930911

CrossRef Full Text | Google Scholar

Amari, S., and Nagaoka, H. (2000). Methods of Information Geometry. Translations of Mathematical Monographs, vol. 191. Oxford: American Mathematical Society.

Google Scholar

Ay, N. (2015). Information geometry on complexity and stochastic interaction. Entropy 17, 2432. doi:10.3390/e17042432

CrossRef Full Text | Google Scholar

Bertschinger, N., Rauh, J., Olbrich, E., Jost, J., and Ay, N. (2014). Quantifying unique information. Entropy 16, 2161. doi:10.3390/e16042161

CrossRef Full Text | Google Scholar

Csiszár, I., and Shields, P. C. (2004). Information theory and statistics: a tutorial. Found. Trends Commun. Inf. Theory 1, 417–528. doi:10.1561/0100000004

CrossRef Full Text | Google Scholar

Griffith, V., and Koch, C. (2014). “Quantifying synergistic mutual information,” in Guided Self-Organization: Inception, Emergence, Complexity and Computation Series, vol. 9. Springer.

Google Scholar

Harder, M., Salge, C., and Polani, D. (2013). Bivariate measure of redundant information. Phys. Rev. E 87, 012130. doi:10.1103/PhysRevE.87.012130

PubMed Abstract | CrossRef Full Text | Google Scholar

Kakihara, Y. (1999). Abstract Methods in Information Theory. Series on Multivariant Analysis, vol. 4. Singapore: World Scientific.

Google Scholar

McGill, W. L. (1954). Multivariate information transmission. Psychometrika 19, 97–116. doi:10.1007/BF02289159

CrossRef Full Text | Google Scholar

Rauh, J., Bertschinger, N., Olbrich, E., and Jost, J. (2014). “Reconsidering unique information: towards a multivariate information decomposition,” in International Symposium on Information Theory (ISIT). IEEE.

Google Scholar

Schneidmann, E., Bialek, W., and Berry, M. J. II (2003a). Synergy, redundancy, and independence in population codes. J. Neurosci. 23, 11539–11553. doi:10.1523/JNEUROSCI.5319-04.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneidmann, E., Still, S., Berry, M. J. II, and Bialek, W. (2003b). Network information and connected correlations. Phys. Rev. Lett. 91, 238701–238704. doi:10.1103/PhysRevLett.91.238701

PubMed Abstract | CrossRef Full Text | Google Scholar

Williams, P. L. and Beer, R. D. (2010). Nonnegative decomposition of multivariate information. Preprint available on arXiv:1004.2151.

Google Scholar

Keywords: synergy, redundancy, hierarchy, projections, divergences, interactions, iterative scaling, information geometry

Citation: Perrone P and Ay N (2016) Hierarchical Quantification of Synergy in Channels. Front. Robot. AI 2:35. doi: 10.3389/frobt.2015.00035

Received: 14 October 2015; Accepted: 08 December 2015;
Published: 08 January 2016

Edited by:

Mikhail Prokopenko, University of Sydney, Australia

Reviewed by:

Daniele Marinazzo, University of Ghent, Belgium
Rick Quax, University of Amsterdam, Netherlands

Copyright: © 2016 Perrone and Ay. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Paolo Perrone, perrone@mis.mpg.de