arxiv: v1 [math.st] 23 Apr 2008

Transcription

1 arxiv: v1 [math.st] 23 Apr 2008 Causal inference using the algorithmic Markov condition Dominik Janzing and Bernhard Schölkopf Max Planck Institute for Biological Cybernetics, Spemannstr. 38, Tübingen, Germany April 23, 2008 Abstract Inferring the causal structure that links n observables is usually based upon detecting statistical dependences and choosing simple graphs that make the joint measure Markovian. Here we argue why causal inference is also possible when only single observations are present. We develop a theory how to generate causal graphs explaining similarities between single objects. To this end, we replace the notion of conditional stochastic independence in the causal Markov condition with the vanishing of conditional algorithmic mutual information and describe the corresponding causal inference rules. We explain why a consistent reformulation of causal inference in terms of algorithmic complexity implies a new inference principle that takes into account also the complexity of conditional probability densities, making it possible to select among Markov equivalent causal graphs. This insight provides a theoretical foundation of a heuristic principle proposed in earlier work. We also discuss how to replace Kolmogorov complexity with decidable complexity criteria. This can be seen as an algorithmic analog of replacing the empirically undecidable question of statistical independence with practical independence tests that are based on implicit or explicit assumptions on the underlying distribution. {dominik.janzing,bernhard.schoelkopf}@tuebingen.mpg.de 1

2 Contents 1 Introduction to causal inference from statistical data Causal Markov condition Seeking for new statistical inference rules Inferring causal relations among individual objects Algorithmic mutual information Markov condition for algorithmic dependences among individual objects Relative causality Novel statistical inference rules from the algorithmic Markov condition Algorithmic independence of Markov kernels Resolving statistical ensembles into individual observations Conditional density estimation on subsamples Plausible Markov kernels in time series Decidable modifications of the inference rule Causal inference using symmetry constraints Resource-bounded complexity Conclusions 44 1 Introduction to causal inference from statistical data Causal inference from statistical data has attracted increasing interest in the past decade. In contrast to traditional statistics where statistical dependences are only taken to prove that some kind of relation between random variables exists, causal inference methods in machine learning are explicitly designed to generate hypotheses on causal directions automatically based upon statistical independence tests [1, 2]. The crucial assumption connecting statistics with causality is the causal Markov condition explained below after we have introduced some notations and terminology. We denote random variables by capitals and their values by the corresponding lowercase letters. Let X 1,...,X n be random variables and G be a directed acyclic graph (DAG) representing the causal structure where an arrow from node X i to node X j indicates a direct causal effect. Here the term direct is understood with respect to the chosen set of variables in the sense that the information flow between the two variables considered is not performed via using one or more of the other variables as intermediate nodes. We will next briefly rephrase the postulates that are required in the statistical theory of inferred causation [2, 1]. 1.1 Causal Markov condition When we consider the causal structure that links n random variables V := {X 1,...,X n } we will implicitly assume that V is causally sufficient in the sense that all common causes of two variables in V are also in V. Then a causal hypothesis G is only acceptable as potential causal structure if the joint distribution P of X 1,...,X n satisfies the Markov condition with respect to G. There are several formulations of the Markov condition that are known to coincide under some technical condition (see Lemma 1). We will first introduce the following version which is sometimes referred to as the parental or the local Markov condition [3]. 2

3 To this end, we introduce the following notations. PA j is the set of parents of X j and ND j the set of non-descendants of X j except itself. If S,T,R are sets of random variables, S T R means S is statistically independent of T, given R. Postulate 1 (statistical causal Markov condition, local) If adirected acyclic graph G formalizes the causal structureamong the random variables X 1,...,X n. Then X j ND j PA j, for all j = 1,...,n. We call this postulate the statistical causal Markov condition because we will later introduce an algorithmic version. The fact that conditional irrelevance not only occurs in the context of statistical dependences has been emphasized in the literature (e.g. [4, 1]) in the context of describing abstract properties (like semi-graphoid axioms) of the relation. We will therefore state the causal Markov condition also in an abstract form that does not refer to any specific notion of conditional informational irrelevance: Postulate 2 (abstract causal Markov condition, local) Given all the direct causes of an observable O, its non-effects provide no additional information on O. Here, observables denote something in the real world that can be observed and the observation of which can be formalized in terms of a mathematical language. In this paper, observables will either be random variables (formalizing statistical quantities) or they will be strings (formalizing the description of objects). Accordingly, information will be statistical or algorithmic mutual information, respectively. The importance of the causal Markov condition lies in the fact that it links causal terms like direct causes and non-effects to informational relevance of observables. The local Markov condition is rather intuitive because it echoes the fact that the information flows from direct causes to their effect and every dependence between a node and its non-descendants involves the direct causes. However, the independences postulated by the local Markov condition imply additional independences. It is therefore hard to decide whether an independence must hold for a Markovian distribution or not, solely on the basis of the local formulation. In contrast, the global Markov condition makes the complete set of independences obvious. To state it we first have to introduce the following graph-theoretical concept. Definition 1 (d-separation) A path p in a DAG is said to be d-separated (or blocked) by a set of nodes Z if and only if 1. p contains a chain i m j or fork i m j such that the middle node m is in Z, or 2. p contains an inverted fork (or collider) i m j such that the middle node m is not in Z and such that no descendant of m is in Z. A set Z is said to d-separate X from Y if and only if Z blocks every (possibly undirected) path from a node in X to a node in Y. The following Lemma shows that d-separation is the correct condition for deciding whether an independence is implied by the local Markov condition [4], Theorem Lemma 1 (equivalent Markov conditions) Let P(X 1,...,X n ) have a density P(x 1,...,x n ) with respect to a product measure. Then the following three statements are equivalent: 3

