Initial Commit

Francisco Coelho
1 parent 8eb857b6
Showing 10 changed files with 6342 additions and 0 deletions Show diff stats
drafts.jl
model_sample-[0.0;0.3;0.2].csv
model_sample-[0.0;0.3;0.5].csv
model_sample-[0.1;0.3;0.2].csv
model_sample-[0.1;0.3;0.5].csv
model_sample-[0.1;0.3;0.8].csv
model_sample.csv
sample.csv
sample_analysis.jl
sampling.jl
@@ -0,0 +1,67 @@
+using Symbolics
+using Latexify
+
+# # γ = 0.5
+# no_s = [
+#     0   ,
+#     23  ,
+#     614 ,
+#     165 ,
+#     169 ,
+#     0   ,
+#     0   ,
+#     4   ,
+#     25  ,
+# ]
+
+# γ = 0.8
+no_s = [
+    0   ,
+    28  ,
+    632 ,
+    246 ,
+    59 ,
+    0   ,
+    0   ,
+    5   ,
+    27  ,
+]
+pr_s = (x -> x // 1000).(no_s)
+
+@variables θ
+num_e = [
+    0           , 
+    0           , 
+    7           , 
+    3 * θ       , 
+    3 * (1 - θ) , 
+    0           , 
+    0           , 
+    3           , 
+    10          ,
+]
+pr_e = (x -> x // 23).(num_e)
+
+target = expand(sum( (x -> x^2).(pr_s - pr_e) ))
+println(latexify(target))
+
+# using Plots
+# g(t) = (20869963/66125000) + (477/52900)*t + (18/529)*(t^2)
+# t = 0:0.1:1
+# plot(t, g.(t))
+
+
+
+function solve2(a, b, c)
+    delta = sqrt(Complex(b^2 - 4 * a * c))
+    return ( (-b - delta)/(2*a), (-b + delta)/(2*a) )
+end
+    
+# # g = target ~ 0
+
+a = float(18//529)
+b = -float(21903//264500)
+c = float(188207311//529000000)
+
+println("a: $a b: $b c: $c")
+println("̂θ = $(-b/a)")
@@ -0,0 +1,1001 @@
+event
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+ac
+ac
+A
+ac
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+ab
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+A
+ab
+A
+A
+A
+ab
+A
+ab
+A
+A
+A
+ab
+A
+A
+A
+ac
+ac
+A
+ab
+A
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+ac
+A
+ac
+A
+A
+ac
+ac
+A
+ac
+ac
+A
+A
+A
+ac
+A
+A
+A
+ac
+A
+ac
+A
+ab
+ac
+A
+ac
+A
+A
+ab
+A
+ac
+A
+ac
+ab
+ac
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+A
+A
+ac
+ac
+ac
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+ac
+ac
+A
+ab
+A
+A
+ac
+ab
+A
+A
+A
+A
+ac
+ac
+ac
+A
+ac
+A
+A
+ac
+A
+A
+ab
+ac
+A
+ac
+ac
+A
+A
+ac
+ac
+A
+ac
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+ac
+A
+A
+ac
+A
+A
+A
+A
+A
+ab
+ac
+A
+ac
+A
+ac
+A
+ac
+A
+A
+ab
+A
+A
+ac
+ac
+A
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+ab
+ac
+A
+A
+ac
+A
+A
+A
+ab
+A
+ac
+ab
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+ac
+ac
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+ac
+A
+ac
+ac
+ac
+A
+A
+ac
+A
+A
+A
+A
+ac
+ac
+A
+A
+A
+A
+A
+A
+ac
+A
+ab
+A
+A
+ab
+ac
+ab
+ac
+ac
+A
+A
+ab
+ac
+A
+ac
+A
+A
+A
+A
+ac
+ac
+A
+A
+ac
+A
+A
+A
+A
+ab
+A
+A
+ac
+ac
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+ac
+A
+A
+A
+A
+ac
+ac
+A
+ac
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+ac
+ac
+A
+ac
+A
+A
+A
+A
+ac
+A
+ac
+A
+ac
+A
+A
+A
+A
+ab
+A
+ac
+A
+A
+ab
+A
+ac
+A
+ac
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+ac
+ab
+ab
+A
+A
+A
+ac
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+ac
+A
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+ab
+A
+ac
+A
+A
+A
+ab
+A
+ac
+ab
+A
+A
+ac
+A
+ab
+A
+ac
+A
+A
+A
+A
+A
+ac
+A
+ac
+ab
+ac
+A
+A
+A
+ac
+A
+A
+ac
+A
+ac
+ab
+A
+ab
+ac
+A
+ac
+ac
+ac
+A
+A
+ab
+A
+ab
+A
+A
+ac
+A
+ac
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+ac
+ab
+A
+A
+A
+A
+ac
+ac
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+ac
+A
+A
+A
+ac
+A
+ac
+ac
+ac
+A
+A
+A
+A
+A
+ab
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+ac
+ac
+A
+A
+ac
+A
+A
+A
+A
+A
+ab
+ac
+A
+A
+ac
+A
+A
+A
+ab
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+ac
+A
+A
+ac
+ac
+ac
+A
+A
+A
+ac
+A
+ac
+A
+A
+A
+A
+ab
+ac
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+ac
+A
+A
+ac
+A
+A
+A
+ab
+A
+ac
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+ab
+ab
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+ac
+ac
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+ac
+ac
+A
+A
+ac
+A
+ab
+A
+ac
+A
+A
+ac
+ac
+A
+A
+ab
+A
+A
+A
+A
+ac
+A
+ac
+ac
+A
+A
+ab
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+ab
+ab
+ac
+A
+A
+ac
+A
+A
+ac
+ac
+A
+ac
+A
+A
+ac
+ab
+A
+A
+A
+ac
+ac
+ac
+A
+A
+A
+ac
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+ac
+A
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+A
+ac
+ac
+ac
+ab
+A
+ac
+A
+A
+A
+A
+A
+ac
+ab
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+ab
+A
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+A
+ac
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+ac
+ac
+ac
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+ac
+A
+ac
+A
+ac
+A
+A
+A
+A
+A
+A
@@ -0,0 +1,1001 @@
+event
+A
+A
+ac
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+ab
+ab
+ab
+A
+A
+A
+ac
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+ac
+ac
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+ab
+A
+A
+A
+ab
+A
+A
+ac
+A
+A
+A
+A
+ab
+A
+A
+ab
+ab
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+ab
+A
+A
+ab
+A
+A
+A
+A
+ab
+A
+ab
+A
+ac
+A
+ab
+A
+A
+ac
+ac
+ab
+A
+ab
+A
+ab
+A
+ab
+ab
+A
+A
+A
+A
+ac
+A
+A
+ab
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+ac
+A
+A
+A
+A
+A
+A
+ab
+ab
+A
+ac
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+ab
+ab
+A
+ac
+ac
+A
+A
+ab
+ac
+A
+A
+A
+A
+ac
+A
+ac
+A
+A
+A
+ab
+ac
+A
+ab
+A
+A
+A
+ac
+A
+ac
+A
+A
+A
+A
+A
+ab
+ab
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+ac
+A
+A
+ab
+ac
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+ab
+ac
+A
+A
+ac
+ac
+A
+ab
+A
+ac
+A
+ab
+ac
+ac
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+ab
+A
+A
+A
+ac
+A
+ac
+A
+ac
+ab
+A
+A
+ac
+A
+A
+ab
+A
+A
+ac
+A
+ab
+ac
+A
+ac
+A
+A
+ab
+ab
+A
+A
+A
+A
+ac
+A
+ac
+A
+A
+A
+A
+A
+ab
+ac
+ac
+A
+A
+ab
+ac
+A
+A
+A
+ab
+A
+A
+ab
+ab
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+ab
+ac
+ac
+ab
+ac
+A
+A
+ac
+A
+A
+ab
+ac
+A
+ab
+A
+A
+A
+ac
+A
+A
+A
+A
+ab
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+ab
+ab
+A
+A
+A
+A
+A
+ab
+ab
+A
+A
+A
+ac
+ab
+A
+A
+A
+ab
+A
+A
+ac
+A
+A
+ab
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+ac
+ab
+A
+A
+A
+ab
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+ab
+A
+A
+ab
+ac
+A
+ab
+A
+ab
+ac
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+ac
+ab
+ac
+A
+A
+A
+ab
+ac
+A
+A
+A
+ab
+A
+ab
+ac
+A
+A
+ac
+ab
+ab
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+ab
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+ab
+A
+A
+A
+ab
+A
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+ac
+ac
+ab
+ab
+A
+ab
+ab
+ab
+A
+ac
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+A
+A
+A
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+ac
+A
+ab
+A
+ac
+A
+A
+ab
+A
+A
+ac
+A
+A
+A
+ab
+ab
+ac
+ab
+A
+A
+A
+A
+A
+A
+A
+ac
+ab
+A
+A
+ab
+A
+A
+A
+A
+ab
+A
+A
+ab
+A
+ab
+A
+A
+A
+ac
+A
+ac
+A
+A
+A
+A
+ab
+ab
+ac
+A
+A
+ab
+ab
+A
+ac
+A
+A
+A
+A
+ab
+ab
+ab
+ab
+ac
+ab
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+ac
+ab
+ac
+A
+A
+ac
+ac
+A
+A
+A
+A
+A
+A
+ab
+A
+ac
+A
+A
+ac
+ab
+ac
+A
+A
+A
+ac
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+ac
+ac
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+ab
+A
+ac
+ab
+A
+ac
+ac
+A
+ab
+A
+A
+A
+ab
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+ac
+ab
+A
+ac
+A
+ab
+A
+A
+A
+A
+A
+A
+ab
+ac
+A
+A
+ac
+ab
+A
+A
+A
+ac
+ab
+A
+ab
+ac
+A
+ab
+ac
+A
+A
+A
+A
+ab
+ab
+A
+A
+A
+ac
+A
+A
+ab
+A
+A
+A
+A
+A
+ab
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+ab
+ab
+ac
+A
+ac
+A
+A
+ac
+ac
+ab
+A
+A
+A
+ac
+A
+A
+A
+ac
+ab
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+ab
+A
+A
+ab
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+ac
+A
+A
+ab
+ac
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+ac
+A
+ab
+A
+ac
+ab
+A
+A
+A
+A
+ab
+ab
+A
+ab
+ab
+A
+ac
+A
+ac
+A
+ac
+A
+A
+ac
+A
+A
+A
+ac
+ab
+A
+ab
+A
+A
+A
+ac
+ab
+A
+ab
+ac
+ac
+ab
+ab
+A
+ab
+A
+ac
+A
+ab
+A
+A
+ab
+ac
+A
+ab
+ab
+ab
+A
+A
+A
+A
+ac
+A
+ab
+A
+ab
+A
+A
+ab
+A
+A
+A
+A
+ac
+ab
+ab
+A
+A
+A
+A
+ab
+ac
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+ac
+A
+ac
@@ -0,0 +1,1001 @@
+event
+ac
+A
+ab
+A
+Cb
+A
+ac
+A
+ac
+A
+b
+A
+A
+A
+ac
+A
+A
+A
+ac
+A
+ac
+A
+A
+A
+A
+Ac
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+ac
+Δ
+A
+A
+A
+A
+A
+ac
+ac
+A
+A
+A
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+Δ
+ac
+A
+A
+A
+ac
+Ca
+a
+A
+A
+A
+A
+A
+ac
+A
+ac
+A
+A
+b
+Δ
+A
+ac
+A
+A
+A
+ac
+AB
+Δ
+Δ
+ac
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+Δ
+A
+A
+A
+A
+A
+A
+A
+C
+A
+A
+A
+A
+A
+c
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+A
+ac
+ab
+A
+A
+A
+ac
+A
+A
+ac
+ac
+ac
+ac
+A
+A
+ac
+ac
+A
+Δ
+A
+A
+A
+ac
+A
+A
+ac
+ab
+ac
+A
+A
+ac
+A
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+A
+C
+A
+A
+ac
+A
+A
+A
+ac
+A
+ac
+C
+A
+A
+A
+A
+ac
+A
+ac
+ac
+ab
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+ac
+A
+ab
+A
+A
+A
+ab
+A
+A
+A
+A
+ac
+A
+ac
+A
+A
+A
+A
+ab
+ac
+Ca
+A
+A
+ab
+A
+ac
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+Δ
+A
+ac
+A
+ac
+A
+A
+ac
+A
+bc
+ABC
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+ab
+ac
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+Δ
+A
+A
+A
+Δ
+ac
+A
+ac
+A
+A
+ab
+ab
+A
+A
+ab
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+ac
+ac
+A
+A
+A
+A
+A
+A
+ac
+A
+ac
+Δ
+A
+A
+A
+A
+ac
+A
+A
+A
+ab
+a
+A
+ac
+A
+A
+A
+ac
+ac
+ab
+ac
+A
+A
+A
+ac
+ac
+b
+A
+A
+A
+ac
+A
+ac
+ab
+A
+A
+B
+ab
+ac
+Δ
+A
+ab
+A
+ab
+A
+A
+A
+A
+A
+a
+A
+A
+Bac
+ac
+ac
+ac
+ac
+ac
+ac
+ac
+A
+A
+ac
+A
+A
+AB
+A
+ac
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+Δ
+A
+A
+A
+abc
+ac
+ac
+A
+A
+A
+A
+A
+ac
+ac
+ABC
+ac
+ac
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+Δ
+A
+A
+A
+A
+ac
+b
+A
+A
+ac
+A
+A
+ac
+ac
+A
+A
+ac
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+ab
+A
+A
+A
+A
+bc
+A
+A
+A
+A
+A
+A
+b
+A
+A
+A
+A
+ac
+ac
+A
+A
+A
+ac
+ac
+ac
+A
+A
+ac
+A
+ac
+A
+ac
+A
+A
+A
+A
+ab
+ac
+ac
+A
+Cab
+A
+ac
+A
+A
+ab
+ab
+b
+A
+BCa
+ac
+A
+A
+ac
+A
+Δ
+A
+A
+ac
+ab
+A
+A
+ab
+ac
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+ac
+ab
+A
+Cab
+A
+A
+A
+ac
+A
+A
+Δ
+ac
+ac
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+Ab
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+ab
+ac
+A
+ac
+ac
+C
+ab
+A
+A
+A
+A
+A
+A
+A
+ac
+ac
+A
+Δ
+A
+A
+A
+ac
+A
+ab
+A
+A
+ab
+A
+ac
+A
+A
+a
+A
+ab
+A
+ac
+A
+A
+A
+ac
+ac
+ac
+A
+ac
+ac
+A
+A
+bc
+ac
+A
+AB
+A
+A
+Bc
+A
+A
+A
+ab
+A
+A
+ac
+ac
+A
+ab
+bc
+A
+ac
+A
+ac
+ac
+ac
+ACb
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+ac
+ab
+ac
+A
+Δ
+A
+A
+A
+A
+A
+c
+ac
+A
+B
+ac
+ab
+A
+A
+ac
+A
+A
+ac
+ac
+A
+A
+A
+A
+ac
+ac
+A
+A
+A
+ac
+A
+A
+A
+ac
+A
+Δ
+A
+A
+ac
+A
+Δ
+ab
+A
+A
+A
+A
+A
+ac
+A
+ac
+A
+ac
+A
+ac
+A
+A
+A
+ab
+A
+A
+ac
+A
+C
+A
+b
+A
+C
+A
+A
+ac
+A
+A
+Δ
+A
+A
+ac
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+Ab
+A
+BCa
+A
+A
+ABc
+A
+A
+A
+ac
+A
+A
+ac
+A
+ac
+A
+A
+Bc
+A
+ac
+ac
+A
+A
+A
+ac
+ac
+A
+ab
+A
+ac
+ac
+ac
+A
+A
+A
+ab
+A
+A
+A
+ac
+A
+ac
+a
+A
+A
+ac
+ac
+A
+ac
+ac
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+ab
+A
+A
+A
+A
+C
+ac
+A
+A
+ac
+A
+ac
+A
+c
+ac
+ac
+A
+A
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+ac
+A
+ac
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+ab
+A
+A
+ac
+c
+A
+A
+ac
+A
+ac
+ac
+A
+ab
+ac
+A
+A
+ac
+ac
+A
+A
+A
+A
+ac
+ab
+Δ
+ac
+A
+A
+A
+A
+A
+ac
+ac
+bc
+A
+A
+A
+ac
+A
+ac
+A
+A
+ac
+A
+ac
+A
+ab
+A
+ac
+A
+A
+Δ
+A
+ac
+ac
+A
+A
+A
+ac
+ac
+A
+ac
+A
+c
+A
+A
+A
+A
+A
+A
+ac
+ac
+A
+A
+A
+ab
+A
+Ac
+BC
+A
+A
+Ca
+A
+ac
+Δ
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+ab
+A
+A
+ac
+ac
+A
+ac
+ac
+A
+A
+ac
+A
+A
+A
+a
+ac
+Δ
+A
+A
+A
+ac
+A
+A
+A
+C
+A
+A
+b
+ac
+A
+ab
+ac
+ac
+c
+A
+A
+A
+A
+A
+ac
+A
+ab
+A
@@ -0,0 +1,1001 @@
+event
+A
+A
+ac
+ac
+ac
+ab
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+ac
+A
+abc
+A
+A
+A
+ab
+A
+ac
+ac
+ac
+A
+A
+ab
+A
+A
+A
+A
+bc
+ab
+A
+A
+ab
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+ac
+AB
+A
+ac
+A
+A
+A
+A
+ab
+A
+ab
+ab
+A
+A
+A
+ac
+ab
+ac
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+Δ
+A
+ac
+ac
+A
+A
+ab
+ac
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+ac
+ac
+A
+A
+A
+ab
+A
+A
+ac
+A
+ab
+A
+A
+ac
+A
+A
+A
+A
+ac
+A
+ab
+Ab
+A
+ac
+ac
+A
+ab
+A
+A
+A
+ac
+A
+A
+A
+Δ
+ab
+ac
+ac
+ab
+A
+ac
+A
+ACb
+ab
+c
+ab
+ac
+ab
+A
+ac
+A
+A
+A
+A
+ac
+ab
+ab
+ac
+A
+ac
+A
+A
+A
+ABC
+A
+A
+A
+A
+ac
+ab
+B
+ab
+A
+A
+Δ
+ab
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+ab
+ab
+A
+A
+ac
+ab
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+ab
+ac
+ac
+A
+ab
+Δ
+A
+B
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+ac
+A
+ab
+A
+Δ
+A
+A
+A
+Abc
+ABC
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+ab
+BC
+A
+A
+ac
+ac
+A
+A
+ac
+A
+C
+A
+ab
+A
+A
+ab
+A
+A
+A
+ac
+A
+A
+A
+ab
+A
+ab
+ac
+A
+A
+A
+A
+ab
+A
+ac
+A
+Bac
+A
+ac
+ab
+A
+A
+A
+A
+A
+Δ
+BCa
+A
+ac
+A
+ac
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+ab
+ab
+A
+ac
+A
+A
+ac
+A
+A
+A
+ab
+A
+ab
+BCa
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+Cb
+A
+A
+A
+A
+A
+ab
+ac
+Δ
+c
+ac
+Cb
+A
+A
+ab
+A
+ac
+A
+ab
+ac
+A
+A
+A
+A
+A
+A
+ab
+ac
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+ab
+A
+ac
+A
+A
+ab
+A
+c
+A
+A
+ab
+ab
+A
+A
+ac
+c
+A
+A
+ac
+A
+AB
+ab
+ac
+A
+A
+b
+A
+A
+A
+Δ
+ab
+A
+A
+A
+A
+Δ
+ab
+b
+A
+A
+ac
+Cb
+c
+A
+ac
+A
+A
+ab
+a
+ab
+A
+ab
+ab
+A
+Δ
+Ab
+A
+ac
+bc
+A
+A
+ab
+ab
+A
+A
+A
+ab
+b
+A
+A
+A
+A
+A
+A
+ab
+A
+Δ
+A
+ab
+A
+ab
+A
+ab
+A
+A
+A
+ac
+ac
+A
+c
+ac
+A
+A
+ab
+A
+A
+A
+A
+ac
+A
+ac
+A
+A
+bc
+ACb
+A
+A
+A
+A
+A
+A
+Δ
+Δ
+c
+A
+ab
+ac
+A
+A
+ac
+A
+A
+A
+Δ
+A
+A
+ab
+A
+A
+ac
+A
+Ca
+ab
+A
+A
+ac
+A
+A
+A
+A
+Ab
+A
+A
+A
+ac
+A
+ab
+A
+ac
+A
+B
+A
+A
+A
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+Ca
+A
+A
+A
+ac
+A
+ab
+ac
+A
+A
+A
+A
+A
+Cab
+A
+ac
+ab
+A
+A
+A
+A
+A
+ac
+ab
+ab
+ab
+A
+A
+ab
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+ab
+ac
+ab
+A
+Δ
+A
+ac
+A
+A
+A
+A
+A
+AC
+A
+ac
+ab
+Ca
+ac
+ABC
+ac
+ac
+A
+ac
+A
+A
+A
+A
+ac
+ac
+A
+A
+ac
+A
+A
+A
+A
+ac
+ab
+A
+A
+A
+A
+A
+A
+C
+A
+ab
+A
+A
+Δ
+A
+A
+A
+A
+A
+A
+ab
+ac
+ac
+A
+A
+ab
+A
+A
+ab
+ab
+A
+A
+A
+A
+ab
+ab
+ac
+ACb
+A
+A
+ac
+A
+A
+ab
+A
+Δ
+ab
+ab
+ab
+A
+A
+ac
+ab
+A
+A
+ab
+A
+A
+A
+A
+A
+ab
+ab
+A
+A
+A
+A
+A
+ab
+A
+Δ
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+ab
+Δ
+ab
+A
+A
+A
+A
+ab
+ac
+ab
+A
+A
+ac
+a
+ab
+A
+A
+Δ
+ac
+A
+A
+ab
+A
+ac
+A
+A
+ac
+ac
+ab
+A
+A
+A
+A
+ac
+ac
+A
+A
+A
+c
+A
+A
+A
+ac
+A
+A
+ab
+a
+ab
+A
+A
+A
+ab
+A
+A
+ab
+ac
+ac
+A
+A
+ac
+ac
+ac
+A
+A
+ab
+A
+A
+A
+A
+A
+ac
+A
+ab
+A
+ab
+A
+A
+A
+A
+ac
+A
+ab
+A
+A
+ac
+ab
+A
+A
+A
+ab
+A
+ab
+Cb
+Bc
+A
+ac
+A
+ab
+A
+ab
+ac
+ab
+A
+ab
+A
+ab
+ac
+ab
+A
+A
+A
+ac
+ab
+A
+Δ
+ab
+ab
+A
+b
+A
+A
+ab
+ac
+A
+A
+ab
+A
+A
+ac
+ab
+A
+A
+A
+ac
+A
+ab
+A
+ab
+A
+A
+ac
+ab
+ac
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+A
+ab
+A
+A
+Ca
+A
+Ab
+A
+A
+A
+A
+A
+ab
+ab
+ac
+A
+A
+A
+A
+A
+Δ
+ab
+A
+ab
+Cb
+ac
+ab
+ac
+A
+ac
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+ab
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+A
+A
+ab
+A
+ac
+c
+ab
+ac
+A
+ab
+A
+ac
+ac
+A
+A
+Ca
+ac
+A
+ab
+A
+A
+ac
+ab
+ab
+A
+A
+ab
+A
+A
+A
+ab
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+ac
+A
+ab
+A
+ac
+ac
+A
+ab
+Cab
+A
+ab
+A
+A
+ac
+A
+Cb
+A
+A
+ac
+A
+A
+ac
+A
+Δ
+Δ
+A
+ac
+A
+ab
+A
+ab
+ab
+A
+A
+ac
+A
+A
+A
+A
+ab
+A
+ab
+ab
+Δ
+ab
+A
+ac
+ab
+A
@@ -0,0 +1,1001 @@
+event
+A
+ab
+ab
+ab
+A
+A
+Δ
+A
+A
+A
+A
+ab
+A
+a
+A
+ab
+A
+A
+A
+A
+A
+ab
+A
+A
+ab
+ac
+A
+A
+ab
+A
+ab
+ab
+A
+ab
+A
+A
+A
+ac
+ab
+A
+ac
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+ab
+A
+ab
+A
+A
+A
+ab
+ab
+A
+A
+ac
+ab
+A
+Ca
+A
+A
+A
+ab
+A
+A
+A
+A
+ab
+Δ
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+abc
+ab
+ab
+ac
+Δ
+ab
+A
+A
+ab
+A
+ab
+ab
+A
+ab
+ab
+A
+ab
+ab
+ab
+A
+A
+A
+Δ
+ab
+ab
+ac
+A
+A
+ab
+A
+ab
+A
+A
+A
+ab
+A
+A
+ab
+ab
+A
+A
+A
+A
+ab
+A
+A
+ab
+ab
+ab
+ab
+A
+Δ
+Δ
+ac
+ab
+ab
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+ab
+AC
+ab
+ab
+Δ
+Ca
+A
+ac
+A
+A
+A
+ab
+ab
+A
+A
+ab
+ab
+A
+ab
+c
+A
+ab
+ab
+A
+A
+A
+A
+A
+ab
+ab
+A
+A
+Δ
+A
+A
+A
+A
+A
+B
+A
+A
+A
+ab
+ab
+A
+A
+A
+A
+A
+A
+Δ
+ab
+ab
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+Bac
+AB
+ab
+A
+Cab
+Ac
+A
+A
+A
+A
+ac
+A
+A
+A
+AC
+ab
+A
+A
+A
+A
+A
+A
+ab
+Ba
+A
+A
+A
+ab
+A
+ab
+A
+Δ
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+ac
+A
+A
+ab
+ab
+ab
+A
+ab
+A
+ab
+Cab
+A
+Abc
+A
+ab
+ac
+A
+ab
+A
+ab
+ab
+A
+A
+A
+A
+A
+ab
+ab
+ab
+ab
+ab
+A
+A
+A
+A
+A
+A
+ac
+C
+ab
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+ab
+A
+A
+A
+A
+A
+ab
+ab
+ab
+A
+Δ
+A
+A
+A
+A
+ab
+ab
+A
+ab
+A
+ab
+ab
+ac
+ab
+ab
+A
+ab
+A
+ac
+A
+A
+A
+A
+B
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+BCa
+A
+ab
+A
+a
+ac
+A
+Δ
+ab
+A
+Ba
+ab
+ab
+A
+A
+A
+A
+ab
+ab
+A
+A
+B
+ab
+ab
+A
+ab
+ab
+A
+A
+A
+A
+c
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+ab
+A
+A
+A
+A
+bc
+A
+A
+A
+ab
+A
+ab
+A
+A
+ac
+ab
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+A
+abc
+A
+A
+ab
+ab
+ab
+Bac
+A
+A
+A
+ab
+A
+A
+A
+ab
+ab
+A
+A
+ab
+ab
+A
+A
+A
+ab
+A
+A
+A
+A
+ac
+ab
+ab
+A
+A
+A
+A
+A
+ab
+ac
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+ab
+A
+A
+A
+ac
+ab
+A
+A
+A
+Δ
+A
+A
+A
+A
+A
+B
+A
+a
+ac
+A
+ab
+A
+A
+A
+ab
+A
+ac
+A
+A
+A
+ab
+ab
+A
+B
+ac
+A
+ab
+Ca
+C
+A
+A
+ac
+A
+Δ
+ab
+A
+ab
+A
+A
+ab
+ab
+A
+A
+ac
+A
+ab
+A
+A
+A
+A
+A
+Ac
+ac
+Δ
+ab
+A
+A
+ab
+A
+ab
+ac
+A
+c
+ab
+Δ
+A
+ac
+A
+A
+A
+A
+A
+ab
+A
+ab
+A
+ac
+Δ
+ab
+A
+A
+Δ
+A
+ab
+A
+ac
+A
+A
+A
+ab
+A
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+Δ
+A
+Bac
+A
+A
+A
+A
+ab
+a
+ab
+b
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+ab
+ac
+ab
+ab
+ab
+AB
+ab
+A
+A
+ab
+A
+A
+ab
+A
+A
+ac
+A
+AC
+Δ
+ac
+A
+Ac
+ab
+A
+ab
+A
+A
+A
+ab
+A
+A
+A
+ab
+Δ
+ab
+Ba
+ab
+A
+A
+A
+ABC
+A
+A
+A
+A
+A
+A
+A
+ab
+ab
+A
+A
+Ba
+Δ
+A
+A
+ab
+ac
+A
+A
+A
+ab
+A
+A
+A
+A
+ab
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+AB
+B
+ab
+C
+A
+A
+A
+Bc
+ab
+A
+ab
+ab
+A
+A
+ab
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+ab
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+A
+ab
+A
+b
+AB
+A
+AB
+A
+A
+ab
+A
+A
+A
+ab
+ab
+A
+ac
+A
+A
+A
+A
+ab
+ac
+A
+A
+c
+ab
+A
+A
+A
+A
+A
+ab
+ab
+ab
+ab
+A
+ab
+A
+A
+A
+A
+ab
+A
+A
+A
+ab
+ab
+A
+A
+A
+A
+ab
+A
+A
+ab
+A
+A
+A
+A
+BCa
+ac
+ab
+ac
+A
+ab
+Δ
+Cb
+A
+A
+A
+ab
+A
+Δ
+A
+ab
+A
+AC
+ab
+ab
+ab
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+Δ
+A
+ABC
+A
+A
+Ac
+A
+A
+A
+ab
+A
+A
+A
+A
+ab
+ab
+bc
+ab
+ac
+ab
+A
+A
+a
+A
+ab
+A
+A
+Δ
+ab
+ab
+ab
+A
+ab
+ab
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+ab
+A
+ab
+A
+A
+ab
+A
+A
+ab
+ab
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+B
+A
+A
+BC
+ab
+ab
+A
+A
+A
+A
+ac
+B
+bc
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+ab
+A
+A
+ac
+ab
+A
+ab
+A
+A
+ab
+A
+A
+ab
+a
+A
+C
+A
+A
+ab
+A
+A
+A
+ab
+A
+ab
+A
+ac
+A
+A
+ab
+ab
+A
+A
+ab
+A
+BC
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+Cab
+A
+ab
+ab
+ab
+A
+A
+c
+A
+A
+A
+ab
+A
+A
+A
+A
+A
+A
+ac
+ac
+Δ
+ab
+A
+ac
+A
@@ -0,0 +1,1001 @@
+event
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+bc
+ac
+c
+A
+ac
+ac
+
+B
+A
+A
+A
+ac
+ac
+A
+A
+ac
+A
+A
+ab
+ab
+A
+A
+a
+A
+b
+ab
+A
+AB
+A
+ab
+ac
+ab
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+ab
+ac
+A
+ac
+ac
+B
+A
+A
+ac
+A
+A
+ac
+A
+A
+ab
+A
+ac
+A
+A
+A
+A
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+
+ac
+A
+ac
+A
+ac
+ac
+A
+A
+A
+A
+ac
+A
+bc
+A
+A
+ac
+ac
+A
+A
+Ba
+A
+ACb
+b
+ac
+A
+A
+ac
+A
+A
+a
+Ca
+ac
+A
+ac
+ac
+A
+A
+A
+ab
+ac
+A
+ACb
+A
+
+A
+A
+A
+
+A
+A
+A
+A
+c
+ac
+A
+A
+A
+A
+A
+ac
+A
+ac
+A
+A
+A
+ab
+ac
+A
+B
+ab
+ac
+ab
+A
+A
+A
+ac
+
+C
+c
+A
+ac
+ac
+A
+ab
+A
+A
+
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+ac
+A
+ac
+A
+A
+A
+
+Cb
+A
+A
+A
+ac
+A
+A
+A
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+ac
+ac
+A
+ac
+bc
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+Cab
+ac
+B
+Ca
+A
+A
+A
+ac
+A
+B
+A
+A
+A
+A
+A
+A
+A
+C
+A
+c
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+A
+A
+ab
+c
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+ac
+A
+ac
+ac
+A
+A
+A
+A
+ac
+A
+A
+
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+ac
+A
+A
+Ac
+A
+A
+A
+ab
+AB
+A
+A
+b
+A
+A
+ac
+ABC
+ac
+ac
+ABC
+b
+A
+ac
+ac
+ac
+A
+A
+A
+A
+A
+A
+ac
+A
+ac
+A
+A
+A
+ac
+A
+A
+ac
+A
+A
+ac
+A
+
+ac
+a
+A
+ab
+A
+ac
+A
+A
+Ac
+A
+A
+ac
+A
+A
+B
+ac
+A
+A
+A
+A
+A
+A
+ac
+A
+ac
+A
+ac
+ac
+ac
+A
+ab
+ac
+A
+A
+A
+ab
+A
+ac
+A
+ab
+A
+
+ac
+A
+bc
+A
+A
+ac
+ac
+ac
+A
+ac
+A
+A
+B
+A
+A
+A
+A
+ac
+A
+A
+A
+b
+A
+A
+A
+A
+A
+A
+A
+A
+A
+bc
+A
+ac
+ac
+A
+Ba
+ac
+A
+A
+A
+A
+A
+ac
+ac
+
+A
+Ca
+A
+ab
+A
+A
+A
+A
+A
+ac
+A
+ac
+A
+ac
+ac
+A
+A
+ac
+A
+A
+ac
+A
+A
+A
+A
+A
+ac
+ac
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+ac
+A
+A
+A
+A
+ac
+c
+A
+ac
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+ac
+ab
+ab
+A
+A
+Ca
+ac
+A
+A
+ac
+A
+ac
+A
+A
+A
+C
+A
+A
+ac
+A
+Cb
+A
+A
+ac
+A
+A
+ac
+A
+A
+ABC
+ab
+A
+ac
+A
+A
+A
+ac
+ac
+A
+A
+A
+ac
+C
+ac
+ac
+A
+ac
+ab
+ac
+A
+A
+ab
+ac
+A
+ab
+A
+A
+ab
+A
+A
+A
+A
+ac
+ac
+ac
+Abc
+A
+A
+A
+ac
+ab
+
+A
+ac
+a
+ab
+A
+A
+ab
+A
+A
+A
+A
+A
+A
+A
+bc
+A
+bc
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+ac
+ac
+A
+A
+ac
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+Ca
+ac
+ac
+ac
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+Ca
+ac
+A
+ab
+A
+ab
+AC
+ab
+A
+A
+A
+A
+A
+A
+
+A
+A
+ac
+A
+A
+ac
+A
+ABC
+A
+A
+A
+ac
+ac
+A
+ac
+A
+A
+A
+A
+A
+A
+ac
+ac
+A
+ac
+A
+A
+ac
+A
+A
+ac
+A
+ab
+A
+ac
+A
+A
+c
+A
+A
+A
+A
+A
+ac
+Cb
+A
+ac
+A
+A
+ac
+A
+A
+A
+A
+ac
+
+A
+ac
+A
+A
+ab
+A
+A
+A
+A
+A
+Bc
+A
+A
+ac
+A
+A
+ac
+A
+A
+Ca
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+ab
+A
+A
+A
+ac
+A
+A
+ac
+A
+ac
+C
+ab
+A
+A
+A
+b
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+AC
+ac
+A
+A
+ac
+A
+A
+A
+A
+ABC
+A
+A
+ac
+
+ac
+A
+A
+ac
+A
+C
+A
+A
+A
+A
+
+ac
+A
+ac
+A
+
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+a
+ab
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+ac
+A
+A
+ab
+A
+ac
+A
+A
+b
+A
+Bac
+A
+ac
+A
+A
+ac
+A
+A
+ac
+A
+A
+Abc
+A
+ac
+ab
+A
+A
+b
+A
+A
+A
+A
+A
+ab
+ac
+
+A
+A
+ac
+A
+A
+A
+A
+ac
+A
+
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+A
+A
+
+ac
+A
+A
+ac
+Ab
+A
+A
+ac
+Ac
+ab
+A
+A
+A
+ac
+ac
+ab
+A
+ac
+A
+A
+A
+A
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+
+A
+A
+ab
+A
+A
+A
+A
+ac
+ac
+AC
+A
+A
+A
+A
+A
+ac
+A
+A
+ac
+A
+ABC
+ac
+A
+A
+A
+A
+A
+A
+ac
+ac
+A
+A
+ac
+ab
+ab
+bc
+A
+ac
+
+ac
+A
+ab
+A
+Ac
+A
+ab
+A
+A
+AC
+A
+A
+ac
+A
+ac
+A
+A
+
+A
+ac
+ac
+ac
+A
+A
+ac
+A
+A
+ab
+A
+A
+ab
+A
+A
+ac
+A
+a
+ac
+A
+A
+A
+ac
+Cb
+A
+ab
+ac
+A
+BCa
+ac
+ac
+A
@@ -0,0 +1,101 @@
+event
+ac
+bA
+bAC
+bAC
+cA
+B
+c
+BC
+acB
+
+cAB
+cA
+abC
+C
+bC
+bcA
+cB
+cAB
+
+bcA
+BC
+C
+c
+abc
+aC
+AB
+bc
+ABC
+ac
+AC
+aBC
+aB
+AB
+A
+cAB
+AB
+B
+B
+bC
+bcA
+A
+a
+cAB
+aC
+AC
+bAC
+bAC
+cA
+bC
+AB
+bc
+abC
+bAC
+B
+AB
+bc
+B
+b
+bcA
+AC
+bA
+a
+A
+ac
+bC
+bA
+bcA
+cB
+A
+a
+bcA
+
+bc
+bAC
+bcA
+aBC
+abc
+b
+a
+cA
+aB
+
+cAB
+b
+aBC
+A
+aC
+ab
+bcA
+bC
+cA
+
+aB
+bcA
+AB
+ac
+aBC
+acB
+C
+A
@@ -0,0 +1,77 @@
+using CSV
+using DataFrames
+using Latexify
+using Plots
+using StatsPlots
+
+function format_event(event)
+    symbols = []
+    if event == "z"
+        symbols = "\\bot"
+    elseif event == "i"
+        symbols = "\\indepclass"
+    else
+        for c in event
+            if isuppercase(c)
+                push!(symbols, "\\co{$(lowercase(c))}")
+            else
+                push!(symbols, c)
+            end
+        end
+    end
+    return join(symbols, "")
+end
+
+cores = Dict(
+    "indepclass" => [ "B", "C", "bc", "Ba", "Bc", "BC", "Ca", "Cb", "BCa" ],
+    "A" => [ "A", "Ab", "Ac", "AB", "AC", "Abc", "ACb", "ABc", "ABC" ],
+    "ab" => [ "b", "ab", "Cab" ],
+    "ac" => [ "c", "ac", "Bac" ],
+    "A,ab" => [  ],
+    "A,ac" => [  ],
+    "ab,ac" => [ "a", "abc" ],
+    "A,ab,ac" => [ "Δ" ],
+)
+
+core_of = Dict{String,String}(
+    "B" => "i", # independent class
+    "C" => "i", 
+    "bc" => "i", 
+    "Ba" => "i", 
+    "Bc" => "i", 
+    "BC" => "i", 
+    "Ca" => "i", 
+    "Cb" => "i", 
+    "BCa" => "i",
+    "A" => "A",
+    "Ab" => "A",
+    "Ac" => "A",
+    "AB" => "A",
+    "AC" => "A",
+    "Abc" => "A",
+    "ACb" => "A",
+    "ABc" => "A",
+    "ABC" => "A",
+    "b" => "ab",
+    "ab" => "ab",
+    "Cab" => "ab",
+    "c" => "ac",
+    "ac" => "ac",
+    "Bac" => "ac",
+    "a" => "ab,ac",
+    "abc" => "ab,ac",
+    "Δ" => "A,ab,ac",
+)
+
+α = 0.1
+β = 0.3
+γ = 0.8
+
+data = CSV.read("model_sample-[$α;$β;$γ].csv", DataFrame)
+transform!(data, :event => ByRow(e -> get(core_of, e, "z")) => :core)
+groups = groupby(data, :core)
+totals = combine(groups, :core => length => :count)
+for row in copy.(eachrow(totals))
+    row_fmt = format_event(row.core)
+    println("$(row_fmt) & $(row.count)")
+end
 \ No newline at end of file
@@ -0,0 +1,91 @@
+using Random
+using Combinatorics
+
+neg(a) = islowercase(a) ? uppercase(a) : lowercase(a)
+
+function random_sample(atoms)
+    result = Char[]
+    for a in atoms
+        if rand(Bool)
+            push!(result, rand(Bool) ? a : neg(a))
+        end
+    end
+    return length(result) > 0 ? join(result) : "Δ"
+end
+
+function random_sample(n::Int, events)
+    result = String[]
+    for _ in 1:n
+        x = rand(events)
+        push!(result, join(x, ""))
+    end
+    return result
+end
+
+function model_sample(n::Int; α=0.01, β=0.3,  γ=0.2)
+    literals = [ 'a', 'b', 'c', 'A', 'B', 'C' ]
+    samples = []
+    for _ in 1:n
+        sample = []
+        if rand() <= α # RCE
+            sample = Random.shuffle(literals)[1:rand(0:3)]
+            while !is_consistent(sample)
+                sample = Random.shuffle(literals)[1:rand(0:3)]
+            end
+            if sample == []
+                push!(sample , 'Δ')
+            else
+                sort!(sample)
+            end
+        elseif rand() <= β
+            # a
+            push!(sample, "a")
+            if rand() <= γ
+                push!(sample, "b")
+            else
+                push!(sample, "c")
+            end
+        else
+            # A
+            push!(sample, "A")
+        end
+        push!(samples, join(sort(sample), ""))
+    end
+    return samples
+end
+
+
+function is_consistent(event)
+    for (i,x) in enumerate(event)
+        for y in event[i:end]
+            if y == neg(x)
+                return false
+            end
+        end
+    end 
+    return true
+end
+
+using DelimitedFiles
+
+literals = [ 'a', 'b', 'c', 'A', 'B', 'C' ]
+omega = collect(powerset(literals))
+consistent = filter(is_consistent, omega)
+consistent = join.(consistent)
+
+ 
+α = 0.1
+β = 0.3
+γ = 0.5
+
+events = model_sample(1000; 
+    α = α,
+    β = β,
+    γ = γ)
+
+println( events )
+
+open("model_sample-[$α;$β;$γ].csv", "w") do io
+    writedlm(io, [ "event" ])
+    writedlm(io,  events)
+end
 \ No newline at end of file
@@ -0,0 +1,67 @@		@@ -0,0 +1,67 @@
	1	+using Symbolics
	2	+using Latexify
	3	+
	4	+# # γ = 0.5
	5	+# no_s = [
	6	+# 0 ,
	7	+# 23 ,
	8	+# 614 ,
	9	+# 165 ,
	10	+# 169 ,
	11	+# 0 ,
	12	+# 0 ,
	13	+# 4 ,
	14	+# 25 ,
	15	+# ]
	16	+
	17	+# γ = 0.8
	18	+no_s = [
	19	+ 0 ,
	20	+ 28 ,
	21	+ 632 ,
	22	+ 246 ,
	23	+ 59 ,
	24	+ 0 ,
	25	+ 0 ,
	26	+ 5 ,
	27	+ 27 ,
	28	+]
	29	+pr_s = (x -> x // 1000).(no_s)
	30	+
	31	+@variables θ
	32	+num_e = [
	33	+ 0 ,
	34	+ 0 ,
	35	+ 7 ,
	36	+ 3 * θ ,
	37	+ 3 * (1 - θ) ,
	38	+ 0 ,
	39	+ 0 ,
	40	+ 3 ,
	41	+ 10 ,
	42	+]
	43	+pr_e = (x -> x // 23).(num_e)
	44	+
	45	+target = expand(sum( (x -> x^2).(pr_s - pr_e) ))
	46	+println(latexify(target))
	47	+
	48	+# using Plots
	49	+# g(t) = (20869963/66125000) + (477/52900)t + (18/529)(t^2)
	50	+# t = 0:0.1:1
	51	+# plot(t, g.(t))
	52	+
	53	+
	54	+
	55	+function solve2(a, b, c)
	56	+ delta = sqrt(Complex(b^2 - 4 * a * c))
	57	+ return ( (-b - delta)/(2a), (-b + delta)/(2a) )
	58	+end
	59	+
	60	+# # g = target ~ 0
	61	+
	62	+a = float(18//529)
	63	+b = -float(21903//264500)
	64	+c = float(188207311//529000000)
	65	+
	66	+println("a: $a b: $b c: $c")
	67	+println("̂θ = $(-b/a)")
@@ -0,0 +1,1001 @@		@@ -0,0 +1,1001 @@
	1	+event
	2	+A
	3	+A
	4	+A
	5	+A
	6	+A
	7	+A
	8	+A
	9	+ab
	10	+A
	11	+A
	12	+A
	13	+A
	14	+A
	15	+ac
	16	+A
	17	+A
	18	+ac
	19	+A
	20	+A
	21	+A
	22	+ab
	23	+A
	24	+A
	25	+A
	26	+A
	27	+A
	28	+A
	29	+ac
	30	+A
	31	+A
	32	+A
	33	+ac
	34	+A
	35	+A
	36	+ac
	37	+ac
	38	+A
	39	+ac
	40	+A
	41	+A
	42	+ac
	43	+A
	44	+A
	45	+A
	46	+A
	47	+A
	48	+A
	49	+A
	50	+A
	51	+A
	52	+ac
	53	+ab
	54	+A
	55	+A
	56	+A
	57	+A
	58	+A
	59	+ac
	60	+A
	61	+A
	62	+A
	63	+ac
	64	+A
	65	+A
	66	+A
	67	+ac
	68	+A
	69	+A
	70	+A
	71	+ac
	72	+A
	73	+A
	74	+A
	75	+A
	76	+ac
	77	+A
	78	+A
	79	+A
	80	+ab
	81	+A
	82	+A
	83	+A
	84	+ab
	85	+A
	86	+ab
	87	+A
	88	+A
	89	+A
	90	+ab
	91	+A
	92	+A
	93	+A
	94	+ac
	95	+ac
	96	+A
	97	+ab
	98	+A
	99	+A
	100	+ac
	101	+A
	102	+A
	103	+A
	104	+A
	105	+ac
	106	+A
	107	+A
	108	+A
	109	+A
	110	+A
	111	+A
	112	+ac
	113	+A
	114	+ac
	115	+A
	116	+A
	117	+ac
	118	+ac
	119	+A
	120	+ac
	121	+ac
	122	+A
	123	+A
	124	+A
	125	+ac
	126	+A
	127	+A
	128	+A
	129	+ac
	130	+A
	131	+ac
	132	+A
	133	+ab
	134	+ac
	135	+A
	136	+ac
	137	+A
	138	+A
	139	+ab
	140	+A
	141	+ac
	142	+A
	143	+ac
	144	+ab
	145	+ac
	146	+A
	147	+A
	148	+ac
	149	+A
	150	+A
	151	+A
	152	+A
	153	+A
	154	+A
	155	+ac
	156	+A
	157	+A
	158	+A
	159	+A
	160	+A
	161	+ac
	162	+A
	163	+A
	164	+ac
	165	+A
	166	+A
	167	+A
	168	+A
	169	+A
	170	+ac
	171	+A
	172	+A
	173	+A
	174	+A
	175	+A
	176	+A
	177	+A
	178	+A
	179	+ac
	180	+A
	181	+A
	182	+ac
	183	+A
	184	+A
	185	+A
	186	+ac
	187	+ac
	188	+ac
	189	+A
	190	+A
	191	+ac
	192	+A
	193	+A
	194	+A
	195	+A
	196	+A
	197	+A
	198	+A
	199	+A
	200	+ab
	201	+A
	202	+A
	203	+A
	204	+A
	205	+ac
	206	+ac
	207	+A
	208	+ab
	209	+A
	210	+A
	211	+ac
	212	+ab
	213	+A
	214	+A
	215	+A
	216	+A
	217	+ac
	218	+ac
	219	+ac
	220	+A
	221	+ac
	222	+A
	223	+A
	224	+ac
	225	+A
	226	+A
	227	+ab
	228	+ac
	229	+A
	230	+ac
	231	+ac
	232	+A
	233	+A
	234	+ac
	235	+ac
	236	+A
	237	+ac
	238	+A
	239	+A
	240	+A
	241	+ac
	242	+A
	243	+A
	244	+A
	245	+A
	246	+A
	247	+A
	248	+A
	249	+A
	250	+A
	251	+A
	252	+A
	253	+A
	254	+A
	255	+A
	256	+ab
	257	+ac
	258	+A
	259	+A
	260	+ac
	261	+A
	262	+A
	263	+A
	264	+A
	265	+A
	266	+ab
	267	+ac
	268	+A
	269	+ac
	270	+A
	271	+ac
	272	+A
	273	+ac
	274	+A
	275	+A
	276	+ab
	277	+A
	278	+A
	279	+ac
	280	+ac
	281	+A
	282	+A
	283	+ac
	284	+A
	285	+A
	286	+A
	287	+A
	288	+ac
	289	+A
	290	+A
	291	+ab
	292	+ac
	293	+A
	294	+A
	295	+ac
	296	+A
	297	+A
	298	+A
	299	+ab
	300	+A
	301	+ac
	302	+ab
	303	+A
	304	+A
	305	+ac
	306	+A
	307	+A
	308	+A
	309	+A
	310	+A
	311	+A
	312	+ac
	313	+ac
	314	+A
	315	+A
	316	+A
	317	+ac
	318	+A
	319	+A
	320	+A
	321	+A
	322	+A
	323	+A
	324	+A
	325	+A
	326	+A
	327	+A
	328	+A
	329	+A
	330	+A
	331	+A
	332	+A
	333	+A
	334	+A
	335	+ac
	336	+A
	337	+ac
	338	+A
	339	+ac
	340	+ac
	341	+ac
	342	+A
	343	+A
	344	+ac
	345	+A
	346	+A
	347	+A
	348	+A
	349	+ac
	350	+ac
	351	+A
	352	+A
	353	+A
	354	+A
	355	+A
	356	+A
	357	+ac
	358	+A
	359	+ab
	360	+A
	361	+A
	362	+ab
	363	+ac
	364	+ab
	365	+ac
	366	+ac
	367	+A
	368	+A
	369	+ab
	370	+ac
	371	+A
	372	+ac
	373	+A
	374	+A
	375	+A
	376	+A
	377	+ac
	378	+ac
	379	+A
	380	+A
	381	+ac
	382	+A
	383	+A
	384	+A
	385	+A
	386	+ab
	387	+A
	388	+A
	389	+ac
	390	+ac
	391	+A
	392	+A
	393	+A
	394	+A
	395	+A
	396	+ac
	397	+A
	398	+A
	399	+A
	400	+A
	401	+A
	402	+A
	403	+A
	404	+ac
	405	+A
	406	+A
	407	+ac
	408	+ac
	409	+A
	410	+A
	411	+A
	412	+A
	413	+ac
	414	+ac
	415	+A
	416	+ac
	417	+A
	418	+ac
	419	+A
	420	+A
	421	+A
	422	+A
	423	+ac
	424	+A
	425	+A
	426	+A
	427	+A
	428	+ac
	429	+A
	430	+A
	431	+A
	432	+A
	433	+A
	434	+A
	435	+A
	436	+A
	437	+A
	438	+ac
	439	+A
	440	+A
	441	+A
	442	+A
	443	+A
	444	+A
	445	+A
	446	+ac
	447	+ac
	448	+A
	449	+ac
	450	+A
	451	+A
	452	+A
	453	+A
	454	+ac
	455	+A
	456	+ac
	457	+A
	458	+ac
	459	+A
	460	+A
	461	+A
	462	+A
	463	+ab
	464	+A
	465	+ac
	466	+A
	467	+A
	468	+ab
	469	+A
	470	+ac
	471	+A
	472	+ac
	473	+ac
	474	+ac
	475	+A
	476	+A
	477	+A
	478	+A
	479	+A
	480	+A
	481	+A
	482	+A
	483	+A
	484	+ac
	485	+A
	486	+A
	487	+A
	488	+A
	489	+A
	490	+ac
	491	+ab
	492	+ab
	493	+A
	494	+A
	495	+A
	496	+ac
	497	+A
	498	+A
	499	+ac
	500	+A
	501	+A
	502	+A
	503	+A
	504	+A
	505	+A
	506	+A
	507	+A
	508	+A
	509	+A
	510	+A
	511	+ab
	512	+A
	513	+ac
	514	+A
	515	+A
	516	+A
	517	+A
	518	+A
	519	+A
	520	+A
	521	+A
	522	+ac
	523	+A
	524	+A
	525	+A
	526	+A
	527	+A
	528	+A
	529	+A
	530	+ac
	531	+A
	532	+A
	533	+A
	534	+A
	535	+ac
	536	+A
	537	+ac
	538	+ac
	539	+A
	540	+A
	541	+A
	542	+A
	543	+A
	544	+A
	545	+A
	546	+A
	547	+A
	548	+ab
	549	+A
	550	+A
	551	+ab
	552	+A
	553	+ac
	554	+A
	555	+A
	556	+A
	557	+ab
	558	+A
	559	+ac
	560	+ab
	561	+A
	562	+A
	563	+ac
	564	+A
	565	+ab
	566	+A
	567	+ac
	568	+A
	569	+A
	570	+A
	571	+A
	572	+A
	573	+ac
	574	+A
	575	+ac
	576	+ab
	577	+ac
	578	+A
	579	+A
	580	+A
	581	+ac
	582	+A
	583	+A
	584	+ac
	585	+A
	586	+ac
	587	+ab
	588	+A
	589	+ab
	590	+ac
	591	+A
	592	+ac
	593	+ac
	594	+ac
	595	+A
	596	+A
	597	+ab
	598	+A
	599	+ab
	600	+A
	601	+A
	602	+ac
	603	+A
	604	+ac
	605	+A
	606	+ac
	607	+A
	608	+A
	609	+A
	610	+A
	611	+A
	612	+A
	613	+A
	614	+ac
	615	+A
	616	+A
	617	+ac
	618	+ac
	619	+ab
	620	+A
	621	+A
	622	+A
	623	+A
	624	+ac
	625	+ac
	626	+A
	627	+A
	628	+A
	629	+A
	630	+ab
	631	+A
	632	+A
	633	+A
	634	+A
	635	+A
	636	+A
	637	+A
	638	+ac
	639	+ac
	640	+A
	641	+A
	642	+A
	643	+A
	644	+A
	645	+A
	646	+A
	647	+A
	648	+ab
	649	+A
	650	+ac
	651	+A
	652	+A
	653	+A
	654	+ac
	655	+A
	656	+ac
	657	+ac
	658	+ac
	659	+A
	660	+A
	661	+A
	662	+A
	663	+A
	664	+ab
	665	+ac
	666	+ac
	667	+A
	668	+A
	669	+A
	670	+A
	671	+A
	672	+A
	673	+A
	674	+ab
	675	+A
	676	+A
	677	+ac
	678	+ac
	679	+A
	680	+A
	681	+ac
	682	+A
	683	+A
	684	+A
	685	+A
	686	+A
	687	+ab
	688	+ac
	689	+A
	690	+A
	691	+ac
	692	+A
	693	+A
	694	+A
	695	+ab
	696	+A
	697	+ab
	698	+A
	699	+A
	700	+A
	701	+A
	702	+A
	703	+A
	704	+A
	705	+ac
	706	+A
	707	+A
	708	+A
	709	+A
	710	+A
	711	+A
	712	+A
	713	+A
	714	+A
	715	+ac
	716	+A
	717	+ac
	718	+A
	719	+A
	720	+ac
	721	+ac
	722	+ac
	723	+A
	724	+A
	725	+A
	726	+ac
	727	+A
	728	+ac
	729	+A
	730	+A
	731	+A
	732	+A
	733	+ab
	734	+ac
	735	+ac
	736	+ac
	737	+A
	738	+A
	739	+A
	740	+A
	741	+A
	742	+A
	743	+A
	744	+A
	745	+ac
	746	+ac
	747	+A
	748	+A
	749	+ac
	750	+A
	751	+A
	752	+A
	753	+ab
	754	+A
	755	+ac
	756	+A
	757	+ac
	758	+A
	759	+A
	760	+A
	761	+A
	762	+A
	763	+A
	764	+A
	765	+ac
	766	+A
	767	+A
	768	+A
	769	+A
	770	+ab
	771	+ab
	772	+A
	773	+ac
	774	+A
	775	+A
	776	+A
	777	+A
	778	+A
	779	+A
	780	+A
	781	+A
	782	+A
	783	+ac
	784	+A
	785	+A
	786	+ac
	787	+A
	788	+A
	789	+A
	790	+ac
	791	+A
	792	+A
	793	+ac
	794	+ac
	795	+A
	796	+ac
	797	+A
	798	+A
	799	+A
	800	+A
	801	+A
	802	+A
	803	+A
	804	+A
	805	+A
	806	+A
	807	+A
	808	+A
	809	+A
	810	+A
	811	+A
	812	+A
	813	+A
	814	+A
	815	+A
	816	+ac
	817	+A
	818	+A
	819	+A
	820	+ac
	821	+ac
	822	+A
	823	+A
	824	+ac
	825	+A
	826	+ab
	827	+A
	828	+ac
	829	+A
	830	+A
	831	+ac
	832	+ac
	833	+A
	834	+A
	835	+ab
	836	+A
	837	+A
	838	+A
	839	+A
	840	+ac
	841	+A
	842	+ac
	843	+ac
	844	+A
	845	+A
	846	+ab
	847	+A
	848	+A
	849	+ac
	850	+A
	851	+A
	852	+A
	853	+ac
	854	+A
	855	+A
	856	+A
	857	+A
	858	+A
	859	+A
	860	+A
	861	+A
	862	+A
	863	+ac
	864	+A
	865	+A
	866	+A
	867	+A
	868	+ac
	869	+A
	870	+A
	871	+ab
	872	+ab
	873	+ac
	874	+A
	875	+A
	876	+ac
	877	+A
	878	+A
	879	+ac
	880	+ac
	881	+A
	882	+ac
	883	+A
	884	+A
	885	+ac
	886	+ab
	887	+A
	888	+A
	889	+A
	890	+ac
	891	+ac
	892	+ac
	893	+A
	894	+A
	895	+A
	896	+ac
	897	+ac
	898	+ac
	899	+A
	900	+A
	901	+A
	902	+A
	903	+A
	904	+A
	905	+A
	906	+A
	907	+A
	908	+A
	909	+A
	910	+A
	911	+ac
	912	+A
	913	+ac
	914	+A
	915	+A
	916	+A
	917	+ac
	918	+A
	919	+A
	920	+A
	921	+ac
	922	+A
	923	+A
	924	+A
	925	+ac
	926	+A
	927	+A
	928	+A
	929	+ac
	930	+ac
	931	+ac
	932	+ab
	933	+A
	934	+ac
	935	+A
	936	+A
	937	+A
	938	+A
	939	+A
	940	+ac
	941	+ab
	942	+A
	943	+A
	944	+A
	945	+A
	946	+ac
	947	+A
	948	+A
	949	+ac
	950	+A
	951	+ab
	952	+A
	953	+A
	954	+ac
	955	+A
	956	+A
	957	+A
	958	+A
	959	+ac
	960	+A
	961	+A
	962	+ac
	963	+A
	964	+A
	965	+ac
	966	+ac
	967	+ac
	968	+A
	969	+A
	970	+A
	971	+A
	972	+A
	973	+A
	974	+A
	975	+A
	976	+ac
	977	+ac
	978	+ac
	979	+ac
	980	+A
	981	+A
	982	+A
	983	+A
	984	+A
	985	+ac
	986	+A
	987	+A
	988	+A
	989	+A
	990	+A
	991	+ac
	992	+A
	993	+ac
	994	+A
	995	+ac
	996	+A
	997	+A
	998	+A
	999	+A
	1000	+A
	1001	+A
@@ -0,0 +1,101 @@		@@ -0,0 +1,101 @@
	1	+event
	2	+ac
	3	+bA
	4	+bAC
	5	+bAC
	6	+cA
	7	+B
	8	+c
	9	+BC
	10	+acB
	11	+
	12	+cAB
	13	+cA
	14	+abC
	15	+C
	16	+bC
	17	+bcA
	18	+cB
	19	+cAB
	20	+
	21	+bcA
	22	+BC
	23	+C
	24	+c
	25	+abc
	26	+aC
	27	+AB
	28	+bc
	29	+ABC
	30	+ac
	31	+AC
	32	+aBC
	33	+aB
	34	+AB
	35	+A
	36	+cAB
	37	+AB
	38	+B
	39	+B
	40	+bC
	41	+bcA
	42	+A
	43	+a
	44	+cAB
	45	+aC
	46	+AC
	47	+bAC
	48	+bAC
	49	+cA
	50	+bC
	51	+AB
	52	+bc
	53	+abC
	54	+bAC
	55	+B
	56	+AB
	57	+bc
	58	+B
	59	+b
	60	+bcA
	61	+AC
	62	+bA
	63	+a
	64	+A
	65	+ac
	66	+bC
	67	+bA
	68	+bcA
	69	+cB
	70	+A
	71	+a
	72	+bcA
	73	+
	74	+bc
	75	+bAC
	76	+bcA
	77	+aBC
	78	+abc
	79	+b
	80	+a
	81	+cA
	82	+aB
	83	+
	84	+cAB
	85	+b
	86	+aBC
	87	+A
	88	+aC
	89	+ab
	90	+bcA
	91	+bC
	92	+cA
	93	+
	94	+aB
	95	+bcA
	96	+AB
	97	+ac
	98	+aBC
	99	+acB
	100	+C
	101	+A
	@@ -0,0 +1,77 @@		@@ -0,0 +1,77 @@
		1	+using CSV
		2	+using DataFrames
		3	+using Latexify
		4	+using Plots
		5	+using StatsPlots
		6	+
		7	+function format_event(event)
		8	+ symbols = []
		9	+ if event == "z"
		10	+ symbols = "\\bot"
		11	+ elseif event == "i"
		12	+ symbols = "\\indepclass"
		13	+ else
		14	+ for c in event
		15	+ if isuppercase(c)
		16	+ push!(symbols, "\\co{$(lowercase(c))}")
		17	+ else
		18	+ push!(symbols, c)
		19	+ end
		20	+ end
		21	+ end
		22	+ return join(symbols, "")
		23	+end
		24	+
		25	+cores = Dict(
		26	+ "indepclass" => [ "B", "C", "bc", "Ba", "Bc", "BC", "Ca", "Cb", "BCa" ],
		27	+ "A" => [ "A", "Ab", "Ac", "AB", "AC", "Abc", "ACb", "ABc", "ABC" ],
		28	+ "ab" => [ "b", "ab", "Cab" ],
		29	+ "ac" => [ "c", "ac", "Bac" ],
		30	+ "A,ab" => [ ],
		31	+ "A,ac" => [ ],
		32	+ "ab,ac" => [ "a", "abc" ],
		33	+ "A,ab,ac" => [ "Δ" ],
		34	+)
		35	+
		36	+core_of = Dict{String,String}(
		37	+ "B" => "i", # independent class
		38	+ "C" => "i",
		39	+ "bc" => "i",
		40	+ "Ba" => "i",
		41	+ "Bc" => "i",
		42	+ "BC" => "i",
		43	+ "Ca" => "i",
		44	+ "Cb" => "i",
		45	+ "BCa" => "i",
		46	+ "A" => "A",
		47	+ "Ab" => "A",
		48	+ "Ac" => "A",
		49	+ "AB" => "A",
		50	+ "AC" => "A",
		51	+ "Abc" => "A",
		52	+ "ACb" => "A",
		53	+ "ABc" => "A",
		54	+ "ABC" => "A",
		55	+ "b" => "ab",
		56	+ "ab" => "ab",
		57	+ "Cab" => "ab",
		58	+ "c" => "ac",
		59	+ "ac" => "ac",
		60	+ "Bac" => "ac",
		61	+ "a" => "ab,ac",
		62	+ "abc" => "ab,ac",
		63	+ "Δ" => "A,ab,ac",
		64	+)
		65	+
		66	+α = 0.1
		67	+β = 0.3
		68	+γ = 0.8
		69	+
		70	+data = CSV.read("model_sample-[$α;$β;$γ].csv", DataFrame)
		71	+transform!(data, :event => ByRow(e -> get(core_of, e, "z")) => :core)
		72	+groups = groupby(data, :core)
		73	+totals = combine(groups, :core => length => :count)
		74	+for row in copy.(eachrow(totals))
		75	+ row_fmt = format_event(row.core)
		76	+ println("$(row_fmt) & $(row.count)")
		77	+end
0	\ No newline at end of file	78	\ No newline at end of file
	@@ -0,0 +1,91 @@		@@ -0,0 +1,91 @@
		1	+using Random
		2	+using Combinatorics
		3	+
		4	+neg(a) = islowercase(a) ? uppercase(a) : lowercase(a)
		5	+
		6	+function random_sample(atoms)
		7	+ result = Char[]
		8	+ for a in atoms
		9	+ if rand(Bool)
		10	+ push!(result, rand(Bool) ? a : neg(a))
		11	+ end
		12	+ end
		13	+ return length(result) > 0 ? join(result) : "Δ"
		14	+end
		15	+
		16	+function random_sample(n::Int, events)
		17	+ result = String[]
		18	+ for _ in 1:n
		19	+ x = rand(events)
		20	+ push!(result, join(x, ""))
		21	+ end
		22	+ return result
		23	+end
		24	+
		25	+function model_sample(n::Int; α=0.01, β=0.3, γ=0.2)
		26	+ literals = [ 'a', 'b', 'c', 'A', 'B', 'C' ]
		27	+ samples = []
		28	+ for _ in 1:n
		29	+ sample = []
		30	+ if rand() <= α # RCE
		31	+ sample = Random.shuffle(literals)[1:rand(0:3)]
		32	+ while !is_consistent(sample)
		33	+ sample = Random.shuffle(literals)[1:rand(0:3)]
		34	+ end
		35	+ if sample == []
		36	+ push!(sample , 'Δ')
		37	+ else
		38	+ sort!(sample)
		39	+ end
		40	+ elseif rand() <= β
		41	+ # a
		42	+ push!(sample, "a")
		43	+ if rand() <= γ
		44	+ push!(sample, "b")
		45	+ else
		46	+ push!(sample, "c")
		47	+ end
		48	+ else
		49	+ # A
		50	+ push!(sample, "A")
		51	+ end
		52	+ push!(samples, join(sort(sample), ""))
		53	+ end
		54	+ return samples
		55	+end
		56	+
		57	+
		58	+function is_consistent(event)
		59	+ for (i,x) in enumerate(event)
		60	+ for y in event[i:end]
		61	+ if y == neg(x)
		62	+ return false
		63	+ end
		64	+ end
		65	+ end
		66	+ return true
		67	+end
		68	+
		69	+using DelimitedFiles
		70	+
		71	+literals = [ 'a', 'b', 'c', 'A', 'B', 'C' ]
		72	+omega = collect(powerset(literals))
		73	+consistent = filter(is_consistent, omega)
		74	+consistent = join.(consistent)
		75	+
		76	+
		77	+α = 0.1
		78	+β = 0.3
		79	+γ = 0.5
		80	+
		81	+events = model_sample(1000;
		82	+ α = α,
		83	+ β = β,
		84	+ γ = γ)
		85	+
		86	+println( events )
		87	+
		88	+open("model_sample-[$α;$β;$γ].csv", "w") do io
		89	+ writedlm(io, [ "event" ])
		90	+ writedlm(io, events)
		91	+end
0	\ No newline at end of file	92	\ No newline at end of file