Skip to content

模型Model

结构因果模型

施工中 · 内容逐步修订

Structural causal model · SCM · Structural equation model for causality

以模块化结构方程和外生变量共同定义观察、干预与反事实分布的因果模型。

形式陈述 ​

观察数据只记录系统自然运行时的结果。结构因果模型额外规定:若外部改变某个变量,其余变量怎样继续生成。形式上,它是四元组 M=(U,V,F,PU)。U 为外生随机变量集合,V=(V1,…,Vk) 为内生变量,PU 给出外生变量联合分布,各 fi 为可测函数,结构方程组为

Vi=fi(pai,Ui),i=1,…,k

其中 pai 表示方程实际依赖的其他内生变量,Ui 可以是外生向量的一部分,各个 Ui 不必独立。方程规定每个内生变量如何由其直接原因和外生扰动生成。递归 SCM 要求存在一个顺序,使 fi 只依赖更早的内生变量;对应箭头 Vj→Vi 构成有向无环图。在这个顺序下,给定 U=u 后可以从第一个方程逐个代入,因而得到唯一解。再让 u 按 PU 随机变化,就得到观测分布 PM(V)。这里的唯一性由递归结构保证;对含反馈的方程组则不能直接断言。

原子干预 do(Vj=v) 以常数方程 Vj:=v 替换原方程,并保持其他方程和 PU 不变。在递归情形下,替换后的模型 Mv 仍有唯一解,从而诱导干预分布 PM(⋅∣do(v))。反事实 Yv(u) 是同一个外生状态 u 在替换方程后的解;同一 u 把现实世界和反事实世界配对。

因果目标量例如 EM(Y1−Y0) 属于干预/反事实分布。识别问题问所有产生同一观测分布且满足图与结构限制的模型是否给它相同值;估计问题才从样本拟合右侧可识别泛函。SCM 提供语义,不保证任意效应由 P(V) 识别。

直觉

每条方程是一个可替换模块:它描述变量在其原因给定时如何响应。观察 A=a 只是筛选自然系统中恰好生成 a 的单位;干预 do(A=a) 则断开生成 A 的旧机制,把所有单位送入同一个赋值模块。两者不同,正是混杂产生时条件分布与干预分布不同的根源。

图只保留方程间的依赖骨架,函数和外生分布保留效应大小、非线性与潜在共同原因。相同图可以容纳很多观测分布;相同观测分布也可来自因果方向不同的模型。因此箭头含义来自模块化因果假设,而非相关矩阵自动画出的方向。

例子与边界

令 UL∼Bernoulli(1/2),UA∼Bernoulli(0.2),E(UY)=0,三者独立,并设

L=UL,A=L⊕UA,Y=A+2L+UY,

其中 ⊕ 为异或。图为 L→A、L→Y、A→Y。由于 P(L=1∣A=1)=0.8、P(L=1∣A=0)=0.2,

E(Y∣A=1)−E(Y∣A=0)=(1+2×0.8)−(2×0.2)=2.2.

把 A 的方程替换为 A:=a 后,L 的分布不变,故 E(Y∣do(1))=1+2E(L)=2,E(Y∣do(0))=1,因果差为 1。关联差 2.2 因而混合了处理作用和 L 的组成差异。这个例子还可以从观测量核验干预均值:由于每个 L 层都有两种处理,且扰动独立,E(Y∣A=a,L=l)=a+2l。按总体中两个 L 层各半平均,得到 a+1,与替换方程的结果一致。第一种计算利用完整机制,第二种计算说明在这些假设下,效应确实可由观测分布识别。

本例的三节点图包含 L→A、L→Y、A→Y,对观测分布没有额外的条件独立约束。因此即使 UA 与 UY 相关,观测分布仍可满足它的因果 Markov 条件;但这不再保证观测条件均值 E(Y∣A,L) 等于只改变处理后的结果机制。相关噪声不使 SCM 失效,却要求在干预解释中保留潜在共同原因。对于确实有缺边的图,遗漏这样的共同原因还可能破坏图要求的条件独立。

SCM 与潜在结果表示可以相互翻译许多具体问题:方程替换生成 Y(a)。但二者不能粗暴称为完全等价。SCM 通常对所有变量的生成机制和跨世界耦合给出更多结构;只给出 Y(0),Y(1) 的潜在结果分析未必指定中介方程或因果图。

循环反馈、多个方程解和随机动态系统超出递归唯一解的简单版本。可以使用带时间展开、平衡选择规则或循环 SCM,但干预语义与可识别结论需要重新声明。用 DAG 强行表示同期反馈会丢失问题的核心。

仍用上面的加性模型,若已观察某人的 L=1,A=1,Y=4,结构方程确定 UY=1。把同一个人的处理改为0后,Y0=0+2×1+1=3;这与总体 E(Y∣do(0))=1 不同。前者保留已观测个体的外生状态,后者对总体所有外生状态平均。这也说明反事实分析需要比总体干预均值更多的条件信息。

推论与应用

当递归 SCM 的外生扰动相互独立时,观测分布按图的父节点分解;这给出 d-separation 可读出的条件独立。干预替换一个方程后,相应入边删除,得到截断分解。由此再推导 backdoor 调整与 do-calculus,而不是把这些公式当回归恒等式。

反事实查询按“溯因—行动—预测”计算:先用已观测证据更新外生状态分布,再替换干预方程,最后在更新后的同一外生状态上求新结果。若只先验模拟一个独立新个体,得到的是总体干预分布,不是对已观察个体的反事实;两种量的条件信息不同。

参考资料
  • Brady Neal, Introduction to Causal Inference from a Machine Learning Perspective, 在线教材,2020,§§2.3、4.5、7.1,作者稿。
  • Judea Pearl, Madelyn Glymour, and Nicholas P. Jewell, Causal Inference in Statistics: A Primer, Wiley, 2016,§§1.5、3.2、4.2–4.3,作者公开稿。
关系图谱9 个相邻概念 · 1 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组