4 I. Recursive form: P admits the factorization P(x 1,...,x n ) = n P(x j pa j ), (1) where P(. pa j ) is shorthand for the conditional probability density, given the values of all parents of X j. II. Local (or parental) Markov condition: for every node X j we have j=1 X j ND j PA j, i.e., it is conditionally independent of its non-descendants (except itself), given its parents. III. Global Markov condition: S T R for all three sets S,T,R of nodes for which S and T are d-separated by R. Moreover, the local and the global Markov condition are equivalent even if P does not have a density with respect to a product measure. The conditional densities P(x j pa j ) are also called the Markov kernels relative to the hypothetical causal graph G. It is important to note that every choice of Markov kernels define a Markovian density P, i.e., the Markov kernels define exactly the set of free parameters remaining after the causal structure has been specified. To select graphs among all those that render P Markovian, we also need an additional postulate: Postulate 3 (causal faithfulness) Among all graphs G for which P is Markovian, prefer the ones for which all the observed conditional independences in the joint measure P(X 1,...,X n ) are imposed by the Markov condition. The idea is that the set of observed independences is typical for the causal structure under consideration rather than being the result of specific choices of the Markov kernels. This becomes even more intuitive when we restrict our attention to random variables with finite value set and observe that the values P(x j pa j ) then define a natural parameterization of the set of Markovian distributions in a finite dimensional space. The non-faithful distributions form a submanifold of lower dimension, i.e., a set of Lebesgue measure zero [5]. They therefore almost surely don t occur if we assume that nature chooses the Markov kernels for the different nodes independently according to some density on the parameter space. The above zero Lebesgue measure argument is close to the spirit of Bayesian approaches [6], where priors on the set of Markov kernels are specified for every possible hypothetical causal DAG and causal inference is performed by maximizing posterior probabilities for hypothetical DAGs, given the observed data. This procedure leads to an implicit preference of faithful structures in the infinite sampling limit given some natural conditions for the priors on the parameter space. The assumption that nature chooses Markov kernels independently, which is also part of the Bayesian approach, will turn out to be closely related to the algorithmic Markov condition postulated in this paper. We now discuss the justification of the statistical causal Markov condition because we will later justify the algorithmic Markov condition in a similar way. To this end, we introduce functional models [1]: 4

5 Postulate 4 (functional model of causality) If a directed acyclic graph G formalizes the causal relation between the random variables X 1,...,X N then every X j can be written as a deterministic function of PA j and a noise variable N j, where all N j are jointly independent. Then we have [1], Theorem 1.4.1: X j = f j (PA j,n j ), Lemma 2 (Markov condition in functional models) Every joint distribution P(X 1,...,X n ) generated according to the functional model in Postulate 4 satisfies the local and the global Markov condition relative to G. We rephrase the proof in [1] because our proof for the algorithmic version will rely on the same idea. Proof of Lemma 2: extend G to a graph G with nodes X 1,...,X n,n 1,...,N n that additionally contains an arrow from each N j to X j. The given joint distribution of noise variables induces a joint distribution P(X 1,...,X n,n 1,...,N n ), that satisfies the local Markovcondition with respect to G: first, every X j is completely determined by its parents making the condition trivial. Second, every N j is parentless and thus we have to check that it is (unconditionally) independent of its non-descendants. The latter are deterministic functions of {N 1,...,N n }\{N j }. Hence the independence follows from the joint independence of all N i. By Lemma 1, P is also globally Markovian w.r.t. G. Then we observe that NDj and X j are d-separated in G (where the parents and non-descendants are defined with respect to G). Hence P satisfies the local Markov condition w.r.t. G and hence also the global Markov condition. Functional models formalize the idea that the outcome of an experiment is completely determined by the values of all relevant parameters where the only uncertainty stems from the fact that some of these parameters are hidden. Even though this kind of determinism is in contrast with the commonly accepted interpretation of quantum mechanics [7], we still consider functional models as a helpful framework for discussing causality in real life since quantum mechanical laws refer mainly to phenomena in micro-physics. Causal inference using the Markov condition and the faithfulness assumption has been implemented as causal learning algorithms [2]. The following fundamental limitations of these methods deserve our further attention: 1. Markov equivalence: There are only few cases where the inference rules provide unique causal graphs. Often one ends up with a large class of Markov equivalent graphs, i.e., graphs that entail the same set of independences. For this reason, additional inference rules are desirable. 2. Dependence on i.i.d. sampling: the whole setting of causal inference relies on the ability to sample repeatedly and independently from the same joint distribution P(X 1,...,X n ). As opposed to this assumption, causal inference in real life also deals with probability distributions that change in time and often one infers causal relations among single observations without referring to statistics at all. The idea of this paper is to develop a theory of probability-free causal inference that helps to construct causal hypotheses based on similarities of single objects. Here, similarities will be defined 5

