Skip to content

结构因果模型

Structural causal model · SCM · Structural equation model for causality

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

领域
统计学
条目类型
模型

形式陈述

结构因果模型是四元组 M=(U,V,F,PU)U 为外生随机变量集合,V=(V1,,Vk) 为内生变量,PU 给出外生变量联合分布,结构方程组

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

规定每个内生变量如何由其直接原因和外生扰动生成。递归 SCM 要求存在一个顺序,使 fi 只依赖更早的内生变量;对应箭头 VjVi 构成有向无环图。给定 U=u 时方程有唯一解,推送 PU 得到观测分布 PM(V)

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

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

直觉

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

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

例子与边界

LBernoulli(1/2)UABernoulli(0.2)E(UY)=0,三者独立,并设

L=UL,A=LUA,Y=A+2L+UY,

其中 为异或。图为 LALYAY。由于 P(L=1A=1)=0.8P(L=1A=0)=0.2

E(YA=1)E(YA=0)=(1+2×0.8)(2×0.2)=2.2.

A 的方程替换为 A:=a 后,L 的分布不变,故 E(Ydo(1))=1+2E(L)=2E(Ydo(0))=1,因果差为 1。这个可复算机制直接展示相关 2.2 不等于干预效应 1。

UAUY 相关,却没有在图中加入表示共同外生原因的双向边语义,观测分布可能不再满足预期的因果 Markov 条件。外生独立不是符号装饰;它决定缺失箭头对应哪些独立约束。

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

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

推论与应用

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

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

参考资料
  • Judea Pearl, Causality: Models, Reasoning, and Inference, 2nd ed., Cambridge University Press, 2009,§§1.3–1.4、3.2 与 7.1。
  • Judea Pearl, Madelyn Glymour, and Nicholas P. Jewell, Causal Inference in Statistics: A Primer, Wiley, 2016,Ch. 1–3。
  • Peter Spirtes, Clark Glymour, and Richard Scheines, Causation, Prediction, and Search, 2nd ed., MIT Press, 2000,Ch. 3。
关系图谱9 个相邻概念 · 1 类关系

拖动节点调整位置。

显示关系

显示:依赖

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