altered spawning model and changed temporal probability on

weight_success

1 files changed, 13 insertions, 13 deletions

diff --git a/experiments/process.py b/experiments/process.py
index 68c3507..a432dab 100644
--- a/experiments/process.py
+++ b/experiments/process.py
@@ -9,18 +9,18 @@ from math import exp
 if __name__ == "__main__":
     if len(sys.argv) < 2:
         sys.exit("usage: {0} <file>".format(sys.argv[0]))
-
     root_victims, victims, non_victims, age = load(open(sys.argv[1]))
-    alphas = 1. / np.arange(0., 1000., 100.)  # parameter of the time component
-    deltas = np.arange(0.01, 0.5, 0.1)  # parameter of the structural component
-    with open("out.log", "w") as fh:
-        for alpha, delta in product(alphas, deltas):
-            beta, roots, ll = ml(root_victims, victims, non_victims, age, alpha, delta)
-            print "\t".join(map(str, [1/alpha, delta, beta, roots, ll]))
-            fh.write("\t".join(map(str, [alpha, delta, beta, roots, ll])) + "\n")
-            fh.flush()
+    
+    # alphas = 1. / np.arange(1000., 10000., 50.)  # parameter of the time component
+    # deltas = np.arange(0.01, 0.5, 0.03)  # parameter of the structural component
+    # with open("out.log", "w") as fh:
+    #     for alpha, delta in product(alphas, deltas):
+    #         beta, roots, ll = ml(root_victims, victims, non_victims, age, alpha, delta)
+    #         print "\t".join(map(str, [1/alpha, delta, beta, roots, ll]))
+    #         fh.write("\t".join(map(str, [alpha, delta, beta, roots, ll])) + "\n")
+    #         fh.flush()
 
-    # alpha = 100.
-    # delta = 0.2
-    # beta, roots, ll = ml(root_victims, victims, non_victims, age, alpha, delta)
-    # print "\t".join(map(str, [1./alpha, delta, beta, roots, ll]))
\ No newline at end of file
+    alpha = 1/10.
+    delta = 0.28
+    beta, roots, ll = ml(root_victims, victims, non_victims, age, alpha, delta)
+    print "\t".join(map(str, [1./alpha, delta, beta, roots, ll]))
\ No newline at end of file