Backward transfer entropy: Informational measure for detecting hidden Markov models and its interpretations in thermodynamics, gambling and causality.

The transfer entropy is a well-established measure of information flow, which quantifies directed influence between two stochastic time series and has been shown to be useful in a variety fields of science. Here we introduce the transfer entropy of the backward time series called the backward transfer entropy, and show that the backward transfer entropy quantifies how far it is from dynamics to a hidden Markov model. Furthermore, we discuss physical interpretations of the backward transfer entropy in completely different settings of thermodynamics for information processing and the gambling with side information. In both settings of thermodynamics and the gambling, the backward transfer entropy characterizes a possible loss of some benefit, where the conventional transfer entropy characterizes a possible benefit. Our result implies the deep connection between thermodynamics and the gambling in the presence of information flow, and that the backward transfer entropy would be useful as a novel measure of information flow in nonequilibrium thermodynamics, biochemical sciences, economics and statistics.

In many scientific problems, we consider directed influence between two component parts of complex system. To extract meaningful influence between component parts, the methods of time series analysis have been widely used123. Especially, time series analysis based on information theory4 provides useful methods for detecting the directed influence between component parts. For example, the transfer entropy (TE)567 is one of the most influential informational methods to detect directed influence between two stochastic time series. The main idea behind TE is that, by conditioning on the history of one time series, informational measure of correlation between two time series represents the information flow that is actually transferred at the present time. Transfer entropy has been well adopted in a variety of research areas such as economics8, neural networks91011, biochemical physics121314 and statistical physics1516171819. Several efforts to improve the measure of TE have also been done202122.

In a variety of fields, a similar concept of TE has been discussed for a long time. In economics, the statistical hypothesis test called as the Granger causality (GC) has been used to detect the causal relationship between two time series2324. Indeed, for Gaussian variables, the statement of GC is equivalent to TE25. In information theory, nearly the same informational measure of information flow called the directed information (DI)2627 has been discussed as a fundamental bound of the noisy channel coding under causal feedback loop. As in the case of GC, DI can be applied to an economic situation2829, that is the gambling with side information430.

In recent studies of a thermodynamic model implementing the Maxwell’s demon3132, which reduces the entropy change in a small subsystem by using information, TE has attracted much attention13141518333435363738. In this context, TE from a small subsystem to other systems generally gives a lower bound of the entropy change in a subsystem151833. As a tighter bound of the entropy change for Markov jump process, another directed informational measure called the dynamic information flow (DIF)34 has also been discussed3334353637383940414243.

In this article, we provide the unified perspective on different measures of information flow, i.e., TE, DI, and DIF. To introduce TE for backward time series1338, called backward transfer entropy (BTE), we clarify the relationship between these informational measures. By considering BTE, we also obtain a tighter bound of the entropy change in a small subsystem even for non Markov process. In the context of time series analysis, this BTE has a proper meaning: an informational measure for detecting a hidden Markov model. From the view point of the statistical hypothesis test, BTE quantifies an anti-causal prediction. These fact implies that BTE would be a useful directed measure of information flow as well as TE.

Furthermore, we also discuss the analogy between thermodynamics for a small system324445 and the gambling with side information430. To considering its analogy, we found that TE and BTE play similar roles in both settings of thermodynamics and gambling: BTE quantifies a loss of some benefit while TE quantifies some benefit. Our result reveals the deep connection between two different fields of science, thermodynamics and gambling.

Results

Setting

We consider stochastic dynamics of interacting systems and , which are not necessarily Markov processes. We consider a discrete time k (=1, …, N), and write the state of () at time k as x (y). Let () be the path of system () from time k − l + 1 to k where l ≥ 1 is the length of the path. The probability distribution of the composite system at time k is represented by p(X = x, Y = y), and that of paths is represented by , where capital letters (e.g., X) represent random variables of its states (e.g., x).

The dynamics of composite system are characterized by the conditional probability such that

where p(A = a|B = b):= p(A = a, B = b)/p(B = b) is the conditional probability of a under the condition of b.

Transfer entropy

Here, we introduce conventional TE as a measure of directed information flow, which is defined as the conditional mutual information4 between two time series under the condition of the one’s past. The mutual information characterizes the static correlation between two systems. The mutual information between X and Y at time k is defined as