6 by comparing the length of the shortest description of single objects to the length of their shortest joint description. Despite the analogy to causal inference from statistical data (which is due to known analogies between statistical and algorithmic information theory) our theory also implies new statistical inference rules. In other words, our approach to address weakness 2 also yields new methods to address 1. The paper is structured as follows. In the remaining part of this Section, i.e., Subsection 1.2, we describe recent approaches from the literature to causal inference from statistical data that address problem 1 above. In Section 2 we develop the general theory on inferring causal relations among individual objects based on algorithmic information. This framework appears, at first sight, as a straightforward adaption of the statistical framework (using well-known correspondences between statistical and algorithmic information theory). However, Section 3 describes that this implies novel causal inference rules for statistical inference because non-statistical algorithmic dependences can even occur in data that were obtained from statistical sampling. In Section 4 we describe how to replace causal inference rules based on the uncomputable algorithmic information with decidable criteria that are still motivated by the uncomputable idealization. The table in fig. 1 summarizes the analogies between the theory of statistical and the theory of algorithmic causal inference described in this paper. The differences, however, which are the main subject of Sections 3 to 4, can hardly be represented in the table. 1.2 Seeking for new statistical inference rules In[8]and[9]wehaveproposedcausalinferencerulesthatarebasedontheideathatthefactorization of P(cause,effect) into P(effect cause) and P(cause) typically leads to simpler terms than the artificial factorization into P(effect)P(cause effect). The generalization of this principle reads: Among all graphs G that render P Markovian prefer the one for which the decomposition in eq. (1) yields the simplest Markov kernels. We have called this vague idea the principle of plausible Markov kernels. Before we describe several options to define simplicity we describe a simple example to illustrate the idea. Assume we have observed that a binary variable X (with values x = 1, 1) and a continuous variable Y (with values in R) are distributed according to a mixture of two Gaussians (see fig. 2). Since this will simplify the further discussion let us assume that the two components are equally weighted, i.e., P(x,y) = 1 1 2σ (y µ xλ) 2 2π e 2σ 2, where λ determines the shift of the mean caused by switching between x = 1 and x = 1. The marginal P(Y) is given by P(y) = 1 ( ) 1 2σ 2 e (y µλ) 2σ 2 e (y µ λ)2 2σ 2. (2) 2π One will prefer the causal structure X Y compared to Y X because the former explains in a natural way why P(Y) is bimodal: the effect of X on Y is simply to shift the Gaussian distribution by 2λ. In the latter model the bimodality ofp(y) remains unexplained. To preferone causalmodel to another one because the corresponding conditionals are simpler seems to be a natural application of Occam s Razor. However, Section 3 will show that such an inference rule also follows from the theory developed in the present paper when simplicity is meant in the sense of low Kolmogorov complexity. In the remaining part of this section we will sketch some approaches to implement the principle of plausible Markov kernels in practical applications. In [8] we have defined a family of plausible Markov kernels by conditionals P(X j PA j ) that are second order exponential models, i.e., logp(x j pa j ) is a polynomial of order two in the variables 6

7 statistical algorithmic observables random variables sequences of strings (vertices of a DAG) observations i.i.d. sampled data strings conditional independence X Y Z x y z.... I(X;Y Z) = 0 I(x : y z) = 0 I. recursion formula P(x 1,...,x n ) K(x 1,...,x n ) = = j p(x j pa j ) j K(x j pa j ) II. local Markov condition X j ND j PA j x j nd j pa j III. global Markov d-separation d-separation condition statistical independence algorithmic independence equivalence of I-III Theorem 3.27 Theorem 3 in [4] functional models Section 1.4 Postulate 6 in [1] functional models Theorem Theorem 4 imply Markov condition in [1] decidable dependence assumptions on Section 4 criteria joint distribution Figure 1: Analogies between statistical and algorithmic causal inference 7

8 Figure 2: Observed joint distribution of X and Y consisting of two Gaussians of equal width shifted against each other. {X j } {PA j } up to some additive partition function (for normalization) that depends only on the variables PA j. For every hypothetical causal graph, one thus obtains a family of plausible joint distributions P(X 1,...,X n ) that are products of the plausible Markov kernels. Then we prefer the causal direction for which the plausible joint distributions provide the best fit for the given observations. In [9] we have proposed the following principle for causal inference: Given a joint distribution of the random variables X 1,...,X n, prefer a causal structure for which n C(P(X j PA j )) (3) j=1 is minimal, where C is some complexity measure on conditional probability densities. There is also another recent proposal for new inference rules that refers to a related simplicity assumption, though formally quite different from the ones above. The authors of [10] observe that there are joint distributions of X 1,...,X n that can be explained by a linear model with additive non-gaussian noise for one causal direction but require non-linear causal influence for the other causal directions. For real data they prefer the causal graph for which the observations are closer to the linear model. To justify the belief that conditionals that correspond to the true causal direction tend to be simpler than non-causal conditionals (which is common to all the approaches above) is one of the main goals of this paper. 2 Inferring causal relations among individual objects It has been emphasized [1] that the application of causal inference principles often benefits from the non-determinism of causal relations between the observed random variables. In contrast, human learning in real-life often is about quite deterministic relations. Apart from that, the most important difference between human causal learning and the inference rules in [2, 1] is that the former is also about causal relations among single objects and does not necessarily require sampling. Assume, for instance, that the comparison of two texts show similarities (see e.g. [11]) such that the author of the text that appeared later is blamed to have copied it from the other one or both are blamed to have copied from a third one. The statement that the texts are similar could be based on a statistical analysis of the occurrences of certain words or letter sequences. However, such kind of 8

