Причинный эффект при регулировке задней двери и передней двери


12

Если мы хотим вычислить причинное влияние на Y на приведенном ниже причинном графике, мы можем использовать как теоремы регулировки задней двери, так и теоремы регулировки передней двери, т. Е. P ( y | do ( X = x ) ) = u P ( y | x , u ) P ( u )XY

P(y|do(X=x))=uP(y|x,u)P(u)

и

P(y|do(X=x))=zP(z|x)xP(y|x,z)P(x).

Легко ли домашнее задание показать, что две корректировки приводят к одинаковому причинному влиянию на Y ?XY

график


Это настоящая домашняя работа? Затем, пожалуйста, добавьте тег самообучения. Тогда люди могут дать вам подсказки, оставив вам мышление (и обучение). Расскажите нам, что вы пытались и где вы застряли. Помните, резюме не для аутсорсинга домашней работы ...
Knarpie

5
Привет Кнарпи, это часть самостоятельной работы, а не домашняя работа. В настоящее время я читаю «Причинный вывод в статистике» Pearl et al. и потратить около 1 часа на размышления над вопросом, который я задал выше, так как это естественный вопрос, но он не может показать равенство. Либо я что-то здесь упускаю, либо два выражения не равны.
Jae

Ответы:


13

Действие соответствует вмешательству в переменную X, которое устанавливает его в значение xdo(x)Xx . Когда мы вмешиваемся в , это означает, что родители X больше не влияют на его значение, что соответствует удалению стрелок, указывающих на X. Итак, давайте представим это вмешательство в новом DAG.XXX

введите описание изображения здесь

Назовем исходное распределение наблюдений и распределение после вмешательства P . Наша цель состоит в том, чтобы выразить Р * в терминах P . Обратите внимание , что в Р * мы имеем , что U X . Кроме того, вероятности до вмешательства и после вмешательства разделяют эти две инвариантности: P ( U ) = P ( U ) и P ( Y | X , U ) = P ( Y |PPPPPUXP(U)=P(U) поскольку мы не касались стрелок, вводящих эти переменные в нашем вмешательстве. Так:P(Y|X,U)=P(Y|X,U)

P(Y|do(X)):=P(Y|X)=UP(Y|X,U)P(U|X)=UP(Y|X,U)P(U)=UP(Y|X,U)P(U)

XZ

P(Z|do(X))=P(Z|X)

P(Y|do(X))XZY

P(Y|do(Z))=XP(Y|X,Z)P(X)

Где я использую штрих для удобства обозначения следующего выражения. Таким образом, эти два выражения уже относятся к распределению перед вмешательством, и мы просто использовали предыдущее обоснование бэкдора для их получения.

XYZYXZP(Y|Z,do(X))=P(Y|do(Z),do(X))=P(Y|do(Z))XYZZYX

P(Y|do(X))=ZP(Y|Z,do(X))P(Z|do(X))=ZP(Y|do(Z))P(Z|do(X))=ZXP(Y|X,Z)P(X)P(Z|X)=ZP(Z|X)XP(Y|X,Z)P(X)

ZP(Y|do(Z))P(Z|do(X))ZYP(Y|do(Z))XZXP(Z|do(X))

Следовательно, две корректировки дают вам такое же распределение после вмешательства на этом графике, как мы показали.


Перечитывая ваш вопрос, я подумал, что вас может заинтересовать прямая демонстрация того, что правая часть двух уравнений одинакова в распределении до вмешательства (которым они должны быть, учитывая наш предыдущий вывод). Это не сложно показать напрямую. Достаточно показать, что в вашем DAG:

XP(Y|Z,X)P(X)=UP(Y|Z,U)P(U)

YX|U,ZUZ|X

XP(Y|Z,X)P(X)=X(UP(Y|Z,X,U)P(U|Z,X))P(X)=X(UP(Y|Z,U)P(U|X))P(X)=UP(Y|Z,U)XP(U|X)P(X)=UP(Y|Z,U)P(U)

Следовательно:

ZP(Z|X)XP(Y|X,Z)P(X)=ZP(Z|X)UP(Y|Z,U)P(U)=UP(U)ZP(Y|Z,U)P(Z|X)=UP(U)ZP(Y|Z,X,U)P(Z|X,U)=UP(Y|X,U)P(U)

2
P(Y|do(X))ZP(Y|do(Z))P(Z|do(X)ZP(Y|Z,do(X))P(Z|do(X))

@JulianSchuessler, поэтому я написал «можно думать как», как способ помочь пониманию, но не буквально говоря, что это так. Что касается деривации парадной двери, не было ясно, как ОП знал, как ее получить, поэтому я и положил ее туда.
Карлос

Отличный ответ. Спасибо, Карлос. Вторая часть вашего ответа была именно то, что я просил. У меня есть два дополнительных вопроса здесь. 1) Какую стратегию поиска вы использовали для алгебраической манипуляции выражениями во втором ответе? (Прищурившись достаточно долго в выражениях?) Поскольку пространство поиска велико, мне интересно, как можно написать алгоритм, позволяющий автоматически прийти к такому же выводу.
Jae

zP(Y|do(Z)P(Z|do(X))do(Z)ZZ

1
P(X,U)P(Y|X,U)MMPPYXU
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.