This mutual information is nonnegative quantity, and vanishes if and only if x and y are statistically independent (i.e., p(X = x, Y = y) = p(X = x)p(Y = y))4. This mutual information quantifies how much the state of y includes the information about x, or equivalently the state of x includes the information about y. In a same way, the mutual information between two paths and is also defined as

While the mutual information is very useful in a variety fields of science4, it only represents statistical correlation between two systems in a symmetric way. In order to characterize the directed information flow from X to Y, Schreiber5 introduced TE defined as

with k ≤ k′. Equation (4) implies that TE is an informational difference about the path of the system that is newly obtained by the path of the system from time k′ to k′ + 1. Thus, TE can be regarded as a directed information flow from to at time k′. This TE can be rewritten as the conditional mutual information4 between the paths of and the state of under the condition of the history of :

which implies that TE is nonnegative quantity, and vanishes if and only if the transition probability in from to y does not depend on the time series , i.e., ). [see also Fig. 1(a)].

Figure 1

Schematics of TE and BTE.

Two graphs (a) and (b) are the Bayesian networks corresponding to the joint probabilities p(X = x, Y = y, Y = y) = p(X = x)p(Y = y|X = x)p(Y = y|X = x, Y = y) and p(X = x, Y = y, Y = y) = p(Y = y)p(Y = y|Y = y)p(X = x|Y = y, Y = y), respectively (see also refs 15, 37 and 60). (a) Transfer entropy corresponds to the edge from X to Y on the Bayesian network. If TE is zero, the edge from X to Y vanishes, i.e., p(X = x, Y = y, Y = y) = p(X = x)p(Y = y|X = x)p(Y = y|Y = y). (b) Backward transfer entropy corresponds to the edge from Y to X on the Bayesian network. If BTE is zero, the edge from Y to X vanishes, i.e., p(X = x, Y = y, Y = y) = p(Y = y)p(Y = y|Y = y)p(X = x|Y = y).

Backward transfer entropy

Here, we introduce BTE as a novel usage of TE for the backward paths. We first consider the backward path of the system (); (), which is the time-reversed trajectories of the system () from time N − k + l to N − k + 1. We now introduce the concept of BTE defined as TE for the backward paths

with m = N − k, m′ = N − k′ and k ≤ k′. In this sense, BTE may represent “the time-reversed directed information flow from the future to the past.” However BTE is well defined as the conditional mutual information, it is nontrivial if such a concept makes any sense information-theoretically or physically where stochastic dynamics of composite system itself do not necessarily have the time-reversal symmetry.

To clarify the proper meaning of BTE, we compare BTE with TE [see Fig. 1]. Transfer entropy quantifies the dependence of X in the transition from time Y to Y [see Fig. 1(a)]. In the same way, BTE quantifies the dependence of Y in the correlation between X and Y [see Fig. 1(b)]. Thus, BTE implies how X depends on Y without the dependence of the past state Y. In other words, BTE is nonnegative and vanishes if and only if a Markov chain Y → Y → X exists, which implies that dynamics of X are given by a hidden Markov model. In general, BTE is nonnegative and vanishes if and only if a Markov chain

exists. Therefore, BTE from to quantifies how far it is from composite dynamics of and to a hidden Markov model in .

Thermodynamics of information

We next discuss a thermodynamic meaning of BTE. To clarify the interpretation of BTE in nonequilibrium stochastic thermodynamics, we consider the following non-Markovian interacting dynamics

where a nonnegative integer n represents the time delay between and . The stochastic entropy change in heat bath attached to the system from time 1 to N in the presence of 15 is defined as

where

is the transition probability of backward dynamics, which satisfies the normalization of the probability . For example, if the system and does not include any odd variable that changes its sign with the time-reversal transformation, the backward probability is given by pB(X = x|X = x, Y = y) = p(X = x|X = x, Y = y) with k ≥ n + 1 (pB(X = x|X = x, Y1 = y1) = p(X = x|X = x, Y1 = y1) with k ≤ n). This definition of the entropy change in the heat bath Eq. (9) is well known as the local detailed balance or the detailed fluctuation theorem45. We define the entropy change in and heat bath as

where is the stochastic Shannon entropy change in .

For the non-Markovian interacting dynamics Eq. (8), we have the following inequality (see Method);

We add that the term vanishes for the Markovian interacting dynamics (n = 0).

