Rel 0.1.0

Author: goedman

examples/SR/SR6.4.2/sr6.4.2.jl

@@ -35,7 +35,7 @@ d = OrderedDict(
 );
 u = [:u]
 
-dag = DAG("sr6_4_2", d, df);
+dag = DAG("sr6_4_2", d, df=df);
 show(dag)
 
 fname = joinpath(mktempdir(), "sr6.4.2.dot")
@@ -46,16 +46,16 @@ to_ggm(dag) |> display
 
 display(dag.s); println()
 
-bs = basis_set(dag, debug=true)
+bs = basis_set(dag)
 display(bs); println()
 
 t = shipley_test(dag)
 display(t); println()
 
 f = [:a]; s = [:b]; conditioning_set = [:u, :c]
 
-e = d_separation(dag, f, s, conditioning_set)
-println("d_separation($(dag.name), $f, $s, $conditioning_set) = $e")
+e = d_separation(dag, f, s, cset=conditioning_set)
+println("d_separation($(dag.name), $f, $s, cset=$conditioning_set) = $e")
 println()
 
 adjustmentsets = adjustment_sets(dag, :x, :y)

examples/SR/SR6.4.3/SR6.4.3.png

@@ -26,7 +26,7 @@ d = OrderedDict(
   :m => [:a]
 );
 
-dag = DAG("sr6_4_3", d, df);
+dag = DAG("sr6_4_3", d, df=df);
 show(dag)
 
 fname = ProjDir * "/sr6.4.3.dot"
@@ -43,10 +43,9 @@ t = shipley_test(dag)
 display(t); println()
 
 f = :w; s = :d;
-#cond = [:m, :a]
 
-e = d_separation(dag, [f], [s], [:m, :a])
-println("d_separation($(dag.name), $f, $s, $cond) = $e\n")
+e = d_separation(dag, [f], [s], cset=[:m, :a])
+println("d_separation($(dag.name), $f, $s, cset=[:m, :a]) = $e\n")
 
 #ap = all_paths(dag, f, s)
 ap = all_paths(dag, f, s)

examples/SR/SR6.4.3/sr6.4.3.jl

@@ -34,7 +34,7 @@ d = OrderedDict(
   :D => [:M]
 );
 
-dag = DAG("chain", d, df);
+dag = DAG("chain", d; df=df);
 show(dag)
 
 display(dag.s); println()
@@ -46,10 +46,10 @@ end
 println()
 
 f = [:A]; s = [:D]; sel = vcat(f, s)
-cond = [:M]
+cset = [:M]
 
 e = d_separation(dag, f, s)
 println("d_separation($(dag.name), $f, $s) = $e")
 
-e = d_separation(dag, f, s, cond)
-println("d_separation($(dag.name), $f, $s, $cond) = $e")
+e = d_separation(dag, f, s; cset=cset)
+println("d_separation($(dag.name), $f, $s; cset=$cset) = $e")

examples/SR/chain/AMD.jl

@@ -32,7 +32,7 @@ d = OrderedDict(
   :D => [:A]
 );
 
-dag = DAG("fork_correlated", d, df);
+dag = DAG("fork_correlated", d, df=df);
 show(dag)
 
 display(dag.s); println()
@@ -44,13 +44,13 @@ end
 println()
 
 f = [:M]; s = [:D]; sel = vcat(f, s)
-cond = [:A]
+cset = [:A]
 
 e = d_separation(dag, f, s)
 println("d_separation($(dag.name), $f, $s) = $e")
 
-e = d_separation(dag, f, s, cond)
-println("d_separation($(dag.name), $f, $s, $cond) = $e")
+e = d_separation(dag, f, s, cset=cset)
+println("d_separation($(dag.name), $f, $s, cset=$cset) = $e")
 println()
 N = 100
 df = DataFrame();
@@ -82,7 +82,7 @@ d = OrderedDict(
   :D => [:A]
 );
 
-dag = DAG("fork_correlated", d, df);
+dag = DAG("fork_correlated", d, df=df);
 show(dag)
 
 display(dag.s); println()
@@ -94,13 +94,13 @@ end
 println()
 
 f = :M; s = :D; sel = vcat(f, s)
-cond = :A
+cset = :A
 
 e = d_separation(dag, f, s)
 println("d_separation($(dag.name), $f, $s) = $e")
 
-e = d_separation(dag, f, s, cond)
-println("d_separation($(dag.name), $f, $s, $cond) = $e")
+e = d_separation(dag, f, s, cset=cset)
+println("d_separation($(dag.name), $f, $s, cset=$cset) = $e")
 
 as = adjustment_sets(dag, :A, :D)
 as |> display