From 9350ee2c6359562a23cf8efdefdd7de80b2a682e Mon Sep 17 00:00:00 2001 From: Ben Green Date: Sat, 1 Aug 2015 12:01:19 -0400 Subject: minor updates --- experiments/process.py | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) (limited to 'experiments/process.py') diff --git a/experiments/process.py b/experiments/process.py index ffa3162..f339fd1 100644 --- a/experiments/process.py +++ b/experiments/process.py @@ -16,12 +16,14 @@ if __name__ == "__main__": root_victims, victims, non_victims, age = load(open(sys.argv[1])) # alphas = np.arange(1, 200, 10) # parameter of the time component - alphas = np.logspace(-3,1,num=16) + alphas = np.logspace(-6,1,num=5) # deltas = np.arange(.00003, .0003, .00002) # parameter of the structural component - deltas = np.logspace(-4,-1,num=12) - lmbda = 1. + deltas = np.logspace(-6,-1,num=5) + lmbdas = np.logspace(-8,-3,num=5) + # lmbda = 1 with open("out.log", "w") as fh: - for alpha, delta in product(alphas, deltas): + for alpha, delta, lmbda in product(alphas, deltas, lmbdas): + # for alpha, delta in product(alphas, deltas): beta, roots, ll = ml(root_victims, victims, non_victims, alpha, delta, lmbda) print "\t".join(map(str, [alpha, delta, beta, roots, ll])) fh.write("\t".join(map(str, [alpha, delta, beta, roots, ll])) + "\n") -- cgit v1.2.3-70-g09d2