These results [Eqs (12) and (13)] can be interpreted as a generalized second law of thermodynamics for the subsystem in the presence of information flow from to . If there is no interaction between and , informational terms vanish, i.e., , , , , I(X;Y) = 0 and I(X1; Y1) = 0. Thus these results reproduce the conventional second law of thermodynamics , which indicates the nonnegativity of the entropy change in and bath45. If there is some interaction between and , can be negative, and its lower bound is given by the sum of TE from X to Y and mutual information between and at initial time;

which is a nonnegative quantity . In information theory, this quantity is known as DI from to 27. Intuitively speaking, quantifies a kind of thermodynamic benefit because its negativity is related to the work extraction in in the presence of 32. Thus, a weaker bound (13) implies that the sum of TE quantifies a possible thermodynamic benefit of in the presence of .

We next consider the sum of TE for the time-reversed trajectories;

which is given by the sum of BTE and the mutual information between and at final time. A tighter bound (12) can be rewritten as the difference between the sum of TE and BTE;

This result implies that a possible benefit should be reduced by up to the sum of BTE . Thus, the sum of BTE means a loss of thermodynamic benefit. We add that a tighter bound is not necessarily nonnegative while a weaker bound is nonnegative.

We here consider the case of Markovian interacting dynamics (n = 0). For Markovian interacting dynamics, we have the following additivity for a tighter bound [see Supplementary information (SI)]

where the sum of TE and BTE for a single time step and are defined as and , respectively. This additivity implies that a tighter bound for multi time steps is equivalent to the sum of a tighter bound for a single time step . We stress that a tighter bound for a single time step has been derived in ref. 13. We next consider the continuous limit x = x(t = kΔt), y = y(t = kΔt), and N = O(Δt−1), where t denotes continuous time, Δt ≪ 1 is an infinitesimal time interval and the symbol O is the Landau notation. Here we clarify the relationship between a tighter bound (16) and DIF34 (or the learning rate18) defined as . For the bipartite Markov jump process18 or two dimensional Langevin dynamics without any correlation between thermal noises in X and Y15, we have the following relationship [see also SI]

Thus a bound by TE and BTE is equivalent to a bound by DIF for such systems in the continuous limit, i.e., .

Gambling with side information

In classical information theory, the formalism of the gambling with side information has been well known as another perspective of information theory based on the data compression over a noisy communication channel430. In the gambling with side information, the mutual information between the result in the gambling and the side information gives a bound of the gambler’s benefit.

This formalism of gambling is similar to the above-mentioned result in thermodynamics of information. In thermodynamics, thermodynamic benefit (e.g., the work extraction) can be obtained by using information. On the other hand, the gambler obtain the benefit by using side information. We here clarify the analogy between gambling and thermodynamics in the presence of information flow. To clarify the analogy between thermodynamics and gambling, BTE plays a crucial role as well as TE.

We introduce the basic concept of the gambling with side information given by the horse race430. Let y be the horse that won the k-th horse race. Let f ≥ 0 and o ≥ 0 be the bet fraction and the odds on the k-th race, respectively. Let m be the gambler’s wealth before the k-th race. Let s be the side information at time k. We consider the set of side information x = {s1, …, s}, which the gambler can access before the k-th race. The bet fraction f is given by the function with k ≥ 2, and f1(y1|x1). The conditional dependence ({x1}) of (f(y1|x1)) implies that the gambler can decide the bet fraction f (f1) by considering the past information ({x1}). We assume normalizations of the bet fractions and , which mean that the gambler bets all one’s money in every race. We also assume that . This condition satisfies if the odds in every race are fair, i.e., 1/o(y) is given by a probability of Y.

The stochastic gambler’s wealth growth rate at k-th race is given by

with k ≥ 2 [], which implies that the gambler’s wealth stochastically changes due to the bet fraction and odds. The information theory of the gambling with side information indicates that the ensemble average of total wealth growth is bounded by the sum of TE (or DI) from X to Y2829 (see Method);

where indicates the ensemble average, and is the Shannon entropy of . This result (21) implies that the sum of TE can be interpreted as a possible benefit of the gambler.

