Goals:
- Show that \(Y(w) \perp\!\!\!\perp
W\) can be extracted from a SWIG
DAG
Define and plot the DAG representing the structural causal model
\[\begin{align*}
W = & f_W(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
## A/B test
## Randomized controlled trial
## Randomized experiment
# DAG
dag = dagitty('dag{
W [exposure,pos = "1,1"]
Y [outcome,pos = "2,1"]
W -> 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*}
W = & f_W(U_W) \\
Y = & f_Y(W,U_Y)
\end{align*}\]
swig_exp = dagitty('dag{
W [exposure,pos = "1,1"]
Yw [outcome,pos = "2,1"]
w [pos = "1.2,1"]
w -> Yw
}')
ggdag(swig_exp) + theme_dag()
and observe that it implies the standard independence \(Y(w) \perp\!\!\!\perp W\) that is created
be randomly allocating \(W\)
impliedConditionalIndependencies(swig_exp)
W _||_ Yw
W _||_ w
LS0tDQp0aXRsZTogIkNhdXNhbCBJbmZlcmVuY2U6IERBRyBhbmQgU1dJRyBmb3IgUkNUcyINCnN1YnRpdGxlOiAiQXBwbGljYXRpb24gbm90ZWJvb2siDQphdXRob3I6ICJNaWNoYWVsIEtuYXVzIg0KZGF0ZTogImByIGZvcm1hdChTeXMudGltZSgpLCAnJW0vJXknKWAiDQpvdXRwdXQ6IA0KICBodG1sX25vdGVib29rOg0KICAgIHRvYzogdHJ1ZQ0KICAgIHRvY19mbG9hdDogdHJ1ZQ0KICAgIGNvZGVfZm9sZGluZzogc2hvdw0KLS0tDQoNCg0KR29hbHM6DQoNCi0gU2hvdyB0aGF0ICRZKHcpIFxwZXJwXCFcIVwhXHBlcnAgVyQgY2FuIGJlIGV4dHJhY3RlZCBmcm9tIGEgU1dJRw0KDQo8YnI+DQoNCiMgREFHDQoNCkRlZmluZSBhbmQgcGxvdCB0aGUgREFHIHJlcHJlc2VudGluZyB0aGUgc3RydWN0dXJhbCBjYXVzYWwgbW9kZWwgDQpcW1xiZWdpbnthbGlnbip9DQogICAgVyA9ICYgZl9XKFVfVykgXFwNCiAgICBZID0gJiBmX1koVyxVX1kpDQpcZW5ke2FsaWduKn1cXQ0KDQp1c2luZyB0aGUgYGRhZ2l0dHlgIGluZnJhc3RydWN0dXJlDQoNCmBgYHtyLCB3YXJuaW5nID0gRkFMU0UsIG1lc3NhZ2U9RkFMU0V9DQpsaWJyYXJ5KHRpZHl2ZXJzZSkgICMgRm9yIGdncGxvdDIgYW5kIGZyaWVuZHMNCmxpYnJhcnkoZGFnaXR0eSkgICAgIyBGb3IgZGVhbGluZyB3aXRoIERBRyBtYXRoDQpsaWJyYXJ5KGdnZGFnKSAgICAgICMgRm9yIG1ha2luZyBEQUdzIHdpdGggZ2dwbG90DQoNCiMjIEEvQiB0ZXN0DQojIyBSYW5kb21pemVkIGNvbnRyb2xsZWQgdHJpYWwNCiMjIFJhbmRvbWl6ZWQgZXhwZXJpbWVudA0KDQojIERBRw0KZGFnID0gZGFnaXR0eSgnZGFnew0KVyBbZXhwb3N1cmUscG9zID0gIjEsMSJdDQpZIFtvdXRjb21lLHBvcyA9ICIyLDEiXQ0KVyAtPiBZDQp9JykNCg0KZ2dkYWcoZGFnKSArIHRoZW1lX2RhZygpDQpgYGANCg0KQ2hlY2sgd2hpY2ggKGNvbmRpdGlvbmFsKSBpbmRlcGVuZGVuY2VzIGJldHdlZW4gb2JzZXJ2ZWQgdmFyaWFibGVzIGFyZSBpbXBsaWVkIGJ5IHRoZSBEQUc6DQoNCmBgYHtyfQ0KaW1wbGllZENvbmRpdGlvbmFsSW5kZXBlbmRlbmNpZXMoZGFnKQ0KYGBgDQoNCk5vdCBzdXJwcmlzaW5nbHkgbm9uZS4NCg0KPGJyPg0KPGJyPg0KDQoNCiMgU1dJRw0KDQpEZWZpbmUgYW5kIHBsb3QgdGhlIFNXSUcgaW1wbGllZCBieSB0aGUgc3RydWN0dXJhbCBjYXVzYWwgbW9kZWwgDQpcW1xiZWdpbnthbGlnbip9DQogICAgVyA9ICYgZl9XKFVfVykgXFwNCiAgICBZID0gJiBmX1koVyxVX1kpDQpcZW5ke2FsaWduKn1cXQ0KDQoNCmBgYHtyfQ0Kc3dpZ19leHAgPSBkYWdpdHR5KCdkYWd7DQpXIFtleHBvc3VyZSxwb3MgPSAiMSwxIl0NCll3IFtvdXRjb21lLHBvcyA9ICIyLDEiXQ0KdyBbcG9zID0gIjEuMiwxIl0NCncgLT4gWXcNCn0nKQ0KDQpnZ2RhZyhzd2lnX2V4cCkgKyB0aGVtZV9kYWcoKQ0KYGBgDQoNCg0KYW5kIG9ic2VydmUgdGhhdCBpdCBpbXBsaWVzIHRoZSBzdGFuZGFyZCBpbmRlcGVuZGVuY2UgJFkodykgXHBlcnBcIVwhXCFccGVycCBXJCB0aGF0IGlzIGNyZWF0ZWQgYmUgcmFuZG9tbHkgYWxsb2NhdGluZyAkVyQNCg0KYGBge3J9DQppbXBsaWVkQ29uZGl0aW9uYWxJbmRlcGVuZGVuY2llcyhzd2lnX2V4cCkNCmBgYA0KDQo=