The do-operator as graph surgery (modularity, ch. 7 p. 154)

do(X = x) deletes every edge into X and leaves every other structural equation untouched.

Usage

morie_dointerv(dag, x)

Arguments

dag

directed edge matrix

x

intervened node(s)

Value

list(edges, removed, nremoved, nkept, nnodes)

Examples

V <- c(1, 2, 3, 4, 5, 6, 7, 8)
morie_dointerv(V, V)
#> $edges
#>      [,1] [,2]
#> 
#> $removed
#> [1] 4
#> 
#> $nremoved
#> [1] 4
#> 
#> $nkept
#> [1] 0
#> 
#> $nnodes
#> [1] 8
#>