Goals:

  • Show that \(Y(w) \perp\!\!\!\perp W \mid X\) can be extracted from a SWIG


DAG

Define and plot the DAG representing the structural causal model \[\begin{align*} X = & f_X(U_X) \\ W = & f_W(X,U_W) \\ Y = & f_Y(W,U_Y) \end{align*}\]

using the dagitty infrastructure

library(tidyverse)  # For ggplot2 and friends
library(dagitty)    # For dealing with DAG math
library(ggdag)      # For making DAGs with ggplot

# DAG
dag = dagitty('dag{
W [exposure,pos = "1,1"]
Y [outcome,pos = "2,1"]
X [pos = "1,0.5"]
W -> Y
X -> W
X -> Y
}')

ggdag(dag) + theme_dag()

Check which (conditional) independences between observed variables are implied by the DAG:

impliedConditionalIndependencies(dag)

Not surprisingly none.



SWIG

Define and plot the SWIG implied by the structural causal model \[\begin{align*} X = & f_X(U_X) \\ W = & f_W(X,U_W) \\ Y = & f_Y(W,U_Y) \end{align*}\]

swig = dagitty('dag{
W [exposure,pos = "1,1"]
Yw [outcome,pos = "2,1"]
w [pos = "1.2,1"]
X [pos = "1,0.5"]
w -> Yw
X -> W
X -> Yw
}')

ggdag(swig) + theme_dag()

and observe that it implies the standard independence \(Y(w) \perp\!\!\!\perp W \mid X\) that is known as - Backdoor adjustment - Conditional independence assumption - Exogeneity - Ignorability - Measured confounding - No unmeasured confounding - Selection-on-observables - Unconfoundedness - …

