tests/testthat/testChainGraphs.R

In dagitty: Graphical Analysis of Structural Causal Models

x <- dagitty("pdag{
	X1 [exposure]
	X2 [exposure]
	X3 [exposure]
	
	Y [outcome]
	
	X1 -> Z1 <- Y
	X2 <- Z2 -> Y
	X3 -- Z3 -> Y
	
	X1 -> Z2
}")

x2 <- dagitty("pdag{
	x [exposure]
	y [outcome]
	z -- x -> y
	x <- c -> y
	z -- c
}")

x3 <- dagitty("pdag{
	x [exposure]
	y [outcome]
	z -- x -> y
	x <- c -> y
}")

test_that("adjustment set construction", {
	expect_equal(length(adjustmentSets(x)), 0)	
	expect_equal(adjustmentSets(x,"X1")[[1]], list())	
	expect_equal(adjustmentSets(x,"X2")[[1]], c("Z2"))	
	expect_equal(length(adjustmentSets(x,"X3")), 0)
	expect_equal(length(adjustmentSets(x2)), 0)
	expect_equal(unlist(adjustmentSets(x3),use.names=FALSE), ("c"))
})