We discuss the analogy between thermodynamics of information and the gambling with side information. A weaker bound in the gambling with side information (21) is similar to a weaker bound in thermodynamics of information (16), where the negative entropy change corresponds to the total wealth growth G. On the other hand, a tighter bound in the gambling with side information (20) is rather different from a tighter bound by the sum of BTE in thermodynamics of information (16). We show that a tighter bound in the gambling is also given by the sum of BTE if we consider the special case that the bookmaker who decides the odds o cheats in the horse race; The odds o can be decided by the unaccessible side information x and information of the future races [see also Fig. 2]. In this special case, the fair odds of the k-th race o can be the conditional probability of the future information with k ≤ N − 1, and 1/o(y) = p(Y = y|X = x). The inequality (20) can be rewritten as

Figure 2

Schematic of the special case of the horse race.

The gambler can only access the past side information x and the past races , and decides the bet fraction f on the k-th race. The bookmaker makes some cheating which can access the future side information x and the future races , and decides the odds on the k-th race.

which implies that the sum of BTE ) represents a loss of the gambler’s benefit because of the cheating by the bookmaker who can access the future information with anti-causality. We stress that Eq. (22) has a same form of the thermodynamic inequality (16) for Markovian interacting dynamics (n = 0). This fact implies that thermodynamics of information can be interpreted as the special case of the gambling with side information; The gambler uses the past information and the bookmaker uses the future information. If we regard thermodynamic dynamics as the gambling, anti-causal effect should be considered.

Causality

We here show that BTE itself is related to anti-causality without considering the gambling. From the view point of the statistical hypothesis test, TE is equivalent to GC for Gaussian variables25. Therefore, it is naturally expected that BTE can be interpreted as a kind of the causality test.

Suppose that we consider two linear regression models

where α (α′) is a constant term, A (A′) is the vector of regression coefficients, ⊕ denotes concatenation of vectors, and () is an error term. The Granger causality of to quantifies how the past time series of in the first model reduces the prediction error of compared to the error in the second model. Performing ordinary mean squares to find the regression coefficients A (A′) and α (α′) that minimize the variance of (), the standard measure of GC is given by

where var() denotes the variance of . Here we assume that the joint probability is Gaussian. Under Gaussian assumption, TE and GC are equivalent up to a factor of 2,

In the same way, we discuss BTE from the view point of GC. Here we assume that the joint probability is Gaussian. Suppose that two linear regression models

where α† (α′†) is a constant term, A† (A′†) is the vector of regression coefficients and () is an error term. These linear regression models give a prediction of the past state of using the future time series of and . Intuitively speaking, we consider GC of to for the rewind playback video of composite dynamics and . We call this causality test the Granger anti-causality of to . Performing ordinary mean squares to find A† (A′†) and α† (α′†) that minimize var() (var()), we define a measure of the Granger anti-causality of to as . The backward transfer entropy is equivalent to the Granger anti-causality up to factor 2,

This fact implies that BTE can be interpreted as a kind of anti-causality test. We stress that composite dynamics of and are not necessarily driven with anti-causality even if a measure of the Granger anti-causality has nonzero value. As GC just finds only the predictive causality2324, the Granger anti-causality also finds only the predictive causality for the backward time series.

Discussion

We proposed that directed measure of information called BTE, which is possibly useful to detect a hidden Markov model (7) and predictive anti-causality (29). In the both setting of thermodynamics and the gambling, the measurement of BTE has a profitable meaning; the detection of a loss of a possible benefit in the inequalities (16) and (22).

The concept of BTE can provide a clear perspective in the studies of the biochemical sensor and thermodynamics of information, because the difference between TE and DIF has attracted attention recently in these fields1435. In ref. 14, Hartich et al. have proposed the novel informational measure for the biochemical sensor called sensory capacity. The sensory capacity is defined as the ratio between TE and DIF . Because DIF can be rewritten by TE and BTE [Eq. (18)] for Markovian interacting dynamics, we have the following expression for the sensory capacity in a stationary state,

where we used I(X; Y) = I(X; Y) in a stationary state. This fact indicates that the ratio between TE and BTE could be useful to quantify the performance of the biochemical sensor. By using this expression (29), we show that the maximum value of the sensory capacity C = 1 can be achieved if a Markov chain of a hidden Markov model Y → Y → X exists. In ref. 35, Horowitz and Sandberg have shown a comparison between two thermodynamic bound by TE and DIF for two dimensional Langevin dynamics. For the Kalman-Bucy filter which is the optimal controller, they have found the fact that DIF is equivalent to TE in a stationary state. This idea can be clarified by the concept of BTE. Because the Kalman-Bucy filter can be interpreted as a hidden Markov model, BTE should be zero, and DIF is equivalent to TE in a stationary state.

