tried generating cascades and recovering parameters, still got the same

result where optimization minimizes alpha, delta as much as possible

1 files changed, 5 insertions, 4 deletions

diff --git a/experiments/process.py b/experiments/process.py
index b5e0d6d..ff722c8 100644
--- a/experiments/process.py
+++ b/experiments/process.py
@@ -11,9 +11,10 @@ if __name__ == "__main__":
         sys.exit("usage: {0} <file>".format(sys.argv[0]))
     root_victims, victims, non_victims, age = load(open(sys.argv[1]))
     
-    # alphas = 1. / np.arange(1000., 10000., 50.)  # parameter of the time component
-    alphas = np.logspace(-2,2,num=5)
-    deltas = np.logspace(-10,-1,num=10)#np.arange(0.000001, 0.005, 0.001)  # parameter of the structural component
+    # alphas = 1. / np.arange(1., 2000., 30.)  # parameter of the time component
+    alphas = np.logspace(-3,-.1,num=50)
+    # deltas = np.arange(0.005, 0.5, 0.005)  # parameter of the structural component
+    deltas = np.logspace(-3,-.5,num=50)
     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)
@@ -21,7 +22,7 @@ if __name__ == "__main__":
             fh.write("\t".join(map(str, [alpha, delta, beta, roots, ll])) + "\n")
             fh.flush()
 
-    # alpha = 20000.
+    # alpha = 1/10.
     # delta = .5
     # 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