From f2891c93b96388442d44512e6c43b092153f8c25 Mon Sep 17 00:00:00 2001 From: Ben Green Date: Sat, 20 Jun 2015 10:13:23 -0400 Subject: tried generating cascades and recovering parameters, still got the same result where optimization minimizes alpha, delta as much as possible --- experiments/process.py | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) (limited to 'experiments/process.py') 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} ".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 -- cgit v1.2.3-70-g09d2