Our results can be interpreted as a generalization of previous works in thermodynamics of information464748. In refs 46 and 47, S. Still et al. discuss the prediction in thermodynamics for Markovian interacting dynamics. In our results, we show the connection between thermodynamics of information and the predictive causality from the view point of GC. Thus, our results give a new insight into these works of the prediction in thermodynamics. In ref. 48, G. Diana and M. Esposito have introduced the time-reversed mutual information for Markovian interacting dynamics. In our results, we introduce BTE, which is TE in the time-reversed way. Thus, our result provides a similar description of thermodynamics by introducing BTE, even for non-Markovian interacting dynamics.

We point out the time symmetry in the generalized second law (12). For Markovian interacting dynamics, the equality in Eq. (12) holds if dynamics of has a local reversibility (see SI). Here we consider a time reversed transformation , and assume a local reversibility such that the backward probability pB(A = a|B = b) equals to the original probability p(A = a|B = b) for any random variables A and B. In a time reversed transformation, we have , and . The generalized second law Eq. (12) changes the sign in a time reversed transformation, . Thus, the generalized second law (12) has the same time symmetry in the conventional second law, i.e., even for non-Markovian interacting dynamics, where ΔStot is the entropy change in total systems. In other words, the generalized second law (12) provides the arrow of time as the conventional second law. This fact may indicate that BTE is useful as well as TE in physical situations where the time symmetry plays a crucial role in physical laws.

We also point out that this paper clarifies the analogy between thermodynamics of information and the gambling. The analogy between the gambling and thermodynamics has been proposed in ref. 49, however, the analogy between Eqs (16) and (22) are different from one in ref. 49. In ref. 49, D. A. Vinkler et al. discuss the particular case of the work extraction in Szilard engine, and consider the work extraction in Szilard engine as the gambling. On the other hand, our result provides the analogy between the general law of thermodynamics of information and the gambling. To clarify this analogy, we may apply the theory of gambling, for example the portfolio theory5051, to thermodynamic situations in general. We also stress that the gambling with side information directly connects with the data compression in information theory4. Therefore, the generalized second law of thermodynamics may directly connect with the data compression in information theory. To consider such applications, BTE would play a tricky role in the theory of the gambling where the odds should be decided with anti-causality.

Finally, we discuss the usage of BTE in time series analysis. In principle, we prepare the backward time series data from the original time series data, and do a calculation of BTE as TE. To calculate BTE, we can estimate how far it is from dynamics of two time series to a hidden Markov model, or detect the predictive causality for the backward time series. In physical situations, we also can detect thermodynamic performance by comparing BTE with TE. If the sum of BTE from the target system to the other systems is larger than the sum of TE from the target system to the other systems, the target system could seem to violate the second law of thermodynamics because of the inequality (16), where the other systems play a similar role of Maxwell’s demon. Therefore, BTE could be useful to detect phenomena of Maxwell’s demon in several settings such as Brownian particles5253, electric devices5455, and biochemical networks135657585960.

Method

The outline of the derivation of inequality (12)

We here show the outline of the derivation of the generalized second law (12) [see also SI for details]. In SI, we show that the quantity , can be rewritten as the Kullbuck-Leiber divergence 4, where and , are nonnegative functions that satisfy the normalizations and . Due to the nonnegativity of the Kullbuck-Leiber divergence, we obtain the inequality (12), i.e., . We add that the integrated fluctuation theorem corresponding to the inequality (12) is also valid, i.e., .

The outline of the derivation of inequality (20)

We here show the outline of the derivation of the gambling inequality (20) [see also SI for details]. The quantity can be rewritten as the Kullbuck-Leiber divergence DK(ρ || π), where and , are nonnegative functions that satisfy the normalizations and . Due to the nonnegativity of the Kullbuck-Leiber divergence, we have the inequality (20), i.e., .

Additional Information

How to cite this article: Ito, S. Backward transfer entropy: Informational measure for detecting hidden Markov models and its interpretations in thermodynamics, gambling and causality. Sci. Rep. 6, 36831; doi: 10.1038/srep36831 (2016).

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.