From a473003961419502b66b5111374de26331bf4fc3 Mon Sep 17 00:00:00 2001 From: Ben Green Date: Sat, 20 Jun 2015 18:33:48 -0400 Subject: More tests with ml2 to find optimal results --- experiments/process.py | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) (limited to 'experiments/process.py') diff --git a/experiments/process.py b/experiments/process.py index d01b3fd..330a49e 100644 --- a/experiments/process.py +++ b/experiments/process.py @@ -12,13 +12,13 @@ if __name__ == "__main__": sys.exit("usage: {0} ".format(sys.argv[0])) root_victims, victims, non_victims, age = load(open(sys.argv[1])) - alphas = np.arange(10., 30., 3.) # parameter of the time component - # alphas = np.logspace(0,2,num=10) - # deltas = np.arange(0.005, 0.5, 0.005) # parameter of the structural component - deltas = np.logspace(-4,-2,num=5) + # alphas = np.arange(5e-4, 3e-3, 1e-4) # parameter of the time component + alphas = np.logspace(-5,-1,num=20) + # deltas = np.arange(0.03, 0.1, 0.005) # parameter of the structural component + deltas = np.logspace(-5,-1,num=20) 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) + beta, roots, ll = ml2(root_victims, victims, non_victims, age, alpha, delta) print "\t".join(map(str, [alpha, delta, beta, roots, ll])) fh.write("\t".join(map(str, [alpha, delta, beta, roots, ll])) + "\n") fh.flush() -- cgit v1.2.3-70-g09d2