9 simple statistical tests can fail in both directions: In Subsection 2.2 (before Theorem 3) we will discuss an example showing that they can erroneously infer causal relations even though they do not exist. This is because parts that are common two both objects, e.g., the two texts, are only suitable to prove a causal link if they are not too straightforward to come up with. On the other hand, causal relations can generate similarities between texts for which every efficient statistical analysis is believed to fail. We will describe an idea from cryptography to show this. A cryptosystem is called ROR-CCA-secure (Real or Random under Chosen Ciphertext Attacks) if there is no efficient method to decide whether a text is random or the encrypted version of some known text without knowing the key [12]. Given that there are ROR-CCA-secure schemes (which is unknown but believed by cryptographers) we have a causal relation leading to similarities that are not detected by any kind of simple counting statistics. However, once an attacker has found the key (maybe by exhaustive search), he recognizes similarities between the encrypted text and the plain text and infers a causal relation. This already suggests two things: (1) detecting similarities involves searching over potential rules how properties of one object can be algorithmically derived from the structure of the other. (2) It is likely that inferring causal relations therefore relies on computationally infeasible decisions (if computable at all) on whether two objects have information in common or not. 2.1 Algorithmic mutual information We will now describe how the information one object provides about the other can be measured in terms of Kolmogorov complexity. We start with some notation and terminology. Below, strings will always be binary strings since every description given in terms of a different alphabet can be converted into a binary word. The set of binary strings of arbitrary length will be denoted by {0,1}. Recall that the Kolmogorov complexity K(s) of a string s {0,1} is defined as the length of the shortest program that generates s using a previously defined universal Turing machine [13, 14, 15, 16, 17, 18, 19]. The conditional Kolmogorov complexity K(t s) [18] of a string t given another string s is the length of the shortest program that can generate t from s. In order to keep our notation simple we use K(x,y) to refer to the complexity of the concatenation of x,y. We will mostly have equations that are valid only up to additive constant terms in the sense that the difference between both sides does not depend on the strings involved in the equation (but it may depend on the Turing machines they refer to). To indicate such constants we denote the corresponding equality by = and likewise for inequalities. In this context it is important to note that the number n of nodes of the causal graph is considered to be a constant. Moreover, for every string s we define s as its shortest description. If the latter is not unique, we consider the first one in an lexicographic order. It is necessary to distinguish between K( s) and K( s ). This is because there is a trivial algorithmic method to generate s from s (just apply the Turing machine to s ), but there is no algorithm of length O(1) that computes the shortest description s from a general input s. One can show [19] that s (s,k(s)). Here, the equivalence symbol means that both sides can be obtained from each other by O(1) programs. The following equation for the joint algorithmic information of two strings x, y will be useful [20]: The conditional version reads [20]: K(x,y) = K(x)K(y x ) = K(x)K(y x,k(x)). (4) K(x,y z) = K(x z)k(y x,k(x z),z) (5) The most important notion in this paper will be the algorithmic mutual information measuring the amount of algorithmic information that two objects have in common. Following [21] we define: 9

10 Definition 2 (algorithmic mutual information) Let x,y be two strings. Then the algorithmic mutual information of x,y is I(x : y) := K(y) K(y x ). The mutual information is the number of bits that can be saved in the description of y when the shortest description of x is already known. The fact that one uses x instead of x ensures that it coincides with the symmetric expression [21]: Lemma 3 (symmetric version of algorithmic mutual information) For two strings x,y we have I(x : y) = K(x)K(y) K(x,y). In the following sections, non-vanishing mutual information will be taken as an indicator for causal relations, but more detailed information on the causal structure will be inferred from conditional mutual information. This is in contrast to approaches from the literature to measure similarity versus differences of single objects that we briefly review now. To measure differences between single objects, e.g. pictures [22, 23], one defines the information distance E(x, y) between the two corresponding strings as the length of the shortest program that computes x from y and y from x. It can be shown [22] that E(x,y) log = max{k(x y),k(y x)}, where log = means equality up to a logarithmic term. However, whether E(x,y) is small or large is not an appropriate condition for the existence and the strength of a causal link. Complex objects can have much information in common even though their distance is large. In order to obtain a measure that relates the amount of information that is disjoint for the two strings to the amount they share, Li et al. [23] and Bennett et al. [11] use the normalized distance measure or d s (x,y) := K(x y ) K(y x ) K(x, y) d(x,y) = max{k(x y),k(y x)} max{k(x),k(y)} I(x : y) = 1 K(x,y), The intuitive meaning of d s (x,y) is obvious from its direct relation to mutual information, and 1 d(x,y) measures the fraction of the information of the more complex string that is shared with the other one. Bennett et al. [11] propose to construct evolutionary histories of chain letters using such kinds of information distance measures. However, like in statistical causal inference, inferring adjacencies on the basis of strongest dependences is only possible for simple causal structures like trees. In the general case, non-adjacent nodes can share more information than adjacent ones when information is propagated via more than one path. Instead of constructing causal neighborhood relations by comparing information distances we will therefore use conditional mutual information. In order to define its algorithmic version, we first observe that Definition 2 can be rewritten into the less concise form I(x : y) = K(y) K(y x,k(x)). This formula generalizes more naturally to the conditional analog [20]: Definition 3 (conditional algorithmic mutual information information) Let x,y,z be three strings. Then the conditional mutual algorithmic information of x,y, given z is I(x : y z) := K(y z) K(y x,k(x z),z).. 10

