在之前的文章因果，因果推断，因果关系是什么？中，我们介绍了什么是因果，什么是因果推断，什么是因果关系。
本文将围绕如下内容介绍如何进行因果关系发现研究：
- 因果发现的三个层级
- 因果发现的重要工具
- 因果发现的具体方法
（本文共4137字，预计阅读时长10分钟）
先从一组数据说起：为了研究某个治疗方案是否对某疾病有疗效，研究人员安排了治疗组和控制组进行实验，每组40人。治疗组将严格按照方案进行治疗，而控制组的人只能得到一些安慰剂。在实验进行一段时间后，得到两组病人的生存率，如图所示。
- 人群总体：50%（治疗组生存率） ＞40% （控制组生存率）
- 男性：60%（治疗组生存率） ＜70% （控制组生存率）
- 女性：20%（治疗组生存率） ＜30% （控制组生存率）
上面的数据给出一个令人迷惑的结论：从人群总体的角度，我们发现治疗组的生存率要高于控制组。但当区分性别来看，会完全颠覆之前的结论，即无论男性还是女性，治疗组的生存率均低于控制组。
你是不是也有一点糊涂？别担心，这就是困扰了统计学家60多年的著名的辛普森悖论，是因为仅通过学习某个条件概率，就去回答因果性问题而产生的。
一、认知因果的三个层次
图灵奖得主Judea Pearl提出，认知因果包含三个层级 [1]。通过观察发现事物之间的关联只是第一个层级；在这之上，还需要对过程进行有目的的干预，才能去回答“如果进行了治疗，生存率是否会提高的问题”。认知因果的第三个层级是反事实推理，回答“假如没有…，那么…“这样的问题。
传统的机器学习最擅长做的是根据数据中呈现出来的相关性学习函数f(Y|X1,…,Xm)去拟合条件概率P(Y|X1,…,Xm)。这种机器学习模型只学习到了认知因果的第一个层次即关联。
干预，是对所有可能影响因果关系的因子进行控制。如果不使用严格的控制实验，仅通过数据进行干预，需要引入由do算子定义的干预分布，直观地描述对Xi进行有意图干预时，其他变量概率分布的变化情况 [2]。例如：P(Y|do(X)=a)描述了当X取值为a时，Y所对应的分布。写成如下形式：
P(X1,…,Xi−1,Xi+1…,Xn|do(Xi)=a)=P(X1,…,Xn)P(Xi|pai)I(Xi=a)
其中，pai表示所有产生Xi的原因；I(Xi=a)为恒等函数，即当Xi的值为a时，其值为1，当Xi的值不为a时，其值为0。等式右边的除数P(Xi|pai)表示在对Xi进行干预时，需要去除pai对Xi的影响。
如果能得到如下情况：
P(Xj|do(Xi=x,X∖ij=c))≠P(Xj|do(Xi=x′,X∖ij=c))
即在保持系统中除Xi和Xj外的其他变量不变的情况下，Xi值的变化引起了Xj的分布变化，那么我们可以说Xi是Xj的直接原因。
我们以辛普森悖论为例，直观地理解干预分布和条件分布之间的差别。
辛普森悖论案例
用条件概率衡量治疗对生存率的影响，计算如下：
同理P(S=survial|T=control)=40%。
用do算子衡量治疗对生存率的影响，计算如下:
同理P(S=survial|do(T)=control)=60%。
对比两组计算：条件概率中Gender对生存率的影响会因为Treatment而变化（见）P(G|T)）；而干预分布中Gender对Treatment的影响（作为除数的）P(T|G)）受到严格控制，其对生存率的影响并不会因为Treatment变化。
通过干预的推理告诉我们通过治疗，人们的存活率为40%，要低于控制组的60%。但这些都是在研究总体效应或平均因果效应，如果针对特定事件或个体层面谈论个性化因果关系还需要进行反事实推理。比如在治疗了一段时间后，某病人的生存率提高，但这真是因为服用了特定的药物吗？还是因为听到了好消息？在反事实层，因为我们不可能观测到未曾发生的事情，只能在想象的世界里，推测现象原因为何。
二、因果发现的工具
因果发现旨在从一堆纷繁复杂的数据中，挖掘出不同变量之间因果影响的网络结构。为了进行因果发现，我们需要认识两个描述系统因果机制的工具：因果图和结构因果模型。
2.1 因果图
因果图是在贝叶斯网络的基础上定义的，它们都通过有向无环图（DAG）的形式并遵从马尔可夫性和忠实性假设，抓住图与数据交互的关键，实现了图的连接性和变量的独立性之间的联结。区别是，贝叶斯网络是由一系列条件概率描述的有向无环图，而因果图引入由do算子定义的干预，突破了条件概率只能学习到相关性的限制，从而到达认知因果的第二个层级，学习到更稳定的结构。下图直观地描述了五个变量之间的因果关系：X1是X2和X3的共同原因, X3和X4又共同产生了X5。
接下来我们识别有向无环图中的几个重要结构：
在
head-to-tail和tail-to-tail中，X和Y关于Z条件独立，即X⊥Y|Z，其中head-to-tail中的Z为中间变量，tail-to-tail中的Z为混淆变量；在head-to-head（v-结构）中，X和Y无条件独立，或者说X和Y关于空集条件独立 X⊥Y|∅，Z被称为对撞点。
我们需要严格控制混淆变量以消除其带来的偏差。然而，长期以来，统计学家时常困惑于具体应该控制哪些变量。如果控制中间变量，会切断X和Y之间的间接因果联系，得出X对Y没有影响的错误结论；如果控制对撞点，则会错误地认为X对Y有因果关系。后门准则和前门准则，可以帮助识别并消除因果图中的混淆变量，将do算子表示的干预分布转化成条件分布，在此基础上可以利用统计方法进行因果推断。
后门准则：后门路径定义为连接X和Y且包含了指向X的箭头的所有路径。阻断X和Y之间的所有后门路径，可以防止X的信息在非因果方向上的传递。存在对撞点的后门路径被认为时天然阻断的。
以研究“吸烟”（原因变量）对“癌症”（结果变量）的影响为例，“吸烟->焦油沉积->癌症”为因果路径，其中，“焦油沉积”为中间变量。“吸烟<-吸烟基因->癌症”为一条从“吸烟”到“癌症”的后门路径，该路径包含指向“吸烟”的箭头，其中“吸烟基因“为混淆变量。此外，整个因果图中只有一条从“吸烟”到“癌症”的后门路径。因此，控制了“吸烟基因”，我们就阻断了从“吸烟”到“癌症”的所有后门路径。
吸烟案例的因果图
为了研究“吸烟”对“癌症”的因果效应，我们对“吸烟”使用do算子：
前门准则：前门路径就是指从X到Y的直接因果路径，即上述：“吸烟->焦油沉积->癌症“的路径。当因为缺乏必要的数据而无法阻断某条后门路径时，就要通过前门准则，将X对Y的因果效应分解为X对Z的因果效应和Z对Y的因果效应。
在吸烟的案例中，假设我们无法对吸烟基因进行测量，但是可以获取”吸烟“，”焦油沉积“，以及“癌症”这三个变量的数据。这时，我们将“吸烟”对“癌症”的平均因果效应癌症吸烟P(癌症|do(吸烟))，转化为焦油沉积吸烟P(焦油沉积|do(吸烟))和癌症焦油沉积P(癌症|do(焦油沉积))的加权。在计算焦油沉积吸烟P(焦油沉积|do(吸烟))时，路径“吸烟<-吸烟基因->癌症<-焦油沉积”中 “癌症”处出现的对撞天然地阻断了这条后门路径。在计算癌症焦油沉积P(癌症|do(焦油沉积))时，存在后门路径“焦油沉积<-吸烟<-吸烟基因->癌症“，通过控制”吸烟“可以阻断。
利用前门准则，最终可以得到：
2.2 结构因果模型
结构因果模型也称作函数因果模型，旨在通过一系列函数方程对由图所描述的因果关系进行定义，将因果发现转化为函数估计问题 [3]。
在结构因果模型中，除了因果变量Xj外，还有一组随机变量Ej，它们只影响所对应的Xj，且相互独立，用于描述Xj受环境影响产生的不确定性。
相比因果图，结构因果模型包含了更多的信息，其不仅蕴含了一个观测分布，还蕴含了干预分布和反事实分布，可以在因果图干预的基础上，进一步支持反事实推理。
结构因果模型不仅蕴含了观测分布，还蕴含了干预分布和反事实
找到一个结构因果模型就意味着找到数据蕴含的唯一的联合分布；根据马尔可夫假设，可以将其按照因果图进行分解，即 P(X1,…,Xd)=∏j=1dP(Xj|Paj)。但当存在head-to-tail以及tail-to-tail结构时，我们观测到的联合分布是一致的。例如：X和Y关于Z相互独立，可以得到如下(X,Y)的联合概率分布：
所以实际上，利用结构因果模型只能够确定一类具有相同骨架（无向图）和v-结构的部分有向无环图（CPDAG），即马尔可夫等价类。
下图中，（a）是真实的因果图，（b）是其骨架，（d）和（e）是（a）的所有马尔可夫等价图；（c）展示了图（a）马尔可夫等价类的CPDAG [4]。
因果图的马尔可夫等价类和CPDAG
当可以通过一个结构因果模型找到一个马尔可夫等价类时，我们就说他是可识别的。这需要对Nj和fj作相应的假设。当假设所有的Nj服从相互独立的高斯分布时，关于fj的形式和结构因果模型的可识别性有如下结果：
三、因果发现的具体方法
基于观测的因果发现，就是要对生成数据的因果结构进行估计。用于确定一个因果图的马尔可夫等价类的方法，主要包括两类：基于约束的方法和基于分数的方法。
基于约束的方法主要通过一系列假设检验，判断在给定一个集合A后，两个节点之间是否条件独立。以该类方法中的PC算法为例 [5]，为了避免对所有可能的子集A进行搜索，在因果图构建过程中，从一个全连接图开始，通过逐渐增加条件集合的大小，来判断两个变量之间是否条件独立。如果可以找到一个集合使得两个变量之间条件独立，就可以去掉两个变量之间的边。当无法删除任何边时，就得到了因果图的骨架。如果两个变量(X,Y)之间条件独立，即它们之间没有直接相连的边，但存在路径X−Z−Y，且Z无法使得(X,Y)条件独立，那么就可以确定一个对撞结构（v-结构）X→Z←Y。在找到图中的所有v-结构后，就确定了CPDAG。
基于分数的方法主要通过对模型类型作相应的假设限制，直接拟合一个结构因果模型结构。拟合的效果通常由一个打分函数进行定义，通过求解
G^:=argmaxG DAG over X S(D,G)
得到最优的图结构。该打分函数通常包括两个部分：一是要最大化对数据的拟合程度；二是对图结构的复杂程度进行惩罚。在对最优图结构进行搜索时，可以使用贪婪搜索方法等局部搜索方法，也可以使用如动态规划、混合整数规划进行精确搜索。
通过有效地结合基于独立性和基于分数的混合方法，节约计算成本，得到更加准确有效的模型估计。
基于独立性和基于分数的方法，只能够确定CPDAG，要画出完整的因果图，可以借助成对比较方法去探索局部信息，主要是利用反映在数据中的因果产生机制的不对称性，确定因果方向。
具体来讲：如果X和Y之间存在直接因果关系，根据概率公理无法区分因果方向，即p(x,y)=p(x|y)p(y)=p(y|x)p(x)。考虑如下结构因果模型：
≔X≔EX≔Y≔f(X)+EY
其中，EX⊥ EY。由EX和EY相互独立，可知X与[Y|X]相互独立；然而在非因果方向上Y与[X|Y]却不一定独立。如下图所示，区别于因果方向上的独立性，即f(X)−Y与EY独立，用线性方程g(Y)去拟合X，得到的残差与Y之间存在明显的相关性。
数据中的不对称性可用于确定因果方向
结语
因果推理是非常强大的解释、分析和建模工具，与机器学习结合，可以提取更稳定和可解释的特征信息，增加模型的泛化能力。
比如，在因果表征学习中，对一个复杂系统的因果环节进行解构，就可以在系统发生变化时，定位变化模块，解释变化原因，进行局部干预，甚至是反事实推理 [6]。以下图为例，系统中一个机械手指的动作发生变化，造成了红色方块掉落。在右图所示的像素空间中，红色方块掉落也遮挡了背景中其他物体，这使得机械手指引起的变动信息与系统中其他互不相干的物体纠缠在一起，无法分解。在左图中，将整个系统的物理机制用因果图进行描述，那么只有机械手指对应的节点和红色方块节点及其子节点发生了变化。由此，在因果表征空间里实现了信息解耦。这种信息解耦，可以帮助探索事物内部的物理机制和逻辑关系，对视频跟踪监控、自动驾驶、飞行控制等自动化系统的精准智能感知、智能辨识与控制预测至关重要。
此外，因果推断在半监督学习、域适应、迁移学习或稳定学习、情景强化学习等多个领域都有极大的研究意义和应用价值。
参考：
系统参考：关河因果分析系统：https://yinguo.grandhoo.com/home
- ^J. Pearl and D. Mackenzie (2018). The book of why. Basic Books.
- ^J. Pearl et al. (2016). Causal Inference in Statistics: A Primer. Wiley.
- ^J. Peters et al. (2017). Elements of Causal Inference: Foundations and Learning Algorithms. MIT Press.
- ^O. Goudet, et al. (2018). Learning Functional Causal Models with Generative Neural Networks. In Explainable and Interpretable Models in Computer Vision and Machine Learning, 39-88.
- ^M. Kalisch and P. Bühlmann (2007), Estimating high-dimensional directed acyclic graphs with the PC-algorithm. Journal of Machine Learning Research.
- ^B. Schölkopf, et al. (2021). Towards Causal Representation Learning. https://arxiv.org/abs/2102.11107.