impliedConditionalIndependencies(swig)
W _||_ Yw | X
W _||_ w
X _||_ w
LS0tDQp0aXRsZTogIkNhdXNhbCBJbmZlcmVuY2U6IERBRyBhbmQgU1dJRyBmb3IgbWVhc3VyZWQgY29uZm91bmRpbmciDQpzdWJ0aXRsZTogIkFwcGxpY2F0aW9uIG5vdGVib29rIg0KYXV0aG9yOiAiTWljaGFlbCBLbmF1cyINCmRhdGU6ICJgciBmb3JtYXQoU3lzLnRpbWUoKSwgJyVtLyV5JylgIg0Kb3V0cHV0OiANCiAgaHRtbF9ub3RlYm9vazoNCiAgICB0b2M6IHRydWUNCiAgICB0b2NfZmxvYXQ6IHRydWUNCiAgICBjb2RlX2ZvbGRpbmc6IHNob3cNCi0tLQ0KDQoNCkdvYWxzOg0KDQotIFNob3cgdGhhdCAkWSh3KSBccGVycFwhXCFcIVxwZXJwIFcgXG1pZCBYJCBjYW4gYmUgZXh0cmFjdGVkIGZyb20gYSBTV0lHDQoNCjxicj4NCg0KIyBEQUcNCg0KRGVmaW5lIGFuZCBwbG90IHRoZSBEQUcgcmVwcmVzZW50aW5nIHRoZSBzdHJ1Y3R1cmFsIGNhdXNhbCBtb2RlbCANClxbXGJlZ2lue2FsaWduKn0NCiAgICBYID0gJiBmX1goVV9YKSBcXA0KICAgIFcgPSAmIGZfVyhYLFVfVykgXFwNCiAgICBZID0gJiBmX1koVyxVX1kpDQpcZW5ke2FsaWduKn1cXQ0KDQp1c2luZyB0aGUgYGRhZ2l0dHlgIGluZnJhc3RydWN0dXJlDQoNCmBgYHtyLCB3YXJuaW5nID0gRkFMU0UsIG1lc3NhZ2U9RkFMU0V9DQpsaWJyYXJ5KHRpZHl2ZXJzZSkgICMgRm9yIGdncGxvdDIgYW5kIGZyaWVuZHMNCmxpYnJhcnkoZGFnaXR0eSkgICAgIyBGb3IgZGVhbGluZyB3aXRoIERBRyBtYXRoDQpsaWJyYXJ5KGdnZGFnKSAgICAgICMgRm9yIG1ha2luZyBEQUdzIHdpdGggZ2dwbG90DQoNCiMgREFHDQpkYWcgPSBkYWdpdHR5KCdkYWd7DQpXIFtleHBvc3VyZSxwb3MgPSAiMSwxIl0NClkgW291dGNvbWUscG9zID0gIjIsMSJdDQpYIFtwb3MgPSAiMSwwLjUiXQ0KVyAtPiBZDQpYIC0+IFcNClggLT4gWQ0KfScpDQoNCmdnZGFnKGRhZykgKyB0aGVtZV9kYWcoKQ0KYGBgDQoNCkNoZWNrIHdoaWNoIChjb25kaXRpb25hbCkgaW5kZXBlbmRlbmNlcyBiZXR3ZWVuIG9ic2VydmVkIHZhcmlhYmxlcyBhcmUgaW1wbGllZCBieSB0aGUgREFHOg0KDQpgYGB7cn0NCmltcGxpZWRDb25kaXRpb25hbEluZGVwZW5kZW5jaWVzKGRhZykNCmBgYA0KDQpOb3Qgc3VycHJpc2luZ2x5IG5vbmUuDQoNCjxicj4NCjxicj4NCg0KDQojIFNXSUcNCg0KRGVmaW5lIGFuZCBwbG90IHRoZSBTV0lHIGltcGxpZWQgYnkgdGhlIHN0cnVjdHVyYWwgY2F1c2FsIG1vZGVsIA0KXFtcYmVnaW57YWxpZ24qfQ0KICAgIFggPSAmIGZfWChVX1gpIFxcDQogICAgVyA9ICYgZl9XKFgsVV9XKSBcXA0KICAgIFkgPSAmIGZfWShXLFVfWSkNClxlbmR7YWxpZ24qfVxdDQoNCg0KYGBge3J9DQpzd2lnID0gZGFnaXR0eSgnZGFnew0KVyBbZXhwb3N1cmUscG9zID0gIjEsMSJdDQpZdyBbb3V0Y29tZSxwb3MgPSAiMiwxIl0NCncgW3BvcyA9ICIxLjIsMSJdDQpYIFtwb3MgPSAiMSwwLjUiXQ0KdyAtPiBZdw0KWCAtPiBXDQpYIC0+IFl3DQp9JykNCg0KZ2dkYWcoc3dpZykgKyB0aGVtZV9kYWcoKQ0KYGBgDQoNCg0KYW5kIG9ic2VydmUgdGhhdCBpdCBpbXBsaWVzIHRoZSBzdGFuZGFyZCBpbmRlcGVuZGVuY2UgJFkodykgXHBlcnBcIVwhXCFccGVycCBXIFxtaWQgWCQgdGhhdCBpcyBrbm93biBhcyANCi0gQmFja2Rvb3IgYWRqdXN0bWVudA0KLSBDb25kaXRpb25hbCBpbmRlcGVuZGVuY2UgYXNzdW1wdGlvbg0KLSBFeG9nZW5laXR5DQotIElnbm9yYWJpbGl0eQ0KLSBNZWFzdXJlZCBjb25mb3VuZGluZw0KLSBObyB1bm1lYXN1cmVkIGNvbmZvdW5kaW5nDQotIFNlbGVjdGlvbi1vbi1vYnNlcnZhYmxlcw0KLSBVbmNvbmZvdW5kZWRuZXNzDQotIC4uLg0KDQoNCmBgYHtyfQ0KaW1wbGllZENvbmRpdGlvbmFsSW5kZXBlbmRlbmNpZXMoc3dpZykNCmBgYA0KDQo=