11 As shown in [20] (Remark II.3), the conditional mutual information also is symmetric up to a constant term: Lemma 4 (symmetric algorithmic conditional mutual information) For three strings x,y,z one has: I(x : y z) = K(x z)k(y z) K(x,y z). Definition 4 (algorithmic conditional independence) Given three strings x,y,z, we call x conditionally independent of y, given z (denoted by x y z) if I(x : y z) 0. In words: Given z, the additional knowledge of y does not allow us a stronger compression of x. This remains true if we are given the Kolmogorov complexity of y, given z. The theory developed below will describe laws where symbols like x, y, z represent arbitrary strings. Then one can always think of sequences of strings of increasing complexity and statements like the equation holds up to constant terms are well-defined. We will then understand conditional independence in the sense of I(x : y z) = 0. However, if we are talking about three fixed strings that represent objects in real-life, this does not make sense and the threshold for considering two strings dependent will heavily depend on the context. For this reason, we will not specify the symbol any further. This is the same arbitrariness as the cutoff rate for statistical dependence tests. The definitions and lemmas presented so far were strongly motivated by the statistical analog. Now we want to focus on a theorem in [21] that provides a mathematical relationship between algorithmic and statistical mutual information. First we rephrase the following theorem Theorem of [18], showing that the Kolmogorov complexity of a random string is approximatively given by the entropy of the underlying probability distribution: Theorem 1 (entropy and Kolmogorov complexity) Let x = x 1,x 2,,x n be a string whose symbols x j A are drawn i.i.d. from a probability distribution P(X) over the finite alphabet A. Slightly overloading notation, set P(x) := P(x 1 ) P(x n ). Let H(.) denote the Shannon entropy of a probability distribution. Then there is a constant c such that H(P(X)) 1 A logn E(K(x n)) H(P(X)) c n n n n, where E(.) is short hand for the expected value with respect to P(x). Hence lim n 1 E(K(X)) = H(P(X)). n However, for our purpose, we need to see the relation between algorithmic and statistical mutual information. If x = x 1,x 2,,x n and y = y 1,y 2,,y n such that each pair (x j,y j ) is drawn i.i.d. from a joint distribution P(X,Y), the theorem already shows that lim n 1 E(I(x : y)) = I(X;Y). n This can be seen by writing statistical mutual information as H(P(X))H(P(Y)) H(P(X,Y)). 11

12 The above translations between entropy and algorithmic information refer to a particular setting and to special limits. The focus of this paper is mainly the situation where the above limits are not justified. Before we rephrase Theorem 5.3 in [21] which provides insights into the general case, we recall that a function f is called recursive if there is a program on a Turing machine that computes f(x) from the input x, and halts on all possible inputs. Theorem 2 (statistical and algorithmic mutual information) Given string-valued random variables X,Y with a recursive probability mass function P(x,y) over pairs (x,y) of strings. We then have I(X;Y) K(P) E(I(x : y)) I(X;Y)2K(P), where K(P) is the length of the shortest prefix-free program that computes P(x,y) from (x,y). We want to provide an intuition about various aspects of this theorem. (1) IfI(X;Y) is largecomparedto K(P)the expected algorithmicmutual informationisdominated by the statistical mutual information. (2) If K(P) is no longer assumed to be small, statistical dependences do not necessarily ensure that the knowledge of x allows us to compress y further than without knowing x. It could be that the description of the statistical dependences requires more memory space than its knowledge would save. (3) On the other hand, knowledge of x could allow us to compress y even in the case of a product measure on x and y. Consider, for instance, the case that we have the point mass distribution on the pair (x,y) with x = y. To describe a more sophisticated example generalizing this case we first have to introduce a family of product probability distributions on {0,1} n that we will need several times throughout the paper. Definition 5 (Defining product distributions by strings) Let P 0,P 1 be two probability distributions on {0,1} and c be a binary string of length n. Then P c := P c1 P c2 P cn defines a distribution on {0,1} n. We will later also needthe following generalization: If P 00,P 01,P 10,P 11 are four distributions on {0, 1}, then P c,d := P c1,d 1 P c2,d 2 P cn,d n defines also a family of product measures on {0,1} n that is labeled by two strings. Denote by P m c the m-fold copy of P c from Definition 5. It describes a distribution on {0,1} nm assigning the probbaility P m c (x) to x {0,1} nm. If Q(x,y) := P m c (x)p m c (y), knowledge of x in the typical case provides knowledge of c, provided m is large enough. Then we can compress y better than without knowing x because we do not have to describe c any more. Hence the algorithmic mutual information is large and the statistical mutual information is zero because Q is by construction a product distribution. In other words, algorithmic dependences in a setting with i.i.d sampling can arise from statistical dependences and from algorithmic dependences between probability distributions. 12

