From 4ff5e31b715a62a2cc62a1a8f8beab6c4c2234d6 Mon Sep 17 00:00:00 2001 From: Ben Green Date: Thu, 18 Jun 2015 00:24:29 -0400 Subject: altered spawning model and changed temporal probability on weight_success --- experiments/process.py | 26 +++++++++++++------------- 1 file changed, 13 insertions(+), 13 deletions(-) (limited to 'experiments/process.py') 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} ".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 -- cgit v1.2.3-70-g09d2