13 2.2 Markov condition for algorithmic dependences among individual objects Now we state the causal Markov condition for individual objects as a postulate that links algorithmic mutual dependences with causal structure: Postulate 5 (algorithmic causal Markov condition) Let x 1,...,x n be n strings representing descriptions of observations whose causal connections are formalized by a directed acyclic graph G with x 1,...,x n as nodes. Let pa j be the concatenation of all parents of x j and nd j the concatenation of all its non-descendants except x j itself. Then x j nd j pa j. As in Definition 4, the appropriate cut-off rate for rejecting G when I(x j : nd j pa j ) > 0 will not be specified here. This formulation is a natural interpretation of Postulate 2 in terms of algorithmic independences. The only point that remains to be justified is why we condition on pa j instead of pa j, i.e., why we aregiven the optimal joint compressionofthe parent strings. The main reasonis that this turns out to yield nice statements on the equivalence of different Markov conditions (in analogy to Lemma 1). Since the differences between I(x j : nd j pa j ) and I(x j : nd j pa j ) can only be logarithmic in the string lengths 1 we will not focus on this issue any further. If we apply Postulate 5 to a trivial graph consisting of two unconnected nodes, we obtain the following statement. Lemma 5 (causal principle for algorithmic information) If the mutual information I(x : y) between two objects x,y is significantly greater than zero they have some kind of common past. Here, common past between two objects means that one has causally influenced the other or there is a third one influencing both. The statistical version of this principle is part of Reichenbach s principle of the common cause [24] stating that statistical dependences between random variables 2 X and Y are always due at least one of the following three types of causal links: (1) X is a cause of Y or (2) vice versa or (3) there is a common cause Z. For objects, the term common past includes all three types of causal relations. For a text, for instance, it reads: similarities of two texts x,y indicate that one author has been influenced by the other or that both have been influenced by a third one. Before we construct a model of causality that makes it possible to prove the causal Markov condition we want to discuss some examples. If one discovers significant similarities in the genome oftwosortsofanimalsonewilltrytoexplainthe similaritiesbyrelatednessinthe senseofevolution. Usually, one would, for instance, assume such a common history if one has identified long substrings that both animals have in common. However, the following scenario shows two observations that superficially look similar, but nevertheless we cannot infer a common past since their algorithmic complexity is low (implying that the algorithmic mutual information is low, too). Assume two persons are instructed to write down a binary string of length 1000 and both decide to write the same string x = It seems straightforward to assume that the persons have communicated and agreed upon this choice. However, after observing that x is just the binary representation of π, one can easily imagine that it was just a coincidence that both wrote the same sequence. In other words, the similarities are no longer significant after observing 1 this is because K(x y) K(x y ) = O(log y ), see [19] 2 The original formulation considers actually dependences between events, i.e., binary variables. 13

14 that they stem from a simple rule. This shows that the length of the pattern that is common to both observations, is not a reasonable criterion on whether the similarities are significant. To understand the algorithmic causal Markov condition we will study its implications as well as its justification. In analogy to Lemma 1 we have Theorem 3 (equivalence of algorithmic Markov conditions) Given the strings x 1,...,x n and a directed acyclic graph G. Then the following conditions are equivalent: I. Recursive form: the joint complexity is given by the sum of complexities of each node, given the optimal compression of its parents: K(x 1,...,x n ) = n K(x j pa j ) (6) II. Local Markov condition: Every node is independent of its non-descendants, given the optimal compression of its parents: j=1 I(x j : nd j pa j) = 0. III. Global Markov condition: if R d-separates S and T. I(S : T R ) = 0 Below we will therefore no longer distinguish between the different versions and just refer to the algorithmic Markov condition. The intuitive meaning of eq. (6) is that the shortest description of all strings is given by describing how to generate every string from its direct causes. A similar kind of modularity of descriptions will also occur later in a different context when we consider description complexity of joint probability distributions. For the proofof Theorem 3 we will need a Lemma that is an analogueof the observationthat for any two random variables X,Y the statistical mutual information satisfies I(f(X);Y) I(X;Y) for every measurable function f. The algorithmic analog is to consider two strings x,y and one string z that is derived from x by a simple rule. Lemma 6 (monotonicity of algorithmic information) Let x,y,z be three strings such that K(z x ) = 0. Then I(z : y) I(x : y). This lemma is a special case of Theorem II.7 in [20]. We will also need the following result: Lemma 7 (monotonicity of conditional information) Let x,y,z be three strings. Then K(z x ) K(z (x,y) ). Note that K(z x ) K(z x,y) and K(z x ) K(z x,y ) is obvious but Lemma 7 is nontrivial because the star operation is jointly applied to x and y. 14

15 Proof of Lemma 7: Clearly the string x can be derived from x,y by a program of length O(1). Lemma 6 therefore implies I(z : x) I(z : x,y), where I(z : x,y) is shorthand for I(z : (x,y)). Hence K(z) K(z x ) = I(z : x) I(z : x,y) = K(z) K(z (x,y) ). Then we obtain the statement by subtracting K(z) and inverting the sign. The following lemma will only be used in Subsection 3.3. We state it here because it is closely related to the ones above. Lemma 8 (generalized data processing inequality) For any three strings x,y,z, I(x : y z ) = 0 implies I(x : y) I(x : z). The name data processing inequality is justified because the assumption x y z may arise from the typical data processing scenario where y is obtained from x via z. Proof of Lemma 8: Using Lemma 7 we have K(x y ) K(x (zy) ) (7) = K(x z,y,k(yz)) = K(x z,y,k(z)k(y z )) K(x z,y,k(z),k(y z ) = K(x z,y,k(y z )), where the second inequality holds because K(z) K(y z ) can obviously be computed from the pair (K(z),K(y z )) by an O(1) program. The last equality uses, again, the equivalence of z and (z, K(z)). Hence we obtain: I(x : y) = K(x) K(x y ) = K(x z )I(x : z) K(x y ) K(x z )I(x : z) K(x y,k(y z ),z ) = I(x : z)i(x : y z ) = I(x : z). The first step is by Definition 2, the second one uses Lemma 7, the third step is a direct application of ineq. (7), the fourth one is due to Definition 3, and the last step is by assumption. Proof of Theorem 3: I III: Define a probability mass function P on ({0,1} ) n, i.e., the set of n-tuples of strings, as follows. Set P(x j pa j ) := 1 z j 2 K(xj pa j ), (8) 15

16 where z j is a normalization factor. In this context, it is important that the symbol pa j refers to conditioning on the k-tuple of strings x i that are parents of x j (in contrast to conditional complexities where we can interpret K(. pa j ) equally well as conditioning on one string given by the concatenation of all those x i ). Note that Kraft s inequality (see [19], Example 3.3.1) implies 2 K(x y) 1, x for every y entailing that the expression is indeed normalizable by z j 1. We have Then we set K(x j pa j) = log 2 P(x j pa j ). P(x 1,...,x n ) := n P(x j pa j ), (9) i.e., P is by construction recursive with respect to G. It is easy to see that K(x 1,...,x n ) can be computed from P: j=1 K(x 1,...,x n ) = n K(x j pa j ) (10) j=1 = n log 2 P(x j pa j ) j=1 = log 2 P(x 1,...,x n ). Remarkably, we can also compute Kolmogorov complexities of subsets of {x 1,...,x n } from the corresponding marginal probabilities. We start by proving K(x 1,...,x n 1 ) = log 2 x n 2 K(x1,...,xn). (11) To this end, we observe x n 2 K(x1,...,xn) = 2 K(x1,...,xn 1) K(xn (x1,...,xn 1) ) x n 2 K(x1,...,xn 1), (12) where = denotes equality up to a multiplicative constant. The equality follows from eq. (4) and the inequality is obtained by applying Kraft s inequality [19] to the conditional complexity K(. (x 1,...,x n 1 ) ). On the other hand we have K(x 1,...,x n 1 ) = K(x 1,...,x n 1,0), since adding the 1-bit string x n = 0 certainly can be performed by a programof length O(1). Hence we have x n 2 K(x1,...,xn) 2 K(x1,...,xn 1,0) = 2 K(x1,...,xn 1). 16

17 Combining this with ineq. (12) yields 2 K(x1,...,xn 1) = x n 2 K(x1,...,xn). Using eq. (10) we obtain K(x 1,...,x n 1 ) = log 2 = log 2 x n 2 K(x1,...,xn) x n P(x 1,...,x n ) = log 2 P(x 1,...,x n 1 ), which proves equation (11). This implies K(x 1,...,x n 1 ) = log 2 P(x 1,...,x n 1 ). Since the same argument holds for marginalizing over any other variable x j we conclude that K(x j1,...,x jk ) = log 2 P(x j1,...,x jk ), (13) for every subset of strings of size k with k n. This follows by induction over n k. Now we can use the relation between marginal probabilities and Kolmogorov complexities to show that conditional complexities are also given by the corresponding conditional probabilities, i.e., for any two subsets S,T {x 1,...,x n } we have K(S T ) = log 2 P(S T). Without loss of generality, set S := {x 1,...,x j } and T := {x j1,...,x k } for j < k n. Using eqs. (4) and (13) we get K(x 1,...,x j (x j1,...,x k ) ) = K(x 1,...,x k ) K(x j1,...,x k ) = log 2 (P(x 1,...,x k )/P(x j1,...,x k ) ) = log 2 P(x 1,...,x j x j1,...,x k ). Let S,T,R be three subsets of {x 1,...,x n } such that R d-separates S and T. Then S T R with respect to P because P satisfies the recursion (9) (see Lemma 1) 3. Hence K(S,T R ) = log 2 P(S,T R) = logp(s R) log 2 P(T R) = K(S R )K(S R ). This proves algorithmic independence of S and T, given R and thus I III. To show that III II it suffices to recall that nd j and x j are d-separated by pa j. Now we show II I in strong analogy to the proof for the statistical version of this statement in [3]: Consider first a terminal node of G. Assume, without loss of generality, that it is x n. Hence all strings 3 Since P is, by construction, a discrete probability function, P the density with respect to a product measure is directly given by the probability function itself. 17

18 x 1,...,x n 1 are non-descendants of x n. We thus have (nd n,pa n ) (x 1,...,x n 1 ) where means that both strings coincide up to a permutation (on one side) and removing those strings that occur twice (on the other side). Due to eq. (4) we have K(x 1,...,x n ) = K(x 1,...,x n 1 )K(x n (nd n,pa n ) ). (14) Using, again, the equivalence of w (w,k(w)) for any string w we have K(x n (nd n,pa n ) ) = K(x n nd n,pa n,k(nd n,pa n )) = K(x n nd n,pa n,k(pa n )K(nd n pa n)) K(x n nd n,pa n,k(nd n pa n )) = K(x n pa n). (15) The second step follows from K(nd n,pa n ) = K(pa n )K(nd n pa n ). The inequality holds because nd n,pa n,k(pa n )K(nd n pa n )canbecomputedfromnd n,pa n,k(nd n pa n )viaaprogramoflength O(1). The last step follows directly from the assumption x n nd n pa n. Combining ineq. (15) with Lemma 7 yields K(x n (nd n,pa n ) ) = K(x n pa n ). (16) Combining eqs. (16) and (14) we obtain Then statement I follows by induction over n. K(x 1,...,x n ) = K(x 1,...,x n 1 )K(x n pa n). (17) To show that the algorithmic Markov condition can be derived from an algorithmic version of the functional model in Postulate 4 we introduce the following model of causal mechanisms. Postulate 6 (algorithmic model of causality) Let G be a DAG formalizing the causal structure among the strings x 1,...,x n. Then every x j is computed by a program q j with length O(1) from its parents pa j and an additional input n j. We write formally x j = q j (pa j,n j ), meaning that the Turing machine computes x j from the input pa j,n j using the additional program q j and halts. The inputs n j are jointly independent in the sense n j n 1,...,n j 1,n j1,n n. By defining new programs that contain n j we can, equivalently, drop the assumption that the programs q j are simple and assume that they are jointly independent instead. We could also have assumed that x j is a function f j of all its parents, but our model is more general since the map defined by the input-output behavior of q j need not be a total function [19], i.e., the Turing machine simulating the process would not necessarily halt on all inputs pa j,n j. The idea to represent causal mechanisms by programs written for some universal Turing machine is basically in the spirit of various interpretations of the Church-Turing thesis. One formulation, given by Deutsch [25], states that every process taking place in the real world can be simulated by 18

19 a Turing machine. Here we assume that the way different systems influence each other by physical signals can be simulated by computation processes that exchange messages of bit strings. 4 Note that mathematics also allows us to construct strings that are linked to each other in an uncomputable way. For instance, let x be an arbitrary binary string and y be defined by y := K(x). However, it is hard to believe that a real causal mechanism could create such kind of relations between objects given that one believes that real processes can always be simulated by algorithms. These remarks are intended to give sufficient motivation for our model. Postulate 6 implies the algorithmic causal Markov condition: Theorem 4 (algorithmic model implies Markov) Let x 1,...,x n be generated by the model in Postulate 6. Then they satisfy the algorithmic Markov condition with respect to G. Proof (straightforward adaption of the proof of Lemma 2): Extend G to a causal structure G with nodes x 1,...,x n,n 1,...,n n. To see that the extended set of nodes satisfy the local Markov condition w.r.t. G, observe first that every node xj is given by its parents via an O(1) program. Second, every n j is parentless and (unconditionally) independent of all its non-descendants because they can be computed from {n 1,...n n }\{n j } via an O(1) program. By Theorem 3 the extended set of nodes is also globally Markovian w.r.t. G. The parents paj d-separate x j and nd j in G (here the parents pa j are still defined with respect to G). This implies the local Markov condition for G. It is trivial to construct examples where the causal Markov condition is violated if the programs q j are mutually dependent (for instance, the trivial graph with two nodes x 1,x 2 and no edge would satisfy I(x 1 : x 2 ) > 0 if the programs q 1,q 2 computing x 1,x 2 from an empty input are dependent). The last sentence of Postulate 6 makes apparent that the mechanisms that generate causal relations are assumed to be independent. This is essential for the general philosophy of this paper. To see that such a mutual independence of mechanisms is a reasonable assumption we recall that the causal graph is meant to formalize all relevant causal links between the objects. If we observe, for instance, that two nodes are generated from their parents by the same complex rule we postulate another causal link between the nodes that explains the similarity of mechanisms Relative causality This subsection explains why it is sensible to define algorithmic dependence and the existence or non-existence of causal links relative to some background information. To this end, we consider genetic sequences s 1,s 2 of two persons that are not relatives. We certainly find high similarity that leads to a significant violation of I(s 1 : s 2 ) = 0 due to the fact that both genes are taken from humans. However, given the background information s 1 is a human genetic sequence, s 1 can be 4 Note, however, that sending quantum systems between the nodes could transmit a kind of information ( quantum information [26]) that cannot be phrased in terms of bits. It is known that this enables completely new communication scenarios, e.g. quantum cryptography. The relevance of quantum information transfer for causal inference is not yet fully understood. It has, for instance, been shown that the violation of Bell s inequality in quantum theory is also relevant for causal inference [27]. This is because some causal inference rules between classical variables break down when the latent factors are represented by quantum states rather than being classical variables. 5 One could argue that this would be just the causal principle implying that similarities of the machines generating x j from pa j has to be explained by a causal relation, i.e., a common past of the machines. However, in the context of this paper, such an argument would be circular. We have argued that the causal principle is a special case of the Markov condition and derived the latter from the algorithmic model above. We will therefore consider the independence of mechanisms as a first principle. 19